The Lapse, Shift, UnitNormalVector and TimeVector are related quantities that, essentially, represent the 4-coordinate degrees of freedom that exist in the representation of a gravitational field (the Lapse and Shift), and identify a particular choice of 3D spacelike hypersurface in the decomposition of spacetime into 3 + 1, its space and time parts (UnitNormalVector and TimeVector). This decomposition is relevant in the context of numerical relativity, and in a Hamiltonian description of gravity.

Regarding the Lapse, Shift, you can set their values using Setup and its lapseandshift keyword; there are three possible values: standard, arbitrary, or a list of algebraic expressions representing them. The value chosen determines the values of the components of all the ThreePlusOne tensors (e.g. ExtrinsicCurvature) and so of tensorial expressions involving them (e.g. the ADMEquations). Those components are always computed first in terms of the Lapse and Shift and the space part $g_{i,j}$ of the 4D metric, then in a second step, if lapseandshift = standard, the Lapse and the Shift are replaced by their expressions in terms of the $g_{0,\mathrm{\mu }}$ part of the 4D metric, according to ${\mathbf{\alpha }}^2={\left(-g_{}^{0,0}\right)}^{-1}$ and ${\mathbf{\beta }}_j=g_{0,j}$, or if lapseandshift was set passing a list to Setup then in terms of the values indicated in that list (these values are not used to set the $g_{0,\mathrm{\mu }}$ part of the 4D metric). When lapseandshift = arbitrary, the second step is not performed and the Lapse and Shift evaluate to themselves, representing an arbitrary value for them. In the three cases, standard, arbitrary, or a list of algebraic expressions, the components of the ThreePlusOne tensors are computed in terms of a metric with line element ${\mathrm{ds}}^2=-{\mathit{\alpha }}^2{\mathrm{dt}}^2+{\mathit{\gamma }}_{i,j}\left({\mathrm{dx}}^i+{\mathit{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}\mathrm{dt}\right)\left({\mathrm{dx}}^j+{\mathit{\beta }}_{\phantom{}\phantom{j}}^{\phantom{}j}\mathrm{dt}\right)$, where $\mathbf{\alpha }$ and ${\mathbf{\beta }}_{}^i$ have the Lapse and Shift values mentioned, and ${\mathbf{\gamma }}_{i,j}=g_{i,j}$ is the ThreePlusOne:-gamma3_ metric of the 3D hypersurface. This design permits working with any 4D metric set and, without changing its value, experimenting with different values of the Lapse and Shift (different values of $g_{0,\mathrm{\mu }}$) for the 3+1 decomposition of Einstein's equations. See also LapseAndShiftConditions.

The presentation in the rest of this page assumes the signature (+ + + -), automatically set when ThreePlusOne is loaded and the denominations and definitions are the standard ones of the ADM formalism (ref [5]). The signature (+ + + -) is physically equivalent to (- + + +), but has the computational advantage of having space indices running from 1 to 3 (instead of from 2 to 4) and the value 0 of a spacetime index pointing to position 4 (instead of position 1). With any of these two signatures, the space components (i, j) of the 4D gamma3_${}_{\mathrm{\mu },\mathrm{\nu }}$ are related to the space components (i, j) of the 4D spacetime metric $g_{\mathrm{\mu },\mathrm{\nu }}$ by

where, as is the case of all the tensors of ThreePlusOne, the 3D metric is displayed in black instead of blue.

From a computational point of view, the Lapse is a scalar, and the Shift is a purely spatial 4D vector, that is, its contravariant timelike component is equal to 0. The UnitNormalVector is a timelike vector (all its covariant space components are zero), perpendicular to the hypersurface determined by the Lapse and Shift, and the TimeVector is neither spacelike nor timelike, but in the adapted coordinates used its contravariant space components are zero. All of them can be defined in terms of a global time function $\mathrm{{\rm T}}\left(x^i\,t\right)=\mathrm{constant}$. In a generic system of references, where a 3+1 split formulates spacetime as a 3D spacelike hypersurface parameterized by $\mathrm{{\rm T}}$ that evolves in time, the Lapse, assumed to be positive, is the lapse of proper time between two 3D neighboring hypersurfaces, measured by observers moving along the direction normal to the hypersurfaces (normal, or Eulerian observers). In the adapted coordinates used in the implementation in Physics, $\mathrm{{\rm T}}\left(x^i\,t\right)=t$. The Lapse is thus defined as

where $▿$ is the 4D covariant derivative operator D_. The UnitNormalVector is a timelike unit vector, normal to the 3D hypersurface, i.e. it is the 4-velocity of a normal observer whose worldline is normal to the 3D hypersurface. The UnitNormalVector is defined in terms of $\mathrm{{\rm T}}\left(x^i\,t\right)$ as

from where $\mathrm{n\_\_\μ}{\mathit{n}}_{}^{\mathrm{\mu }}=\mathrm{-1}$. The 3D metric gamma3_, that contains information on the space part of the 4D metric g_ and that projects into this 3D spacelike hypersurface determined by $\mathrm{{\rm T}}\left(x^i\,t\right)=\mathrm{constant}$, is also expressed in terms of g_ and the UnitNormalVector as

The TimeVector is tangent to the time lines (lines of constant spatial coordinates), so it can also be expressed in terms of the global function $\mathrm{{\rm T}}\left(x^i\,t\right)$ identifying the 3D hypersurface as

from where ${\mathit{t}}_{}^{\mathrm{\mu }}{\mathit{n}}_{\mathrm{\mu }}=-\mathbf{\alpha }$. The Shift can be defined as the projection of the TimeVector onto the 3D hypersurface, i.e.

from where it follows the relationship between these three vectors

From a physical point of view, choosing $\mathrm{{\rm T}}\left(x^i\,t\right)=t$ as the time coordinate, the Lapse and Shift arise naturally when defining simultaneity in the general theory of relativity as the possibility of synchronizing clocks located at different points in space (ref [1], sec. (84)). For that synchronization to be possible, all the components $g_{0,j}$ of the spacetime metric must be equal to zero, and that condition can always be achieved, in any gravitational field, by an appropriate choice of coordinates (ref [1], sec. (84) and (97)). If, in addition, $g_{0,0}=1$, the reference system is synchronous. In this referential the time lines (the TimeVector is tangent to them) are geodesics in the 4D spacetime and normal to the hypersurface $t=\mathrm{constant}$ (because from $g_{0,j}=0$ follows that ${\mathbf{\beta }}_{}^{\mathrm{\mu }}=0$ and so  ${\mathit{t}}_{}^{\mathrm{\mu }}$ and  ${\mathit{n}}_{}^{\mathrm{\mu }}$ are parallel), and the 4D infinitesimal interval is written as.

Under a time rescaling $t\text{'}=G\left(t\right)$ and a general transformation of the space coordinates ${\left(x\text{'}\right)}^i=F^i\left(x^j\,t\right)$, and taking into account ${\mathbf{\gamma }}_{i,j}=g_{i,j}$, this infinitesimal interval changes into

