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

Online Help

All Products    Maple    MapleSim


Physics[Weyl] - The Weyl tensor and the Weyl scalars for the tetrad set

Physics[Tetrads][WeylScalars] - The Weyl scalars for any given tetrad

Calling Sequence

Weyl[alpha, beta, mu, nu]

Weyl[keyword]

WeylScalars(tetrad)

Parameters

alpha, beta, mu, nu

-

the indices, as names representing integer numbers between 0 and the spacetime dimension, they can also be the numbers themselves

keyword

-

can be any of: scalars, definition, scalarsdefinition, nonzero, ~, to respectively get the Weyl scalars, definition of the tensor or of its scalars, a set with the nonzero components, or an Array with the contravariant components

tetrad

-

optional, related to Tetrads:-WeylScalar, tetrad is any 4 x 4 matrix for which Tetrads:-IsTetrad returns true.

Description

• 

The Weyl[alpha, beta, mu, nu], displayed as C__α,β,μ,ν, is a computational representation for the Weyl tensor, defined in terms of the Riemann, Ricci and the spacetime metric g_ tensors as

Cα,β,μ,ν=Rα,β,μ,ν+1n2Rα,μgβ,ν+Rα,νgβ,μ+Rβ,μgα,νRβ,νgα,μ+1n1gα,μgβ,νgα,νgβ,μRρρ

  

where n is the dimension of spacetime, by default equal to 4. From this definition, the Weyl tensor has all the symmetries properties of the Riemann tensor; i.e. it is antisymmetric with respect to interchanging the position of its 1st and 2nd indices, or 3rd and 4th indices, and symmetric with respect to interchanging the positions of the 1st and 2nd pair of indices. In addition, it vanishes when contracted on any pair of indices.

• 

The Tetrads:-WeylScalars commands is part of the Tetrads package and computes the Weyl scalars that correspond to a given tetrad, passed as a 4 x 4 matrix that satisfies the tetrad definition. For such tetrads, Tetrads:-IsTetrad returns true. The Weyl command can also compute the Weyl scalars, explained further below, by entering Weyl[scalars], but only for the tetrad set with Setup. In this sense Tetrads:-WeylScalars allows for easy experimentation with different tetrads without having to set the tetrad at every step (see related tetrads example in what is new in Maple 2021).

• 

When the indices of Weyl 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 Weyl is returned. When working with a signature that has the timelike component in position 4, the values 0 and 4, or for the case any dimension instead of 4, set for the spacetime, represent the same object. When the indices have symbolic values Weyl returns unevaluated after normalizing its indices taking into account their symmetry properties.

• 

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 C__α,β,μ,ν 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 Weyl 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 the Christoffel symbols, used to construct the tensors entering the definition of the Weyl tensor. For that purpose see Coordinates or Setup.

• 

Besides being indexed with four indices, Weyl accepts the following keywords:

– 

array: (synonym: Array, Matrix, matrix, or no indices whatsoever, as in Weyl[]) 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 Weyl. If this keyword is passed preceded by the tensor indices, that can be covariant or contravariant, the values in the resulting array are computed taking into account the character of the given indices. If only two of these indices are names, the rest have numerical values (if contravariant then preceded by ~), the returned object is the corresponding 2 x 2 Matrix.

– 

~: returns an Array with the all-contravariant components of the Weyl tensor

– 

definition: returns the definition of the Weyl tensor in terms of the Riemann and Ricci tensors.

– 

nonzero: returns a set of equations, with the left-hand side as a sequence of four positive numbers identifying the element of C__α,β,μ,ν 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.

– 

scalars: returns the five complex Weyl scalars Psi[0], Psi[1], Psi[2], Psi[3], Psi[4] of the Newman-Penrose formalism, relevant in the classification of the PetrovType of a spacetime.

– 

scalarsdefinition: returns the definition of the Weyl scalars; note that a general sign in this definition depends on the signature, that you can query entering Setup(signature) and change it using Setup to any of the four possible signatures predefined.

• 

Some automatic checking and normalization are carried out each time you enter Weyl[...]. The checking is concerned with possible syntax errors. The automatic normalization takes into account the symmetry properties of the indices of Weyl[mu,nu,alpha] as described in the first paragraph.

• 

The %Weyl command is the inert form of Weyl, 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 querying Setup, as in Setup(metric);, or directly entering the metric as g_[], with no indices

g_

gμ,ν=−10000−10000−100001

(3)

Check the nonzero components of Christoffel, used to construct the Weyl tensor entering the definition of Weyl: because the default spacetime is of Minkowski type, there are none

Christoffelnonzero

Γα,μ,ν=

