 LapseAndShiftConditions - Maple Help

Home : Support : Online Help : Physics : ThreePlusOne : LapseAndShiftConditions

Physics[ThreePlusOne][LapseAndShiftConditions] Calling Sequence LapseAndShiftConditions() LapseAndShiftConditions(keyword, options) Parameters

 keyword - optional : any of the valid keywords associated to slicing and spatial gauge conditions, returned when no arguments are given additionalconditions = ... - optional : the right-hand side is a set of equations to be added to the set of equations related to keyword, to be returned or solved evaluate - optional : when given, the set of equations to be returned if first evaluated using value, then an array of equations is constructed with the tensor components output = ... - optional, the right-hand side can be equations (default) or solution, in which case an attempt is made to solve the equations for the Lapse and Shift parameters = ... - optional, the right-hand side is a set to be considered additional solving parameters (typically entering the components of the spacetime metric) when solving for the Lapse and Shift quiet - optional, when given, verbosity is turned OFF. setconditions - optional, when given, an attempt is made to solve the slicing and gauge equations for the Lapse and Shift, and when successful, they are set using Setup and its lapseandshift keyword. simplifier - optional, the indication of a simplifier to be applied to each of the condition equations being returned or solved. Description

 • The LapseAndShiftConditions() command returns a list of keywords associated to the most common slicing and spatial gauge conditions for the Lapse and Shift that fix the freedom of choice of coordinates of the 4D formulation of Einstein's equations describing gravity.
 • The LapseAndShiftConditions(keyword) returns a Vector column with the equations for the Lapse and Shift associated to that keyword, expressed in terms of the inert forms of the tensors involved; these are mainly the ExtrinsicCurvature, the gamma3_ metric of the spacelike hypersurface of the 3+1 splitting, the NormalUnitVector, and the EnergyMomentum tensor.
 • When the option evaluate is passed, all the inert forms of tensors are first transformed into active by calling value, then the Vector column of equations is transformed into an array of equations calling TensorArray. If in addition the option simplifier = ... is given, the right-hand side of this option is applied to each of the elements of the resulting array.
 • The evaluated equations returned for the Lapse and Shift can also then be expressed in terms of the components of the 4D metric ${g}_{\mathrm{\mu },\mathrm{\nu }}$; for that purpose use $\mathrm{convert}\left(\mathrm{equations},\mathrm{g_}\right)$. This rewriting in terms of the metric can also be done for the non-evaluated equations returned by default (when no options are passed), and can be applied selectively to one or several of the tensors appearing in the equations, e.g. via $\mathrm{convert}\left(\mathrm{equations},\mathrm{g_},\mathrm{only}=\left\{\mathrm{ExtrinsicCurvature}\right\}\right)$ where the set can contain one or several tensors entering equations.
 • Depending on the keyword and the 4D metric set, in some cases it is possible to solve these slicing and gauge conditions for the Lapse and Shift exactly. To attempt that solution you can either convert the array of equations returned into a set of equations, then use PDEtools:-Solve to tackle them, or just pass the option output = solution.
 • When the equations for the Lapse and Shift can be solved, by passing the optional argument setconditions, LapseAndShiftConditions will automatically set these values using Setup and its lapseandshift keyword.
 • You can also set the value of the lapseandshift directly using Setup; there are three possible values: standard, arbitrary, or a list of algebraic expressions representing the Lapse and Shift. 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}=-{\mathbit{\alpha }}^{2}{\mathrm{dt}}^{2}+{\mathbit{\gamma }}_{i,j}\left({\mathrm{dx}}^{i}+{\mathbit{\beta }}_{\phantom{}\phantom{i}}^{\phantom{}i}\mathrm{dt}\right)\left({\mathrm{dx}}^{j}+{\mathbit{\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.
 • Depending on the manipulations to be done with the equations for the Lapse and Shift it is frequently convenient to set the value of lapseandshift via $\mathrm{Setup}\left(\mathrm{lapseandshift}=\mathrm{arbitrary}\right)$ in order to avoid their spontaneous evaluation in terms of the components of the metric ${g}_{0,0}$ and ${g}_{0,j}$ during intermediate steps. Examples

 > $\mathrm{with}\left(\mathrm{Physics}\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 }}}{,}{\mathbit{▿}}_{{\mathrm{\mu }}}{,}{\mathbf{\Gamma }}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}}{,}{\mathbit{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{\mathbit{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{\mathbf{\beta }}_{{\mathrm{\mu }}}{,}{\mathbit{n}}_{{\mathrm{\mu }}}{,}{\mathbit{t}}_{{\mathrm{\mu }}}{,}{\mathbf{Κ}}_{{\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{Systems of spacetime coordinates are:}{}{}{}\left\{X=\left(x{,}y{,}z{,}t\right)\right\}$
 ${}{}\mathrm{_______________________________________________________}$
 $\left[{\mathrm{ADMEquations}}{,}{\mathrm{Christoffel3}}{,}{\mathrm{D3_}}{,}{\mathrm{ExtrinsicCurvature}}{,}{\mathrm{Lapse}}{,}{\mathrm{LapseAndShiftConditions}}{,}{\mathrm{Ricci3}}{,}{\mathrm{Riemann3}}{,}{\mathrm{Shift}}{,}{\mathrm{TimeVector}}{,}{\mathrm{UnitNormalVector}}{,}{\mathrm{gamma3_}}\right]$ (1)
 > $\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]$ (2)

When Physics is initialized, the default spacetime metric is of Minkowski type. You can see the metric querying Setup, as in Setup(metric), or directly entering the metric as g_[], with no indices

 > $\mathrm{g_}\left[\right]$
 ${{\mathrm{g_}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{rrrr}1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\\ 0& 0& 0& -1\end{array}\right]\right)$ (3)

The list of slicing and spatial gauge conditions implemented are visible via

 > $\mathrm{LapseAndShiftConditions}\left(\right)$
 $\left[{\mathrm{advectiveonepluslog}}{,}{\mathrm{generalizedharmonicslicing}}{,}{\mathrm{geodesicslicing}}{,}{\mathrm{harmoniccoordinates}}{,}{\mathrm{harmonicslicing}}{,}{\mathrm{kdriver}}{,}{\mathrm{maximalslicing}}{,}{\mathrm{onepluslog}}\right]$ (4)

To see all the equations associated to these keywords, expressed in terms of the inert forms of the ThreePlusOne tensors, use

 > $\mathbf{for}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{key}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{in}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{LapseAndShiftConditions}\left(\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{do}\phantom{\rule[-0.0ex]{0.0em}{0.0ex}}\phantom{\rule[-0.0ex]{2.0em}{0.0ex}}\mathrm{key}=\mathrm{LapseAndShiftConditions}\left(\mathrm{key}\right)\phantom{\rule[-0.0ex]{0.0em}{0.0ex}}\mathbf{end}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{do}$
 ${\mathrm{advectiveonepluslog}}{=}\left(\left[\begin{array}{c}{\mathrm{%LieDerivative}}_{\mathrm{UnitNormalVector}}{}\left(\mathrm{Lapse}\right)=-2{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ {\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)-{\mathrm{Shift}}_{\mathrm{~j}}{}{\mathrm{%d_}}_{j}{}\left(\mathrm{Lapse}\right)=-2{}\mathrm{Lapse}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\end{array}\right]\right)$
 ${\mathrm{generalizedharmonicslicing}}{=}\left(\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)=-{\mathrm{Lapse}}^{2}{}\mathrm{%f}{}\left(\mathrm{Lapse}\right){}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ {\mathrm{Shift}}_{\mathrm{~i}}=0\end{array}\right]\right)$
 ${\mathrm{geodesicslicing}}{=}\left(\left[\begin{array}{c}\mathrm{Lapse}=1\\ {\mathrm{Shift}}_{\mathrm{~i}}=0\end{array}\right]\right)$
 ${\mathrm{harmoniccoordinates}}{=}\left(\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)-{\mathrm{Shift}}_{\mathrm{~j}}{}{\mathrm{%d_}}_{j}{}\left(\mathrm{Lapse}\right)=-{\mathrm{Lapse}}^{2}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ {\mathrm{%d_}}_{0}{}\left({\mathrm{Shift}}_{\mathrm{~i}}\right)-{\mathrm{Shift}}_{\mathrm{~j}}{}{\mathrm{%d_}}_{j}{}\left({\mathrm{Shift}}_{\mathrm{~i}}\right)=-{\mathrm{Lapse}}^{2}{}\left({\mathrm{%gamma3_}}_{\mathrm{~i},\mathrm{~j}}{}{\mathrm{%d_}}_{j}{}\left(\mathrm{%ln}{}\left(\mathrm{Lapse}\right)\right)+{\mathrm{%gamma3_}}_{\mathrm{~j},\mathrm{~k}}{}{\mathrm{%Christoffel3}}_{\mathrm{~i},j,k}\right)\end{array}\right]\right)$
 ${\mathrm{harmonicslicing}}{=}\left(\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)=-{\mathrm{Lapse}}^{2}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ \mathrm{Lapse}=\mathrm{%C}{}\left(x,y,z,t\right){}\sqrt{{\mathrm{%gamma3_}}_{\mathrm{determinant}}}\\ {\mathrm{Shift}}_{\mathrm{~i}}=0\end{array}\right]\right)$
 ${\mathrm{kdriver}}{=}\left(\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)=\mathrm{%c}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ {\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)=-\mathrm{%epsilon}{}\left({\mathrm{%d_}}_{0}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)+\mathrm{%c}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)\\ 0<\mathrm{%c},0<\mathrm{%epsilon}\end{array}\right]\right)$
 ${\mathrm{maximalslicing}}{=}\left(\left[\begin{array}{c}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}=0,{\mathrm{%d_}}_{0}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)=0\\ {\mathrm{%D_}}_{\mathrm{~mu}}{}\left({\mathrm{%D_}}_{\mathrm{μ}}{}\left(\mathrm{Lapse}\right)\right)=\mathrm{Lapse}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{μ},\mathrm{ν}}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{~mu},\mathrm{~nu}}+4{}\mathrm{π}{}\left({\mathrm{%EnergyMomentum}}_{\mathrm{ν},\mathrm{α}}{}{\mathrm{%gamma3_}}_{\mathrm{μ},\mathrm{~nu}}{}{\mathrm{%gamma3_}}_{\mathrm{~mu},\mathrm{~alpha}}+{\mathrm{%EnergyMomentum}}_{\mathrm{~mu},\mathrm{~nu}}{}{\mathrm{%UnitNormalVector}}_{\mathrm{μ}}{}{\mathrm{%UnitNormalVector}}_{\mathrm{ν}}\right)\right)\end{array}\right]\right)$
 ${\mathrm{onepluslog}}{=}\left(\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)=-2{}\mathrm{Lapse}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ \mathrm{Lapse}=1+\mathrm{ln}{}\left({\mathrm{%gamma3_}}_{\mathrm{determinant}}\right)\\ {\mathrm{Shift}}_{\mathrm{~i}}=0\end{array}\right]\right)$ (5)