where now ${\mathbf{\beta }}_{}^{\mathrm{\mu }}\ne 0$ and so the new time lines $t=\mathrm{constant}$ (and hence ${\mathit{t}}_{}^{\mathrm{\mu }}$) are no longer orthogonal to the 3D hypersurface whose geometry is described by $\mathrm{\γ\_\_i,j}$. In this decomposition of the infinitesimal interval, the 4-coordinate degrees of freedom are represented by the freedom in choosing the Lapse and the Shift (corresponding to different choices of the system of references; conversely, one can depart from this general form of the infinitesimal interval and construct a transformation to a synchronous reference system, where $g_{0,0}=1$ and $g_{0,j}=0$, by solving a related Hamilton-Jacobi equation see [1] sec.(97)). In this system of references (so-called adapted coordinates system of references), the Lapse and Shift are expressed in terms of the contravariant spatial components of the 4D spacetime metric g_ as

Besides being indexed with a spacetime or a space index, the Shift, TimeVector and UnitNormalVector accept the following keywords as an index:

definition: when passed alone, gamma3_ returns its 4D definition in terms of the UnitNormalVector and the spacetime metric g_. When passed with space indices, it returns its expression in terms of the spatial components of g_.

matrix: (synonym: Matrix, array, Array, or no indices whatsoever, as in Shift[]) returns an Array that when indexed with numerical values from 1 to the dimension of spacetime it returns the value of each of the components of Shift[mu]. If this keyword is passed together with indices, that can be covariant or contravariant, for example Shift[~mu, matrix], the resulting array takes into account the character of the indices.