(4)

Hence

Christoffelμ,α,β

0

(5)

Weylμ,ν,α,β

0

(6)

and the same is valid for all the general relativity tensors defined in terms of Christoffel and derivatives of the metric g_. To set the scenario as a curved spacetime set the metric using Setup, for instance indicating the square of the spacetime interval. In this example, we also choose to work in spherical coordinates, so consider for instance the metric defined by

ds2r2dtheta2+r2sinθ2dphi22dtdr2kr,t2dt2

ds2r2dθ2+r2sinθ2dφ22dtdr2kr,t2dt2

(7)

Setupcoordinates=spherical,metric=ds2

Systems of spacetime coordinates are:X=r,θ,φ,t

_______________________________________________________

Coordinates: r,θ,φ,t. Signature: - - - +

_______________________________________________________

gμ,ν=000−10r20000r2sinθ20−1002kr,t2

_______________________________________________________

Setting lowercaselatin_is letters to represent space indices

_______________________________________________________

coordinatesystems=X,metric=1,4=−1,2,2=r2,3,3=r2sinθ2,4,4=2kr,t2,spaceindices=lowercaselatin_is

(8)

To avoid the repeated display of the functionality in kr,t in the following examples and have the display of derivatives in jet notation, indexed, use CompactDisplay

CompactDisplay

kr,twill now be displayed ask

(9)

Now when the indices are not numerical, Weyl returns itself after normalizing its indices taking advantage of their symmetry properties, so that different forms of the same tensor enter computations in the same manner, for example, if you interchange the positions as in

Weylα,β,μ,ν

Cα,β,μ,ν

(10)

Weylμ,ν,β,α

Cα,β,μ,ν

(11)

+

0

(12)

The definition of the Weyl scalars in terms of the null vectors l_, n_, m_ and mb_ of the Newman-Penrose formalism (see Tetrads), and the scalars themselves for the spacetime metric (9)

Weylscalarsdefinition

ψ__0=Cμ,ν,α,βμ,ν,α,βlμmνlαmβ,ψ__1=Cμ,ν,α,βμ,ν,α,βlμnνlαmβ,ψ__2=Cμ,ν,α,βμ,ν,α,βlμmνm&conjugate0;αnβ,ψ__3=Cμ,ν,α,βμ,ν,α,βlμnνm&conjugate0;αnβ,ψ__4=Cμ,ν,α,βμ,ν,α,βnμm&conjugate0;νnαm&conjugate0;β

(13)

Weylscalars

ψ__0=0,ψ__1=0,ψ__2=2kr2r22kr,rkr2+4krkr2k2+16r2,ψ__3=0,ψ__4=0

(14)

Note that these scalars depend on the choice of orientation for the axis of the tetrad system of references. To change the tetrad or that orientation of axis see Tetrads:-TransformTetrad. These same scalars can be computed with the Tetrads:-WeylScalars, which, when called with no arguments, it computes the scalars for the tetrad set, so the same way as Weyl[scalars] above

withTetrads

_______________________________________________________

Setting lowercaselatin_ah letters to represent tetrad indices

Defined as tetrad tensors see ?Physics,tetrads, 𝔢a,μ , ηa,b , γa,b,c , λa,b,c

Defined as spacetime tensors representing the NP null vectors of the tetrad formalism see ?Physics,tetrads, lμ , nμ , mμ , m&conjugate0;μ

_______________________________________________________

IsTetrad,NullTetrad,OrthonormalTetrad,PetrovType,SegreType,TransformTetrad,WeylScalars,e_,eta_,gamma_,l_,lambda_,m_,mb_,n_

(15)

WeylScalars

ψ__0=0,ψ__1=0,ψ__2=2kr2r22kr,rkr2+4krkr2k2+16r2,ψ__3=0,ψ__4=0

(16)

The scalars above are constructed using the scalars definition shown lines above departing from the tetrad

e_

𝔢a,μ=I22k0000−Ir0000−Irsinθ0I22k00Ik2

(17)

With WeylScalars, however, you can also compute the scalars for any given tetrad. Construct, for example, one possible different form of a tetrad, say one for which ψ40, via

Tetrads:-TransformTetradnullrotationwithfixedl_

02r2I22rsinθE&conjugate0;k000k1k2rE+E&conjugate0;2I2sinθ2rEE&conjugate0;kEE&conjugate0;+102r2I22rsinθEk

(18)

To verify that a 4 x 4 matrix is a valid tetrad, i.e., satisfies the tetrad definition,

e_definition

𝔢a,μ𝔢bμbμ=ηa,b

(19)

you can use IsTetrad

IsTetrad

