computes and solves the Killing equations - Maple Help

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Physics : Physics/KillingVectors

Physics[KillingVectors] - computes and solves the Killing equations

Calling Sequence

KillingVectors(V, options = ...)

Parameters

V

-

a name to be used for the vector field (rank 1 spacetime tensor) representing the Killing vector

integrabilityconditions = ...

-

optional - can be true (default) or false, to reduce or not the determining PDE output system taking into account its integrability conditions

output = ...

-

optional - the right-hand-side can be equations, solutions (synonym: componentsolutions) or vectorsolutions (default), to return the Killing equations or attempt returning their solutions for each of the Killing vector components, or a vector solution with the Killing vector on the left-hand side and a list with its components on the right-hand side

parameters = ...

-

optional - the right-hand-side can be a set or a list of names or functions, that parameterize the Killing vectors; the system of equations or solutions will be split according to the different cases of these parameters

specializeconstants = ...

-

optional - can be true (default) or false, to specialize or not the integration constants _Cn that appear when solving the Killing differential equations

Description

• 

KillingVectors computes and attempts to solve the Killing equations, that is, the PDE system

𝒟νVα+𝒟αVν=Vgα,ν=0

  

where 𝒟 is the covariant derivative operator D_, and V entering Vα is the name representing the Killing vectors, provided, with or without a spacetime index, and Vgα,ν is the LieDerivative of the spacetime metric with respect to Vμ. In a spacetime of dimension n there are at most 12nn+1 independent Killing vectors, where the maximum happens for a Minkowski spacetime. As it is the case of the other commands in Physics, KillingVectors operates in tensorial notation. To perform the same operation using differential forms notation see DifferentialGeometry[Tensor][KillingVectors].

• 

When the name V is provided without an index, the equations computed are satisfied by the contravariant components Vα. To request the equations or solutions for the covariant components Vα, pass it with a covariant index, say as in V[alpha]. Note that when the spacetime is non-galilean, while the system may succeed in computing solutions for the covariant components, depending on the metric g_ it may fail in computing solutions for the contravariant ones, or vice-versa.

• 

By default, KillingVectors automatically computes the equations and attempts solving them in one step. When successful, KillingVectors returns a list of Killing vectors, obtained after specializing the integration constants that appear when solving the system of differential equations (that you can see using the option output = equations). When no solution is found, KillingVectors returns NULL, the same way the PDE solver pdsolve does.

• 

A more compact solution can be obtained by passing the optional argument specializeconstants = false, in which case KillingVectors returns an equation with the Killing vector on the left-hand side and a list with its components, depending on integration constants, on the right-hand side.

• 

When passing the optional argument output = solutions, instead of returning a vector equation, the output is a set of solution equations, with each Killing vector component on the left-hand side and the corresponding solution on the right-hand side. This is useful, for instance, for verification purposes: you can test the solutions computed against the output of KillingVectors using the option output = equations by using the pdetest command. -When passing the optional argument output = equations, the Killing equations themselves are returned without attempting solving them. This is useful when KillingVectors fails in finding a solution, or when particular solutions of different forms may be of interest; these solutions can frequently be computed using the symmetry commands of PDEtools like FunctionFieldSolutions and PolynomialSolutions, or mainly InvariantSolutions used together with its various options for particularizing solutions.

• 

The PDE system returned by KillingVectors with the option output = equations is by default returned as a list of equations and automatically reduced taking into account its integrability conditions. The order of the equations in the list corresponds to a total degree, obtained with an orderly ranking (see details in casesplit). Simplifying the PDE system taking into account integrability conditions is a frequently desired but in some cases expensive mathematical computation. To avoid it, pass the optional argument integrabilityconditions = false, in which case KillingVectors will return faster, a set of equations, not a list, and with no differential elimination simplifications.

• 

When the metric g_ depends on parameters, either symbols or functions of spacetime variables, the Killing vectors computed using KillingVectors 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 right-hand-side is a set or a list with the parameters, and the system of equations will be split into cases with respect to these parameters using differential algebra techniques, and the Killing vectors will be expressed as a piecewise function according to the different cases for these parameters. As frequently happens when splitting into cases, the resulting lists of equations may involve inequations as well. Note: passing parameters automatically forces reducing the PDE system using integrability conditions, regardless of the value of the keyword integrabilityconditions.

Examples

withPhysics:

Setupmathematicalnotation=true

mathematicalnotation=true

(1)

Set a system of Coordinates

Setupcoordinates=X

* 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

(2)

Define a tensor with one index to represent the Killing vector

DefineV

Defined objects with tensor properties

V,γμ,σμ,Xμ,μ,gμ,ν,δμ,ν,εα,β,μ,ν

(3)

The contravariant components  Vμ

KillingVectorsV

Vμμ=x2,x1,0,0,Vμμ=0,0,1,0,Vμμ=x3,0,x1,0,Vμμ=x4,0,0,x1,Vμμ=1,0,0,0,Vμμ=0,x4,0,x2,Vμμ=0,0,x4,x3,Vμμ=0,0,0,1,Vμμ=0,x3,x2,0,Vμμ=0,1,0,0

(4)

The same result but directly for the covariant components and without specializing the integration constants,

KillingVectorsVμ,specialize=false

* Partial match of 'specialize' against keyword 'specializeconstants'