nonzero: returns a set of equations, with the left-hand side as a sequence of two positive numbers identifying the element of $\mathrm{\γ\_\_i,j}$ and the corresponding value on the right-hand side. Note that this set is actually the output of the ArrayElems command when passing to it the Matrix obtained with the keyword matrix.

You can change the value of the 3D metric gamma3_ by changing the value of the 4D metric g_, using Setup or g_. Regarding the geometry in this 3D hypersurface, a 3D covariant derivative, related Christoffel symbols, and the Ricci and Riemann tensors are represented respectively by the D3_, Christoffel3, Ricci3 and Riemann3 commands.

The %Lapse, %Shift, %TimeVector and %UnitNormalVector command are the inert forms of the corresponding active commands, that is represent the same object but entering them does not result in performing any computation. To perform the related computations, e.g. as if %Lapse were Lapse, use value.

$\mathrm{with}\left(\mathrm{Physics}\right)\:$

$\mathrm{Setup}\left(\mathrm{mathematicalnotation}=\mathrm{true}\,\mathrm{coordinates}=\mathrm{cartesian}\right)$

$\mathrm{Systems\; of\; spacetime\; coordinates\; are:}\left\{X=\left(x\,y\,z\,t\right)\right\}$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$\left[\mathrm{coordinatesystems}=\left\{X\right\}\,\mathrm{mathematicalnotation}=\mathrm{true}\right]$

$\mathrm{with}\left(\mathrm{ThreePlusOne}\right)\:$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$\mathrm{Setting}\mathrm{lowercaselatin\_is}\mathrm{letters\; to\; represent}\mathrm{space}\mathrm{indices}$

$\mathrm{Defined\; as\; 4D\; spacetime\; tensors}\left(\mathrm{see\; ?Physics,ThreePlusOne}\right),{\mathbf{\gamma }}_{\mathrm{\mu },\mathrm{\nu }},{\mathit{▿}}_{\mathrm{\mu }},{\mathbf{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }},{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }},{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }},{\mathbf{\beta }}_{\mathrm{\mu }},{\mathit{n}}_{\mathrm{\mu }},{\mathit{t}}_{\mathrm{\mu }},{\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}$

$\mathrm{Changing\; the\; signature\; of\; spacetime\; from}\left(\mathrm{-\; -\; -\; +}\right)\mathrm{to}\left(\mathrm{+\; +\; +\; -}\right)\mathrm{in\; order\; to\; match\; the\; signature\; customarily\; used\; in\; the\; ADM\; formalism}$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$\mathrm{g\_}\left[\right]$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{Setup}\left(\mathrm{spacetimeindices}\,\mathrm{spaceindices}\right)$

$\left[\mathrm{spaceindices}=\mathrm{lowercaselatin\_is}\,\mathrm{spacetimeindices}=\mathrm{greek}\right]$

$\mathrm{gamma3\_}\left[\mathrm{definition}\right]$

${\mathbf{\gamma }}_{\mathrm{\mu },\mathrm{\nu }}={\mathit{n}}_{\mathrm{\mu }}{\mathit{n}}_{\mathrm{\nu }}+g_{\mathrm{\mu },\mathrm{\nu }}$

$\mathrm{gamma3\_}\left[i,j,\mathrm{definition}\right]$

${\mathbf{\gamma }}_{i,j}=g_{i,j}$

$\mathrm{UnitNormalVector}\left[\mathrm{definition}\right]$

${\mathit{n}}_{\mathrm{\mu }}=-\mathbf{\alpha }{▿}_{\mathrm{\mu }}\left(t\right)$

$\mathrm{TimeVector}\left[\mathrm{definition}\right]$

${\mathit{t}}_{\mathrm{\mu }}=\mathbf{\alpha }{\mathit{n}}_{\mathrm{\mu }}+{\mathbf{\beta }}_{\mathrm{\mu }}$

$\mathrm{ExtrinsicCurvature}\left[\mathrm{definition}\right]$

${\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}=-\frac{{ℒ}_{\mathit{n}}\left({\mathbf{\gamma }}_{\mathrm{\mu },\mathrm{\nu }}\right)}{2}$

$\mathrm{g\_}\left[\mathrm{sc}\right]$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$\mathrm{Systems\; of\; spacetime\; coordinates\; are:}\left\{X=\left(r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right)\right\}$