Type of tetrad: null

true

(20)

Now, without setting this tetrad using Setup, compute the corresponding Weyl scalars:

WeylScalars

ψ__0=0,ψ__1=0,ψ__2=2kr2r2+2kr,rkr24krkr+2k216r2,ψ__3=E&conjugate0;2kr2r2+2kr,rkr24krkr+2k212r2,ψ__4=E&conjugate0;22kr2r2+2kr,rkr24krkr+2k21r2

(21)

We see in this result the expected ψ40. When the tetrad is what you want or produces the Weyl scalars you expect, you can set it the standard way using Setup, as in Setup(tetrad ...) where the right-hand side is the 4x4 matrix.

Back to the Weyl tensor, to express it in terms of the Riemann tensor, the Christoffel symbols and its derivatives, or the metric g_ itself, use convert

=convert,Riemann

Cα,β,μ,ν=Rα,β,μ,νRκα,κ,μκα,κ,μgβ,ν2+Rκα,κ,νκα,κ,νgβ,μ2+Rκβ,κ,μκβ,κ,μgα,ν2Rκβ,κ,νκβ,κ,νgα,μ2+Rλκ,λκλκ,λκgα,μgβ,νgα,νgβ,μ6

(22)

=convert,Christoffel

Cα,β,μ,ν=gα,ωμΓωβ,νωβ,ννΓωβ,μωβ,μ+Γωμ,ζωμ,ζΓζβ,νζβ,νΓων,ζων,ζΓζβ,μζβ,μλΓλα,μλα,μμΓλα,λλα,λ+Γσα,μσα,μΓλλ,σλλ,σΓσα,λσα,λΓλμ,σλμ,σgβ,ν2+λΓλα,νλα,ννΓλα,λλα,λ+Γσα,νσα,νΓλλ,σλλ,σΓσα,λσα,λΓλν,σλν,σgβ,μ2+λΓλβ,μλβ,μμΓλβ,λλβ,λ+Γσβ,μσβ,μΓλλ,σλλ,σΓσβ,λσβ,λΓλμ,σλμ,σgα,ν2λΓλβ,νλβ,ννΓλβ,λλβ,λ+Γσβ,νσβ,νΓλλ,σλλ,σΓσβ,λσβ,λΓλν,σλν,σgα,μ2+gχ,κχ,κσΓσχ,κσχ,κχΓσκ,σσκ,σ+Γτχ,κτχ,κΓσσ,τσσ,τΓτκ,στκ,σΓσχ,τσχ,τgα,μgβ,νgα,νgβ,μ6

(23)

=convert,g_

Cα,β,μ,ν=gα,ωμgλ,ωλ,ωνgβ,λ+βgλ,νλgβ,ν2+gλ,ωλ,ωμνgβ,λ+βμgλ,νλμgβ,ν2νgκ,ωκ,ωμgβ,κ+βgκ,μκgβ,μ2gκ,ωκ,ωμνgβ,κ+βνgκ,μκνgβ,μ2+gω,σω,σζgμ,σ+μgσ,ζσgμ,ζgχ,ζχ,ζνgβ,χ+βgχ,νχgβ,ν4gω,τω,τζgν,τ+νgτ,ζτgν,ζgυ,ζυ,ζμgβ,υ+βgμ,υυgβ,μ4λgλ,τλ,τμgα,τ+αgμ,ττgα,μ2+gλ,τλ,τλμgα,τ+αλgμ,τλτgα,μ2μgκ,λκ,λαgκ,λ2gκ,λκ,λαμgκ,λ2+gω,σω,σμgα,ω+αgμ,ωωgα,μgλ,υλ,υσgλ,υ4gψ,σψ,σλgα,ψ+αgλ,ψψgα,λgχ,λχ,λσgχ,μ+μgχ,σχgμ,σ4gβ,ν2+λgλ,τλ,τνgα,τ+αgν,ττgα,ν2+gλ,τλ,τλνgα,τ+αλgν,τλτgα,ν2νgκ,λκ,λαgκ,λ2gκ,λκ,λανgκ,λ2+gω,σω,σνgα,ω+αgν,ωωgα,νgλ,υλ,υσgλ,υ4gψ,σψ,σλgα,ψ+αgλ,ψψgα,λgχ,λχ,λσgχ,ν+νgχ,σχgν,σ4gβ,μ2+λgλ,τλ,τμgβ,τ+βgμ,ττgβ,μ2+gλ,τλ,τλμgβ,τ+βλgμ,τλτgβ,μ2μgκ,λκ,λβgκ,λ2gκ,λκ,λβμgκ,λ2+gω,σω,σμgβ,ω+βgμ,ωωgβ,μgλ,υλ,υσgλ,υ4gψ,σψ,σλgβ,ψ+βgλ,ψψgβ,λgχ,λχ,λσgχ,μ+μgχ,σχgμ,σ4gα,ν2λgλ,τλ,τνgβ,τ+βgν,ττgβ,ν2+gλ,τλ,τλνgβ,τ+βλgν,τλτgβ,ν2νgκ,λκ,λβgκ,λ2gκ,λκ,λβνgκ,λ2+gω,σω,σνgβ,ω+βgν,ωωgβ,νgλ,υλ,υσgλ,υ4gψ,σψ,σλgβ,ψ+βgλ,ψψgβ,λgχ,λχ,λσgχ,ν+νgχ,σχgν,σ4gα,μ2+gχ,κχ,κσgλ,σλ,σκgχ,λ+χgκ,λλgχ,κ2+gλ,σλ,σκσgχ,λ+χσgκ,λλσgχ,κ2χgψ,σψ,σκgψ,σ2gψ,σψ,σχκgψ,σ2+gτ,ζτ,ζκgχ,ζ+χgκ,ζζgχ,κgω,σω,στgω,σ4gι,τι,τσgι,κ+κgι,σιgκ,σgσ,υσ,υτgχ,υ+χgτ,υυgχ,τ4gα,μgβ,νgα,νgβ,μ6

