Physics[Geodesics]  computes and solves the geodesic equations for the spacetime metric

Calling Sequence


Geodesics(tau, options = ...)


Parameters


tau



optional  a name to represent the affine parameter that parametrizes the geodesic equations; if not given, the symbol is used

integrabilityconditions = ...



optional  can be true (default) or false, to reduce, or not reduce, the resulting ODE output system to a canonical form (see RIF)

output = ...



optional  the righthandside can be equations (default) or solutions, to return the geodesic equations or also attempt solving them

parameters = ...



optional  the righthandside can be a set or a list of names or functions, that parameterize the spacetime metric; the system of geodesic equations (or solutions) will be split according to the different cases of these parameters

tensornotation = ...



optional  can be true or false (default), to return the geodesic equations as an algebraic equation in tensor notation instead of as a set or list of equations for the components





Description


•

Geodesics computes and attempts to solve the geodesic equations for the spacetime metric g_. These equations are parametrized by an affine parameter typically representing the proper time for a timelike curve, or distance for a spacelike curve. These geodesic equations computed are

•

When a name  say tau is provided as first argument, it will be used to represent the parameter for the geodesic equations; otherwise, the symbol itself is used.

•

By default, Geodesics computes only the equations, without attempting solving them. To request compute and solve the equations in one step, pass the option output = solutions. When Geodesics or dsolve fail in finding a solution for the equations, you can still look for particular solutions of different forms that may be of interest using the symmetry commands of PDEtools like FunctionFieldSolutions and PolynomialSolutions, or InvariantSolutions used together with its various options for particularizing solutions.

•

The ODE system returned by default by Geodesics is returned as a list of equations and automatically simplified to a canonical form (reduced involutive form  RIF and DifferentialAlgebra) using differential algebra techniques. The order of the equations in the list corresponds to a total degree, obtained with an orderly ranking (see details in casesplit). Simplifying the ODE system this way is a frequently desired but in some cases expensive mathematical computation. To avoid it, pass the optional argument integrabilityconditions = false, in which case Geodesics will return faster, a set of equations, not a list, and with no differential simplification.

•

When computing just the geodesic equations, it also possible to request them in tensor notation, as an algebraic equation, instead of as a set or list of equations for the tensorial components. Although this kind of output always looks just like the definition for the geodesic equations shown above, the algebraic equation returned actually contains the computable objects like the SpaceTimeVector and Christoffel symbols, and can be manipulated in different ways, including getting from it all the equations for its components.

•

When the metric g_ depends on parameters, either symbols or functions of spacetime variables, the equations computed with Geodesics are expected to be valid for arbitrary values of these parameters. It is sometimes of interest, however, to investigate the different kinds of solutions that may exist for different particular values of these parameters. To perform such an investigation, use the optional argument parameters = ... where the righthandside is a set or a list with the parameters. The system of geodesic equations will then be split into cases with respect to these parameters using differential algebra techniques. As frequently happens when splitting into cases, the resulting lists of equations may involve inequations as well. Note: passing parameters automatically forces simplifying the ODE system to a canonical (reduced involutive RIF) form regardless of the value of the keyword integrabilityconditions.



Compatibility


•

The Physics[Geodesics] command was introduced in Maple 17.



Examples


>


>


 (1) 
Set a system of Coordinates
>


 (2) 
The spacetime metric at this point is the one that loads with Physics, i.e. a Minkowski type metric
>


 (3) 
The geodesics equations for this metric, first as an algebraic equation in tensor notation, then as a list of equations for the tensorial components
>


 (4) 
>


 (5) 
The straightlines solution computed directly
>


 (6) 
Consider now a nongalilean spacetime, for instance set the Schwarzschild metric as the current metric (see g_):
>


 (7) 
Compute the geodesic equations again.
>


 (8) 
>


 (9) 
This system of ODEs, as is, it is out of reach of the DE solvers of the system mainly due to the presence of nonrational objects like sin and cos having for arguments one of the unknowns of the system, . On the other hand, we know the geodesics for the Schwarzschild metric describe the motion of particles of infinitesimal mass in the gravitational field of a central fixed large mass. So to investigate the solvability of these equations one can assume is constant an due to the rotational symmetry choose a value for it that simplifies the equations, for example, .
>


 (10) 
And this system is already solvable with dsolve; the first two of six solutions are:
>


 (11) 
>




See Also


casesplit, Coordinates, DifferentialAlgebra, DifferentialGeometry[Tensor]Geodesics], FunctionFieldSolutions, g_, InvariantSolutions, Library, pdsolve, Physics, Physics conventions, Physics examples, PolynomialSolutions, Setup


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.


