Tensor[NPDirectionalDerivatives] - create procedures for the Newman-Penrose directional derivatives
Calling Sequences
NPDirectionalDerivatives(NTetrad, output)
NPDirectionalDerivatives(Fr, output)
Parameters
NTetrad - a list of 4 vectors defining a null tetrad
Fr - the name of an initialized anholonomic frame, created from a null tetrad
output - (optional) keyword argument output = "sequence"
Description
Examples
See Also
Let g be a metric with signature 1,−1,−1,−1 and L,N,M,M‾ be a null tetrad for g. The Newman-Penrose directional derivative operators take derivatives of functions in the directions of the 4 vectors of the tetrad. See Newman and Penrose, Stewart.
The NPDirectionalDerivatives command returns a table with 4 entries "D", "Delta", "delta", "barDelta" corresponding to the customary labels assigned to the NP directional derivative operators (D,Δ,δ, δ‾). With optional keyword argument output = "sequence", the directional derivative operators are returned as a sequence of 4 Maple procedures.
This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form NPDirectionalDerivatives(...) only after executing the commands with(DifferentialGeometry); with(Tensor); in that order. It can always be used in the long form DifferentialGeometry:-Tensor:-NPDirectionalDerivatives.
with⁡DifferentialGeometry:with⁡Tensor:
Example 1.
Define a manifold S with coordinates t,x,y,z.
DGsetup⁡t,x,y,z,S
frame name: S
Define a metric g of signature (1,-1,-1,-1).
g≔evalDG⁡x2⁢dt&tdt−y2⁢dx&tdx−z2⁢dy&tdy−t2⁢dz&tdz
g:=x2⁢dt⁢dt−y2⁢dx⁢dx−z2⁢dy⁢dy−t2⁢dz⁢dz
Define a null tetrad NTetrad for the metric g. Use GRQuery to verify that NTetrad is a null tetrad.
NTetrad≔evalDG⁡12⁢212x⁢D_t+12⁢212t⁢D_z,12⁢212x⁢D_t−12⁢212t⁢D_z,12⁢212y⁢D_x+12⁢I⁢212z⁢D_y,12⁢212y⁢D_x−12⁢I⁢212z⁢D_y
NTetrad:=22⁢x⁢D_t+22⁢t⁢D_z,22⁢x⁢D_t−22⁢t⁢D_z,22⁢y⁢D_x+I2⁢2z⁢D_y,22⁢y⁢D_x−I2⁢2z⁢D_y
GRQuery⁡NTetrad,g,NullTetrad
true
Calculate the NP directional derivative operators defined by the null tetrad NTetrad. The result is a table with 4 indices.
NPDirDer≔NPDirectionalDerivatives⁡NTetrad
NPDirDer:=tableDelta=f→LieDerivative⁡E2,f,delta=f→LieDerivative⁡E3,f,D=f→LieDerivative⁡E1,f,bardelta=f→LieDerivative⁡E4,f
indices⁡NPDirDer
Delta,delta,D,bardelta
The individual directional derivatives operators can be extracted from the table NPDirDer.
X≔NPDirDerD
X:=NPDirDerD
X⁡f⁡t,x,y,z
12⁢2⁢∂∂t⁢f⁡t,x,y,zx+12⁢2⁢∂∂z⁢f⁡t,x,y,zt
Note that this coincides with the Lie derivative of ft,x,y,z with respect to the first vector in NTetrad.
LieDerivative⁡NTetrad1,f⁡t,x,y,z
Example 2.
With the keyword argument output = "sequence", the command NPDirectionalDerivatives will return the directional derivatives operators as a sequence. (Note that the name D is protected by Maple.)
d,Δ,δ,bardela≔NPDirectionalDerivatives⁡NTetrad,output=Sequence
d,Δ,δ,bardela:=f→DifferentialGeometry:-LieDerivative⁡E1,f,f→DifferentialGeometry:-LieDerivative⁡E2,f,f→DifferentialGeometry:-LieDerivative⁡E3,f,f→DifferentialGeometry:-LieDerivative⁡E4,f
Δ⁡f⁡t,x,y,z
12⁢2⁢∂∂t⁢f⁡t,x,y,zx−12⁢2⁢∂∂z⁢f⁡t,x,y,zt
Example 3.
When working with the NP formalism, it is usually advantageous to work with the anholonomic frame defined by the null tetrad. To create anholonomic frames in DifferentialGeometry, see FrameData.
FD≔FrameData⁡NTetrad,NP
FD:=E1,E2=12⁢2⁢E1t⁢x−12⁢2⁢E2t⁢x,E1,E3=14⁢2⁢E1x⁢y+14⁢2⁢E2x⁢y−14⁢2⁢E3t⁢z+14⁢2⁢E4t⁢z,E1,E4=14⁢2⁢E1x⁢y+14⁢2⁢E2x⁢y+14⁢2⁢E3t⁢z−14⁢2⁢E4t⁢z,E2,E3=14⁢2⁢E1x⁢y+14⁢2⁢E2x⁢y+14⁢2⁢E3t⁢z−14⁢2⁢E4t⁢z,E2,E4=14⁢2⁢E1x⁢y+14⁢2⁢E2x⁢y−14⁢2⁢E3t⁢z+14⁢2⁢E4t⁢z,E3,E4=−12⁢I⁢2⁢E3z⁢y−12⁢I⁢2⁢E4z⁢y
DGsetup⁡FD
frame name: NP
We can now calculate the directional derivatives operators for the null tetrad with the second calling sequence.
NPDirectionalDerivatives⁡NP
tableDelta=f→DifferentialGeometry:-LieDerivative⁡E2,f,delta=f→DifferentialGeometry:-LieDerivative⁡E3,f,D=f→DifferentialGeometry:-LieDerivative⁡E1,f,bardelta=f→DifferentialGeometry:-LieDerivative⁡E4,f
DifferentialGeometry, Tensor, LieDerivative, GRQuery, NullTetrad
Download Help Document