(24)

Check the value of C__α,β,μ,ν, say for α=1, β = 2, μ = 2, ν = 4

Weyl1,2,2,4

kr2r23+kr,rkr232krkr3+k2316

(25)

Check now the value of Cβ,μ,να , with the 1st index contravariant, for the same values of the indices (note you enter the value of the contravariant index prefixed by ~)

Weyl~1,2,2,4

0

(26)

To compute with a representation for Weyl without actually performing the operation, use the inert form %Weyl. To afterwards perform the operation use value

%Weyl~1,2,2,4

C12,2,412,2,4

(27)

value

0

(28)

The all-covariant and all-contravariant components of the Weyl tensor as Arrays

Weyl

Weyl`~`

The (nonzero) Array elements of any these arrays can be computed using ArrayElems.

The nonzero values of C__α,β,μ,ν and of Cβ,μ,να (note you enter the character of the contravariant indices prefixed by ~)

Weylnonzero

Cα,β,μ,ν=1,2,2,4=kr2r23+kr,rkr232krkr3+k2316,1,2,4,2=kr2r23kr,rkr23+2krkr3k23+16,1,3,3,4=sinθ212+kr2r2+kr,rkr22krkr+k23,1,3,4,3=sinθ212+kr2r2+kr,rkr22krkr+k23,1,4,1,4=2kr2r2+2kr,rkr24krkr+2k213r2,1,4,4,1=2kr2r22kr,rkr2+4krkr2k2+13r2,2,1,2,4=kr2r23kr,rkr23+2krkr3k23+16,2,1,4,2=kr2r23+kr,rkr232krkr3+k2316,2,3,2,3=212+kr2r2+kr,rkr22krkr+k2r2sinθ23,2,3,3,2=212+kr2r2+kr,rkr22krkr+k2r2sinθ23,2,4,1,2=kr2r23+kr,rkr232krkr3+k2316,2,4,2,1=kr2r23kr,rkr23+2krkr3k23+16,2,4,2,4=212+kr2r2+kr,rkr22krkr+k2k23,2,4,4,2=212+kr2r2+kr,rkr22krkr+k2k23,3,1,3,4=sinθ212+kr2r2+kr,rkr22krkr+k23,3,1,4,3=sinθ212+kr2r2+kr,rkr22krkr+k23,3,2,2,3=212+kr2r2+kr,rkr22krkr+k2r2sinθ23,3,2,3,2=212+kr2r2+kr,rkr22krkr+k2r2sinθ23,3,4,1,3=sinθ212+kr2r2+kr,rkr22krkr+k23,3,4,3,1=sinθ212+kr2r2+kr,rkr22krkr+k23,3,4,3,4=2k212+kr2r2+kr,rkr22krkr+k2sinθ23,3,4,4,3=2k212+kr2r2+kr,rkr22krkr+k2sinθ23,4,1,1,4=2kr2r22kr,rkr2+4krkr2k2+13r2,4,1,4,1=2kr2r2+2kr,rkr24krkr+2k213r2,4,2,1,2=kr2r23kr,rkr23+2krkr3k23+16,4,2,2,1=kr2r23+kr,rkr232krkr3+k2316,4,2,2,4=212+kr2r2+kr,rkr22krkr+k2k23,4,2,4,2=212+kr2r2+kr,rkr22krkr+k2k23,4,3,1,3=sinθ212+kr2r2+kr,rkr22krkr+k23,4,3,3,1=sinθ212+kr2r2+kr,rkr22krkr+k23,4,3,3,4=2k212+kr2r2+kr,rkr22krkr+k2sinθ23,4,3,4,3=2k212+kr2r2+kr,rkr22krkr+k2sinθ23