$\mathrm{Default\; differentiation\; variables\; for\; d\_,\; D\_\; and\; dAlembertian\; are:}\left\{X=\left(r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right)\right\}$

$\mathrm{The\; Schwarzschild\; metric\; in\; coordinates}\left[r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right]$

$\mathrm{Parameters:}\left[m\right]$

$\mathrm{Signature:}\left(\mathrm{+\; +\; +\; -}\right)$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{Shift}\left[\right]$

${\mathbf{\beta }}_{\mathrm{\mu }}=\left[\right]$

$\mathrm{UnitNormalVector}\left[\right]$

${\mathit{n}}_{\mathrm{\mu }}=\left[\right]$

$\mathrm{TimeVector}\left[\right]$

${\mathit{t}}_{\mathrm{\mu }}=\left[\right]$

$N\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{UnitNormalVector}\left[\mathrm{\mu }\right]\mathrm{UnitNormalVector}\left[\mathrm{\nu }\right]$

$N_{\mathrm{\mu },\mathrm{\nu }}={\mathit{n}}_{\mathrm{\mu }}{\mathit{n}}_{\mathrm{\nu }}$

$\mathrm{Define}\left(\right)$

$\mathrm{Defined\; objects\; with\; tensor\; properties}$

$\left\{{\mathit{▿}}_{\mathrm{\mu }}\,{▿}_{\mathrm{\mu }}\,{\mathrm{\gamma }}_{\mathrm{\mu }}\,N_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\sigma }}_{\mathrm{\mu }}\,R_{\mathrm{\mu },\mathrm{\nu }}\,{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }}\,R_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\mathbf{\beta }}_{\mathrm{\mu }}\,C_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\partial }_{\mathrm{\mu }}\,g_{\mathrm{\mu },\mathrm{\nu }}\,{\mathbf{\gamma }}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\gamma }}_{i,j}\,{\mathrm{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,{\mathbf{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,G_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\,{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\mathit{t}}_{\mathrm{\mu }}\,{\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathit{n}}_{\mathrm{\mu }}\,X_{\mathrm{\mu }}\right\}$

$N\left[\right]$

$N_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$N\left[\mathrm{`~`}\right]$

$N_{\phantom{}\phantom{\mathrm{\mu },\mathrm{\nu }}}^{\phantom{}\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{g\_}\left[\mathrm{\mu },\mathrm{\nu }\right]=\left(\mathrm{gamma3\_}\left[\mathrm{\mu },\mathrm{\alpha }\right]\mathrm{gamma3\_}\left[\mathrm{\nu },\mathrm{\beta }\right]+N\left[\mathrm{\mu },\mathrm{\alpha }\right]N\left[\mathrm{\nu },\mathrm{\beta }\right]\right)\mathrm{g\_}\left[\mathrm{~alpha},\mathrm{~beta}\right]$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left(N_{\mathrm{\alpha },\mathrm{\mu }}{N}_{\mathrm{\beta },\mathrm{\nu }}+{\mathbf{\gamma }}_{\mathrm{\alpha },\mathrm{\mu }}{\mathbf{\gamma }}_{\mathrm{\beta },\mathrm{\nu }}\right)g_{\phantom{}\phantom{\mathrm{\alpha },\mathrm{\beta }}}^{\phantom{}\mathrm{\alpha },\mathrm{\beta }}$

$\mathrm{TensorArray}\left(\right)$

$\mathrm{Ricci}\left[\right]$

$R_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{Ricci3}\left[\right]$

${\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{Define}\left(\mathrm{dx}\left[\mathrm{\mu }\right]\right)$

$\mathrm{Defined\; objects\; with\; tensor\; properties}$

$\left\{{\mathit{▿}}_{\mathrm{\mu }}\,{▿}_{\mathrm{\mu }}\,{\mathrm{\gamma }}_{\mathrm{\mu }}\,N_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\sigma }}_{\mathrm{\mu }}\,R_{\mathrm{\mu },\mathrm{\nu }}\,{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }}\,R_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\mathbf{\beta }}_{\mathrm{\mu }}\,C_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\partial }_{\mathrm{\mu }}\,{\mathrm{dx}}_{\mathrm{\mu }}\,g_{\mathrm{\mu },\mathrm{\nu }}\,{\mathbf{\gamma }}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\gamma }}_{i,j}\,{\mathrm{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,{\mathbf{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,G_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\,{\mathit{R}}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\mathit{t}}_{\mathrm{\mu }}\,{\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathit{n}}_{\mathrm{\mu }}\,X_{\mathrm{\mu }}\right\}$

$\mathrm{alias}\left(\mathrm{dt}=\mathrm{dx}\left[\mathrm{~0}\right]\right)\:$

${\mathrm{ds}}^2=\left(-{\mathrm{Lapse}}^2+{\mathrm{Shift}\left[i\right]}^2\right){\mathrm{dt}}^2+2\mathrm{Shift}\left[i\right]\mathrm{dt}\mathrm{dx}\left[\mathrm{~i}\right]+\mathrm{gamma3\_}\left[i,j\right]\mathrm{dx}\left[\mathrm{~i}\right]\mathrm{dx}\left[\mathrm{~j}\right]$

${\mathrm{ds}}^2={\mathrm{dt}}^2\left(-{\mathbf{\alpha }}^2+{\mathbf{\beta }}_i{\mathbf{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}\right)+2{\mathbf{\beta }}_i\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}+{\mathbf{\gamma }}_{i,j}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}$

${\mathrm{ds}}^2=\mathrm{\%g\_}\left[\mathrm{\mu },\mathrm{\nu }\right]\mathrm{dx}\left[\mathrm{~mu}\right]\mathrm{dx}\left[\mathrm{~nu}\right]$

${\mathrm{ds}}^2=g_{\mathrm{\mu },\mathrm{\nu }}{\mathrm{dx}}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}{\mathrm{dx}}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}$

