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

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Physics : Christoffel

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

withPhysics:

Setupmathematicalnotation=true

mathematicalnotation=true

(1)

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

Setupcoordinatesystems=cartesian

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

coordinatesystems=X

(2)

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μ,ν=

(3)

Check the nonzero components of Christoffel: because the default spacetime is of Minkowski type, there are none

Christoffelnonzero

Γα,μ,ν=

(4)

Hence

Christoffelμ,α,β

0

(5)

Consequently, all the general relativity tensors defined in terms of Christoffel and derivatives of the metric g_ are also automatically zero

Einsteinμ,ν

0

(6)

Ricciμ,ν

0

(7)

Riemannμ,ν,ρ,σ

0

(8)

Weylμ,ν,ρ,σ

0

(9)

For the same reason the covariant derivative ▿__mu=∂__mu; for illustration purposes Define an arbitrary spacetime tensor A

DefineA

Defined objects with tensor properties

A,γμ,σμ,μ,gμ,ν,εα,β,μ,ν,Xμ

(10)

D_μA~nuX

μAννX

(11)

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)

Setupmetric=arbitrary

_______________________________________________________

Setting lowercaselatin_is letters to represent space indices

The arbitrary metric in coordinates x,y,z,t

Signature: - - - +

_______________________________________________________

gμ,ν=

metric=1,1=f__1X,1,2=f__2X,1,3=f__3X,1,4=f__4X,2,2=f__5X,2,3=f__6X,2,4=f__7X,3,3=f__8X,3,4=f__9X,4,4=f__10X,spaceindices=lowercaselatin_is

(12)

g_

gμ,ν=

(13)

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

Γμ,α,β

(14)

Christoffelμ,β,α

Γμ,α,β

(15)

0

(16)

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

(17)

To rewrite expressions like this one back in terms of Christoffel and its derivatives use convert(... Christoffel)

convert,Christoffel

Γμ,α,β

(18)

Riemann~alpha,β,μ,ν

Rαβ,μ,ναβ,μ,ν

(19)

convert,Christoffel

μΓαβ,ναβ,ννΓαβ,μαβ,μ+Γαμ,υαμ,υΓυβ,νυβ,νΓαν,υαν,υΓυβ,μυβ,μ

(20)

Christoffel~mu,α,β

Γμα,βμα,β

(21)

convert,g_

gμ,νμ,νβgα,ν+αgβ,ννgα,β2

(22)

convert,Christoffel

gμ,νμ,νΓν,α,β

(23)

Simplify

Γμα,βμα,β

(24)

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__1x,y,z,twill now be displayed asf__1

f__10x,y,z,twill now be displayed asf__10

f__2x,y,z,twill now be displayed asf__2

f__3x,y,z,twill now be displayed asf__3

f__4x,y,z,twill now be displayed asf__4

f__5x,y,z,twill now be displayed asf__5

f__6x,y,z,twill now be displayed asf__6

f__7x,y,z,twill now be displayed asf__7

f__8x,y,z,twill now be displayed asf__8

f__9x,y,z,twill now be displayed asf__9

(25)

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

(26)

Now of Γα,βμ for α = μ = ν = 1 (note you enter the contravariant index prefixed by ~)

Christoffel~1,1,1

2f__4xf__8f__7f__22f__4xf__9f__6f__22f__4xf__6f__7f__3+f__1yf__8f__10f__2f__1yf__8f__7f__4f__1yf__10f__6f__3+f__1yf__9f__6f__4+f__1yf__9f__7f__3+f1.f__5f__8f__4f1.f__5f__9f__3f1.f__8f__7f__2+f1.f__9f__6f__2+f1.f__6f__7f__3+f__1zf__5f__10f__3f__1zf__5f__9f__4f__1zf__10f__6f__2+f__1zf__9f__7f__2+f__1zf__6f__7f__42f__2xf__8f__10f__2+2f__2xf__8f__7f__4+2f__2xf__10f__6f__32f__2xf__9f__6f__42f__2xf__9f__7f__32f__3xf__5f__10f__3+2f__3xf__5f__9f__4+2f__3xf__10f__6f__22f__3xf__9f__7f__22f__3xf__6f__7f__42f__4xf__5f__8f__4+2f__4xf__5f__9f__3+f__1xf__5f__8f__10+2f__1xf__9f__6f__7f1.f__62f__4f__1xf__5f__92f__1xf__8f__72f__1xf__10f__62+2f__2xf__92f__2+2f__3xf__72f__3+2f__4xf__62f__4f__1zf__72f__3f__1yf__92f__22f__5f__8f__10f__1f__5f__8f__42f__5f__10f__32f__5f__92f__1+2f__5f__9f__3f__4f__8f__10f__22f__8f__72f__1+2f__8f__7f__2f__4f__10f__62f__1+2f__10f__6f__2f__3+f__92f__22+2f__9f__6f__7f__12f__9f__6f__2f__42f__9f__7f__2f__3+f__62f__422f__6f__7f__3f__4+f__72f__32

(27)

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

(28)

value

f__1x2

(29)

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

f__1x2

(30)

and this is the value of Γα,βμ for α = μ = ν = 1

C21,1,1

2f__4xf__8f__7f__22f__4xf__9f__6f__22f__4xf__6f__7f__3+f__1yf__8f__10f__2f__1yf__8f__7f__4f__1yf__10f__6f__3+f__1yf__9f__6f__4+f__1yf__9f__7f__3+f1.f__5f__8f__4f1.f__5f__9f__3f1.f__8f__7f__2+f1.f__9f__6f__2+f1.f__6f__7f__3+f__1zf__5f__10f__3f__1zf__5f__9f__4f__1zf__10f__6f__2+f__1zf__9f__7f__2+f__1zf__6f__7f__42f__2xf__8f__10f__2+2f__2xf__8f__7f__4+2f__2xf__10f__6f__32f__2xf__9f__6f__42f__2xf__9f__7f__32f__3xf__5f__10f__3+2f__3xf__5f__9f__4+2f__3xf__10f__6f__22f__3xf__9f__7f__22f__3xf__6f__7f__42f__4xf__5f__8f__4+2f__4xf__5f__9f__3+f__1xf__5f__8f__10+2f__1xf__9f__6f__7f1.f__62f__4f__1xf__5f__92f__1xf__8f__72f__1xf__10f__62+2f__2xf__92f__2+2f__3xf__72f__3+2f__4xf__62f__4f__1zf__72f__3f__1yf__92f__22f__5f__8f__10f__1f__5f__8f__42f__5f__10f__32f__5f__92f__1+2f__5f__9f__3f__4f__8f__10f__22f__8f__72f__1+2f__8f__7f__2f__4f__10f__62f__1+2f__10f__6f__2f__3+f__92f__22+2f__9f__6f__7f__12f__9f__6f__2f__42f__9f__7f__2f__3+f__62f__422f__6f__7f__3f__4+f__72f__32

(31)

Compare with the output of Γ1,11 in eq (27)

normal

0

(32)

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

Signature: - - - +

_______________________________________________________

gμ,ν=

(33)

seqChristoffelj,α,β,matrix,j=~1,~2,~3,~4

Γ1α,β1α,β=?,Γ2α,β2α,β=?,Γ3α,β3α,β=?,Γ4α,β4α,β=?

(34)

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.