These equations expressed in terms of inert tensors can be transformed into explicit equations for the Lapse and Shift using the option evaluate. For illustration purposes set the spacetime to curved, for instance by setting the Schwarzschild metric in spherical coordinates

 > $\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{_______________________________________________________}$
 ${{\mathrm{g_}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{cccc}-\frac{r}{2{}m-r}& 0& 0& 0\\ 0& {r}^{2}& 0& 0\\ 0& 0& {r}^{2}{}{\mathrm{sin}{}\left(\mathrm{θ}\right)}^{2}& 0\\ 0& 0& 0& \frac{2{}m-r}{r}\end{array}\right]\right)$ (6)

The kdriver equations as explicit differential equations for the Lapse

 > $\mathrm{LapseAndShiftConditions}\left(\mathrm{kdriver}\right)$
 $\left[\begin{array}{c}{\mathrm{%d_}}_{0}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)=\mathrm{%c}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\\ {\mathrm{%d_}}_{0}{}\left(\mathrm{Lapse}\right)=-\mathrm{%epsilon}{}\left({\mathrm{%d_}}_{0}{}\left({\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)+\mathrm{%c}{}{\mathrm{%ExtrinsicCurvature}}_{\mathrm{trace}}\right)\\ 0<\mathrm{%c},0<\mathrm{%epsilon}\end{array}\right]$ (7)
 > $\mathrm{LapseAndShiftConditions}\left(\mathrm{kdriver},\mathrm{evaluate}\right)$
 $\left[\begin{array}{c}\frac{-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}{\mathrm{Shift}}_{\mathrm{~2}}{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},t\right),\mathrm{θ}\right)-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~3}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right),t\right)-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},r\right),t\right)-2{}\left(m-\frac{1}{2}{}r\right){}\left(-\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},t\right),t\right){}r{}\mathrm{Lapse}+\mathrm{diff}{}\left(\mathrm{Lapse},t\right){}\left(2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},t\right){}r-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}-2{}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right){}\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~3}}-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},r\right){}r{}{\mathrm{Shift}}_{\mathrm{~1}}+\mathrm{Lapse}{}\left(\mathrm{cos}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}\right)\right)\right)}{r{}{\mathrm{Lapse}}^{3}{}\left(2{}m-r\right){}\mathrm{sin}{}\left(\mathrm{θ}\right)}=-\frac{2{}c{}\left(\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{θ}\right)+\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~3}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right)+\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~1}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},r\right)-\left(m-\frac{1}{2}{}r\right){}\left(\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},t\right){}r+\mathrm{Lapse}{}\left(\mathrm{cos}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}\right)\right)\right)}{r{}{\mathrm{Lapse}}^{2}{}\left(2{}m-r\right){}\mathrm{sin}{}\left(\mathrm{θ}\right)}\\ \mathrm{diff}{}\left(\mathrm{Lapse},t\right)=-\frac{2{}\left(-\mathrm{sin}{}\left(\mathrm{θ}\right){}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}{\mathrm{Shift}}_{\mathrm{~2}}{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},t\right),\mathrm{θ}\right)-\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~3}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right),t\right)-\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}r{}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},r\right),t\right)+\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}\mathrm{Lapse}{}\left(m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{diff}{}\left(\mathrm{Lapse},t\right),t\right)-2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\left(m-\frac{1}{2}{}r\right){}r{}{\mathrm{diff}{}\left(\mathrm{Lapse},t\right)}^{2}+\left(m-\frac{1}{2}{}r\right){}\left(2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}+2{}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right){}\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~3}}+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},r\right){}r{}{\mathrm{Shift}}_{\mathrm{~1}}+\left(\left(c{}r-2{}{\mathrm{Shift}}_{\mathrm{~1}}\right){}\mathrm{sin}{}\left(\mathrm{θ}\right)-\mathrm{cos}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}\right){}\mathrm{Lapse}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},t\right)+c{}\left(-\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{θ}\right)-\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~3}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right)-\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~1}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},r\right)+\left(m-\frac{1}{2}{}r\right){}\mathrm{Lapse}{}\left(\mathrm{cos}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}\right)\right){}\mathrm{Lapse}\right){}\mathrm{ε}}{r{}{\mathrm{Lapse}}^{3}{}\left(2{}m-r\right){}\mathrm{sin}{}\left(\mathrm{θ}\right)}\\ 0 (8)