$\mathrm{Decompose}\left(\right)$

${\mathrm{ds}}^2=g_{0,0}{\mathrm{dt}}^2+g_{0,j}\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}+g_{i,0}\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}+g_{i,j}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}$

$\mathrm{op}\left(2\,\mathrm{rhs}\left(\right)\right)=\mathrm{op}\left(3\,\mathrm{rhs}\left(\right)\right)$

$g_{0,j}\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}=g_{i,0}\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}$

$\mathrm{subs}\left(\,\right)$

${\mathrm{ds}}^2=g_{0,0}{\mathrm{dt}}^2+2g_{i,0}\mathrm{dt}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}+g_{i,j}{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}$

$\mathrm{simplify}\left(-\right)$

$0=\left({\mathbf{\alpha }}^2-{\mathbf{\beta }}_i{\mathbf{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}+g_{0,0}\right){\mathrm{dt}}^2+2{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}\left(g_{i,0}-{\mathbf{\beta }}_i\right)\mathrm{dt}-{\mathrm{dx}}_{\phantom{}\phantom{i}}^{\phantom{}i}{\mathrm{dx}}_{\phantom{}\phantom{j}}^{\phantom{}j}\left({\mathbf{\gamma }}_{i,j}-g_{i,j}\right)$

$\mathrm{eq}\left[1\right]≔\mathrm{coeff}\left(\mathrm{coeff}\left(\mathrm{rhs}\left(\right)\,\mathrm{dt}\right)\,\mathrm{dx}\left[\mathrm{~i}\right]\right)=0$

${\mathrm{eq}}_1≔2g_{i,0}-2{\mathbf{\beta }}_i=0$

$\mathrm{eq}\left[2\right]≔\mathrm{coeff}\left(\mathrm{rhs}\left(\right)\,{\mathrm{dt}}^2\right)=0$

${\mathrm{eq}}_2≔{\mathbf{\alpha }}^2-{\mathbf{\beta }}_i{\mathbf{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}+g_{0,0}=0$

$\mathrm{eq}\left[3\right]≔\mathrm{coeff}\left(\mathrm{coeff}\left(\mathrm{rhs}\left(\right)\,\mathrm{dx}\left[\mathrm{~i}\right]\right)\,\mathrm{dx}\left[\mathrm{~j}\right]\right)=0$

${\mathrm{eq}}_3≔-{\mathbf{\gamma }}_{i,j}+g_{i,j}=0$

$\mathrm{isolate}\left(\mathrm{eq}\left[1\right]\,\mathrm{Shift}\left[i\right]\right)$

${\mathbf{\beta }}_i=g_{i,0}$

$\mathrm{isolate}\left(\mathrm{eq}\left[2\right]\,{\mathrm{Lapse}}^2\right)$

${\mathbf{\alpha }}^2={\mathbf{\beta }}_i{\mathbf{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}-g_{0,0}$

$\mathrm{isolate}\left(\mathrm{eq}\left[3\right]\,\mathrm{gamma3\_}\left[i,j\right]\right)$

${\mathbf{\gamma }}_{i,j}=g_{i,j}$

$\mathrm{Setup}\left(\mathrm{lapseandshift}\right)$

$\left[\mathrm{lapseandshift}=\mathrm{standard}\right]$

$\mathrm{ExtrinsicCurvature}\left[\right]$

${\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{Setup}\left(\mathrm{lapseandshift}=\mathrm{arbitrary}\right)$

$\left[\mathrm{lapseandshift}=\mathrm{arbitrary}\right]$

$\mathrm{ExtrinsicCurvature}\left[\right]$

${\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{\%UnitNormalVector}\left[0\right]$

${\mathit{n}}_0$

$\mathrm{value}\left(\right)$

$-\mathbf{\alpha }$

$\mathrm{Setup}\left(\mathrm{lapseandshift}=\left[1\,\left[0\,0\,0\right]\right]\right)$

$\left[\mathrm{lapseandshift}=\left[1\,\left[0\,0\,0\right]\right]\right]$

$\mathrm{g\_}\left[\right]$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left[\right]$

$\mathrm{UnitNormalVector}\left[\right]$

${\mathit{n}}_{\mathrm{\mu }}=\left[\right]$

$\mathrm{Christoffel}\left[\mathrm{definition}\right]$

${\mathrm{\Gamma }}_{\mathrm{\alpha },\mathrm{\mu },\mathrm{\nu }}=\frac{{\partial }_{\mathrm{\nu }}\left(g_{\mathrm{\alpha },\mathrm{\mu }}\right)}{2}+\frac{{\partial }_{\mathrm{\mu }}\left(g_{\mathrm{\alpha },\mathrm{\nu }}\right)}{2}-\frac{{\partial }_{\mathrm{\alpha }}\left(g_{\mathrm{\mu },\mathrm{\nu }}\right)}{2}$

$\mathrm{Christoffel3}\left[\mathrm{definition}\right]$

${\mathbf{\Gamma }}_{\mathrm{\alpha },\mathrm{\mu },\mathrm{\nu }}=\frac{{\mathbf{\gamma }}_{\mathrm{\alpha }\phantom{\mathrm{\beta }}}^{\phantom{\mathrm{\alpha }}\mathrm{\beta }}{\mathbf{\gamma }}_{\mathrm{\mu }\phantom{\mathrm{\kappa }}}^{\phantom{\mathrm{\mu }}\mathrm{\kappa }}{\mathbf{\gamma }}_{\mathrm{\nu }\phantom{\mathrm{\lambda }}}^{\phantom{\mathrm{\nu }}\mathrm{\lambda }}\left({\partial }_{\mathrm{\lambda }}\left({\mathbf{\gamma }}_{\mathrm{\beta },\mathrm{\kappa }}\right)+{\partial }_{\mathrm{\kappa }}\left({\mathbf{\gamma }}_{\mathrm{\beta },\mathrm{\lambda }}\right)-{\partial }_{\mathrm{\beta }}\left({\mathbf{\gamma }}_{\mathrm{\kappa },\mathrm{\lambda }}\right)\right)}{2}$

$\mathrm{Ricci}\left[\mathrm{definition}\right]$

$R_{\mathrm{\mu },\mathrm{\nu }}={\partial }_{\mathrm{\alpha }}\left({\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\mu },\mathrm{\nu }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\mu },\mathrm{\nu }}}\right)-{\partial }_{\mathrm{\nu }}\left({\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\mu },\mathrm{\alpha }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\mu },\mathrm{\alpha }}}\right)+{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\beta }}\mathrm{\mu },\mathrm{\nu }}^{\phantom{}\mathrm{\beta }\phantom{\mathrm{\mu },\mathrm{\nu }}}{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\beta },\mathrm{\alpha }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\beta },\mathrm{\alpha }}}-{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\beta }}\mathrm{\mu },\mathrm{\alpha }}^{\phantom{}\mathrm{\beta }\phantom{\mathrm{\mu },\mathrm{\alpha }}}{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\nu },\mathrm{\beta }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\nu },\mathrm{\beta }}}$

