Function Repository Resource:

ADMDecomposition

Source Notebook

Represent a canonical decomposition of the metric for a Riemannian or pseudo-Riemannian manifold via the ADM formalism

Contributed by: Jonathan Gorard

ResourceFunction["ADMDecomposition"][MetricTensor[],t,α,βi]

represents an ADM decomposition with the given (initial) spatial MetricTensor, distinguished "time" coordinate t, lapse function α and shift vector (field) βi.

ResourceFunction["ADMDecomposition"]["name",]

represents a named ADM decomposition "name".

ResourceFunction["ADMDecomposition"][{"name",},]

represents a named parameterized ADM decomposition "name", with additional parameter(s) specified within a list.

ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],t]

transforms a specified ResourceFunction["ADMDecomposition"] into one with new distinguished "time" coordinate t.

ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],α,βi]

transforms a specified ResourceFunction["ADMDecomposition"] into one with new lapse function α and new shift vector (field) βi.

ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],t,α,βi]

transforms a specified ResourceFunction["ADMDecomposition"] into one with new distinguished "time" coordinate t, new lapse function α and new shift vector (field) βi.

Details

Formally, the ADM decomposition (due to Arnowitt, Deser and Misner) represents a foliation of a Riemannian or pseudo-Riemannian manifold into a sequence of submanifolds, i.e. it designates an equivalence relation on the manifold whose equivalence classes are injectively-immersed submanifolds of uniform dimension. Typically, these submanifolds are known as the "leaves", "slices" or "hypersurfaces" of the foliation.
Typically, the ADM formalism is applied in the case of general relativity, in which the ambient manifold is a Lorentzian spacetime and the submanifolds are usually spacelike hypersurfaces of codimension-1, with the equivalence classes being the level surfaces of a smooth universal "time function" whose gradient is everywhere timelike. Such a decomposition thus allows one to formulate the Einstein field equations as an explicit time evolution/initial value problem.
The dynamical variables of the ADM formalism are the components of the spatial metric tensor (i.e. the metric tensor field associated to each submanifold/spacelike hypersurface), along with their respective conjugate momenta, otherwise known as the components of the extrinsic curvature tensor (i.e. the rank-2 tensor field characterizing how the submanifolds/spacelike hypersurfaces are embedded into the ambient manifold/spacetime). Together, these variables allow one to define a Hamiltonian, thus permitting one to write the equations of motion for general relativity as a special case of Hamilton's equations, describing how the components of the spatial metric tensor and the extrinsic curvature tensor "evolve" from one submanifold/spacelike hypersurface to the next.
There is, in addition, a family of Lagrange multipliers associated with the ADM formalism, namely the lapse function α and the shift vector (field) βi; these correspond to the degrees of gauge/coordinate freedom inherent to the ADM decomposition. Geometrically, the lapse function α represents the "timelike" coordinate distance between neighboring submanifolds/spacelike hypersurfaces, and therefore determines the "slicing" of the ambient manifold/spacetime, while the shift vector (field) βi represents how the spatial coordinates are relabeled between neighboring submanifolds/spacelike hypersurfaces.
If the manifold/spacetime is globally hyperbolic, then there exists a foliation (and therefore a gauge choice) in which the submanifolds/spacelike hypersurfaces do not intersect, and thus in which the initial value problem formulation of the Einstein field equations is guaranteed to be well-posed.
The Einstein field equations lead to the evolution equations/equations of motion for the spatial metric tensor and the extrinsic curvature tensor, while the (contracted) Bianchi identities lead to the so-called Hamiltonian constraint and momentum constraint equations. The evolution equations are typically hyperbolic in character, while the Hamiltonian and momentum constraint equations are typically elliptic in character (albeit with many notable exceptions, such as harmonic and generalized harmonic gauge, in which the constraint equations are intentionally formulated so as to be hyperbolic). The resulting system can then be closed by specifying appropriate governing equations for the lapse and shift gauge variables.
ResourceFunction["ADMDecomposition"] comes equipped with several in-built slicing conditions (i.e. governing equations for the lapse function α) as well as coordinate conditions (i.e. governing equations for the shift vector field βi), including the geodesic, maximal, harmonic and 1+log slicing conditions for the lapse, and the normal, harmonic, minimal distortion (also known as the generalized Smarr-York conditions) and pseudo-minimal distortion coordinate conditions for the shift.
Furthermore, any valid ADM decomposition must additionally satisfy the Gauss equations (relating projections of the Riemann curvature tensor over the ambient manifold/spacetime to components of the Riemann curvature tensor over the submanifolds/spacelike hypersurfaces) and the Codazzi-Mainardi equations (relating projections of the Ricci curvature tensor over the ambient manifold/spacetime to covariant derivatives of the extrinsic curvature tensor over the submanifolds/spacelike hypersurfaces), which specify the compatibility conditions necessary for the submanifolds/spacelike hypersurfaces in the decomposition to "embed" correctly into the ambient manifold/spacetime.
By default, ResourceFunction["ADMDecomposition"][MetricTensor[],t,α,βi] represents an ADM decomposition by a given (initial) spatial metric tensor, a distinguished "time" coordinate t, a lapse function α and a shift vector (field) βi, although it can be applied to construct decompositions of arbitrary Riemannian and pseudo-Riemannian manifolds that do not necessarily have distinguished "time" and "space" directions in the physical sense. In this way, ResourceFunction["ADMDecomposition"] can be used to set up a Cauchy problem for the Einstein field equations in general relativity.
ResourceFunction["ADMDecomposition"] does not assume any particular number of dimensions, nor any particular metric convention, for the ambient manifold (although some properties, such as "LorentzianConditions", implicitly assume a (-,+,+,+,) signature). ResourceFunction["ADMDecomposition"] also does not require the ambient metric to be strictly symmetric (i.e. spin and torsion connections are also supported), although the computation of the ADM decomposition implicitly assumes the torsion-free Levi-Civita connection by default.
ResourceFunction["ADMDecomposition"] implicitly keeps track of all known algebraic equivalences between tensor expressions, and can apply them for simplification purposes where necessary. Requesting any property with "Reduced" in its name has the effect of applying all known tensor equivalences and simplifying (note that this can have the effect of increasing computation time significantly).
By default, ResourceFunction["ADMDecomposition"] evaluates all partial derivatives of the metric tensor(s) and gauge variables automatically. In certain cases, however, these partial derivatives may be difficult or even impossible to compute, in which case the evaluation may not terminate in a reasonable time. Requesting any property with "Symbolic" in its name has the effect of leaving all partial derivative operators unevaluated instead (note that this can have the effect of increasing expression length significantly).
In ResourceFunction["ADMDecomposition"]["name",] or ResourceFunction["ADMDecomposition"][{"name",},], the following named (in-built) ADM decompositions are supported:
"Minkowski"ADM decomposition of the metric for 4-dimensional (or 1+3-dimensional) flat/Minkowski spacetime in Cartesian coordinates
{"Minkowski",d}ADM decomposition of the metric for d-dimensional (or 1+(d-1)-dimensional) flat/Minkowski spacetime in Cartesian coordinates
"Schwarzschild"ADM decomposition of the metric describing the exterior spacetime geometry surrounding an uncharged, spherically-symmetric, non-rotating mass distribution (e.g. an uncharged, non-rotating black hole) in Schwarzschild/spherical polar coordinates, with purely symbolic mass "M"
{"Schwarzschild",M}ADM decomposition of the metric describing the exterior spacetime geometry surrounding an uncharged, spherically-symmetric, non-rotating mass distribution (e.g. an uncharged, non-rotating black hole) in Schwarzschild/spherical polar coordinates, with numerical mass M
"Kerr"ADM decomposition of the metric describing the exterior spacetime geometry surrounding an uncharged, axially-symmetric, rotating mass distribution (e.g. an uncharged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with purely symbolic mass "M" and purely symbolic angular momentum "J"
{"Kerr",M}ADM decomposition of the metric describing the exterior spacetime geometry surrounding an uncharged, axially-symmetric, rotating mass distribution (e.g. an uncharged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with numerical mass M and purely symbolic angular momentum "J"
{"Kerr",M,J}ADM decomposition of the metric describing the exterior spacetime geometry surrounding an uncharged, axially-symmetric, rotating mass distribution (e.g. an uncharged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with numerical mass M and numerical angular momentum J
"ReissnerNordstrom"ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, spherically-symmetric, non-rotating mass distribution (e.g. a charged, non-rotating black hole) in Schwarzschild/spherical polar coordinates, with purely symbolic (total) mass "M" and purely symbolic electric charge "Q"
{"ReissnerNordstrom",M}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, spherically-symmetric, non-rotating mass distribution (e.g. a charged, non-rotating black hole) in Schwarzschild/spherical polar coordinates, with numerical (total) mass M and purely symbolic electric charge "Q"
{"ReissnerNordstrom",M,Q}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, spherically-symmetric, non-rotating mass distribution (e.g. a charged, non-rotating black hole) in Schwarzschild/spherical polar coordinates, with numerical (total) mass M and numerical electric charge Q
"KerrNewman"ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, axially-symmetric, rotating mass distribution (e.g. a charged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with purely symbolic (total) mass "M", purely symbolic angular momentum "J" and purely symbolic electric charge "Q"
{"KerrNewman",M}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, axially-symmetric, rotating mass distribution (e.g. a charged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with numerical (total) mass M, purely symbolic angular momentum "J" and purely symbolic electric charge "Q"
{"KerrNewman",M,J}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, axially-symmetric, rotating mass distribution (e.g. a charged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with numerical (total) mass M, numerical angular momentum J and purely symbolic electric charge "Q"
{"KerrNewman",M,J,Q}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a charged, axially-symmetric, rotating mass distribution (e.g. a charged, spinning black hole) in Boyer-Lindquist/oblate spheroidal coordinates, with numerical (total) mass M, numerical angular momentum J and numerical electric charge Q
"BrillLindquist"ADM decomposition of the metric describing the exterior spacetime geometry surrounding a pair of uncharged, spherically-symmetric, non-rotating mass distributions (e.g. an uncharged, non-rotating binary black hole system) in Schwarzschild/spherical polar coordinates, with purely symbolic mass "M" and purely symbolic separation distance "z0"
{"BrillLindquist",M}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a pair of uncharged, spherically-symmetric, non-rotating mass distributions (e.g. an uncharged, non-rotating binary black hole system) in Schwarzschild/spherical polar coordinates, with numerical mass M and purely symbolic separation distance "z0"
{"BrillLindquist",M,z0}ADM decomposition of the metric describing the exterior spacetime geometry surrounding a pair of uncharged, spherically-symmetric, non-rotating mass distributions (e.g. an uncharged, non-rotating binary black hole system) in Schwarzschild/spherical polar coordinates, with numerical mass M and numerical separation distance z0
"FLRW"ADM decomposition of the metric describing the spacetime geometry of a homogeneous, isotropic universe which is either contracting or expanding in spherical polar coordinates, with purely symbolic global curvature "k" and purely symbolic scale factor "a" (treated as a function of coordinate time)
{"FLRW",k}ADM decomposition of the metric describing the spacetime geometry of a homogeneous, isotropic universe which is either contracting or expanding in spherical polar coordinates, with numerical global curvature k and purely symbolic scale factor "a" (treated as a function of coordinate time)
{"FLRW",k,a}ADM decomposition of the metric describing the spacetime geometry of a homogeneous, isotropic universe which is either contracting or expanding in spherical polar coordinates, with numerical global curvature k and numerical scale factor a (treated as a function of coordinate time)
Calling ResourceFunction["ADMDecomposition"][] returns a list of all named (built-in) ADM decompositions.
If the function succeeds in performing the specified ADM decomposition, it will return an ResourceFunction["ADMDecomposition"] expression.
Calling ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],t] has the effect of transforming the ADM decomposition (as well as its underlying spatial metric tensor) to use the new distinguished "time" coordinate t. Calling ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],α,βi] has the effect of transforming the ADM decomposition to use the new lapse function α and shift vector (field) βi. Calling ResourceFunction["ADMDecomposition"][ResourceFunction["ADMDecomposition"][],t,α,βi] has the effect of performing both transformations simultaneously.
Based on the eigenvalues of the matrix representation of the ambient metric tensor in covariant form, the ambient manifold/spacetime will be classified by ResourceFunction["ADMDecomposition"] as either Riemannian (all eigenvalues positive or all eigenvalues negative), pseudo-Riemannian (some eigenvalues positive and some eigenvalues negative), Lorentzian (all eigenvalues positive except for one negative, or all eigenvalues negative except for one positive) or Indeterminate.
The property "RiemannianConditions" returns the conditions necessary to guarantee that the eigenvalues of the matrix representation of the ambient metric tensor are strictly positive; "PseudoRiemannianConditions" returns the conditions necessary to guarantee that the eigenvalues of the matrix representation of the ambient metric tensor are all non-zero; "LorentzianConditions" returns the conditions necessary to guarantee that the eigenvalue corresponding to the {1,0,0,} eigenvector of the ambient metric tensor (if it exists) if negative, with all other eigenvalues being positive. Note that this is strictly less general than the behavior of properties such as "RiemannianQ" (described in the point above), since it assumes certain features of the metric signature.
In ResourceFunction["ADMDecomposition"], the following properties are supported:
"SpatialMetricTensor"underlying (spatial) metric tensor on each submanifold/spacelike hypersurface of the ADM decomposition
"SpacetimeMetricTensor"metric tensor on the resulting ambient manifold/spacetime of the ADM decomposition
"NormalVector"(future-pointing, timelike) unit vector normal to each submanifold/spacelike hypersurface of the ADM decomposition
"ReducedNormalVector"(future-pointing, timelike) unit vector normal to each submanifold/spacelike hypersurface of the ADM decomposition, modulo all tensor equivalences
"SymbolicNormalVector"(future-pointing, timelike) unit vector normal to each submanifold/spacelike hypersurface of the ADM decomposition, with purely symbolic partial derivative operators
"TimeVector"(future-pointing, timelike) "time vector" connecting corresponding points on neighboring submanifolds/spacelike hypersurfaces of the ADM decomposition
"SymbolicTimeVector"(future-pointing, timelike) "time vector" connecting corresponding points on neighboring submanifolds/spacelike hypersurfaces of the ADM decomposition, with purely symbolic partial derivative operators
"TimeCoordinate"distinguished "time" coordinate symbol for the ADM decomposition
"SpatialCoordinates"list of distinguished "spatial" coordinate symbols for the ADM decomposition
"CoordinateOneForms"list of differential 1-form symbols for the ambient/spacetime coordinates of the ADM decomposition
"LapseFunction"lapse function determining "timelike" coordinate distance between neighboring submanifolds/spacelike hypersurfaces of the ADM decomposition
"ShiftVector"shift vector (field) determining how "spatial" coordinates are relabeled between neighboring submanifolds/spacelike hypersurfaces of the ADM decomposition
"GaussEquations"list of Gauss equations relating projections of the Riemann curvature tensor over the ambient manifold/spacetime to components of the Riemann curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition
"SymbolicGaussEquations"list of Gauss equations relating projections of the Riemann curvature tensor over the ambient manifold/spacetime to components of the Riemann curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition, with purely symbolic partial derivative operators
"CodazziMainardiEquations"list of Codazzi-Mainardi equations relating projections of the Ricci curvature tensor over the ambient manifold/spacetime to covariant derivatives of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition
"SymbolicCodazziMainardiEquations"list of Codazzi-Mainardi equations relating projections of the Ricci curvature tensor over the ambient manifold/spacetime to covariant derivatives of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition, with purely symbolic partial derivative operators
"GeodesicSlicingCondition"condition on the lapse function required to guarantee geodesic slicing (i.e. unit lapse over the entire ambient manifold/spacetime) for the ADM decomposition
"ReducedGeodesicSlicingCondition"condition on the lapse function required to guarantee geodesic slicing (i.e. unit lapse over the entire ambient manifold/spacetime) for the ADM decomposition, modulo all tensor equivalences
"MaximalSlicingCondition"condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition
"ReducedMaximalSlicingCondition"condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition, modulo all tensor equivalences
"SymbolicMaximalSlicingCondition"condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition, with purely symbolic partial derivative operators
"HarmonicSlicingCondition"condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the "time" coordinate) for the ADM decomposition
"ReducedHarmonicSlicingCondition"condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the "time" coordinate) for the ADM decomposition, modulo all tensor equivalences
"SymbolicHarmonicSlicingCondition"condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the "time" coordinate) for the ADM decomposition, with purely symbolic partial derivative operators
"OnePlusLogSlicingCondition"condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition
"ReducedOnePlusLogSlicingCondition"condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition, modulo all tensor equivalences
"SymbolicOnePlusLogSlicingCondition"condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition, with purely symbolic partial derivative operators
"NormalCoordinateConditions"list of conditions on the shift vector (field) required to guarantee normal coordinates (i.e. vanishing shift over the entire ambient manifold/spacetime) for the ADM decomposition
"ReducedNormalCoordinateConditions"list of conditions on the shift vector (field) required to guarantee normal coordinates (i.e. vanishing shift over the entire ambient manifold/spacetime) for the ADM decomposition, modulo all tensor equivalences
"HarmonicCoordinateConditions"list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the "spatial" coordinates) for the ADM decomposition
"ReducedHarmonicCoordinateConditions"list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the "spatial" coordinates) for the ADM decomposition, modulo all tensor equivalences
"SymbolicHarmonicCoordinateConditions"list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the "spatial" coordinates) for the ADM decomposition, with purely symbolic partial derivative operators
"MinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition
"ReducedMinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition, modulo all tensor equivalences
"SymbolicMinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition, with purely symbolic partial derivative operators
"PseudoMinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition
"ReducedPseudoMinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition, modulo all tensor equivalences
"SymbolicPseudoMinimalDistortionConditions"list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition, with purely symbolic partial derivative operators
"Dimensions"number of dimensions of the ambient manifold/spacetime described by the ADM decomposition
"Signature"list of +1s and -1s designating the signature of the ambient manifold/spacetime described by the ADM decomposition (+1 for each positive eigenvalue of the metric, -1 for each negative eigenvalue of the metric)
"RiemannianQ"whether the ambient manifold/spacetime described by the ADM decomposition is Riemannian (i.e. all eigenvalues of the metric have the same sign)
"PseudoRiemannianQ"whether the ambient manifold/spacetime described by the ADM decomposition is pseudo-Riemannian (i.e. all eigenvalues of the metric are non-zero, but not all have the same sign)
"LorentzianQ"whether the ambient manifold/spacetime described by the ADM decomposition is Lorentzian (i.e. all eigenvalues of the metric have the same sign, except for one eigenvalue which has the opposite sign)
"RiemannianConditions"list of conditions required to guarantee that the ambient manifold/spacetime described by the ADM decomposition is Riemannian (i.e. all eigenvalues of the metric are positive)
"PseudoRiemannianConditions"list of conditions required to guarantee that the ambient manifold/spacetime described by the ADM decomposition is pseudo-Riemannian (i.e. all eigenvalues of the metric are non-zero)
"LorentzianConditions"list of conditions required to guarantee that the ambient manifold/spacetime described by the ADM decomposition is Lorentzian (i.e. the "time" eigenvalue of the metric is negative and all other eigenvalues are positive)
"EvolutionEquations"list of equations characterizing the "time" evolution of components of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition
"SymbolicEvolutionEquations"list of equations characterizing the "time" evolution of components of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition, with purely symbolic partial derivative operators
"HamiltonianConstraint"value of the Hamiltonian constraint (derived from the "timelike" component of the contracted Bianchi identities) of the ADM decomposition
"SymbolicHamiltonianConstraint"value of the Hamiltonian constraint (derived from the "timelike" component of the contracted Bianchi identities) of the ADM decomposition, with purely symbolic partial derivative operators
"MomentumConstraints"list of values of the momentum constraints (derived from the "spacelike" components of the contracted Bianchi identities) of the ADM decomposition
"SymbolicMomentumConstraints"list of values of the momentum constraints (derived from the "spacelike" components of the contracted Bianchi identities) of the ADM decomposition, with purely symbolic partial derivative operators
"HamiltonianConstraintEquation"condition required to guarantee that the Hamiltonian constraint (derived from the "timelike" component of the contracted Bianchi identities) of the ADM decomposition vanishes identically
"SymbolicHamiltonianConstraintEquation"condition required to guarantee that the Hamiltonian constraint (derived from the "timelike" component of the contracted Bianchi identities) of the ADM decomposition vanishes identically, with purely symbolic partial derivative operators
"MomentumConstraintEquations"list of conditions required to guarantee that the momentum constraints (derived from the "spacelike" components of the contracted Bianchi identities) of the ADM decomposition vanish identically
"SymbolicMomentumConstraintEquations"list of conditions required to guarantee that the momentum constraints (derived from the "spacelike" components of the contracted Bianchi identities) of the ADM decomposition vanish identically, with purely symbolic partial derivative operators
"Properties"list of properties