The conditions for the harmonicslicing in the context of the Schwarzschild metric

 > $\mathrm{LapseAndShiftConditions}\left(\mathrm{harmonicslicing},\mathrm{evaluate}\right)$
 $\left[\begin{array}{c}0=0\\ \mathrm{diff}{}\left(\mathrm{Lapse},t\right)=\frac{2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{θ}\right)+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~3}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},\mathrm{φ}\right)+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~1}}{}\left(-\frac{1}{2}{}{\mathrm{Lapse}}^{2}{}r+m-\frac{1}{2}{}r\right){}\mathrm{diff}{}\left(\mathrm{Lapse},r\right)-2{}\left(m-\frac{1}{2}{}r\right){}\left(\mathrm{sin}{}\left(\mathrm{θ}\right){}\mathrm{diff}{}\left(\mathrm{Lapse},t\right){}r+\mathrm{Lapse}{}\left(\mathrm{cos}{}\left(\mathrm{θ}\right){}r{}{\mathrm{Shift}}_{\mathrm{~2}}+2{}\mathrm{sin}{}\left(\mathrm{θ}\right){}{\mathrm{Shift}}_{\mathrm{~1}}\right)\right)}{r{}\left(2{}m-r\right){}\mathrm{sin}{}\left(\mathrm{θ}\right)}\\ \mathrm{Lapse}=\frac{C{}\left(r,\mathrm{θ},\mathrm{φ},t\right){}{r}^{2}{}\mathrm{sin}{}\left(\mathrm{θ}\right)}{\mathrm{Lapse}}\\ \left[\begin{array}{ccc}{\mathrm{Shift}}_{\mathrm{~1}}=0& {\mathrm{Shift}}_{\mathrm{~2}}=0& {\mathrm{Shift}}_{\mathrm{~3}}=0\end{array}\right]\end{array}\right]$ (9)