$\mathrm{Ricci3}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{convert}\left(\mathrm{Ricci3}\left[\mathrm{\mu },\mathrm{\nu }\right]\,\mathrm{Christoffel}\right)$

${\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }}={\mathbf{\gamma }}_{\mathrm{\mu }\phantom{\mathrm{\sigma }}}^{\phantom{\mathrm{\mu }}\mathrm{\sigma }}{\mathbf{\gamma }}_{\mathrm{\nu }\phantom{\mathrm{\tau }}}^{\phantom{\mathrm{\nu }}\mathrm{\tau }}\left({\partial }_{\mathrm{\alpha }}\left({\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\sigma },\mathrm{\tau }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\sigma },\mathrm{\tau }}}\right)-{\partial }_{\mathrm{\tau }}\left({\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\alpha },\mathrm{\sigma }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\alpha },\mathrm{\sigma }}}\right)+{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\beta }}\mathrm{\sigma },\mathrm{\tau }}^{\phantom{}\mathrm{\beta }\phantom{\mathrm{\sigma },\mathrm{\tau }}}{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\alpha },\mathrm{\beta }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\alpha },\mathrm{\beta }}}-{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\beta }}\mathrm{\alpha },\mathrm{\sigma }}^{\phantom{}\mathrm{\beta }\phantom{\mathrm{\alpha },\mathrm{\sigma }}}{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\beta },\mathrm{\tau }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\beta },\mathrm{\tau }}}\right)$