Examples

Basic Examples (4) 

Perform an ADM decomposition for the Schwarzschild metric (e.g. for an uncharged, non-rotating black hole with symbolic mass "M") in standard spherical polar coordinates, using the most general/maximally-unconstrained choice of gauge:

In[1]:=
adm = ResourceFunction["ADMDecomposition"]["Schwarzschild"]
Out[1]=

Show the spatial metric tensor for the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[2]:=
adm["SpatialMetricTensor"]
Out[2]=
In[3]:=
%["MatrixRepresentation"] // MatrixForm
Out[3]=

Show the spacetime metric tensor for the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[4]:=
adm["SpacetimeMetricTensor"]
Out[4]=
In[5]:=
%["MatrixRepresentation"] // MatrixForm
Out[5]=

Show the lapse function for the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[6]:=
adm["LapseFunction"]
Out[6]=

Show the shift vector for the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[7]:=
adm["ShiftVector"]
Out[7]=

Show the future-pointing, timelike unit vector normal to spacelike hypersurfaces in the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[8]:=
adm["NormalVector"]
Out[8]=

Show the condition on the lapse function that must hold in order to achieve geodesic slicing in the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[9]:=
adm["GeodesicSlicingCondition"]
Out[9]=

Show the list of conditions on the shift vector that must hold in order to achieve normal coordinates in the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[10]:=
adm["NormalCoordinateConditions"]
Out[10]=

