Physics[ThreePlusOne][D3_] - indexed covariant differential operator with respect to curvilinear space coordinates acting on a 3D hypersurface
Calling Sequence
D3_[mu](A)
D3_[mu](A, [X])
D3_(A)
D3_(A, [X])
Parameters
mu
-
a name representing an integer number between 1 and the spacetime dimension or a space index running from 1 to the dimension - 1, can also be the number itself or preceded by ~ to represent a contravariant index
A
any mathematical expression or relation between expressions, or a set or list of them, or an rtable
[X]
a list of differentiation variables to which the index j refers
Description
As is the case of all the tensors of the ThreePlusOne package, D3_[mu] is a 4D spacetime tensor with special meaning (acting) on a 3D hypersurface specified by the values of the Lapse and Shift. The D3_[mu] command is a computational representation for the covariant differential operator in that hypersurface. Its space components also form a tensor in the 3D hypersurface, so D3_ can be used with 3D space indices, for instances as in D_[i](A[j]).
On the screen, D3_, is displayed as the corresponding D_ command, but in black ▿__μ, instead of in blue ▿__μ.
Note that the space components of - say - D3_[mu](A[nu](x, y, z, t)) are not necessarily the same as the space components of the 4D tensor D_[mu](A[nu](x, y, z, t)). Instead, the relationship between these two is given by
▿μ⁡Aν=γμα⁢γνβ⁢▿α⁡Aβ
where γμα is the gamma3_, a 4D tensor, the spatial metric induced in the 3D hypersurface by the 4D g_ metric, also a projection operator that projects 4D tensors into the 3D hypersurface, resulting in purely spatial 4D tensors, all of whose components are equal to zero when their indices are contravariant and any of them has a timelike value (i.e. contravariant 0).
D3_ can be used as well without an index, as in D3_(A) displayed as ▿⁡A, in which case it represents the total differential in curvilinear coordinates, and the output comes automatically expanded as D3_[mu](A) * D3_(X[~mu]), where X[~mu] is the corresponding spacetime vector.
The %D3_ command is the inert form of D3_, so it represents the same mathematical operation but without performing it. To perform the operation, use value.
Examples
with⁡Physics:with⁡ThreePlusOne
⁢_______________________________________________________
Setting lowercaselatin_is letters to represent space indices
Defined as 4D spacetime tensors see ?Physics,ThreePlusOne,γμ,ν,▿μ,Γμ,ν,α,Rμ,ν,Rμ,ν,α,β,βμ,nμ,tμ,Κμ,ν
Changing the signature of spacetime from - - - + to + + + - in order to match the signature customarily used in the ADM formalism
⁢Systems of spacetime coordinates are:⁢X=x,y,z,t
ADMEquations,Christoffel3,D3_,ExtrinsicCurvature,Lapse,LapseAndShiftConditions,Ricci3,Riemann3,Shift,TimeVector,UnitNormalVector,gamma3_
Setup⁡mathematicalnotation=true
mathematicalnotation=true
Set up a coordinate system to work with - the first one to be set is automatically taken as the differentiation variables for D_ and so also for D3_
Setup⁡coordinatesystems=cartesian
coordinatesystems=X
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 g_ with no indices
g_
gμ,ν=100001000010000−1
Check the nonzero components of Christoffel: since this is a Minkowski spacetime, there are none, then ▿__j=∂__j
Christoffelnonzero
Γα,μ,ν=∅
The nonzero components of Christoffel3 are also none
Christoffel3nonzero
Note the difference in color: the 4D tensors are displayed in blue while the corresponding tensors of ThreePlusOne are displayed in black.
Because the Christoffel symbols are all equal to zero, the covariant derivative is equal to the standard derivative expressed using the d_ operator, both for D_ and D3_
D3_⁡X~nu
ⅆ⁡X⁢ν⁢ν
Set the scenario as a Schwarzschild spacetime in spherical coordinates; you can do this entering Setup(metric = Schwarzschild) or in the simpler form taking advantage of abbreviations and directly using the spacetime metric g_ command
g_sc
⁢Systems of spacetime coordinates are:⁢X=r,θ,φ,t
⁢Default differentiation variables for d_, D_ and dAlembertian are:⁢X=r,θ,φ,t
⁢The Schwarzschild metric in coordinates ⁢r,θ,φ,t
Parameters: m
Signature: + + + -
gμ,ν=−r2⁢m−r0000r20000r2⁢sin⁡θ200002⁢m−rr
The covariant derivative of a scalar is always equal to the ∂ derivative
D3_j⁡Φ⁡X
γjμjμ⁢∂μ⁡Φ⁡X
The covariant differential of a scalar function
D3_⁡Φ⁡X
∂μ⁡Φ⁡X⁢ⅆ⁡X⁢μ⁢μ
For illustration purposes Define an arbitrary 3D space tensor A; to indicate that this is a 3D tensor, define it with its space indices explicit, for example
Define⁡Aj
Defined objects with tensor properties
Aj,▿μ,▿μ,γμ,σμ,Rμ,ν,Rμ,ν,Rμ,ν,α,β,βμ,Cμ,ν,α,β,∂μ,gμ,ν,γμ,ν,γi,j,Γμ,ν,α,Γμ,ν,α,Gμ,ν,εα,β,μ,ν,Rμ,ν,α,β,tμ,Κμ,ν,nμ,Xμ
Use a compact display for this function to avoid redundant repeated display of its functionality plus having derivatives displayed indexed by the differentiation variables
CompactDisplay⁡A⁡X
A⁡r,θ,φ,t⁢will now be displayed as⁢A
Now we have, for the total 3D covariant differential and the 3D covariant derivative, respectively,
CompactDisplay⁡Aj⁡X
D3_⁡A~j⁡X
▿μ⁡A⁢j⁢j⁢ⅆ⁡X⁢μ⁢μ
D3_j⁡A~j⁡X
▿j⁡A⁢j⁢j
This D3_ covariant derivative is expressed in terms of the Christoffel3 symbols in the same way the covariant derivative D_ is expressed in terms of the 4D Christoffel. To see that you can use expand or convert to d_
convert⁡,d_
γjαjα⁢γβjβj⁢∂α⁡A⁢β⁢β+Γ⁢βα,i⁢βα,i⁢A⁢i⁢i
The generalized divergence of a vector in curvilinear coordinates
expand⁡
TensorArray⁡
A⁢1⁢1r⁢2⁢m−r⁢r+m⁢A⁢1⁢1+A⁢2⁢2θ⁢2⁢m−r⁢r+2⁢A⁢1⁢1⁢2⁢m−r+A⁢3⁢3φ⁢2⁢m−r⁢r+cot⁡θ⁢A⁢2⁢2⁢2⁢m−r⁢r2⁢m−r⁢r
To compute with a representation for D3_ without actually performing the operation, use the inert form %D3_. To afterwards perform the operation use value. For example, when the space part of the 3+1 metric gamma3_ is diagonal, its covariant derivative is equal to 0
%D3_i⁡gamma3_j,k
▿i⁡γj,k
value⁡
0
The general form of such a 3+1 metric happens when the space part of the 4D metric g_ is diagonal. For example, using arbitrary functions of X for the other components of the 4D metric, we have
Setup⁡metric=1,1=f__1⁡X,1,4=f__14⁡X,2,2=f__2⁡X,2,4=f__24⁡X,3,3=f__3⁡X,3,4=f__34⁡X,4,4=f__4⁡X
Coordinates: r,θ,φ,t. Signature: + + + -
gμ,ν=f__1⁡X00f__14⁡X0f__2⁡X0f__24⁡X00f__3⁡Xf__34⁡Xf__14⁡Xf__24⁡Xf__34⁡Xf__4⁡X
metric=1,1=f__1⁡X,1,4=f__14⁡X,2,2=f__2⁡X,2,4=f__24⁡X,3,3=f__3⁡X,3,4=f__34⁡X,4,4=f__4⁡X
D3_μ⁡gamma3_~mu,~nu
Set the spacetime metric g_ by giving the square of the spacetime interval such that g1,2≠0
ds2≔x2⁢dx2+y2⁢dy2+z2⁢dz2+x⁢y⁢dx⁢dy−dt2
ds2≔x2⁢dx2+x⁢y⁢dx⁢dy+y2⁢dy2+z2⁢dz2−dt2
Setup⁡coordinates=cartesian,metric=ds2
Coordinates: x,y,z,t. Signature: + + + -
gμ,ν=x2x⁢y200x⁢y2y20000z20000−1
_______________________________________________________
coordinatesystems=X,metric=1,1=x2,1,2=x⁢y2,2,2=y2,3,3=z2,4,4=−1
The divergence is now not zero
Compute for instance its first component, for ν=1
simplify⁡SumOverRepeatedIndices⁡D3_μ⁡gamma3_~mu,~1
See Also
Christoffel, CompactDisplay, Coordinates, D_, Define, g_, Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Setup, ThreePlusOne:-gamma3_, value
References
[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.
Compatibility
The Physics[ThreePlusOne][D3_] command was introduced in Maple 2017.
For more information on Maple 2017 changes, see Updates in Maple 2017.
Download Help Document