$\mathrm{Riemann}\left[\mathrm{definition}\right]$

$R_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}=g_{\mathrm{\alpha },\mathrm{\lambda }}\left({\partial }_{\mathrm{\mu }}\left({\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\beta },\mathrm{\nu }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\beta },\mathrm{\nu }}}\right)-{\partial }_{\mathrm{\nu }}\left({\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\beta },\mathrm{\mu }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\beta },\mathrm{\mu }}}\right)+{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\upsilon },\mathrm{\mu }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\upsilon },\mathrm{\mu }}}{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\upsilon }}\mathrm{\beta },\mathrm{\nu }}^{\phantom{}\mathrm{\upsilon }\phantom{\mathrm{\beta },\mathrm{\nu }}}-{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\upsilon },\mathrm{\nu }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\upsilon },\mathrm{\nu }}}{\mathrm{\Gamma }}_{\phantom{}\phantom{\mathrm{\upsilon }}\mathrm{\beta },\mathrm{\mu }}^{\phantom{}\mathrm{\upsilon }\phantom{\mathrm{\beta },\mathrm{\mu }}}\right)$

$\mathrm{Riemann3}\left[\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{convert}\left(\mathrm{Riemann3}\left[\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }\right]\,\mathrm{Christoffel}\right)$

${\mathit{R}}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}={\mathbf{\gamma }}_{\mathrm{\alpha },\mathrm{\lambda }}{\mathbf{\gamma }}_{\mathrm{\beta }\phantom{\mathrm{\sigma }}}^{\phantom{\mathrm{\beta }}\mathrm{\sigma }}{\mathbf{\gamma }}_{\mathrm{\mu }\phantom{\mathrm{\tau }}}^{\phantom{\mathrm{\mu }}\mathrm{\tau }}{\mathbf{\gamma }}_{\mathrm{\nu }\phantom{\mathrm{\kappa }}}^{\phantom{\mathrm{\nu }}\mathrm{\kappa }}\left({\partial }_{\mathrm{\tau }}\left({\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\kappa },\mathrm{\sigma }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\kappa },\mathrm{\sigma }}}\right)-{\partial }_{\mathrm{\kappa }}\left({\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\sigma },\mathrm{\tau }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\sigma },\mathrm{\tau }}}\right)+{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\tau },\mathrm{\upsilon }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\tau },\mathrm{\upsilon }}}{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\upsilon }}\mathrm{\kappa },\mathrm{\sigma }}^{\phantom{}\mathrm{\upsilon }\phantom{\mathrm{\kappa },\mathrm{\sigma }}}-{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\lambda }}\mathrm{\kappa },\mathrm{\upsilon }}^{\phantom{}\mathrm{\lambda }\phantom{\mathrm{\kappa },\mathrm{\upsilon }}}{\mathbf{\Gamma }}_{\phantom{}\phantom{\mathrm{\upsilon }}\mathrm{\sigma },\mathrm{\tau }}^{\phantom{}\mathrm{\upsilon }\phantom{\mathrm{\sigma },\mathrm{\tau }}}\right)$