Show the condition on the lapse function that must hold in order to achieve 1+log slicing in the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric:

In[11]:=
adm["OnePlusLogSlicingCondition"]
Out[11]=

Show the condition on the lapse function that must hold in order to achieve 1+log slicing in the (maximally-unconstrained) ADM decomposition of the Schwarzschild metric, with all algebraic equivalences imposed:

In[12]:=
adm["ReducedOnePlusLogSlicingCondition"]
Out[12]=

Perform the same ADM decomposition for the Schwarzschild metric in spherical polar coordinates, but with numerical mass 1, time coordinate symbol t, spatial coordinate symbols r, a1 and a2, using a partially-constrained choice of gauge (defined in terms of scalar functions a and b):

In[13]:=
adm2 = ResourceFunction["ADMDecomposition"][{"Schwarzschild", 1}, t, {r, a1, a2}, a[t, r, a1, a2], {b[t, r, a1, a2], 0, 0}]
Out[13]=

Show the spacetime metric tensor:

In[14]:=
adm2["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[14]=

Show the list of equations of motion describing the evolution of the extrinsic curvature tensor with respect to coordinate time t for the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[15]:=
adm2["EvolutionEquations"]
Out[15]=

Show that the evolution equations for the (partially-constrained) ADM decomposition of the Schwarzschild metric all hold identically:

In[16]:=
FullSimplify[%]
Out[16]=

Show the list of Gauss equations relating projections of the spacetime Riemann curvature tensor to components of the Riemann curvature tensor over spacelike hypersurfaces for the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[17]:=
adm2["GaussEquations"]
Out[17]=

Show that the Gauss equations for the (partially-constrained) ADM decomposition of the Schwarzschild metric all hold identically:

In[18]:=
FullSimplify[%]
Out[18]=

Show the Codazzi-Mainardi equations relating projections of the spacetime Ricci curvature tensor to covariant derivatives of the extrinsic curvature tensor over spacelike hypersurfaces for the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[19]:=
Short[adm2["CodazziMainardiEquations"], 10]
Out[19]=

Show that the Codazzi-Mainardi equations for the (partially-constrained) ADM decomposition of the Schwarzschild metric all hold identically:

In[20]:=
FullSimplify[%]
Out[20]=

Show the equation requiring that the Hamiltonian constraint (derived from the timelike component of the contracted Bianchi identities) vanishes identically for the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[21]:=
adm2["HamiltonianConstraintEquation"]
Out[21]=

Show that the Hamiltonian constraint for the (partially-constrained) ADM decomposition of the Schwarzschild metric vanishes identically:

In[22]:=
FullSimplify[%]
Out[22]=

Show the list of equations requiring that the momentum constraints (derived from the spacelike components of the contracted Bianchi identities) vanish identically for the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[23]:=
adm2["MomentumConstraintEquations"]
Out[23]=

Show that the momentum constraints for the (partially-constrained) ADM decomposition of the Schwarzschild metric vanish identically:

In[24]:=
FullSimplify[%]
Out[24]=

Show the condition on the lapse function that must hold in order to achieve maximal slicing in the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[25]:=
adm2["MaximalSlicingCondition"]
Out[25]=

Show that condition on the lapse function that must hold in order to achieve maximal slicing in the (partially-constrained) ADM decomposition of the Schwarzschild metric, with all algebraic equivalences imposed:

In[26]:=
adm2["ReducedMaximalSlicingCondition"]
Out[26]=

Show the list of conditions on the shift vector that must hold in order to achieve harmonic coordinates in the (partially-constrained) ADM decomposition of the Schwarzschild metric:

In[27]:=
Short[adm2["HarmonicCoordinateConditions"], 10]
Out[27]=

Show the list of conditions on the shift vector that must hold in order to achieve harmonic coordinates in the (partially-constrained) ADM decomposition of the Schwarzschild metric, with all algebraic equivalences imposed:

In[28]:=
adm2["ReducedHarmonicCoordinateConditions"]
Out[28]=

Perform an ADM decomposition for the Kerr metric (e.g. for an uncharged, spinning black hole with symbolic mass "M" and symbolic angular momentum "J") in Boyer-Lindquist/oblate spheroidal coordinates, using a partially-constrained choice of gauge (defined in terms of scalar functions a and b):

In[29]:=
adm = ResourceFunction["ADMDecomposition"]["Kerr", t, {r, a1, a2}, a[t, r, a1, a2], {0, 0, b[t, r, a1, a2]}]
Out[29]=
In[30]:=
adm["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[30]=

Show the list of evolution equations for the (partially-constrained) ADM decomposition of the Kerr metric:

In[31]:=
adm["EvolutionEquations"]
Out[31]=

Show that the evolution equations for the (partially-constrained) ADM decomposition of the Kerr metric all hold identically:

In[32]:=
FullSimplify[%]
Out[32]=

Compute the Hamiltonian constraint (derived from the timelike component of the contracted Bianchi identities) for the (partially-constrained) ADM decomposition of the Kerr metric:

In[33]:=
adm["HamiltonianConstraint"]
Out[33]=

Show that the Hamiltonian constraint for the (partially-constrained) ADM decomposition of the Kerr metric vanishes identically:

In[34]:=
FullSimplify[%]
Out[34]=

Compute the list of momentum constraints (derived from the spacelike components of the contracted Bianchi identities) for the (partially-constrained) ADM decomposition of the Kerr metric:

In[35]:=
adm["MomentumConstraints"]
Out[35]=

Show that the momentum constraints for the (partially-constrained) ADM decomposition of the Kerr metric vanish identically:

In[36]:=
FullSimplify[%]
Out[36]=

Show the condition on the lapse function that must hold in order to achieve harmonic slicing in the (partially-constrained) ADM decomposition of the Kerr metric:

In[37]:=
Short[adm["HarmonicSlicingCondition"], 10]
Out[37]=

Show the condition on the lapse function that must hold in order to achieve harmonic slicing in the (partially-constrained) ADM decomposition of the Kerr metric, with all algebraic equivalences imposed:

In[38]:=
adm["ReducedHarmonicSlicingCondition"]
Out[38]=

Show the list of conditions on the shift vector that must hold in order to achieve minimal distortion coordinates in the (partially-constrained) ADM decomposition of the Kerr metric:

In[39]:=
Short[adm["MinimalDistortionConditions"], 10]
Out[39]=

Show the list of conditions on the shift vector that must hold in order to achieve minimal distortion coordinates in the (partially-constrained) ADM decomposition of the Kerr metric, with all algebraic equivalences imposed:

In[40]:=
adm["ReducedMinimalDistortionConditions"]
Out[40]=

Show the list of conditions on the shift vector that must hold in order to achieve pseudo-minimal distortion coordinates in the (partially-constrained) ADM decomposition of the Kerr metric:

In[41]:=
Short[adm["PseudoMinimalDistortionConditions"], 10]
Out[41]=

Show the list of conditions on the shift vector that must hold in order to achieve pseudo-minimal distortion coordinates in the (partially-constrained) ADM decomposition of the Kerr metric, with all algebraic equivalences imposed:

In[42]:=
adm["ReducedPseudoMinimalDistortionConditions"]
Out[42]=

Perform an ADM decomposition from an (initial) spatial metric tensor directly, using the most general/maximally-unconstrained choice of gauge:

In[43]:=
metric = ResourceFunction["MetricTensor"]["Euclidean"]
Out[43]=
In[44]:=
adm = ResourceFunction["ADMDecomposition"][metric]
Out[44]=
In[45]:=
adm["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[45]=

Show the list of gauge conditions that must hold for the resulting ambient/spacetime metric to be Riemannian:

In[46]:=
adm["RiemannianConditions"]
Out[46]=

Show the list of gauge conditions that must hold for the resulting ambient/spacetime metric to be pseudo-Riemannian:

In[47]:=
adm["PseudoRiemannianConditions"]
Out[47]=

Perform an ADM decomposition for the Brill-Lindquist metric (e.g. for an uncharged, non-rotating binary black hole system with symbolic mass "M" and initial separation distance "z0") in standard spherical polar coordinates, using a partially constrained choice of gauge (defined in terms of a single scalar function a):

In[48]:=
adm = ResourceFunction["ADMDecomposition"]["BrillLindquist", t, {r, a1, a2}, a[t, r, a1, a2], {0, 0, 0}]
Out[48]=

Show the list of evolution equations for the (partially-constrained) ADM decomposition of the Brill-Lindquist metric:

In[49]:=
adm["EvolutionEquations"]
Out[49]=

Show that the evolution equations for the (partially-constrained) ADM decomposition of the Brill-Lindquist metric all hold identically:

In[50]:=
FullSimplify[%]
Out[50]=

Show the Hamiltonian constraint equation for the (partially-constrained) ADM decomposition of the Brill-Lindquist metric:

In[51]:=
adm["HamiltonianConstraintEquation"]
Out[51]=

Show that the Hamiltonian constraint equation for the (partially-constrained) ADM decomposition of the Brill-Lindquist metric holds identically:

In[52]:=
FullSimplify[%]
Out[52]=

Scope (5) 

ADM decompositions can be performed directly from an (initial) spatial MetricTensor expression:

In[53]:=
adm = ResourceFunction["ADMDecomposition"][
  ResourceFunction["MetricTensor"]["Euclidean"]]
Out[53]=
In[54]:=
adm["SpacetimeMetricTensor"]["LineElement"]
Out[54]=

Additional arguments can be used to specify the distinguished "time" coordinate symbol (otherwise the default symbol "t" will be chosen automatically):

In[55]:=
adm2 = ResourceFunction["ADMDecomposition"][
  ResourceFunction["MetricTensor"]["Euclidean"], T]
Out[55]=
In[56]:=
adm2["SpacetimeMetricTensor"]["LineElement"]
Out[56]=

Or the lapse/shift gauge conditions:

In[57]:=
adm3 = ResourceFunction["ADMDecomposition"][
  ResourceFunction["MetricTensor"]["Euclidean"], 1, {0, 0, 0}]
Out[57]=
In[58]:=
adm3["SpacetimeMetricTensor"]["LineElement"]
Out[58]=

Or both simultaneously:

In[59]:=
adm4 = ResourceFunction["ADMDecomposition"][
  ResourceFunction["MetricTensor"]["Euclidean"], T, 1, {0, 0, 0}]
Out[59]=
In[60]:=
adm4["SpacetimeMetricTensor"]["LineElement"]
Out[60]=

Common ADM decompositions can also be performed using an in-built name:

In[61]:=
adm = ResourceFunction["ADMDecomposition"]["ReissnerNordstrom"]
Out[61]=
In[62]:=
adm["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[62]=

When an in-built/named ADM decomposition has one or more parameters, those parameters can be left unspecified (in which case they are filled with purely symbolic defaults, such as "M" in the above), or can be specified explicitly in list form:

In[63]:=
adm2 = ResourceFunction[
  "ADMDecomposition"][{"ReissnerNordstrom", 1, 1/3}]
Out[63]=
In[64]:=
adm2["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[64]=

If only some parameters are explicitly specified, then the remainder are filled with symbolic defaults (e.g. if one specifies only a numerical mass for the Reissner-Nordström metric, then ADMDecomposition will use a purely symbolic electric charge, namely "Q"):

In[65]:=
adm3 = ResourceFunction["ADMDecomposition"][{"ReissnerNordstrom", 1}]
Out[65]=
In[66]:=
adm3["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[66]=

Coordinate and lapse/shift gauge information can also be specified for in-built/named ADM decompositions:

In[67]:=
adm4 = ResourceFunction["ADMDecomposition"][{"ReissnerNordstrom", 1}, t, {r, a1, a2}, a[t, r, a1, a2], {b[t, r, a1, a2], 0, 0}]
Out[67]=
In[68]:=
adm4["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[68]=

A new distinguished "time" coordinate symbol can be specified for any ADM decomposition:

In[69]:=
adm = ResourceFunction["ADMDecomposition"]["FLRW"]
Out[69]=
In[70]:=
adm["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[70]=
In[71]:=
adm2 = ResourceFunction["ADMDecomposition"][adm, T]
Out[71]=
In[72]:=
adm2["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[72]=

New lapse/shift gauge conditions can also be specified for any ADM decomposition:

In[73]:=
adm3 = ResourceFunction["ADMDecomposition"][adm, 1, {0, 0, 0}]
Out[73]=
In[74]:=
adm3["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[74]=

New coordinate symbols and new lapse/shift gauge conditions can also be specified simultaneously:

In[75]:=
adm4 = ResourceFunction["ADMDecomposition"][adm, T, 1, {0, 0, 0}]
Out[75]=
In[76]:=
adm4["SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[76]=

Show the list of all in-built/named ADM decompositions:

In[77]:=
ResourceFunction["ADMDecomposition"][]
Out[77]=

Perform an ADM decomposition for the 4-dimensional (i.e. 1+3-dimensional) Minkowski metric (default) in Cartesian coordinates:

In[78]:=
ResourceFunction["ADMDecomposition"]["Minkowski"][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[78]=

Perform an ADM decomposition for the 7-dimensional (i.e. 1+6-dimensional) Minkowski metric in Cartesian coordinates:

In[79]:=
ResourceFunction["ADMDecomposition"][{"Minkowski", 7}][
  "SpacetimeMetricTensor"]["MatrixRepresentation"]
Out[79]=

Perform an ADM decomposition for the Schwarzschild metric, e.g. the exterior spacetime of an uncharged, non-rotating black hole, with symbolic mass "M" (default) in spherical polar coordinates:

In[80]:=
ResourceFunction["ADMDecomposition"]["Schwarzschild"][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[80]=

Perform an ADM decomposition for the Schwarzschild metric with numerical mass 1 in spherical polar coordinates:

In[81]:=
ResourceFunction["ADMDecomposition"][{"Schwarzschild", 1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[81]=

Perform an ADM decomposition for the Kerr metric, e.g. the exterior spacetime of an uncharged, spinning black hole, with symbolic mass "M" and symbolic angular momentum "J" (default) in Boyer-Lindquist/oblate spheroidal coordinates:

In[82]:=
ResourceFunction["ADMDecomposition"]["Kerr"]["SpacetimeMetricTensor"][
  "MatrixRepresentation"] // MatrixForm
Out[82]=

Perform an ADM decomposition for the Kerr metric with numerical mass 1 and symbolic angular momentum "J" in Boyer-Lindquist/oblate spheroidal coordinates:

In[83]:=
ResourceFunction["ADMDecomposition"][{"Kerr", 1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[83]=

Perform an ADM decomposition for the Kerr metric with numerical mass 1 and numerical angular momentum 1/2 in Boyer-Lindquist/oblate spheroidal coordinates:

In[84]:=
ResourceFunction["ADMDecomposition"][{"Kerr", 1, 1/2}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[84]=

Perform an ADM decomposition for the Reissner-Nordström metric, e.g. the exterior spacetime of a charged, non-rotating black hole, with symbolic mass "M" and symbolic electric charge "Q" (default) in spherical polar coordinates:

In[85]:=
ResourceFunction["ADMDecomposition"]["ReissnerNordstrom"][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[85]=

Perform an ADM decomposition for the Reissner-Nordström metric with numerical mass 1 and symbolic electric charge "Q" in spherical polar coordinates:

In[86]:=
ResourceFunction["ADMDecomposition"][{"ReissnerNordstrom", 1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[86]=

Perform an ADM decomposition for the Reissner-Nordström metric with numerical mass 1 and numerical electric charge 1/3 in spherical polar coordinates:

In[87]:=
ResourceFunction["ADMDecomposition"][{"ReissnerNordstrom", 1, 1/3}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[87]=

Perform an ADM decomposition for the Kerr-Newman metric, e.g. the exterior spacetime of a charged, spinning black hole, with symbolic mass "M", symbolic angular momentum "J" and symbolic electric charge "Q" (default) in Boyer-Lindquist/oblate spheroidal coordinates:

In[88]:=
ResourceFunction["ADMDecomposition"]["KerrNewman"][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[88]=

Perform an ADM decomposition for the Kerr-Newman metric with numerical mass 1, symbolic angular momentum "J" and symbolic electric charge "Q" in Boyer-Lindquist/oblate spheroidal coordinates:

In[89]:=
ResourceFunction["ADMDecomposition"][{"KerrNewman", 1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[89]=

Perform an ADM decomposition for the Kerr-Newman metric with numerical mass 1, numerical angular momentum 1/2 and symbolic electric charge "Q" in Boyer-Lindquist/oblate spheroidal coordinates:

In[90]:=
ResourceFunction["ADMDecomposition"][{"KerrNewman", 1, 1/2}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[90]=

Perform an ADM decomposition for the Kerr-Newman metric with numerical mass 1, numerical angular momentum 1/2 and numerical electric charge 1/3 in Boyer-Lindquist/oblate spheroidal coordinates:

In[91]:=
ResourceFunction["ADMDecomposition"][{"KerrNewman", 1, 1/2, 1/3}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[91]=

Perform an ADM decomposition for the Brill-Lindquist metric, e.g. the exterior spacetime of a pair of uncharged, non-rotating black holes, with symbolic mass "M" and symbolic separation distance "z0" (default) in spherical polar coordinates:

In[92]:=
ResourceFunction["ADMDecomposition"]["BrillLindquist"][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[92]=

Perform an ADM decomposition for the Brill-Lindquist metric with numerical mass 1 and symbolic separation distance "z0" in spherical polar coordinates:

In[93]:=
ResourceFunction["ADMDecomposition"][{"BrillLindquist", 1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[93]=

Perform an ADM decomposition for the Brill-Lindquist metric with numerical mass 1 and numerical separation distance 2 in spherical polar coordinates:

In[94]:=
ResourceFunction["ADMDecomposition"][{"BrillLindquist", 1, 2}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[94]=

Perform an ADM decomposition for the Friedmann-Lemaître-Robertson-Walker/FLRW metric, i.e. the metric for a homogeneous, isotropic and uniformly expanding/contracting universe, with symbolic global curvature "k" and symbolic scale function "a" (default) in spherical polar coordinates:

In[95]:=
ResourceFunction["ADMDecomposition"]["FLRW"]["SpacetimeMetricTensor"][
  "MatrixRepresentation"] // MatrixForm
Out[95]=

Perform an ADM decomposition for the FLRW metric with numerical global curvature -1 and symbolic scale function "a" in spherical polar coordinates:

In[96]:=
ResourceFunction["ADMDecomposition"][{"FLRW", -1}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[96]=

Perform an ADM decomposition for the FLRW metric with numerical global curvature -1 and numerical scale function (#*3)& in spherical polar coordinates:

In[97]:=
ResourceFunction["ADMDecomposition"][{"FLRW", -1, (#*3) &}][
   "SpacetimeMetricTensor"]["MatrixRepresentation"] // MatrixForm
Out[97]=

Perform an ADM decomposition for the Schwarzschild metric, with symbolic mass "M":

In[98]:=
adm = ResourceFunction["ADMDecomposition"]["Schwarzschild"]
Out[98]=

Show the list of properties:

In[99]:=
adm["Properties"]
Out[99]=

Show the spatial metric tensor (i.e. the metric tensor on submanifolds/spacelike hypersurfaces) for the ADM decomposition:

In[100]:=
adm["SpatialMetricTensor"]
Out[100]=

Show the spacetime metric tensor (i.e. the metric tensor on the ambient manifold/spacetime) for the ADM decomposition:

In[101]:=
adm["SpacetimeMetricTensor"]
Out[101]=

Show the future-pointing, timelike unit vector normal to submanifolds/spacelike hypersurfaces for the ADM decomposition:

In[102]:=
adm["NormalVector"]
Out[102]=

Show the future-pointing, timelike unit vector normal to submanifolds/spacelike hypersurfaces for the ADM decomposition, with all algebraic equivalences imposed:

In[103]:=
adm["ReducedNormalVector"]
Out[103]=

Show the future-pointing, timelike unit vector normal to submanifolds/spacelike hypersurfaces for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[104]:=
adm["SymbolicNormalVector"]
Out[104]=

Show the future-pointing, timelike "time vector" for the ADM decomposition:

In[105]:=
adm["TimeVector"]
Out[105]=

Show the future-pointing, timelike "time vector" for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[106]:=
Short[adm["SymbolicTimeVector"], 10]
Out[106]=

Show the distinguished time coordinate symbol for the ADM decomposition:

In[107]:=
adm["TimeCoordinate"]
Out[107]=

Show the list of distinguished spatial coordinate symbols for the ADM decomposition:

In[108]:=
adm["SpatialCoordinates"]
Out[108]=

Show the list of differential 1-form symbols for the (ambient/spacetime) coordinates of the ADM decomposition:

In[109]:=
adm["CoordinateOneForms"]
Out[109]=

Show the lapse function for the ADM decomposition:

In[110]:=
adm["LapseFunction"]
Out[110]=

Show the shift vector (field) for the ADM decomposition:

In[111]:=
adm["ShiftVector"]
Out[111]=

Show the list of Gauss equations relating projections of the Riemann curvature tensor over the ambient manifold/spacetime to components of the Riemann curvature tensor over the submanifolds/spacelike hypersurfaces for the ADM decomposition:

In[112]:=
adm["GaussEquations"]
Out[112]=

Show the list of Gauss equations relating projections of the Riemann curvature tensor over the ambient manifold/spacetime to components of the Riemann curvature tensor over the submanifolds/spacelike hypersurfaces for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[113]:=
adm["SymbolicGaussEquations"]
Out[113]=

Show the list of Codazzi-Mainardi equations relating projections of the Ricci curvature tensor over the ambient manifold/spacetime to covariant derivatives of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition:

In[114]:=
adm["CodazziMainardiEquations"]
Out[114]=

Show the list of Codazzi-Mainardi equations relating projections of the Ricci curvature tensor over the ambient manifold/spacetime to covariant derivatives of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces of the ADM decomposition, with all partial derivative operators left purely symbolic:

In[115]:=
adm["SymbolicCodazziMainardiEquations"]
Out[115]=

Show the condition on the lapse function required to guarantee geodesic slicing (i.e. unit lapse over the entire ambient manifold/spacetime) for the ADM decomposition:

In[116]:=
adm["GeodesicSlicingCondition"]
Out[116]=

Show the condition on the lapse function required to guarantee geodesic slicing (i.e. unit lapse over the entire ambient manifold/spacetime) for the ADM decomposition, with all algebraic equivalences imposed:

In[117]:=
adm["ReducedGeodesicSlicingCondition"]
Out[117]=

Show the condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition:

In[118]:=
Short[adm["MaximalSlicingCondition"], 10]
Out[118]=

Show the condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition, with all algebraic equivalences imposed:

In[119]:=
adm["ReducedMaximalSlicingCondition"]
Out[119]=

Show the condition on the lapse function required to guarantee maximal slicing (i.e. maximum spatial volume of each submanifold/spacelike hypersurface) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[120]:=
Short[adm["SymbolicMaximalSlicingCondition"], 10]
Out[120]=

Show the condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the time coordinate) for the ADM decomposition:

In[121]:=
Short[adm["HarmonicSlicingCondition"], 10]
Out[121]=

Show the condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the time coordinate) for the ADM decomposition, with all algebraic equivalences imposed:

In[122]:=
adm["ReducedHarmonicSlicingCondition"]
Out[122]=

Show the condition on the lapse function required to guarantee harmonic slicing (i.e. vanishing of the curved spacetime d'Alembertian of the time coordinate) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[123]:=
adm["SymbolicHarmonicSlicingCondition"]
Out[123]=

Show the condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition:

In[124]:=
adm["OnePlusLogSlicingCondition"]
Out[124]=

Show the condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition, with all algebraic equivalences imposed:

In[125]:=
adm["ReducedOnePlusLogSlicingCondition"]
Out[125]=

Show the condition on the lapse function required to guarantee 1+log slicing (i.e. generalized harmonic slicing) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[126]:=
adm["SymbolicOnePlusLogSlicingCondition"]
Out[126]=

Show the list of conditions on the shift vector (field) required to guarantee normal coordinates (i.e. vanishing shift over the entire ambient manifold/spacetime) for the ADM decomposition:

In[127]:=
adm["NormalCoordinateConditions"]
Out[127]=

Show the list of conditions on the shift vector (field) required to guarantee normal coordinates (i.e. vanishing shift over the entire ambient manifold/spacetime) for the ADM decomposition, with all algebraic equivalences imposed:

In[128]:=
adm["ReducedNormalCoordinateConditions"]
Out[128]=

Show the list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the spatial coordinates) for the ADM decomposition:

In[129]:=
Short[adm["HarmonicCoordinateConditions"], 10]
Out[129]=

Show the list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the spatial coordinates) for the ADM decomposition, with all algebraic equivalences imposed:

In[130]:=
adm["ReducedHarmonicCoordinateConditions"]
Out[130]=

Show the list of conditions on the shift vector (field) required to guarantee harmonic coordinates (i.e. vanishing of the curved spacetime d'Alembertian of the spatial coordinates) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[131]:=
adm["SymbolicHarmonicCoordinateConditions"]
Out[131]=

Show the list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition:

In[132]:=
Short[adm["MinimalDistortionConditions"], 10]
Out[132]=

Show the list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition, with all algebraic equivalences imposed:

In[133]:=
adm["ReducedMinimalDistortionConditions"]
Out[133]=

Show the list of conditions on the shift vector (field) required to guarantee minimal distortion coordinates (i.e. minimum strain on each submanifold/spacelike hypersurface) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[134]:=
adm["SymbolicMinimalDistortionConditions"]
Out[134]=

Show the list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition:

In[135]:=
Short[adm["PseudoMinimalDistortionConditions"], 10]
Out[135]=

Show the list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition, with all algebraic equivalences imposed:

In[136]:=
adm["ReducedPseudoMinimalDistortionConditions"]
Out[136]=

Show the list of conditions on the shift vector (field) required to guarantee pseudo-minimal distortion coordinates (i.e. minimal distortion coordinates, but with all covariant derivatives replaced with partial derivatives) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[137]:=
Short[adm["SymbolicPseudoMinimalDistortionConditions"], 10]
Out[137]=

Show the number of dimensions of the ambient manifold/spacetime represented by the ADM decomposition:

In[138]:=
adm["Dimensions"]
Out[138]=

Show the signature of the ambient manifold/spacetime represented by the ADM decomposition (with +1s representing positive eigenvalues and -1s representing negative eigenvalues of the metric tensor):

In[139]:=
adm["Signature"]
Out[139]=

Determine whether the ambient manifold/spacetime represented by the ADM decomposition is Riemannian (i.e. all eigenvalues of the metric tensor have the same sign):

In[140]:=
adm["RiemannianQ"]
Out[140]=

Determine whether the ambient manifold/spacetime represented by the ADM decomposition is pseudo-Riemannian (i.e. all eigenvalues are non-zero, but not all have the same sign):

In[141]:=
adm["PseudoRiemannianQ"]
Out[141]=

Determine whether the ambient manifold/spacetime represented by the ADM decomposition is Lorentzian (i.e. all eigenvalues of the metric tensor have the same sign, except for one eigenvalue which has the opposite sign):

In[142]:=
adm["LorentzianQ"]
Out[142]=

Show the list of conditions on the coordinates required to guarantee that the ambient manifold/spacetime represented by the ADM decomposition is Riemannian (i.e. all eigenvalues of the metric tensor are positive):

In[143]:=
Short[adm["RiemannianConditions"], 20]
Out[143]=

Show the list of conditions on the coordinates required to guarantee that the ambient manifold/spacetime represented by the ADM decomposition is pseudo-Riemannian (i.e. all eigenvalues of the metric tensor are non-zero):

In[144]:=
Short[adm["PseudoRiemannianConditions"], 20]
Out[144]=

Show the list of conditions on the coordinates required to guarantee that the ambient manifold/spacetime represented by the ADM decomposition is Lorentzian (i.e. the "time" eigenvalue is negative, and all other eigenvalues are positive):

In[145]:=
adm["LorentzianConditions"]
Out[145]=

Show the list of equations characterizing the time evolution of components of the extrinsic curvature tensor over submanifold/spacelike hypersurfaces for the ADM decomposition:

In[146]:=
adm["EvolutionEquations"]
Out[146]=

Show the list of equations characterizing the time evolution of components of the extrinsic curvature tensor over submanifolds/spacelike hypersurfaces for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[147]:=
adm["SymbolicEvolutionEquations"]
Out[147]=

Show the value of the Hamiltonian constraint (i.e. the timelike component of the contracted Bianchi identities) for the ADM decomposition:

In[148]:=
adm["HamiltonianConstraint"]
Out[148]=

Show the value of the Hamiltonian constraint (i.e. the timelike component of the contracted Bianchi identities) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[149]:=
adm["SymbolicHamiltonianConstraint"]
Out[149]=

Show the list of values of the momentum constraints (i.e. the spacelike components of the contracted Bianchi identities) for the ADM decomposition:

In[150]:=
adm["MomentumConstraints"]
Out[150]=

Show the list of values of the momentum constraints (i.e. the spacelike components of the contracted Bianchi identities) for the ADM decomposition, with all partial derivative operators left purely symbolic:

In[151]:=
adm["SymbolicMomentumConstraints"]
Out[151]=

Show the condition required to guarantee that the Hamiltonian constraint (i.e. the timelike component of the contracted Bianchi identities) for the ADM decomposition vanishes identically:

In[152]:=
adm["HamiltonianConstraintEquation"]
Out[152]=

Show the condition required to guarantee that the Hamiltonian constraint (i.e. the timelike component of the contracted Bianchi identities) for the ADM decomposition vanishes identically, with all partial derivative operators left purely symbolic:

In[153]:=
adm["SymbolicHamiltonianConstraintEquation"]
Out[153]=

Show the list of conditions required to guarantee that the momentum constraints (i.e. the spacelike components of the contracted Bianchi identities) for the ADM decomposition vanish identically:

In[154]:=
adm["MomentumConstraintEquations"]
Out[154]=

Show the list of conditions required to guarantee that the momentum constraints (i.e. the spacelike components of the contracted Bianchi identities) for the ADM decomposition vanish identically, with all partial derivative operators left purely symbolic:

In[155]:=
adm["SymbolicMomentumConstraintEquations"]
Out[155]=

Publisher

Jonathan Gorard

Version History

  • 1.0.0 – 22 March 2023

Source Metadata

Related Resources

License Information