Physics[Christoffel] - The Christoffel symbols
Calling Sequence
Christoffel[μ,α,β]
Christoffel[keyword]
convert(expression, Christoffel)
Parameters
μ,α,β
-
the indices, as names representing integer numbers between 0 and the spacetime dimension, they can also be the numbers themselves
keyword
optional, it can be definition, array, nonzero
Description
The all covariant Christoffel[mu, alpha, beta], displayed as Γ__μ,α,β, is a computational representation for the Christoffel symbols, defined in terms of derivatives of the spacetime metric as
Γμ,α,β=⁢12∂β⁡gα,μ+∂α⁡gβ,μ−∂μ⁡gα,β
where the spacetime differentiation operator ∂__mu is implemented in Maple as d_[mu]. From this definition, Christoffel is symmetric with respect to interchanging the second and third indices.
Christoffel can also have contravariant indices, typically the first one. When any of the indices of Christoffel is contravariant, the definition for its components is according to Γ⁢⁢α,βμ=gμ,ρ⁢Γρ,α,β,⁢⁢⁢Γμ⁢β⁢α=gα,ρ⁢Γμ,ρ,β and Γμ,α⁢⁢⁢⁢⁢⁢⁢β=gβ,ρ⁢Γμ,α,ρ and not by raising the indices in the terms on the right-hand side of the definition above for Γ__μ,α,β. Note also that, despite computing with them as tensors, the Christoffel symbols are not actual tensors in that, when you transform the coordinates, they do not transform like the coordinate differentials (regarding their contravariant indices) nor like the derivatives of a scalar (regarding their covariant indices).
To express Γ__μ,α,β or Γ⁢⁢α,βμ using this definition in terms of derivatives of the spacetime metric use convert to g_. Sometimes it is also convenient to rewrite tensorial expressions the other way around, in terms of the Christoffel symbols and its derivatives. For example: non-covariant derivatives of the metric or any of the general relativity tensors that can be expressed in terms of them. For this purpose you can use convert(expression, Christoffel) - see the Examples section.
When the indices of Christoffel assume integer values they are expected to be between 0 and the spacetime dimension, prefixed by ~ when they are contravariant, and the corresponding value of Christoffel is returned. The values 0 and 4, or for the case any dimension set for the spacetime, represent the same object. When the indices have symbolic values Christoffel returns unevaluated after normalizing its indices taking into account the symmetry property of the 2nd and 3rd indices.
Computations performed with the Physics package commands take into account Einstein's sum rule for repeated indices - see `.` and Simplify. The distinction between covariant and contravariant indices in the input of tensors is done by prefixing contravariant ones with ~, say as in ~mu; in the output, contravariant indices are displayed as superscripts. For contracted indices, you can enter them one covariant and one contravariant. Note however that - provided that the spacetime metric is galilean (Euclidean or Minkowski), or the object is a tensor also in curvilinear coordinates - this distinction in the input is not relevant, and so contracted indices can be entered as both covariant or both contravariant, in which case they will be automatically rewritten as one covariant and one contravariant. Tensors can have spacetime and space indices at the same time. To change the type of letter used to represent spacetime or space indices see Setup.
During a Maple session, the value of any component of Γ__μ,α,β is automatically determined by the value of the spacetime metric at the moment. When Physics is loaded, the spacetime is set to Minkowski type, and so all the elements of Christoffel are automatically zero. To set the spacetime metric to something different use Setup. Also, at least one system of coordinates must be set in order to compute the derivatives entering the definition of Christoffel. For that purpose see Coordinates or Setup.
Besides being indexed with three indices, Christoffel accepts three keywords:
definition: returns the definition of the Christoffel symbols in terms of derivatives of the metric g_.
array: (synonym: Array, or no indices whatsoever, as in Christoffel[]) returns an Array that when indexed with numerical values from 1 to the dimension of spacetime it returns the value of each of the components of Christoffel. If this keyword is passed together with indices, that can be covariant or contravariant, the resulting array takes into account the character of the indices. If one of these indices has a numerical value, if contravariant then preceded by ~, the returned object is the corresponding 2 x 2 Matrix.
nonzero: returns a set of equations, with the left-hand side as a sequence of three positive numbers identifying the element of Γ__μ,α,β and the corresponding value on the right-hand side. Note that this set is actually the output of the ArrayElems command when passing to it the Array obtained with the keyword array.
Some automatic checking and normalization are carried out each time you enter Christoffel[...]. The checking is concerned with possible syntax errors. The automatic normalization takes into account the symmetry of Christoffel[mu, alpha, beta] with respect to interchanging the positions of the indices alpha and beta.
The %Christoffel command is the inert form of Christoffel, so it represents the same mathematical operation but without performing it. To perform the operation, use value.
Examples
with⁡Physics:
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_, the covariant derivative D_ and the dAlembertian
Setup⁡coordinatesystems=cartesian
⁢Systems of spacetime coordinates are:⁢X=x,y,z,t
coordinatesystems=X
When Physics is initialized, the default spacetime metric is of Minkowski type. You can see the metric by querying Setup, as in Setup(metric);, or by directly entering the metric g_ with no indices
g_
gμ,ν=
Check the nonzero components of Christoffel: because the default spacetime is of Minkowski type, there are none
Christoffelnonzero
Γα,μ,ν=∅
Hence
Christoffelμ,α,β
0
Consequently, all the general relativity tensors defined in terms of Christoffel and derivatives of the metric g_ are also automatically zero
Einsteinμ,ν
Ricciμ,ν
Riemannμ,ν,ρ,σ
Weylμ,ν,ρ,σ
For the same reason the covariant derivative ▿__mu=∂__mu; for illustration purposes Define an arbitrary spacetime tensor A
Define⁡A
Defined objects with tensor properties
A,γμ,σμ,∂μ,gμ,ν,εα,β,μ,ν,Xμ
D_μ⁡A~nu⁡X
∂μ⁡A⁢ν⁢ν⁡X
To set the scenario to be an arbitrary non-galilean spacetime, set the metric to be arbitrary (that is, depending on 10 arbitrary functions of the spacetime coordinates X)
Setup⁡metric=arbitrary
⁢_______________________________________________________
Setting lowercaselatin_is letters to represent space indices
⁢The arbitrary metric in coordinates ⁢x,y,z,t
Signature: - - - +
metric=1,1=f__1⁡X,1,2=f__2⁡X,1,3=f__3⁡X,1,4=f__4⁡X,2,2=f__5⁡X,2,3=f__6⁡X,2,4=f__7⁡X,3,3=f__8⁡X,3,4=f__9⁡X,4,4=f__10⁡X,spaceindices=lowercaselatin_is
When the indices are not numerical, Christoffel returns itself after normalizing its second and third indices taking advantage of the symmetry under exchange of them
Γμ,α,β
Christoffelμ,β,α
−
To express Christoffel in terms of its definition in terms of the spacetime metric g_, use convert to g_
convert⁡,g_
∂β⁡gα,μ2+∂α⁡gβ,μ2−∂μ⁡gα,β2
To rewrite expressions like this one back in terms of Christoffel and its derivatives use convert(... Christoffel)
convert⁡,Christoffel
Riemann~alpha,β,μ,ν
R⁢αβ,μ,ν⁢αβ,μ,ν
∂μ⁡Γ⁢αβ,ν⁢αβ,ν−∂ν⁡Γ⁢αβ,μ⁢αβ,μ+Γ⁢αμ,υ⁢αμ,υ⁢Γ⁢υβ,ν⁢υβ,ν−Γ⁢αν,υ⁢αν,υ⁢Γ⁢υβ,μ⁢υβ,μ
Christoffel~mu,α,β
Γ⁢μα,β⁢μα,β
g⁢μ,ν⁢μ,ν⁢∂β⁡gα,ν+∂α⁡gβ,ν−∂ν⁡gα,β2
g⁢μ,ν⁢μ,ν⁢Γν,α,β
Simplify⁡
Related to displaying the components of Christoffel, it is convenient to have a compact display of the 10 arbitrary functions entering the metric as well as of their derivatives; for that purpose use the enhanced display scheme of the differential equation packages.
PDEtools:-declare⁡
f__1⁡x,y,z,t⁢will now be displayed as⁢f__1
f__10⁡x,y,z,t⁢will now be displayed as⁢f__10
f__2⁡x,y,z,t⁢will now be displayed as⁢f__2
f__3⁡x,y,z,t⁢will now be displayed as⁢f__3
f__4⁡x,y,z,t⁢will now be displayed as⁢f__4
f__5⁡x,y,z,t⁢will now be displayed as⁢f__5
f__6⁡x,y,z,t⁢will now be displayed as⁢f__6
f__7⁡x,y,z,t⁢will now be displayed as⁢f__7
f__8⁡x,y,z,t⁢will now be displayed as⁢f__8
f__9⁡x,y,z,t⁢will now be displayed as⁢f__9
So from herein, in the display of expressions, derivatives are displayed as indexed and the functionality of the arbitrary functions is suppressed.
Check the value of Γ__α,μ,ν for α = μ = ν = 1
Christoffel1,1,1
f__1x2
Now of Γα,βμ for α = μ = ν = 1 (note you enter the contravariant index prefixed by ~)
Christoffel~1,1,1
2⁢f__4x⁢f__8⁢f__7⁢f__2−2⁢f__4x⁢f__9⁢f__6⁢f__2−2⁢f__4x⁢f__6⁢f__7⁢f__3+f__1y⁢f__8⁢f__10⁢f__2−f__1y⁢f__8⁢f__7⁢f__4−f__1y⁢f__10⁢f__6⁢f__3+f__1y⁢f__9⁢f__6⁢f__4+f__1y⁢f__9⁢f__7⁢f__3+f1.⁢f__5⁢f__8⁢f__4−f1.⁢f__5⁢f__9⁢f__3−f1.⁢f__8⁢f__7⁢f__2+f1.⁢f__9⁢f__6⁢f__2+f1.⁢f__6⁢f__7⁢f__3+f__1z⁢f__5⁢f__10⁢f__3−f__1z⁢f__5⁢f__9⁢f__4−f__1z⁢f__10⁢f__6⁢f__2+f__1z⁢f__9⁢f__7⁢f__2+f__1z⁢f__6⁢f__7⁢f__4−2⁢f__2x⁢f__8⁢f__10⁢f__2+2⁢f__2x⁢f__8⁢f__7⁢f__4+2⁢f__2x⁢f__10⁢f__6⁢f__3−2⁢f__2x⁢f__9⁢f__6⁢f__4−2⁢f__2x⁢f__9⁢f__7⁢f__3−2⁢f__3x⁢f__5⁢f__10⁢f__3+2⁢f__3x⁢f__5⁢f__9⁢f__4+2⁢f__3x⁢f__10⁢f__6⁢f__2−2⁢f__3x⁢f__9⁢f__7⁢f__2−2⁢f__3x⁢f__6⁢f__7⁢f__4−2⁢f__4x⁢f__5⁢f__8⁢f__4+2⁢f__4x⁢f__5⁢f__9⁢f__3+f__1x⁢f__5⁢f__8⁢f__10+2⁢f__1x⁢f__9⁢f__6⁢f__7−f1.⁢f__62⁢f__4−f__1x⁢f__5⁢f__92−f__1x⁢f__8⁢f__72−f__1x⁢f__10⁢f__62+2⁢f__2x⁢f__92⁢f__2+2⁢f__3x⁢f__72⁢f__3+2⁢f__4x⁢f__62⁢f__4−f__1z⁢f__72⁢f__3−f__1y⁢f__92⁢f__22⁢f__5⁢f__8⁢f__10⁢f__1−f__5⁢f__8⁢f__42−f__5⁢f__10⁢f__32−f__5⁢f__92⁢f__1+2⁢f__5⁢f__9⁢f__3⁢f__4−f__8⁢f__10⁢f__22−f__8⁢f__72⁢f__1+2⁢f__8⁢f__7⁢f__2⁢f__4−f__10⁢f__62⁢f__1+2⁢f__10⁢f__6⁢f__2⁢f__3+f__92⁢f__22+2⁢f__9⁢f__6⁢f__7⁢f__1−2⁢f__9⁢f__6⁢f__2⁢f__4−2⁢f__9⁢f__7⁢f__2⁢f__3+f__62⁢f__42−2⁢f__6⁢f__7⁢f__3⁢f__4+f__72⁢f__32
To compute with a representation for Christoffel without actually performing the operation, use the inert form %Christoffel. To afterwards perform the operation use value
%Christoffel1,1,1
Γ1,1,1
value⁡
Compute Arrays with all the components of Γ__μ,α,β and Γα,βμ
Christoffelμ,α,β,array
Christoffel~mu,α,β,array
Because these arrays of values are computed by already taking into account the character (covariant/contravariant) of each index, you do not need to indicate furthermore that character. Thus, this is the all covariant Γ__α,μ,ν with α = μ = ν = 1
C11,1,1
and this is the value of Γα,βμ for α = μ = ν = 1
C21,1,1
Compare with the output of Γ1,11 in eq (27)
normal⁡−
The Schwarzschild metric in spherical coordinates (see g_) and the four 2 x 2 matrices for Γ⁢⁢α,βμ by giving values to its first index
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
seq⁡Christoffelj,α,β,matrix,j=~1,~2,~3,~4
Γ⁢1α,β⁢1α,β=?,Γ⁢2α,β⁢2α,β=?,Γ⁢3α,β⁢3α,β=?,Γ⁢4α,β⁢4α,β=?
See Also
`.`, Array, ArrayElems, Coordinates, D_, d_, Define, DifferentialGeometry[Tensor][Christoffel], Einstein, PDEtools[declare], Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Ricci, Riemann, Setup, value, Weyl
References
Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2, fourth revised English edition. Elsevier, 1975.
Compatibility
The Physics[Christoffel] command was introduced in Maple 16.
For more information on Maple 16 changes, see Updates in Maple 16.
The Physics[Christoffel] command was updated in Maple 2017.
Download Help Document