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 2015 has been Vector Analysis, symbolic Tensor manipulations, Quantum Mechanics, and General Relativity. With more than 400 enhancements throughout the entire package to increase robustness and versatility, two new commands, Assume and SubstituteTensor, a new Tetrads subpackage with 13 commands, as well as 26 new Physics:-Library commands to support further explorations and extensions, and an enlargement of the database of solutions to Einstein's equations with more than 100 new metrics, Maple 2015 extends again the range of Physics-related algebraic computations that can be done, in a natural way, using computer algebra software.
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 2015. Examples illustrating the use of the new capabilities in the context of more general problems are found in the MaplePrimes post Computer Algebra for Theoretical Physics.
Simplification
Tensors
Tetrads in General Relativity
More Metrics in the Database of Solutions to Einstein's Equations
Commutators, AntiCommutators, and Dirac notation in quantum mechanics
New Assume command and new enhanced Mode: automaticsimplification
Vectors Package
The Physics Library
Miscellaneous
Simplification is perhaps the most common operation performed in a computer algebra system. In Physics, this typically entails simplifying tensorial expressions, or expressions involving noncommutative operators that satisfy certain commutator/anticommutator rules, or sums and integrals involving quantum operators and Dirac delta functions in the summands and integrands. Relevant enhancements were introduced in Maple 2015 for all these cases, including enhancements in the simplification of:
Products of LeviCivita tensors in curved spacetimes when LeviCivita represents the Galilean pseudo-tensor (related to Setup(levicivita = Galilean)), instead of its generalization to curved spaces (related to Setup(levicivita = nongalilean)).
Tensorial expressions in general that have spacetime, space, and/or tetrad contracted indices, possibly at the same time.
New option tryhard, that resolves zero recognition in an important number of nontrivial situations.
Expressions involving the Dirac function.
Vectorial expressions involving cylindrical or spherical coordinates and related unit vectors.
Expressions simplified with respect to side relations (equations) in the presence of quantum vectorial equations.
Expressions involving products of quantum operators entering parameterized algebra rules.
Expressions involving vectorial quantum operators simplified with respect to other vectorial equations.
Add support for the simplification and integration of spherical harmonics (SphericalY ) relevant in quantum mechanics.
Examples
Enhancements in the simplification of tensorial expressions.
restart;withPhysics: Setupnotation = true, coordinates = X
* Partial match of 'notation' against keyword 'mathematicalnotation'
* Partial match of 'coordinates' against keyword 'coordinatesystems'
Default differentiation variables for d_, D_ and dAlembertian are: X=x1,x2,x3,x4
Systems of spacetime Coordinates are: X=x1,x2,x3,x4
coordinatesystems=X,mathematicalnotation=true
Define a tensor lμX and define the metric in terms of the components of this tensor so that gμ,ν=ημ,ν+lμ⁢lν and ημ,ν is a Minkowski metric.
lμ,ημ,ν=−rhs⁡g_Minkowski
The Minkowski metric in cartesian coordinates
Defined objects with tensor properties
γμ,σμ,Xμ,∂μ,ημ,ν,gμ,ν,lμ,δμ,ν,εα,β,μ,ν
Avoid redundant display:
PDEtools:-declarelX
l⁡x1,x2,x3,x4⁢will now be displayed as⁢l
The metric:
g_mu,nu = ημ,ν+lμ⁡X⁢lν⁡X
gμ,ν=ημ,ν+lμ⁢lν
New: you can define it directly using a tensorial equation like (4), using either Define or Setup.
Define
𝒟μ,γμ,σμ,Rμ,ν,Rμ,ν,α,β,Cμ,ν,α,β,Xμ,∂μ,ημ,ν,gμ,ν,lμ,Γμ,ν,α,Gμ,ν,δμ,ν,εα,β,μ,ν
Verify the resulting gμ,ν:
g_
gμ,ν=1+l12l1⁢l2l1⁢l3l1⁢l4l1⁢l21+l22l2⁢l3l2⁢l4l1⁢l3l2⁢l31+l32l3⁢l4l1⁢l4l2⁢l4l3⁢l4−1+l42
New: you can query about the definition you gave for any tensor directly from the tensor itself.
g_definition
Enhanced simplification capabilities: show that this expression involving derivatives of lμ⁡X is equal to zero.
−d_α⁡lλ⁡X⁢η⁢α,β⁢α,β⁢d_β⁡lκ⁡X⁢l⁢κ⁢κ⁡X⁢l⁢λ⁢λ⁡X+η⁢α,β⁢α,β⁢d_β⁡lκ⁡X⁢d_λ⁡lα⁡X⁢l⁢κ⁢κ⁡X⁢l⁢λ⁢λ⁡X+d_λ⁡lκ⁡X⁢d_ρ⁡l⁢ρ⁢ρ⁡X⁢l⁢κ⁢κ⁡X⁢l⁢λ⁢λ⁡X+d_β⁡d_λ⁡lκ⁡X⁢l⁢β⁢β⁡X⁢l⁢κ⁢κ⁡X⁢l⁢λ⁢λ⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢d_ρ⁡l⁢ρ⁢ρ⁡X⁢d_μ⁡lν⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢ρ⁢ρ⁡X⁢d_ρ⁡d_μ⁡lν⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢d_ρ⁡l⁢ρ⁢ρ⁡X⁢d_ν⁡lμ⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢ρ⁢ρ⁡X⁢d_ρ⁡d_ν⁡lμ⁡X+12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢η⁢α,ρ⁢α,ρ⁢d_ρ⁡lμ⁡X⁢d_α⁡lν⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢η⁢α,ρ⁢α,ρ⁢d_μ⁡lα⁡X⁢d_ρ⁡lν⁡X+12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢η⁢α,ρ⁢α,ρ⁢d_α⁡lμ⁡X⁢d_ρ⁡lν⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢η⁢α,ρ⁢α,ρ⁢d_ν⁡lα⁡X⁢d_ρ⁡lμ⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢α⁢α⁡X⁢l⁢ρ⁢ρ⁡X⁢d_ρ⁡lν⁡X⁢d_α⁡lμ⁡X+12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢α⁢α⁡X⁢l⁢ρ⁢ρ⁡X⁢d_μ⁡lα⁡X⁢d_ρ⁡lν⁡X−12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢α⁢α⁡X⁢l⁢ρ⁢ρ⁡X⁢d_ρ⁡lμ⁡X⁢d_α⁡lν⁡X+12⁢l⁢μ⁢μ⁡X⁢l⁢ν⁢ν⁡X⁢l⁢α⁢α⁡X⁢l⁢ρ⁢ρ⁡X⁢d_ν⁡lα⁡X⁢d_ρ⁡lμ⁡X
−∂α⁡lλ⁢η⁢α,β⁢α,β⁢∂β⁡lκ⁢l⁢κ⁢κ⁢l⁢λ⁢λ+η⁢α,β⁢α,β⁢∂β⁡lκ⁢∂λ⁡lα⁢l⁢κ⁢κ⁢l⁢λ⁢λ+∂λ⁡lκ⁢∂ρ⁡l⁢ρ⁢ρ⁢l⁢κ⁢κ⁢l⁢λ⁢λ+∂β⁡∂λ⁡lκ⁢l⁢β⁢β⁢l⁢κ⁢κ⁢l⁢λ⁢λ−l⁢μ⁢μ⁢l⁢ν⁢ν⁢∂ρ⁡l⁢ρ⁢ρ⁢∂μ⁡lν2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢ρ⁢ρ⁢∂ρ⁡∂μ⁡lν2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢∂ρ⁡l⁢ρ⁢ρ⁢∂ν⁡lμ2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢ρ⁢ρ⁢∂ρ⁡∂ν⁡lμ2+l⁢μ⁢μ⁢l⁢ν⁢ν⁢η⁢α,ρ⁢α,ρ⁢∂ρ⁡lμ⁢∂α⁡lν2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢η⁢α,ρ⁢α,ρ⁢∂μ⁡lα⁢∂ρ⁡lν2+l⁢μ⁢μ⁢l⁢ν⁢ν⁢η⁢α,ρ⁢α,ρ⁢∂α⁡lμ⁢∂ρ⁡lν2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢η⁢α,ρ⁢α,ρ⁢∂ν⁡lα⁢∂ρ⁡lμ2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢α⁢α⁢l⁢ρ⁢ρ⁢∂ρ⁡lν⁢∂α⁡lμ2+l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢α⁢α⁢l⁢ρ⁢ρ⁢∂μ⁡lα⁢∂ρ⁡lν2−l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢α⁢α⁢l⁢ρ⁢ρ⁢∂ρ⁡lμ⁢∂α⁡lν2+l⁢μ⁢μ⁢l⁢ν⁢ν⁢l⁢α⁢α⁢l⁢ρ⁢ρ⁢∂ν⁡lα⁢∂ρ⁡lμ2
Simplify
0
A new tryhard option significantly improves zero recognition, which also got improved regardless of the new option.
Set spacetimeindices to be represented with lowercase Latin letters to make the input simpler.
Setupspacetimeindices = lowercaselatin
spacetimeindices=lowercaselatin
DefineAa, Ba, b, Fa, b, c, Ha, b, c, d, Ja, b, c, d, e
Aa,Ba,b,𝒟μ,γμ,Fa,b,c,Ha,b,c,d,Ja,b,c,d,e,σμ,Rμ,ν,Rμ,ν,α,β,Cμ,ν,α,β,Xμ,∂μ,ημ,ν,gμ,ν,lμ,Γμ,ν,α,Gμ,ν,δμ,ν,εα,β,μ,ν
Consider the following tensorial expressions.
e__1≔Ac⁢Hd,h,c,b⁢Fh,d,f⁢Ha,e,e,f−Ah⁢Ha,c,c,f⁢Fe,d,f⁢Hd,e,h,b
e__1≔Ha,eefa,eef⁢Hd,h,c,b⁢F⁢h,d,f⁢h,d,f⁢A⁢c⁢c−Ha,ccfa,ccf⁢Hd,e,h,b⁢F⁢e,d,f⁢e,d,f⁢A⁢h⁢h
e__2≔Jf,a,b,g,c⁢Fe,f,d⁢Fe,d,g−Je,a,b,f,c⁢Fd,e,g⁢Fd,g,f
e__2≔−Je,a,b,f,c⁢Fdegdeg⁢F⁢d,g,f⁢d,g,f+Jf,a,b,g,c⁢Fe,dge,dg⁢F⁢e,f,d⁢e,f,d
e__3≔Jg,c,k,j,e⁢Hi,h,a,h⁢He,b,i,g⁢Fd,k,j−Hh,b,g,j⁢Fd,k,e⁢Jj,c,k,e,h⁢Hg,i,a,i
e__3≔−Jj,c,k,e,h⁢Hg,i,aig,i,ai⁢H⁢hbg,j⁢hbg,j⁢Fdk,edk,e+Jg,c,k,j,e⁢H⁢eb,ig⁢eb,ig⁢H⁢ih,ah⁢ih,ah⁢Fdk,jdk,j
There are two, three and four free indices, and five, four and six repeated indices respectively in each of these three expressions, e__1,e__2,e__3 (see Check).
mapu → Checku,all,quiet, e__1,e__2,e__3
c,d,e,f,h,c,d,e,f,h,a,b,d,e,f,g,d,e,f,g,a,b,c,e,g,h,i,j,k,e,g,h,i,j,k,a,b,c,d
Simplifye__1,tryhard
Simplifye__2,tryhard
Simplifye__3
Improve the simplification with respect to side relations (equations) in the presence of quantum vectorial equations.
withVectors:
Setupredo, hermitianoperators = L,r_, p_, x,y,z, p, quantumoperators = A, B
hermitianoperators=L,p,p→,r→,x,y,z,quantumoperators=A,B,L,p,p→,r→,x,y,z
The identification of the vectorial character of expressions got improved. For example, this commutator is actually a non-projected vector.
Commutatorp_·p_,r_
p→2,r→−
Identify
5
Consider an expression involving vectorial quantum operators.
interfaceimaginaryunit = i
I
_i⋅ⅈ⁢x⁢pz−ⅈ⁢z⁢px+_j⋅ⅈ⁢y⁢pz−ⅈ⁢z⁢py
i∧⁢ⅈ⁢x⁢pz−ⅈ⁢z⁢px+j∧⁢ⅈ⁢y⁢pz−ⅈ⁢z⁢py
Simplify (23) taking into account the definition of angular momentum.
y⁢pz−z⁢py=Lx,z⁢px−x⁢pz=Ly
simplify,
−ⅈ⁢i∧⁢Ly+ⅈ⁢j∧⁢Lx
Enhancements in the simplification and integration of spherical harmonics (SphericalY) relevant in quantum mechanics.
18648640⁢7293⁢ⅇ6⁢I⁢theta⁢1−cos⁡phi23⁢20270252⁢cos⁡phi2−1351352π+SphericalY8,6,phi,theta
7293⁢ⅇ6⁢ⅈ⁢θ⁢1−cos⁡φ23⁢2027025⁢cos⁡φ22−13513528648640⁢π+Y86⁡φ,θ
simplify
conjugateSphericalY8,6,phi,theta*1+1/5*sin6*theta*sin5*phi*sinphi
Y86⁡φ,θ&conjugate0;⁢1+sin⁡6⁢θ⁢sin⁡5⁢φ5⁢sin⁡φ
IntInt,theta=0..2 Pi, phi=0..Pi
∫0π∫02⁢πY86⁡φ,θ&conjugate0;⁢1+sin⁡6⁢θ⁢sin⁡5⁢φ5⁢sin⁡φⅆθⅆφ
value
−37⁢ⅈ81920⁢π32⁢7293
A number of relevant changes happened in the tensor routines of the Physics package, towards making the routines pack more functionality, the simplification more powerful, and the handling of symmetries, substitutions, and other operations more flexible and natural.
Physics now works with four kinds of Minkowski spaces (different signatures) to accommodate the typical situations seen in textbooks; to these, correspond the signatures +---, ---+, -+++ and +++-.
Allow setting the metric by specifying the signature directly, as in g_[`-`] or g_[`+---`], or Setup⁡metric=`---+` or Setup⁡gμ,ν=`---+`.
The signature keyword of the Physics Setup is now in use, to set the metric and to indicate the form of the orthonormal tetrad, in turn used to derive the form of a null tetrad.
Automatic detection of the position of t as the time variable when you set the coordinates automatically sets the signature of the default Minkowski spacetime metric accordingly to ---+ or +---.
New keywords with special meaning when indexing the Physics (also the user defined) tensors: · `~`; for example g_[`~`] returns the all-contravariant matrix form of the metric. · definition; for example Ricci[definition] returns the definition of the Ricci tensor; works also with user-defined tensors. · scalars; for example Weyl[scalars] and Ricci[scalars] return the five Weyl and seven Ricci scalars used to perform a Petrov classification and in the Newman-Penrose formalism. · scalarsdefinition, and invariantsdefinition; for example Weyl[scalarsdefinition] or Riemann[invariantsdefinition] return the corresponding definitions for the scalars and invariants. · nullvectors; for example, when the new Tetrads subpackage is loaded, e_[nullvectors] returns a sequence of null vectors with their products normalized according to the Newman-Penrose formalism. · matrix; this keyword was introduced in previous releases, and in Maple 2015 it can appear after a space index (not spacetime), in which case a matrix with only the space components is returned.
Tensorial expressions can now have spacetime indices (related to a global system of references) and tetrad indices (related to a local system of references) at the same time, or they be rewritten in one (spacetime) or the other (tetrad) frames.
The matrix keyword can be used with spacetime, space, or tetrad indices, resulting in the corresponding matrix
Implement automatic determination of symmetry under permutation of tensor indices when the tensor is defined as a matrix.
New conversions from the Weyl to the Ricci tensors, and from Weyl to the Christoffel symbols.
New option evaluatetrace = true or false within convert/Ricci, to avoid automatically evaluating the Ricci trace when performing conversions that involve this trace.
New option 'evaluate' to convert/g_, convert/Christoffel and convert/Ricci. With this option set to false, it is possible to see the algebraic form of the result (that is, of the tensors involved) before evaluating it.
The Maple 18 Library:-SubstituteTensor command, got enhanced and transformed into one of the main Physics commands, that substitutes tensorial equation(s) Eqs into an expression, taking care of the free and repeated indices, such that: 1) equations in Eqs are interpreted as mappings having the free indices as parameters, 2) repeated indices in Eqs do not clash with repeated indices in the expression and 3) spacetime, space, and tetrad indices are handled independently, so they can all be present in Eqs and in the expression at the same time. This new command can also substitute algebraic sub-expressions of type product or sum within the expression, generalizing and unifying the functionality of the subs and algsubs commands for algebraic tensor expressions.
Minkowski spacetime: four different conventions.
restart; withPhysics:CoordinatesX
X
By default, the metric has a signature - - - +, with time in the last position.
x0
x4
You can now set the signature of a Minkowski spacetime directly from the metric, or using the Setup keyword signature to use any of the four conventions frequently found in textbooks (- - - +), (+ - - -), (+++ -), and (- +++), and to indicate an Euclidean spacetime you can indicate (++++) or just +.
g_`+`
Changing the signature of the tensor spacetime to: + + + +
The Euclidean metric in coordinates x1,x2,x3,x4
g_`+++-`
The Minkowski metric, with signature + + + -, in coordinates x1,x2,x3,x4
Using Setup's signature keyword:
Setupsignature = `- +++`
signature=- + + +
The position of the different sign in a Minkowski spacetime always, either the first or the last symbol in the signature, refers to the time variable. So if you indicate t as one of the coordinates, say in the first position, and the signature is not in agreement with time in the first position, it is now automatically corrected and the metric set. For example, in this moment, due to (33), time is implicitly assumed to be in the first position; set coordinates with t in the last (fourth) position.
CoordinatesX=x,y,z,t
Detected `t`, the time variable, in position 4. Changing the signature of the spacetime metric accordingly, to: + + + -
Default differentiation variables for d_, D_ and dAlembertian are: X=x,y,z,t
Systems of spacetime Coordinates are: X=x,y,z,t
New shortcut for the matrix form of the all contravariant components of a tensor.
Since previous releases, you can always request the matrix (or array in the case of more than 2 indices) form of a tensor by adding the keyword matrix as a last index. For example, g_μ,ν,matrix . A shortcut notation for the all covariant components is to omit the indices, as in . A new shortcut notation for the all contravariant components is to pass only they tilde `~`. To see the difference, consider a non-Minkowski spacetime, for instance set it in one go via:
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
The all contravariant components of the metric are:
g_`~`
It is now possible to also use the keyword matrix passing only space (not spacetime) indices, resulting in a matrix with only the space components.
Setupspaceindices=lowercase_is
spaceindices=lowercaselatin_is
g_~i,~j, matrix
g⁢i,j⁢i,j=−r+2⁢mr000−1r2000−1r2⁢sin⁡θ2
All the tensor predefined in the Physics package, or defined using defining equations and the Define command, accept a new keyword, definition, that returns their definition. For example:
Riccidefinition
Rμ,ν=∂α⁡Γ⁢αμ,ν⁢αμ,ν−∂ν⁡Γ⁢αμ,α⁢αμ,α+Γ⁢βμ,ν⁢βμ,ν⁢Γ⁢αβ,α⁢αβ,α−Γ⁢βμ,α⁢βμ,α⁢Γ⁢αν,β⁢αν,β
Einsteindefinition
Gμ,ν=Rμ,ν−gμ,ν⁢Rαααα2
Define a tensor using a tensorial equation:
PDEtools:-declareAX
A⁡r,θ,φ,t⁢will now be displayed as⁢A
Fmu,nu = d_muAnuX − d_nuAmuX
Fμ,ν=∂μ⁡Aν−∂ν⁡Aμ
Note that for consistency of free indices, if you now define the left-hand side Fμ,ν as a tensor, in the right-hand side Aμ is also a tensor and so is automatically defined as well.
Aμ,𝒟μ,γμ,Fμ,ν,σμ,Rμ,ν,Rμ,ν,α,β,Cμ,ν,α,β,Xμ,∂μ,gμ,ν,Γμ,ν,α,Gμ,ν,δμ,ν,εα,β,μ,ν
Check the definition of Fμ,ν:
Fdefinition
A new keyword scalars for the Weyl and Ricci tensors, generates the Weyl and Ricci scalars of the Newman-Penrose formalism. Correspondingly, another new keyword scalarsdefinition return the definition of these scalars.
Weylscalars
ψ__0=0,ψ__1=0,ψ__2=mr3,ψ__3=0,ψ__4=0
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;β
In the definition earlier, the tensors lμ,nμ,mμ,m&conjugate0;μ are the null tensors of the Newman-Penrose formalism, implemented in Maple 2015 within the new Tetrads package.
withTetrads
Setting lowercaselatin_ah 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,SimplifyTetrad,TransformTetrad,e_,eta_,gamma_,l_,lambda_,m_,mb_,n_
In this new Tetrads package, 𝔢a,μ is the tetrad (vierbein), by default an orthonormal tetrad that can be set to be a null tetrad using Setup⁡tetrad=null, γa,b,c,λa,b,c are respectively the Ricci rotation coefficients and the lambda tensor defined in the "Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2" (definitions (98.9) and (98.10))
For instance, the form of the orthonormal tetrad for the Schwarzschild metric set in is:
e_
𝔢a,μ=r−r+2⁢m0000−I⁢r0000−I⁢r⁢sin⁡θ0000−r+2⁢mr
A new keyword for the tetrad 𝔢a,μ is nullvectors:
e_nullvectors
lμ=2⁢r2⁢−r+2⁢m002⁢−r+2⁢m2⁢r,nμ=−2⁢r2⁢−r+2⁢m002⁢−r+2⁢m2⁢r,mμ=0−I2⁢2⁢r2⁢r⁢sin⁡θ20,m&conjugate0;μ=0−I2⁢2⁢r−2⁢r⁢sin⁡θ20
The matrix keyword can be used with different kinds of indices, representing different objects. For example, the spacetime components of the electromagnetic tensor defined in equation (40) are:
F
Fμ,ν=0A2r−A1θA3r−A1φA4r−A1tA1θ−A2r0A3θ−A2φA4θ−A2tA1φ−A3rA2φ−A3θ0A4φ−A3tA1t−A4rA2t−A4θA3t−A4φ0
The components of this tensor in the local inertial (tetrad) system of references are:
Fa,b,matrix
Fa,b=0I⁢−r+2⁢m⁢A2r−A1θr32I⁢−r+2⁢m⁢A3r−A1φr32⁢sin⁡θA1t−A4r−I⁢−r+2⁢m⁢A2r−A1θr320A2φ−A3θr2⁢sin⁡θI⁢A2t−A4θ−r+2⁢m⁢r−I⁢−r+2⁢m⁢A3r−A1φr32⁢sin⁡θA3θ−A2φr2⁢sin⁡θ0I⁢A3t−A4φ−r+2⁢m⁢r⁢sin⁡θA4r−A1tI⁢A4θ−A2tr⁢−r+2⁢mI⁢A4φ−A3tr⁢sin⁡θ⁢−r+2⁢m0
Likewise, the nonzero components of the Riemann tensor in the global (spacetime, Greek indices) and local (tetrad, lowercase roman indices from a to h) system of references are:
Riemannnonzero
Rμ,ν,α,β=1,2,1,2=−m−r+2⁢m,1,2,2,1=m−r+2⁢m,1,3,1,3=−m⁢sin⁡θ2−r+2⁢m,1,3,3,1=m⁢sin⁡θ2−r+2⁢m,1,4,1,4=2⁢mr3,1,4,4,1=−2⁢mr3,2,1,1,2=m−r+2⁢m,2,1,2,1=−m−r+2⁢m,2,3,2,3=−2⁢r⁢m⁢sin⁡θ2,2,3,3,2=2⁢r⁢m⁢sin⁡θ2,2,4,2,4=−r+2⁢m⁢mr2,2,4,4,2=−2⁢m2+m⁢rr2,3,1,1,3=m⁢sin⁡θ2−r+2⁢m,3,1,3,1=−m⁢sin⁡θ2−r+2⁢m,3,2,2,3=2⁢r⁢m⁢sin⁡θ2,3,2,3,2=−2⁢r⁢m⁢sin⁡θ2,3,4,3,4=sin⁡θ2⁢−r+2⁢m⁢mr2,3,4,4,3=−sin⁡θ2⁢−r+2⁢m⁢mr2,4,1,1,4=−2⁢mr3,4,1,4,1=2⁢mr3,4,2,2,4=−2⁢m2+m⁢rr2,4,2,4,2=−r+2⁢m⁢mr2,4,3,3,4=−sin⁡θ2⁢−r+2⁢m⁢mr2,4,3,4,3=sin⁡θ2⁢−r+2⁢m⁢mr2
Riemanna,b,c,d,nonzero
Ra,b,c,d=1,2,1,2=mr3,1,2,2,1=−mr3,1,3,1,3=mr3,1,3,3,1=−mr3,1,4,1,4=2⁢mr3,1,4,4,1=−2⁢mr3,2,1,1,2=−mr3,2,1,2,1=mr3,2,3,2,3=−2⁢mr3,2,3,3,2=2⁢mr3,2,4,2,4=−mr3,2,4,4,2=mr3,3,1,1,3=−mr3,3,1,3,1=mr3,3,2,2,3=2⁢mr3,3,2,3,2=−2⁢mr3,3,4,3,4=−mr3,3,4,4,3=mr3,4,1,1,4=−2⁢mr3,4,1,4,1=2⁢mr3,4,2,2,4=mr3,4,2,4,2=−mr3,4,3,3,4=mr3,4,3,4,3=−mr3
Define a tensor using a symmetric matrix on the right-hand-side:
Mmu, nu = Matrix4, symbol=m, shape=symmetric
Aμ,𝒟μ,γμ,Fμ,ν,Mμ,ν,σμ,Rμ,ν,Rμ,ν,α,β,Cμ,ν,α,β,Xμ,∂μ,𝔢a,μ,ηa,b,gμ,ν,γa,b,c,lμ,λa,b,c,mμ,m&conjugate0;μ,nμ,Γμ,ν,α,Gμ,ν,δμ,ν,εα,β,μ,ν
Check whether the system detected that ηα,β is symmetric under permutation of α,β.
Library:−IsTensorialSymmetric⁡Mα,β
true
Hence, the indices of this tensor are automatically normalized taking this symmetry into account, so that:
Mα,β
Mβ,α
−
The formalism of tetrads in general relativity got implemented within Physics as a new package, Physics:-Tetrads, with 13 commands, mainly the null vectors of the Newman-Penrose formalism, the tetrad tensors 𝔢a,μ,ηa,b,γa,b,c,λa,b,c, respectively: the tetrad, the tetrad metric, the Ricci rotation coefficients, and the lambda tensor, plus five algebraic manipulation commands: IsTetrad, NullTetrad, OrthonormalTetrad, SimplifyTetrad, and TransformTetrad to construct orthonormal and null tetrads of different forms and using different methods.
The new Tetrads package contains 13 commands for computing in a local (tetrad) frame. Tensor components in the local frame are represented with tetrad indices using a type of letter different than the one representing global spacetime indices. You can set the type of letter using Setup, or just load the new package and the type of letter will be set automatically.
restart;withPhysics:withTetrads;
Setting lowercaselatin letters to represent tetrad indices
The most relevant commands are 𝔢a,μ and ηa,b, respectively representing the tetrad (also vierbein; by default, this is an orthonormal tetrad) and the tetrad metric (that is, the metric of the local - by default inertial - frame)
e_definition
𝔢a,μ⁢𝔢bμbμ=ηa,b
Then γa,b,c are the Ricci rotation coefficients and λa,b,c is a linear combination of them, according to the definitions in the "Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2" (definitions (98.9) and (98.10)).
gamma_definition
γa,b,c=𝒟ν⁡𝔢a,μ⁢𝔢bμbμ⁢𝔢cνcν
lambda_definition
λa,b,c=𝒟ν⁡𝔢a,μ−𝒟μ⁡𝔢a,ν⁢𝔢bμbμ⁢𝔢cνcν
The lμ,nμ,mμ,m&conjugate0;μ tensors are the null tensors of the Newman-Penrose formalism and the commands IsTetrad,NullTetrad,OrthonormalTetrad,SimplifyTetrad,TransformTetrad are for manipulating and exploring different forms of tetrads.
In a flat space, the spacetime and tetrad metrics are the same, so the orthonormal tetrad is just the identity.
eta_
𝔢a,μ=1000010000100001
In a curved spacetime, for instance, set a "Local Rotational Symmetry metric" metric:
g_13,7,5
The metric in coordinates x,y,z,t
Parameters: ε,A⁡t,B⁡t
Comments: _ⅇpsⅈlon=1 or _ⅇpsⅈlon=-1
PDEtools:-declareAt, Bt
A⁡t⁢will now be displayed as⁢A
B⁡t⁢will now be displayed as⁢B
The default orthonormal tetrad is:
𝔢a,μ=−I⁢ε⁢A0000−I⁢cosh⁡2⁢x⁢B⁢ⅇA⁢xI⁢sinh⁡2⁢x⁢B⁢ⅇA⁢xcosh⁡2⁢x000−I⁢B⁢ⅇA⁢xcosh⁡2⁢x0000I⁢ε
The following null vectors correspond to this tetrad:
lμ=−I2⁢2⁢ε⁢A00I2⁢2⁢ε,nμ=I2⁢2⁢ε⁢A00I2⁢2⁢ε,mμ=0−I2⁢2⁢cosh⁡2⁢x⁢B⁢ⅇA⁢x2⁢B⁢ⅇA⁢x⁢I⁢sinh⁡2⁢x+12⁢cosh⁡2⁢x0,m&conjugate0;μ=0−I2⁢2⁢cosh⁡2⁢x⁢B⁢ⅇA⁢x2⁢B⁢ⅇA⁢x⁢I⁢sinh⁡2⁢x−12⁢cosh⁡2⁢x0
You can compute these null vectors directly since these are also part of the Tetrads package:
l_mu2,l_mu n_mu,l_mu m_mu,l_mu mb_mu
lμ⁢l⁢μ⁢μ,lμ⁢n⁢μ⁢μ,lμ⁢m⁢μ⁢μ,m&conjugate0;μ⁢l⁢μ⁢μ
mapu → u = SumOverRepeatedIndicesu,
lμ⁢l⁢μ⁢μ=0,lμ⁢n⁢μ⁢μ=1,lμ⁢m⁢μ⁢μ=0,m&conjugate0;μ⁢l⁢μ⁢μ=0
You can query about their definition in the same way you can now query any other tensor:
m_definition
mμ⁢l⁢μ⁢μ=0,mμ⁢n⁢μ⁢μ=0,mμ⁢m⁢μ⁢μ=0,mμ⁢m&conjugate0;⁢μ⁢μ=−1,gμ,ν=lμ⁢nν+lν⁢nμ−mμ⁢m&conjugate0;ν−mν⁢m&conjugate0;μ
Now that this definition depends on the signature:
Setupsignature
signature=- - - +
Change the signature, both in sign and placing time in position 1:
Setupsignature=`- +++`
So, now mμ⁢m&conjugate0;⁢μ⁢μ=1 instead of −1 and there is a change in the sign of the definition in terms of the metric gμ,ν if compared with (67).
mμ⁢l⁢μ⁢μ=0,mμ⁢n⁢μ⁢μ=0,mμ⁢m⁢μ⁢μ=0,mμ⁢m&conjugate0;⁢μ⁢μ=1,gμ,ν=−lμ⁢nν−lν⁢nμ+mμ⁢m&conjugate0;ν+mν⁢m&conjugate0;μ
You can verify these tensorial identities using TensorArray; for example, for the last equation:
−1
gμ,ν=−lμ⁢nν−lν⁢nμ+mμ⁢m&conjugate0;ν+mν⁢m&conjugate0;μ
TensorArray,simplifier=simplify@expand
ε⁢A2=ε⁢A20=00=00=00=0B2⁢ⅇ2⁢A⁢x⁢2⁢cosh⁡x2−1=B2⁢ⅇ2⁢A⁢x⁢2⁢cosh⁡x2−1−2⁢B2⁢ⅇ2⁢A⁢x⁢cosh⁡x⁢sinh⁡x=−2⁢B2⁢ⅇ2⁢A⁢x⁢cosh⁡x⁢sinh⁡x0=00=0−2⁢B2⁢ⅇ2⁢A⁢x⁢cosh⁡x⁢sinh⁡x=−2⁢B2⁢ⅇ2⁢A⁢x⁢cosh⁡x⁢sinh⁡xB2⁢ⅇ2⁢A⁢x⁢2⁢cosh⁡x2−1=B2⁢ⅇ2⁢A⁢x⁢2⁢cosh⁡x2−10=00=00=00=0−ε=−ε
To compute with a null tetrad instead of an orthonormal tetrad, set the tetrad or tetradmetric to null using Setup.
Setuptetradmetric = null
tetradmetric=1,2=−1,3,4=1
So, now the metric of the local (tetrad) system of references is:
Note: this form of the null tetrad metric is consistent with time in position 1, a change done in (69)- and not with time in position 4 (default). You can in any case redefine the tetrad metric in any particular way also using Setup⁡tetradmetric=...a symmetric matrix here....
The null tetrad is now:
𝔢a,μ=−I2⁢2⁢ε⁢A2⁢cosh⁡2⁢x⁢B⁢ⅇA⁢x2−2⁢sinh⁡2⁢x⁢B⁢ⅇA⁢x2⁢cosh⁡2⁢x0−I2⁢2⁢ε⁢A−2⁢cosh⁡2⁢x⁢B⁢ⅇA⁢x22⁢sinh⁡2⁢x⁢B⁢ⅇA⁢x2⁢cosh⁡2⁢x0002⁢B⁢ⅇA⁢x2⁢cosh⁡2⁢x−2⁢ε2002⁢B⁢ⅇA⁢x2⁢cosh⁡2⁢x2⁢ε2
Compare this result with the orthonormal tetrad (63).
It is possible to test whether these tetrads satisfy the tetrad definition:
Using the manipulation commands of the package:
IsTetrad