D3_ - Maple Help
For the best experience, we recommend viewing online help using Google Chrome or Microsoft Edge.

Online Help

All Products    Maple    MapleSim


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

withPhysics:withThreePlusOne

_______________________________________________________

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_

(1)

Setupmathematicalnotation=true

mathematicalnotation=true

(2)

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_

Setupcoordinatesystems=cartesian

Systems of spacetime coordinates are:X=x,y,z,t

coordinatesystems=X

(3)

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

(4)

Check the nonzero components of Christoffel: since this is a Minkowski spacetime, there are none, then ▿__j=∂__j

Christoffelnonzero

Γα,μ,ν=

(5)

The nonzero components of Christoffel3 are also none

Christoffel3nonzero

Γα,μ,ν=

(6)

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νν

(7)

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μ,ν=r2mr0000r20000r2sinθ200002mrr

(8)

The covariant derivative of a scalar is always equal to the  derivative

D3_jΦX

γjμjμμΦX

(9)

The covariant differential of a scalar function

D3_ΦX

μΦXⅆXμμ

(10)

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

DefineAj

Defined objects with tensor properties

Aj,μ,μ,γμ,σμ,Rμ,ν,Rμ,ν,Rμ,ν,α,β,βμ,Cμ,ν,α,β,μ,gμ,ν,γμ,ν,γi,j,Γμ,ν,α,Γμ,ν,α,Gμ,ν,εα,β,μ,ν,Rμ,ν,α,β,tμ,Κμ,ν,nμ,Xμ

(11)

Use a compact display for this function to avoid redundant repeated display of its functionality plus having derivatives displayed indexed by the differentiation variables

CompactDisplayAX

Ar,θ,φ,twill now be displayed asA

(12)

Now we have, for the total 3D covariant differential and the 3D covariant derivative, respectively,

CompactDisplayAjX

Ar,θ,φ,twill now be displayed asA

(13)

D3_A~jX

μAjjⅆXμμ

(14)

D3_jA~jX

jAjj

(15)

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βα,iAii

(16)

The generalized divergence of a vector in curvilinear coordinates

D3_jA~jX

jAjj

(17)

expand

γjαjαγβjβjαAββ+Γβα,iβα,iAii

(18)

TensorArray

A11r2mrr+mA11+A22θ2mrr+2A112mr+A33φ2mrr+cotθA222mrr2mrr

(19)

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_igamma3_j,k

iγj,k

(20)

value

0

(21)

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

Setupmetric=1,1=f__1X,1,4=f__14X,2,2=f__2X,2,4=f__24X,3,3=f__3X,3,4=f__34X,4,4=f__4X

_______________________________________________________

Coordinates: r,θ,φ,t. Signature: + + + -

_______________________________________________________

gμ,ν=f__1X00f__14X0f__2X0f__24X00f__3Xf__34Xf__14Xf__24Xf__34Xf__4X

_______________________________________________________

metric=1,1=f__1X,1,4=f__14X,2,2=f__2X,2,4=f__24X,3,3=f__3X,3,4=f__34X,4,4=f__4X

(22)

D3_μgamma3_~mu,~nu

0

(23)

Set the spacetime metric g_ by giving the square of the spacetime interval such that g1,20

ds2x2dx2+y2dy2+z2dz2+xydxdydt2

ds2x2dx2+xydxdy+y2dy2+z2dz2dt2

(24)

Setupcoordinates=cartesian,metric=ds2

Systems of spacetime coordinates are:X=x,y,z,t

_______________________________________________________

Coordinates: x,y,z,t. Signature: + + + -

_______________________________________________________

gμ,ν=x2xy200xy2y20000z20000−1

_______________________________________________________

_______________________________________________________

coordinatesystems=X,metric=1,1=x2,1,2=xy2,2,2=y2,3,3=z2,4,4=−1

(25)

The divergence is now not zero

D3_μgamma3_~mu,~nu

0

(26)

Compute for instance its first component, for ν=1

simplifySumOverRepeatedIndicesD3_μgamma3_~mu,~1

0

(27)

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.