Physics
Maple provides a state-of-the-art environment for algebraic computations in Physics, with emphasis on ensuring that the computational experience is as natural as possible. The theme of the Physics project for Maple 2017 has been the consolidation of the functionality introduced in previous releases, together with significant enhancements and new functionality in General Relativity, in connection with classification of solutions to Einstein's equations and tensor representations to work in the embedded 3D curved space, relevant in numerical relativity, and Field Theory, in connection with computational representations for the objects entering the Standard Model in particle physics.
Taking all together, there are more than 300 enhancements throughout the entire package increasing robustness, versatility and functionality, extending once more the range of Physics-related algebraic computations that can be done using computer algebra software and in a natural way.
As part of its commitment to providing the best possible environment for algebraic computations in Physics, Maplesoft launched a Maple Physics: Research and Development web site with Maple 18, which enabled users to download research versions, ask questions, and provide feedback. The results from this accelerated exchange with people around the world have been incorporated into the Physics package in Maple 2017.
General Relativity: classification of solutions to Einstein's equations and the Tetrads package
The 3D metric and the ThreePlusOne (3 + 1) new Physics subpackage
Tensors in Special and General Relativity
The StandardModel new Physics subpackage
In Maple 2016, the digitizing of the database of solutions to Einstein's equations was finished, added to the standard Maple library, with all the metrics from "Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C.; and Herlt, E., Exact Solutions to Einstein's Field Equations". These metrics can be loaded to work with them, or change them, or searched using g_ (the Physics command representing the spacetime metric that also sets the metric to your choice in one go) or using the command DifferentialGeometry:-Library:-MetricSearch.
Related to these developments, in Maple 2017, the Physics:-Tetrads package has been vastly improved and extended, now including new commands like PetrovType and SegreType to classify these metrics, and the TransformTetrad now has an option canonicalform to automatically derive a transformation and put the tetrad in canonical form (reorientation of the axis of the local system of references), a relevant step in resolving the equivalence between two metrics.
In the PDEtools package, you have the mathematical tools - including a complete symmetry approach - to work with the underlying partial differential equations. By combining the functionality of the Physics:-Tetrads package, the Physics:-TransformCoordinates command, and the ability to compute Riemann invariants and Weyl scalars, you can also formulate and, depending on the metrics also resolve, the equivalence problem; that is: to answer whether or not, given two metrics, they can be obtained from each other by a transformation of coordinates, as well as compute the transformation.
Depending on the context, or textbook, two different signatures are used in the description of spacetime: (− + + +), used for instance in the database of solutions of Einstein's equations and related textbook, and (+ − − −) typically used in particle Physics and in Landau's Course for Theoretical Physics, the main reference for the Maple Physics project. Two other computational representations of those same signatures are (+ + + −) and (− − − +), with time in position 4; for the computer, this difference in representation is relevant since it indicates the computational position where the component 0 is to be found (either 1 or 4). This difference also changes the position of the time coordinate in the ordered list of the coordinates. Changing the signature to follow a textbook or just for convenience thus entails redefining the metric rearranging lines and columns in its matrix representation and/or changing the sign, and reordering the coordinates, operations prone to mistakes. To handle these redefinitions of the metric and coordinates according to a change in the signature, a new command Redefine (options fromsignature, tosignature), got added to the Physics package.
Examples
Petrov and Segre types, tetrads in canonical form
restart;withPhysics:withTetrads;
Setting lowercaselatin letters to represent tetrad indices
Defined as tetrad tensors see ?Physics,tetrads,𝔢a,μ,ηa,b,γa,b,c,λa,b,c
Defined as spacetime tensors representing the NP null vectors of the tetrad formalism see ?Physics,tetrads,lμ,nμ,mμ,m&conjugate0;μ
IsTetrad,NullTetrad,OrthonormalTetrad,PetrovType,SegreType,TransformTetrad,e_,eta_,gamma_,l_,lambda_,m_,mb_,n_
There are six Petrov types: I, II, III, D, N and O. Start with a spacetime metric of Petrov type "I" (the numbers always refer to the equation number in the "Exact solutions to Einstein's field equations" textbook)
g_12,21,1
Systems of spacetime Coordinates are: X=t,x,y,φ
Default differentiation variables for d_, D_ and dAlembertian are: X=t,x,y,φ
The McLenaghan, Tariq (1975), Tupper (1976) metric in coordinates t,x,y,φ
Parameters: a,k,κ0
Comments: k paramⅇtrⅈzⅇs thⅇ most gⅇnⅇral ⅇlⅇctromagnⅇtⅈc ⅈnvarⅈant wⅈth rⅇspⅇct to thⅇ last 3 Kⅈllⅈng vⅇctors
Resetting the signature of spacetime from "- - - +" to `- + + +` in order to match the signature in the database of metrics:
The Weyl scalars are derived from the Weyl tensor and the null vectors of the Newman-Penrose formalism (l_, n_, m_ and mb_)
Weylscalarsdefinition
ψ__0=C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢mν⁢lα⁢mβ,ψ__1=C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢nν⁢lα⁢mβ,ψ__2=C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢mν⁢m&conjugate0;α⁢nβ,ψ__3=C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢nν⁢m&conjugate0;α⁢nβ,ψ__4=C⁢μ,ν,α,β⁢μ,ν,α,β⁢nμ⁢m&conjugate0;ν⁢nα⁢m&conjugate0;β
Weylscalars
ψ__0=I⁢xa2⁢x,ψ__1=0,ψ__2=−1a2,ψ__3=0,ψ__4=I⁢xa2⁢x
For presentation simplicity, without absolute values, assume that x is positive or 0 and recompute the Weyl scalars
Assume⁡x ≥ 0
x::0,∞
ψ__0=Ia2,ψ__1=0,ψ__2=−1a2,ψ__3=0,ψ__4=Ia2
Relevant in connection with the equivalence problem between two metrics, the Petrov type of this metric is
PetrovType
I
and the Plebanski-Segre type is
SegreType;
D,[(1,1)(11)]
Two spacetime metrics with different classification cannot be different coordinate representations of the same metric (if the classification is the same the question is still undecided). The same computation but tracking the principal roots behind this classification
PetrovType⁡principalroots=ROOTS
ROOTS
12+I2⁢2⁢5−3⁢2,−12−I2⁢2⁢5−3⁢2,12−I2⁢2⁢5+3⁢2,−12+I2⁢2⁢5+3⁢2
The principal polynomial only, having four different roots corresponding to Petrov type "I"
PetrovType⁡principalpolynomial=only
Z4+1+6⁢I⁢Z2
An example of Petrov type III and a canonical form for the corresponding tetrad.
g_12, 35, 1
Systems of spacetime Coordinates are: X=u,x,y,z
Default differentiation variables for d_, D_ and dAlembertian are: X=u,x,y,z
The Kaigorodov (1962), Cahen (1964), Siklos (1981), Ozsvath (1987) metric in coordinates u,x,y,z
Parameters: Λ
For presentation simplicity, assume that z and the parameter Λ are positive
Assume⁡0<z,0<Λ
Λ::0,∞,z::0,∞
The Petrov and Plebanski-Segre types, the default tetrad and the corresponding Weyl scalars:
PetrovType⁡
III
SegreType
O,[(1,111)]
e_
𝔢a,μ=−I⁢ⅇ−4⁢z−I⁢ⅇ2⁢z−2⁢I⁢ⅇ−z0I⁢ⅇ−4⁢z⁢20I⁢2⁢ⅇ−z00003Λⅇ−4⁢z000
ψ__0=11⁢Λ4−2⁢Λ⁢2,ψ__1=3⁢Λ4−Λ⁢22,ψ__2=Λ4,ψ__3=−3⁢Λ4−Λ⁢22,ψ__4=11⁢Λ4+2⁢Λ⁢2
These scalars are not in canonical form. Compute then a canonical form for this tetrad (i.e. rotate the tetrad such that the Weyl scalars match the canonical form shown in the Petrov classification table) and set this rotated tetrad (using Setup), then recompute the Weyl scalars
Setup⁡e_=TransformTetrad⁡canonicalform:
ψ__0=0,ψ__1=0,ψ__2=0,ψ__3=1,ψ__4=0
The scalars now match the requirement for canonical form shown in the Petrov classification table for type III. In this example, the canonical form computed for the tetrad
𝔢a,μ=0−I2⁢Λ⁢ⅇ2⁢z⁢2−I2⁢Λ⁢ⅇ−z⁢2Λ⁢32I2⁢ⅇ−4⁢z⁢2Λ−13⁢I8⁢ⅇ2⁢z⁢2Λ−9⁢I8⁢ⅇ−z⁢2Λ−7⁢38⁢Λ3203⁢I4⁢ⅇ2⁢z⁢2I4⁢ⅇ−z⁢2−34⁢Λ−I⁢ⅇ−4⁢z⁢23⁢I4⁢ⅇ2⁢z⁢2I4⁢ⅇ−z⁢2−34⁢Λ
is not much more complicated than the default tetrad (14) that is not in canonical form, but depending on the example the canonical form can be significantly more complicated.
Examples of transforming the tetrad into canonical form for spacetimes of Petrov types II, N and D can be seen by performing the same steps above departing from the metrics of the database with numbers g_[[24, 37, 7]], g_[[12, 6, 1]] and g_[[12, 8, 4]].
Together with the developments above, the TransformTetrad got thoroughly reviewed and its transformation functionality extended in order to perform the six traditional rotations that leave different null vectors of the Newman-Penrose formalism unchanged in direction or length,
Load the Schwarzschild metric in spherical coordinates; you can input g_[Schwarzschild] or the simpler
g_sc
Systems of spacetime Coordinates are: X=t,r,θ,φ
Default differentiation variables for d_, D_ and dAlembertian are: X=t,r,θ,φ
The Schwarzschild metric in coordinates t,r,θ,φ
Parameters: m
An orthonormal tetrad (this is the one used by default) is
Setuptetrad = orthonormal: e_
𝔢a,μ=−I⁢−r+2⁢mr0000−I⁢r−r+2⁢m0000r0000r⁢sin⁡θ
IsTetrad%
Type of tetrad: orthonormal
true
The key observation here is that the orientation of the axis of the local system of references (the tetrad system) is arbitrary and with them the values of the Weyl scalars too. So, you can perform a Lorentz transformation on the tetrad, without changing the spacetime metric g_ or its Petrov classification while fixing the values of the Weyl scalars in one of the possible forms (classes) shown in the Petrov classification table.
The transformations that you can perform with TransformTetrad are either the standard ones, of classes I, II or III specified in the description of TransformTetrad, or arbitrary. For example, transform the orthonormal tetrad above using a null rotation with fixed l_, where l_ and the related NULL vectors of the Newman-Penrose formalism (commands of the Tetrads package) satisfy
l_definition
lμ⁢l⁢μ⁢μ=0,lμ⁢n⁢μ⁢μ=−1,lμ⁢m⁢μ⁢μ=0,lμ⁢m&conjugate0;⁢μ⁢μ=0,gμ,ν=−lμ⁢nν−lν⁢nμ+mμ⁢m&conjugate0;ν+mν⁢m&conjugate0;μ
TransformTetrad⁡nullrotationwithfixedl_
−I2⁢2⁢−r+2⁢mr−I2⁢2⁢r−r+2⁢m00−I2⁢2⁢−r+2⁢m⁢E⁢E&conjugate0;+1r−I2⁢2⁢r⁢E⁢E&conjugate0;−1−r+2⁢m2⁢r⁢E+E&conjugate0;2−I2⁢2⁢r⁢sin⁡θ⁢E−E&conjugate0;−I2⁢E⁢2⁢−r+2⁢mr−I2⁢E⁢2⁢r−r+2⁢m2⁢r2I2⁢2⁢r⁢sin⁡θ−I2⁢E&conjugate0;⁢2⁢−r+2⁢mr−I2⁢E&conjugate0;⁢2⁢r−r+2⁢m2⁢r2−I2⁢2⁢r⁢sin⁡θ
Note the transformation parameter E introduced by TransformTetrad. You can replace it by a value (see subs and eval) or assign it, or you may want to Assume that this parameter is, for instance, real
Assume⁡E::real
E::real
map⁡eval,
−I2⁢2⁢−r+2⁢mr−I2⁢2⁢r−r+2⁢m00−I2⁢2⁢−r+2⁢m⁢E2+1r−I2⁢2⁢r⁢E2−1−r+2⁢mE⁢2⁢r0−I2⁢E⁢2⁢−r+2⁢mr−I2⁢E⁢2⁢r−r+2⁢m2⁢r2I2⁢2⁢r⁢sin⁡θ−I2⁢E⁢2⁢−r+2⁢mr−I2⁢E⁢2⁢r−r+2⁢m2⁢r2−I2⁢2⁢r⁢sin⁡θ
Now is this a tetrad? Of what type?
IsTetrad⁡
Type of tetrad: null
TransformTetrad applies a transformation on a tetrad but does not set the value of e_ to this result. For that purpose, i.e. to effectively rotate the tetrad system of references to match the orientation/tetrad you want, use Setup. Following this reorientation of the tetrad, all the components of tensors in the local (tetrad) system of references, as for instance gamma_, lambda_ and the tetrad components of spacetime tensors, automatically change. For example, the Ricci rotation coefficients gamma_ follow this change of value of the components of 𝔢__a,μ according to their definition
gamma_definition
γa,b,c=𝒟ν⁡𝔢a,μ⁢𝔢bμbμ⁢𝔢cνcν
So before setting (23) as the new tetrad, for comparison purposes below, track first the value of the Ricci rotation coefficient with covariant indices 1, 2, 1
gamma_1,2,1
−I⁢m−r+2⁢m⁢r32
Set now (23) as the new tetrad
Setup⁡e_=
tetrad=1,1=−I2⁢2⁢−r+2⁢mr,1,2=−I2⁢2⁢r−r+2⁢m,2,1=−I2⁢2⁢−r+2⁢m⁢E2+1r,2,2=−I2⁢2⁢r⁢E2−1−r+2⁢m,2,3=E⁢2⁢r,3,1=−I2⁢E⁢2⁢−r+2⁢mr,3,2=−I2⁢E⁢2⁢r−r+2⁢m,3,3=2⁢r2,3,4=I2⁢2⁢r⁢sin⁡θ,4,1=−I2⁢E⁢2⁢−r+2⁢mr,4,2=−I2⁢E⁢2⁢r−r+2⁢m,4,3=2⁢r2,4,4=−I2⁢2⁢r⁢sin⁡θ
The new value of γ1,2,1
I2⁢2⁢m−r+2⁢m⁢r32
Equivalence for Schwarzschild metric (spherical and Kruskal coordinates)
Formulation of the problem (remove mixed coordinates)
restart;
withPhysics: withTetrads: Setupauto=true, mathematicalnotation = true;
* Partial match of 'auto' against keyword 'automaticsimplification'
automaticsimplification=true,mathematicalnotation=true
The departure point, Schwarzschild metric in spherical coordinates
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
gμ,ν=r−r+2⁢m0000−r20000−r2⁢sin⁡θ20000r−2⁢mr
Introduce now Kruskal coordinates following the literature (see Wikipedia) and the corresponding line element involving "mixed" coordinates
CoordinatesK=u,ϑ,ϕ,v
Systems of spacetime Coordinates are: K=u,ϑ,ϕ,v,X=r,θ,φ,t
K,X
ds2≔16⁢∂⁡v⁢m2⁢ⅇ−r2⁢m⁢∂⁡u−1−cos⁡ϑ2⁢∂⁡ϕ2+∂⁡ϑ2⁢r3r
The mixing of variables is visible: in the line element above is in Kruskal coordinates but you also see r , which belongs to the X (not K) coordinates.
For the purpose of formulating problem free of this mixing of coordinates, set the metric now to be (32)
Setupdiff = K,metric=, quiet:
g_
gμ,ν=0008⁢m2⁢ⅇ−r2⁢mr0−r20000−r2⁢sin⁡ϑ208⁢m2⁢ⅇ−r2⁢mr000
To remove the mix of coordinates, introduce a transformation with unknown transformation functions f,h, change variables, and resolve for the transformation functions f,g (this in itself is resolving a form of equivalence problem).
tr__0 ≔ u=fr,t, v=hr,t
tr__0≔u=f⁡r,t,v=h⁡r,t
CompactDisplaytr__0
f⁡r,t⁢will now be displayed as⁢f
h⁡r,t⁢will now be displayed as⁢h
TransformCoordinatestr__0,g_mu,nu
16⁢fr⁢hr⁢m2⁢ⅇ−r2⁢mr008⁢m2⁢ⅇ−r2⁢m⁢fr⁢ht+hr⁢ftr0−r20000−r2⁢sin⁡θ208⁢m2⁢ⅇ−r2⁢m⁢fr⁢ht+hr⁢ftr0016⁢ft⁢ht⁢m2⁢ⅇ−r2⁢mr
Equate to (30) and solve
convert=rhs,setofequations
0=0,−r2=−r2,−r2⁢sin⁡θ2=−r2⁢sin⁡θ2,8⁢m2⁢ⅇ−r2⁢m⁢fr⁢ht+hr⁢ftr=0,16⁢fr⁢hr⁢m2⁢ⅇ−r2⁢mr=r−r+2⁢m,16⁢ft⁢ht⁢m2⁢ⅇ−r2⁢mr=r−2⁢mr
pdsolve
f=_C1+_C2⁢r−2⁢m⁢ⅇr4⁢m⁢ⅇ−t4⁢m,h=−ⅇr+t4⁢m⁢r−2⁢m_C2+_C3,f=_C1+_C2⁢r−2⁢m⁢ⅇr4⁢m⁢ⅇt4⁢m,h=−r−2⁢m⁢ⅇr−t4⁢m_C2+_C3
Without loss of generality, set _C1=0,_C2=1,_C3=0
tr ≔ combinesubs_C1=0,_C2=1,_C3=0, eval,1
tr≔u=r−2⁢m⁢ⅇr−t4⁢m,v=−ⅇr+t4⁢m⁢r−2⁢m
Check it out:
TransformCoordinatestr,g_mu,nu,X,K
r−r+2⁢m0000−r20000−r2⁢sin⁡θ20000r−2⁢mr
Compute the inverse of the transformation (39)
itr ≔ simplifynormalsolve,r,t,expanded
itr≔r=2⁢W⁡−u⁢v⁢ⅇ−12⁢m+1⁢m,t=2⁢ln⁡−vu⁢m
This inverse transformation involves the LambertW function. Set now the metric to be the standard Schwarzschild's metric in spherical coordinates (30) and use itr to get he form of the metric entirely in Kruskal coordinates
So this is Schwarzschild's solution all in Kruskal coordinates
TransformCoordinatesitr,g_mu,nu,K,X
000−8⁢W⁡−u⁢v⁢ⅇ−12⁢m⁢m2W⁡−u⁢v⁢ⅇ−12⁢m+1⁢u⁢v0−4⁢W⁡−u⁢v⁢ⅇ−12⁢m+12⁢m20000−4⁢W⁡−u⁢v⁢ⅇ−12⁢m+12⁢m2⁢sin⁡ϑ20−8⁢W⁡−u⁢v⁢ⅇ−12⁢m⁢m2W⁡−u⁢v⁢ⅇ−12⁢m+1⁢u⁢v000
This metric involves the LambertW function in a non-simplifiable form.
Solving the Equivalence
We now have the two forms: (41) in spherical and (44) in Kruskal coordinates, so we can formulate the equivalence problem from one coordinate system to the other one.
The transformation to be resolved does not need to involve φ because neither φ nor ϕ enter either of the two metrics.
The transformation does not need to involve θ or ϑ because they enter the metrics in exactly the same position and with the same dependence.
In addition the Weyl scalars of both metrics are in canonical form and the only scalar different from zero, that is Ψ__2 does not depend on any of φ,θ, ϕ,ϑ)
So we look for a generic transformation from spherical to Kruskal of the form
r=RK, t=TK
r=R⁡K,t=T⁡K
CompactDisplay
R⁡u,ϑ,ϕ,v⁢will now be displayed as⁢R
T⁡u,ϑ,ϕ,v⁢will now be displayed as⁢T
The metric set in this moment is in spherical coordinates, (43), so change using (45) and equate to (44) in Kruskal coordinates
convertTransformCoordinates,g_mu,nu,K,X=, setofequations
−4⁢−R2+m2⁢Tu2+Ru2⁢R2−R+2⁢m⁢R=0,−4⁢−R2+m2⁢Tv2+Rv2⁢R2−R+2⁢m⁢R=0,−4⁢−R2+m2⁢Tϕ2+2⁢R2⁢Rϕ22+cos⁡ϑ+1⁢−R2+m⁢cos⁡ϑ−1⁢R−R+2⁢m⁢R=−4⁢W⁡−u⁢v⁢ⅇ−12⁢m+12⁢m2⁢sin⁡ϑ2,−4⁢Tv⁢−R2+m2⁢Tu+Rv⁢Ru⁢R2−R+2⁢m⁢R=−8⁢W⁡−u⁢v⁢ⅇ−12⁢m⁢m2W⁡−u⁢v⁢ⅇ−12⁢m+1⁢u⁢v,−4⁢Tv⁢−R2+m2⁢Tϕ+Rv⁢Rϕ⁢R2−R+2⁢m⁢R=0,−4⁢Tv⁢−R2+m2⁢Tϑ+Rv⁢Rϑ⁢R2−R+2⁢m⁢R=0,−4⁢Tϕ⁢−R2+m2⁢Tu+Rϕ⁢Ru⁢R2−R+2⁢m⁢R=0,−4⁢Tϕ⁢−R2+m2⁢Tϑ+Rϕ⁢Rϑ⁢R2−R+2⁢m⁢R=0,−4⁢Tϑ⁢−R2+m2⁢Tu+Rϑ⁢Ru⁢R2−R+2⁢m⁢R=0,−Rϑ2⁢RR−2⁢m−R2+Tϑ2⁢R−2⁢mR=−4⁢W⁡−u⁢v⁢ⅇ−12⁢m+12⁢m2
This is a nonlinear, non-rational PDE system in two unknowns depending on two independent variables (see (45)). You can now either call pdsolve on (47), solving the problem in one step, or first split into cases without solving any differential equation, just doing differential elimination, to see the cases
PDEtools:-casesplit
Tu=2⁢mu,Tϑ=0,Tϕ=0,Tv=−2⁢mv,R=2⁢LambertW⁡−12⁢u⁢v⁢ⅇ−1m+1⁢m &where ,Tu=−2⁢mu,Tϑ=0,Tϕ=0,Tv=2⁢mv,R=2⁢LambertW⁡−12⁢u⁢v⁢ⅇ−1m+1⁢m &where
By using differential elimination we removed all nonlinearities, the problem is actually easy for the differential equation routines
R=2⁢W⁡−u⁢v⁢ⅇ−12⁢m+1⁢m,T=−2⁢m⁢ln⁡u+2⁢m⁢ln⁡v+_C1,R=2⁢W⁡−u⁢v⁢ⅇ−12⁢m+1⁢m,T=2⁢m⁢ln⁡u−2⁢m⁢ln⁡v+_C1
So the transformation of coordinates resolving the equivalence between (41) and (44) is
eval,1
r=2⁢W⁡−u⁢v⁢ⅇ−12⁢m+1⁢m,t=−2⁢m⁢ln⁡u+2⁢m⁢ln⁡v+_C1
Check this result transforming (41) fully written in spherical coordinates into (44) fully written in Kruskal coordinates
TransformCoordinates,g_mu,nu,K,X
Redefining the spacetime metric and coordinates according to the signature
restart; withPhysics:
CoordinatesX=cartesian
Default differentiation variables for d_, D_ and dAlembertian are: X=x,y,z,t
Systems of spacetime Coordinates are: X=x,y,z,t
X
Consider the generic form of a spacetime metric, for the purpose of illustrating its redefinition (reordering of lines and columns) when changing the signature, it doesn't need to depend on the coordinates
G ≔ Matrix4,symbol=g, shape=symmetric
Set the metric, and in the same call check the value of signature: it is (- - - +)
Setupsignature, metric = G
metric=1,1=g1,1,1,2=g1,2,1,3=g1,3,1,4=g1,4,2,2=g2,2,2,3=g2,3,2,4=g2,4,3,3=g3,3,3,4=g3,4,4,4=g4,4,signature=- - - +
Track the line element
g_lineelement
g3,3⁢∂⁡z2+2⁢∂⁡x⁢g1,3+2⁢g2,3⁢∂⁡y⁢∂⁡z+g1,1⁢∂⁡x2+2⁢g1,2⁢∂⁡x⁢∂⁡y+g2,2⁢∂⁡y2+∂⁡t⁢∂⁡t⁢g4,4+2⁢∂⁡x⁢g1,4+2⁢∂⁡y⁢g2,4+2⁢∂⁡z⁢g3,4
Change now the signature by reversing the position of the time-like component, from 4 to 1
Setupsignature=+---
signature=+ - - -
Neither the metric nor the coordinates change: t is still in position 4 and the component (1, 1) of the metric is still g__11
Coordinates
This design gives you freedom to set things as necessary. However, after changing the signature we may also want to redefine the coordinates - place t in position 1 - and possibly also the metric reordering its lines and rows accordingly. A new command for that purpose is Redefine, to which you need to indicate the previous signature (or a new signature, to explore the effect of a change before doing it) and an indication of whether you want to redefine the metric, the coordinates or all:
Redefinecoordinates, fromsignature = ---+
t,x,y,z
Redefinemetric, fromsignature = ---+
These redefinitions however were not set, the keywords coordinates, metric, or all only trigger the change, t is still in position 4 and the component (1,1) of the metric is still g__1,1
You can now either use the output of this routine to redefine things using the Setup command, or do all in one go using any of the keywords setcoordinates, setmetric, or setall, for example
Redefinesetall, fromsignature = ---+
Now t in position 1 and the component (1,1) of the metric equal to g__4,4
Systems of spacetime Coordinates are: X=t,x,y,z
Note that, despite the reordering of lines and columns in the metric, because we also reordered the variables, the line element has not changed:
g2,2⁢∂⁡y2+2⁢∂⁡t⁢g2,4+2⁢∂⁡x⁢g1,2⁢∂⁡y+g4,4⁢∂⁡t2+2⁢g1,4⁢∂⁡x⁢∂⁡t+g1,1⁢∂⁡x2+∂⁡z⁢2⁢∂⁡x⁢g1,3+2⁢g2,3⁢∂⁡y+2⁢∂⁡t⁢g3,4+∂⁡z⁢g3,3
normal−
0
This routine is particularly useful when working with metrics from the database of solutions to Einstein's equations, all of which reset the signature to (-+++) when loaded. For example:
Resetting the signature of spacetime from "+ - - -" to `- + + +` in order to match the signature in the database of metrics:
Note above the message about resetting the signature; query about:
Setupsignature
signature=- + + +
How would this metric and coordinates look with the original signature (---+) ?
Redefineall, tosignature = ---+
By entering the command above replacing all by setall not only the list of coordinates and metric matrix form are returned but they are also set in one go.
In the general theory of relativity, x1, x2, x3 can be any quantities defining the position of bodies in space and the time coordinate x0 can be defined by an arbitrarily running clock. In order to define simultaneity (synchronizing clocks located at different points in space) as well as determine actual space distances and time intervals in terms of these quantities x0, x1, x2, x3, it is relevant to split the spacetime mathematical description of gravity into its space and time parts. For this purpose, a new package, ThreePlusOne, and three new commands, were added to the Physics package:
Decompose, to decompose 4D tensorial expressions (free and/or contracted indices) into the space and time parts.
gamma_, representing the three-dimensional metric tensor, with which the element of spatial distance is defined as dl 2=γi,j⁢dxi⁢dxj.
Redefine, to redefine the coordinates and the spacetime metric according to changes in the signature from any of the four possible signatures(− + + +), (+ − − −), (+ + + −) and ((− + + +) to any of the other ones.
restart;with⁡Physics:
Setup⁡mathematicalnotation=true
mathematicalnotation=true
Define now an arbitrary tensor A
Define⁡A
Defined objects with tensor properties
A,γμ,σμ,∂μ,gμ,ν,δμ,ν,εα,β,μ,ν
So Aμ is a 4D tensor with only one free index, where the position of the time-like component is the position of the different sign in the signature, that you can query about via
Setup⁡signature
signature=- - - +
To perform a decomposition into space and time, set - for instance - the lowercase latin letters from i to s to represent spaceindices and
Setupspaceindices = lowercase_is
spaceindices=lowercaselatin_is
Accordingly, the 3+1 decomposition of Aμ is
DecomposeA~mu;
A⁢i⁢iA⁢0⁢0
The 3+1 decomposition of the inert representation %g_[mu,nu] of the 4D spacetime metric; use the inert representation when you do not want the actual components of the metric appearing in the output
Decompose%g_mu, nu;
gi,jgi,0g0,jg0,0
Note the position of the component %g_[0, 0], related to the trailing position of the time-like component in the signature − − − +.
Compare the decomposition of the 4D inert with the decomposition of the 4D active spacetime metric
g
Decomposeg_mu, nu;
gi,j001
Note that in general the 3D space part of gμ,ν is not equal to the 3D metric γi,j whose definition includes another term (see [1] Landau & Lifshitz, eq.(84.7)).
γi,j=−gi,j+g0,i⁢g0,jg0,0
The 3D space part of −g⁢μ,ν⁢μ,ν is actually equal to the 3D metric γ⁢i,j⁢i,j
gamma_~,definition;
γ⁢i,j⁢i,j=−g⁢i,j⁢i,j
To derive the formula (71) for the covariant components of the 3D metric, Decompose into 3+1 the identity
%g_~alpha, ~mu %g_mu, beta = KroneckerDelta~alpha, beta;
g⁢α,μ⁢α,μ⁢gμ,β=δβαβα
To the side, for illustration purposes, these are the 3 + 1 decompositions, first excluding the repeated indices, then excluding the free indices
Eq≔Decompose⁡,repeatedindices=false
Eq≔g⁢i,μ⁢i,μ⁢gμ,j=δjijig⁢i,μ⁢i,μ⁢gμ,0=0g⁢0,μ⁢0,μ⁢gμ,j=0g⁢0,μ⁢0,μ⁢gμ,0=1
Eq≔Decompose⁡,freeindices=false
Eq≔g0,β⁢g⁢α,0⁢α,0+gi,β⁢g⁢α,i⁢α,i=δβαβα
Compare with a full decomposition
Eq≔Decompose⁡
Eq≔g0,k⁢g⁢j,0⁢j,0+gi,k⁢g⁢j,i⁢j,i=δkjkjg0,0⁢g⁢j,0⁢j,0+gi,0⁢g⁢j,i⁢j,i=0g0,k⁢g⁢0,0⁢0,0+gi,k⁢g⁢0,i⁢0,i=0g0,0⁢g⁢0,0⁢0,0+gi,0⁢g⁢0,i⁢0,i=1
Eq is a symmetric matrix of equations involving non-contracted occurrences of g0,0, gj,0 and gj,i. Isolate, in Eq1,2, gj,0, that you input as %g_[~j, ~0], and substitute into Eq1,1
isolateEq1, 2, %g_~j, ~0;
g⁢j,0⁢j,0=−gi,0⁢g⁢j,i⁢j,ig0,0
subs⁡,Eq1,1
−g0,k⁢gi,0⁢g⁢j,i⁢j,ig0,0+gi,k⁢g⁢j,i⁢j,i=δkjkj
Collect gj,i, that you input as %g_[~j, ~i]
collect⁡,%g_~j, ~i
−g0,k⁢gi,0g0,0+gi,k⁢g⁢j,i⁢j,i=δkjkj
Since the right-hand side is the identity matrix and, from (72), gi,j=− γi,j, the expression between parenthesis, multiplied by -1, is the reciprocal of the contravariant 3D metric γi,j, that is the covariant 3D metric γi,j, in accordance to its definition for the signature - - - +
gamma_definition;
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.
Maple 2017 introduces a new package, ThreePlusOne, to cast Einstein's equations in a 3+1 form, that is, representing spacetime as a stack of nonintersecting 3-hypersurfaces Σ (not necessarily actual space). This 3+1 description is key in the Hamiltonian formulation of gravity as well as in the study of gravitational waves, black holes, neutron stars, and in general to study the evolution of physical system in general relativity by running numerical simulations as traditional initial value (Cauchy) problems. ThreePlusOne includes:
Computational representations for the spatial metric γi,j that is induced by gμ,ν on the 3-dimensional hypersurfaces, and the related covariant derivative, Christoffel symbols and Ricci and Riemann tensors.
Computational representations for the Lapse, Shift, Unit normal and Time vectors and Extrinsic curvature related to the ADM equations.
restart; withPhysics: Setupcoordinatesystems=cartesian;
coordinatesystems=X
withThreePlusOne
Setting lowercaselatin_is letters to represent space indices
Defined as 4D, spacetime tensors that are purely spatialsee ?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
ADMEquations,Christoffel3,D3_,ExtrinsicCurvature,Lapse,Ricci3,Riemann3,Shift,TimeVector,UnitNormalVector,gamma3_
Note the different color for γμ,ν, now a 4D tensor representing the metric of a generic 3-dimensional hypersurface (not necessarily the 3D space) induced by the 4D spacetime metric gμ,ν. All the ThreePlusOne tensors are displayed in black to distinguish them of the corresponding 4D or 3D tensors. The particular hypersurface γμ,ν operates is parameterized by the Lapse α and the Shift βμ.
The induced metric γμ,νis defined in terms of the UnitNormalVector nμ and the 4D metric gμ,ν as
gamma3_definition
γμ,ν=nμ⁢nν+gμ,ν
where nμ is defined in terms of the Lapse α and the derivative of a scalar function t that can be interpreted as a global time function
UnitNormalVectordefinition
nμ=−α⁢𝒟μ⁡t
The TimeVector is defined in terms of the Lapse α and the Shift βμ and this vector nμ as
TimeVectordefinition
tμ=α⁢nμ+βμ
The ExtrinsicCurvature is defined in terms of the LieDerivative of γμ,ν
ExtrinsicCurvaturedefinition
Κμ,ν=−ℒn⁡γμ,ν2
The metric γμ,νis also a projection tensor in that it projects 4D tensors into the 3D hypersurface Σ. The definition for any 4D tensor that is also a 3D tensor in Σ, can thus be written directly by contracting their indices with γμ,ν. In the case of Christoffel3, Ricci3 and Riemann3, these tensors can be defined by replacing the 4D metric gμ,ν by γμ,ν and the 4D Christoffel symbols Γμ,ν,α by the ThreePlusOne Γμ,ν,α in the definitions of the corresponding 4D tensors. So, for instance
Christoffel3definition
Γμ,ν,α=∂α⁡γμ,ν2+∂ν⁡γμ,α2−∂μ⁡γν,α2
Ricci3definition
Rμ,ν=∂α⁡Γ⁢αμ,ν⁢αμ,ν−∂ν⁡Γ⁢αμ,α⁢αμ,α+Γ⁢βμ,ν⁢βμ,ν⁢Γ⁢αβ,α⁢αβ,α−Γ⁢βμ,α⁢βμ,α⁢Γ⁢αν,β⁢αν,β
Riemann3definition
Rμ,ν,α,β=γμ,λ⁢∂α⁡Γ⁢λν,β⁢λν,β−∂β⁡Γ⁢λν,α⁢λν,α+Γ⁢λκ,α⁢λκ,α⁢Γ⁢κν,β⁢κν,β−Γ⁢λκ,β⁢λκ,β⁢Γ⁢κν,α⁢κν,α
When working with the ADM formalism, the line element of an arbitrary spacetime metric can be expressed in terms of the differentials of the coordinates dxμ, the Lapse, the Shift and the spatial components of the 3D metric gamma3_. From this line element one can derive the relation between the Lapse, the spatial part of the Shift, the spatial part of the gamma3_ metric and the g0,j components of the 4D spacetime metric.
For this purpose, define a tensor representing the differentials of the coordinates and an alias dt=dx0
Definedxmu;
𝒟μ,γμ,σμ,Rμ,ν,βμ,Xμ,∂μ,dxμ,gμ,ν,γμ,ν,γi,j,Γμ,ν,α,δμ,ν,εα,β,μ,ν,Rμ,ν,α,β,tμ,Κμ,ν,nμ
aliasdt = dx~0:
The expression for the line element in terms of the Lapse and Shift is (see [2], eq.(2.123))
ds2 = Shifti2−Lapse2 dt2+2 Shifti dt dx~i+gamma_i, j dx