Vμ=_C1x2+_C2x3+_C3x4+_C4,_C1x1+_C5x3+_C6x4+_C7,_C2x1_C5x2+_C8x4+_C9,_C3x1_C6x2_C8x3+_C10

(5)

Note that there are 10 integration constants, according to the fact that the dimension of spacetime is n=4 and there are at most 12nn+1 independent Killing vectors, as it is the case of a Minkowski spacetime. To understand this result, recall that the Killing vectors generate transformations that leave the metric g_ invariant in form; these are isometries. In this case, the spacetime loaded with Physics by default is galilean

g_[]

gμ,ν=1000010000100001

(6)

The result (4) generates translations along the four axis and rotations that leave invariant the distance between two points and thus invariant in form the metric.

Consider now a non-galilean spacetime, for instance set the Schwarzschild metric as the current metric (see g_):

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

gμ,ν=rr+2m0000r20000r2sinθ20000r2mr

(7)

Compute the Killing vectors again. For illustration purposes, compute the covariant components and compute first only the equations, with and without simplifying using integrability conditions, then compute everything automatically. First using integrability conditions. To avoid displaying VX repeatedly use the compact notation available of the PDEtools package

PDEtools:-declareVX

Vr,θ,φ,twill now be displayed asV

(8)

KillingVectorsVα,output=equations

V3r=2V3r,V3θ=2cosθV3V2φsinθsinθ,V3φ=sinθcosθV2,V3t=0,V2φ,φ=V2sinθ2+cosθ2,V2r=2V2r,V2θ=0,V2t=0,V4r=2mV4rr+2m,V4θ=0,V4φ=0,V4t=0,V1=0

(9)

where derivatives are displayed indexed due to the use of declare. Compare with

KillingVectorsVα,output=equations,integrabilityconditions=false

2V4tr3+4V1mr2+m=0,2V2θ+4m+2rV1=0,V4φ+V3t=0,V4θ+V2t=0,2mrr2V1t+2mrr2V4r+2mV4=0,V2rr+V1θr2V2=0,V3θsinθ+V2φsinθ2cosθV3=0,V3rr+V1φr2V3=0,4V1rrm2V1rr22mV1=0,2V3φ+4m2rcosθ24m+2rV1+2sinθcosθV2=0

(10)

The contravariant components of the Killing vectors all in one step

KillingVectorsV

Vμμ=0,0,0,1,Vμμ=0,sinφ,cosφtanθ,0,Vμμ=0,cosφ,sinφtanθ,0,Vμμ=0,0,1,0

(11)

The same solution but with the constants not specialized

KillingVectorsV,specializeconstants=false

Vμμ=0,_C2sinφ+_C3cosφ,_C2cosφ_C3sinφtanθ+_C4,_C1

(12)

Taking m as parameter, there is another solution (from the point of view of differential algebra this other solution is a singular solution). These solutions can also be expressed as solution equations with the vector components on the left-hand sides by using the optional argument output = solutions:

KillingVectorsV,parameters=m

Vμμ=cosθt,sinθtr,0,cosθr,Vμμ=0,0,1,0,Vμμ=sinφtsinθ,sinφ+θ+sinφθt2r,cosφtsinθr,rsinφsinθ,Vμμ=cosφtsinθ,cosφθ+cosφ+θt2r,sinφtsinθr,rcosφsinθ,Vμμ=cosθ,sinθr,0,0,Vμμ=sinφsinθ,sinφcosθr,cosφsinθr,0,Vμμ=cosφsinθ,cosφcosθr,sinφsinθr,0,Vμμ=0,0,0,1,Vμμ=0,sinφ,cosφtanθ,0,Vμμ=0,cosφ,sinφtanθ,0m=0Vμμ=0,0,0,1,Vμμ=0,sinφ,cosφtanθ,0,Vμμ=0,cosφ,sinφtanθ,0,Vμμ=0,0,1,0otherwise

(13)

When parameters are indicated together with the specializeconstants option, the result is a sequence of piecewise functions, with each conditional on the parameters expressed as a branch of the piecewise solution

KillingVectorsV,parameters=m,specializeconstants

Vμμ=cosθt,sinθtr,0,cosθr,Vμμ=0,0,1,0,Vμμ=sinφtsinθ,sinφ+θ+sinφθt2r,cosφtsinθr,rsinφsinθ,Vμμ=cosφtsinθ,cosφθ+cosφ+θt2r,sinφtsinθr,rcosφsinθ,Vμμ=cosθ,sinθr,0,0,Vμμ=sinφsinθ,sinφcosθr,cosφsinθr,0,Vμμ=cosφsinθ,cosφcosθr,sinφsinθr,0,Vμμ=0,0,0,1,Vμμ=0,sinφ,cosφtanθ,0,Vμμ=0,cosφ,sinφtanθ,0m=0Vμμ=0,0,0,1,Vμμ=0,sinφ,cosφtanθ,0,Vμμ=0,cosφ,sinφtanθ,0,Vμμ=0,0,1,0otherwise

(14)

See Also

casesplit, Coordinates, D_, declare, Define, DifferentialAlgebra, DifferentialGeometry[Tensor]KillingVectors], FunctionFieldSolutions, g_, InvariantSolutions, Library, LieBracket, LieDerivative, pdsolve, Physics, Physics conventions, Physics examples, PolynomialSolutions, Riemann, 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.

  

Weinberg, S. Gravitation and Cosmology: Principles and Applications of The General Theory of Relativity, John Wiley & Sons, Inc, 1972.


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam