Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
248 views

Mathematica PDF

This document describes a Mathematica program to calculate Christoffel symbols and the geodesic equation from a given metric tensor. It defines the Christoffel symbols in terms of the metric and its inverse, and calculates the geodesic equation as the derivative of the four-velocity with respect to proper time using the Christoffel symbols and four-velocity. The program is demonstrated on a wormhole metric.

Uploaded by

Asanka Amarasinghe
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
248 views

Mathematica PDF

This document describes a Mathematica program to calculate Christoffel symbols and the geodesic equation from a given metric tensor. It defines the Christoffel symbols in terms of the metric and its inverse, and calculates the geodesic equation as the derivative of the four-velocity with respect to proper time using the Christoffel symbols and four-velocity. The program is demonstrated on a wormhole metric.

Uploaded by

Asanka Amarasinghe
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 3

Christoffel Symbols and Geodesic Equation

This is a Mathematica program to compute the Christoffel and the geodesic equations, starting from a
given metric g . The Christoffel symbols are calculated from the formula

1
2

where g is the matrix inverse of g called the inverse metric. This is the solution of the relation (8.19)
and the notation for the inverse metric is standard [cf (20.17)]. The components of the geodesic equation are
du /d = -

u u .

You must input the covariant components of the metric tensor g by editing the relevant input line in
this Mathematica notebook. You may also wish to change the names of the coordinates. The nonzero
components of the above quantities are displayed as the output.

Clearing the values of symbols:


First clear any values that may already have been assigned to the names of the various objects to be
calculated. The names of the coordinates that you will use are also cleared.
Clear coord, metric, inversemetric, affine, r, , , t

Setting The Dimension


The dimension n of the spacetime (or space) must be set:
n

Defining a list of coordinates:


The example given here is the wormhole metric (7.40). Note that for convenience t is denoted by x4
rather than x0 and summations run from 1 to 4 rather than 0 to 3.
coord

r, , , t

r, , , t

You can change the names of the coordinates by simply editing the definition of coord, for example, to
coord = {x, y, z, t}, when another set of coordinate names is more appropriate.

Defining the metric:


Input the metric as a list of lists, i.e., as a matrix. You can input the components of any metric here, but
you must specify them as explicit functions of the coordinates.

Printed by Wolfram Mathematica Student Edition

christoffel.nb

metric

1, 0, 0, 0 , 0, r ^ 2

1, 0, 0, 0 , 0, b2

b ^ 2, 0, 0 , 0, 0, r ^ 2

r2 , 0, 0, 0, 0, b2

r2 Sin

b ^ 2 Sin ^ 2, 0 , 0, 0, 0,

, 0, 0, 0, 0,

You can also display this in matrix form:


metric
1

MatrixForm
0

0 b2

0
r2

0
b

r Sin
2

0
1

Note:
It is important not to use the symbols, i, j, k, l, n, or s as constants or coordinates in the metric that you
specify above. The reason is that the first four of those symbols are used as summation or table
indices in the calculations done below. The last is the dimension of the space.

Calculating the inverse metric:


The inverse metric is obtained through matrix inversion.
inversemetric

Simplify Inverse metric

1, 0, 0, 0 , 0,

1
b2

r2

, 0, 0, 0, 0,

Csc
b2

, 0, 0, 0, 0,

r2

The inverse metric can also be displayed in matrix form:


inversemetric
1

1
b2 r2

MatrixForm

0
0

Csc 2
b2 r2

Calculating the affine connection:


The calculation of the components of the affine connection is done by transcribing the definition given
earlier into the notation of Mathematica and using the Mathematica functions D for taking partial derivatives, Sum for summing over repeated indices, Table for forming a list of components, and Simplify for
simplifying the result.
affine : affine Simplify Table 1 2
D metric s, j , coord k
D metric s, k , coord j
i, 1, n , j, 1, n , k, 1, n

Sum

inversemetric

D metric

j, k

Displaying the affine connection:

Printed by Wolfram Mathematica Student Edition

i, s

, coord

, s, 1, n

christoffel.nb

The nonzero components of the affine connections are displayed below. You need not follow the details
of constructing the functions that we use for that purpose. Because the affine connection is symmetric
under interchange of the last two indices, only the independent components are displayed.
listaffine : Table If UnsameQ affine i, j, k , 0 ,
ToString
i, j, k , affine i, j, k
, i, 1, n , j, 1, n , k, 1, j
TableForm Partition DeleteCases Flatten listaffine , Null , 2 , TableSpacing
1, 2, 2

1, 3, 3

r Sin

2, 2, 1

2, 2

r
b2 r2

Cos Sin

2, 3, 3
3, 3, 1

r
b2 r2

3, 3, 2

Cot

Calculating the geodesic equations:


The geodesic equations are calculated by asking Mathematica to carry out the sum

u , where

u are the components of the four-velocity. (This gives the derivitive of u with respect to proper time
. (This is replaced by s if the geodesics are spacelike.)
geodesic : geodesic Simplify Table
k, 1, n , i, 1, n

Sum affine

i, j, k

u j u k , j, 1, n ,

Displaying the geodesic equations:


listgeodesic : Table

"d d" ToString u i

TableForm listgeodesic, TableSpacing


d d u 1

r u 2

Sin

u 3

, i, 1, n

d d u 2

2ru 1 u 2
b2 r2

Cos Sin u 3

d d u 3

2ru 1 u 3
b2 r2

2 Cot u 2 u 3

d d u 4

, " ", geodesic

Acknowledgment
This program was adapted from the notebook Curvature and the Einstein equation kindly written by
Leonard Parker especially for this text.

Printed by Wolfram Mathematica Student Edition

You might also like