These equations are actually simpler than what they seem to be and can be solved exactly

 > $\mathrm{LapseAndShiftConditions}\left(\mathrm{harmonicslicing},\mathrm{output}=\mathrm{solution}\right)$
 $\left\{\mathbf{\alpha }{=}{-}\sqrt{{\mathrm{sin}}{}\left({\mathrm{\theta }}\right)}{}\sqrt{{\mathrm{_F1}}{}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}\right)}{}{r}{,}{m}{=}{m}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{1}}}^{\phantom{{}}{1}}{=}{0}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{2}}}^{\phantom{{}}{2}}{=}{0}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{3}}}^{\phantom{{}}{3}}{=}{0}{,}{C}{}\left({X}\right){=}{\mathrm{_F1}}{}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}\right)\right\}{,}\left\{\mathbf{\alpha }{=}\sqrt{{\mathrm{sin}}{}\left({\mathrm{\theta }}\right)}{}\sqrt{{\mathrm{_F1}}{}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}\right)}{}{r}{,}{m}{=}{m}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{1}}}^{\phantom{{}}{1}}{=}{0}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{2}}}^{\phantom{{}}{2}}{=}{0}{,}{\mathbf{\beta }}_{\phantom{{}}\phantom{{3}}}^{\phantom{{}}{3}}{=}{0}{,}{C}{}\left({X}\right){=}{\mathrm{_F1}}{}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}\right)\right\}$ (10)
 > See Also References

  Alcubierre, M., Introduction to 3+1 Numerical Relativity, International Series of Monographs on Physics 140, Oxford University Press, 2008.
  Baumgarte, T.W., Shapiro, S.L., Numerical Relativity, Solving Einstein's Equations on a Computer, Cambridge University Press, 2010.
  Gourgoulhon, E., 3+1 Formalism and Bases of Numerical Relativity, Lecture notes, 2007, https://arxiv.org/pdf/gr-qc/0703035v1.pdf. Compatibility

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