$\mathrm{Ricci3}\left[\mathrm{definition}\right]$

${\mathit{R}}_{\mathrm{\mu },\mathrm{\nu }}={\mathbf{\gamma }}_{\phantom{}\phantom{\mathrm{\alpha },\mathrm{\beta }}}^{\phantom{}\mathrm{\alpha },\mathrm{\beta }}{\mathbf{\gamma }}_{\mathrm{\mu }\phantom{\mathrm{\kappa }}}^{\phantom{\mathrm{\mu }}\mathrm{\kappa }}{\mathbf{\gamma }}_{\mathrm{\alpha }\phantom{\mathrm{\lambda }}}^{\phantom{\mathrm{\alpha }}\mathrm{\lambda }}{\mathbf{\gamma }}_{\mathrm{\nu }\phantom{\mathrm{\sigma }}}^{\phantom{\mathrm{\nu }}\mathrm{\sigma }}{R}_{\mathrm{\beta },\mathrm{\kappa },\mathrm{\lambda },\mathrm{\sigma }}+{\mathbf{{\rm K}}}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\nu }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\nu }}}{\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\alpha }}-{\mathbf{{\rm K}}}_{\phantom{}\phantom{\mathrm{\alpha }}\mathrm{\alpha }}^{\phantom{}\mathrm{\alpha }\phantom{\mathrm{\alpha }}}{\mathbf{{\rm K}}}_{\mathrm{\mu },\mathrm{\nu }}$

$\mathrm{Riemann3}\left[\mathrm{definition}\right]$

${\mathit{R}}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}={\mathbf{\gamma }}_{\mathrm{\alpha }\phantom{\mathrm{\kappa }}}^{\phantom{\mathrm{\alpha }}\mathrm{\kappa }}{\mathbf{\gamma }}_{\mathrm{\beta }\phantom{\mathrm{\lambda }}}^{\phantom{\mathrm{\beta }}\mathrm{\lambda }}{\mathbf{\gamma }}_{\mathrm{\mu }\phantom{\mathrm{\sigma }}}^{\phantom{\mathrm{\mu }}\mathrm{\sigma }}{\mathbf{\gamma }}_{\mathrm{\nu }\phantom{\mathrm{\tau }}}^{\phantom{\mathrm{\nu }}\mathrm{\tau }}{R}_{\mathrm{\kappa },\mathrm{\lambda },\mathrm{\sigma },\mathrm{\tau }}-{\mathbf{{\rm K}}}_{\mathrm{\alpha },\mathrm{\mu }}{\mathbf{{\rm K}}}_{\mathrm{\beta },\mathrm{\nu }}+{\mathbf{{\rm K}}}_{\mathrm{\alpha },\mathrm{\nu }}{\mathbf{{\rm K}}}_{\mathrm{\beta },\mathrm{\mu }}$

[1] Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2, fourth revised English edition. Elsevier, 1975.

[2] Alcubierre, M., Introduction to 3+1 Numerical Relativity, International Series of Monographs on Physics 140, Oxford University Press, 2008.

[3] Baumgarte, T.W., Shapiro, S.L., Numerical Relativity, Solving Einstein's Equations on a Computer, Cambridge University Press, 2010.

[4] Gourgoulhon, E., 3+1 Formalism and Bases of Numerical Relativity, Lecture notes, 2007, https://arxiv.org/pdf/gr-qc/0703035v1.pdf.

[5] Arnowitt, R., Dese, S., Misner, C.W., The Dynamics of General Relativity, Chapter 7 in Gravitation: an introduction to current research (Wiley, 1962), https://arxiv.org/pdf/gr-qc/0405109v1.pdf

The Physics[ThreePlusOne][Lapse], Physics[ThreePlusOne][Shift], Physics[ThreePlusOne][UnitNormalVector] and Physics[ThreePlusOne][TimeVector] commands were introduced in Maple 2017.

For more information on Maple 2017 changes, see Updates in Maple 2017.

The Physics[ThreePlusOne][Lapse] command was introduced in Maple 2020.

For more information on Maple 2020 changes, see Updates in Maple 2020.