(29)

Weyl~alpha,β,μ,ν,nonzero

Cαβ,μ,ναβ,μ,ν=1,1,1,4=2kr2r2+2kr,rkr24krkr+2k213r2,1,1,4,1=2kr2r22kr,rkr2+4krkr2k2+13r2,1,2,1,2=kr2r23+kr,rkr232krkr3+k2316,1,2,2,1=kr2r23kr,rkr23+2krkr3k23+16,1,3,1,3=sinθ212+kr2r2+kr,rkr22krkr+k23,1,3,3,1=sinθ212+kr2r2+kr,rkr22krkr+k23,1,4,1,4=412+kr2r2+kr,rkr22krkr+k2k23r2,1,4,4,1=412+kr2r2+kr,rkr22krkr+k2k23r2,2,1,2,4=2kr2r22kr,rkr2+4krkr2k2+16r2,2,1,4,2=2kr2r2+2kr,rkr24krkr+2k216r2,2,3,2,3=2sinθ212+kr2r2+kr,rkr22krkr+k23,2,3,3,2=2sinθ212+kr2r2+kr,rkr22krkr+k23,2,4,1,2=2kr2r2+2kr,rkr24krkr+2k216r2,2,4,2,1=2kr2r22kr,rkr2+4krkr2k2+16r2,2,4,2,4=212+kr2r2+kr,rkr22krkr+k2k23r2,2,4,4,2=212+kr2r2+kr,rkr22krkr+k2k23r2,3,1,3,4=2kr2r22kr,rkr2+4krkr2k2+16r2,3,1,4,3=2kr2r2+2kr,rkr24krkr+2k216r2,3,2,2,3=2kr2r23+2kr,rkr234krkr3+2k2313,3,2,3,2=2kr2r232kr,rkr23+4krkr32k23+13,3,4,1,3=2kr2r2+2kr,rkr24krkr+2k216r2,3,4,3,1=2kr2r22kr,rkr2+4krkr2k2+16r2,3,4,3,4=212+kr2r2+kr,rkr22krkr+k2k23r2,3,4,4,3=212+kr2r2+kr,rkr22krkr+k2k23r2,4,2,2,4=kr2r23kr,rkr23+2krkr3k23+16,4,2,4,2=kr2r23+kr,rkr232krkr3+k2316,4,3,3,4=sinθ212+kr2r2+kr,rkr22krkr+k23,4,3,4,3=sinθ212+kr2r2+kr,rkr22krkr+k23,4,4,1,4=2kr2r22kr,rkr2+4krkr2k2+13r2,4,4,4,1=2kr2r2+2kr,rkr24krkr+2k213r2

(30)

This is the array form of Cβ,μ,να (note you enter the character of the contravariant indices prefixed by ~)

Weyl~alpha,β,μ,ν,array

To use this array, because its components were computed already taking into account the (covariant/contravariant) character of its indices, you do not need to indicate furthermore that character. So for C11 we have

C1,2,2,4

0

(31)

Compare with the output obtained entering Weyl[~1, 2, 2, 4] in eq (14)

0

(32)

Verify that the Weyl tensor vanishes when contracted on a pair of indices; take for instance the 1st and 3rd indices, you can see the corresponding matrix (note the first index as contravariant)

Weyl~1,β,1,ν,matrix

C1β,1,ν1β,1,ν=0002kr2r2+2kr,rkr24krkr+2k213r20kr2r23+kr,rkr232krkr3+k23160000sinθ212+kr2r2+kr,rkr22krkr+k230000412+kr2r2+kr,rkr22krkr+k2k23r2

(33)

Add the four matrices completing the contraction of indices

Weyl~1,β,1,ν,matrix+Weyl~2,β,2,ν,matrix+Weyl~3,β,3,ν,matrix+Weyl~4,β,4,ν,matrix

C1β,1,ν1β,1,ν+C2β,2,ν2β,2,ν+C3β,3,ν3β,3,ν+C4β,4,ν4β,4,ν=0002kr2r2+2kr,rkr24krkr+2k213r2+2kr2r22kr,rkr2+4krkr2k2+13r2000000002kr2r2+2kr,rkr24krkr+2k213r2+2kr2r22kr,r