Weyl - Maple Help

All Products Maple MapleSim

Home : Support : Online Help : Physics : Weyl

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_{α, β, μ, ν} + \frac{1}{n - 2} (- R_{α, μ} g_{β, ν} + R_{α, ν} g_{β, μ} + R_{β, μ} g_{α, ν} - R_{β, ν} g_{α, μ} + \frac{1}{n - 1} (g_{α, μ} 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

>	$with (Physics) &colon;$

>	$Setup (mathematicalnotation = 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

>	$Setup (coordinatesystems = 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_{μ, ν} = [\begin{array}{c} −1 & 0 & 0 & 0 \\ 0 & −1 & 0 & 0 \\ 0 & 0 & −1 & 0 \\ 0 & 0 & 0 & 1 \end{array}]$

(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

>	$Christoffel [nonzero]$

$Γ_{α, μ, ν} = \emptyset$

(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

>	$ds2 ≔ r^{2} {dtheta}^{2} + r^{2} {\sin (θ)}^{2} {dphi}^{2} - 2 dt dr - 2 {k (r, t)}^{2} {dt}^{2}$

$ds2 ≔ r^{2} {dθ}^{2} + r^{2} {\sin (θ)}^{2} {dφ}^{2} - 2 dt dr - 2 {k (r, t)}^{2} {dt}^{2}$

(7)

>	$Setup (coordinates = spherical, metric = ds2)$

$Systems of spacetime coordinates are: \{X = (r, θ, φ, t)\}$

$_______________________________________________________$

$Coordinates: [r, θ, φ, t] . Signature: (- - - +)$

$_______________________________________________________$

$g_{μ, ν} = [\begin{array}{c} 0 & 0 & 0 & −1 \\ 0 & r^{2} & 0 & 0 \\ 0 & 0 & r^{2} {\sin (θ)}^{2} & 0 \\ −1 & 0 & 0 & - 2 {k (r, t)}^{2} \end{array}]$

$_______________________________________________________$

$Setting lowercaselatin_is letters to represent space indices$

$_______________________________________________________$

$[coordinatesystems = \{X\}, metric = \{(1, 4) = −1, (2, 2) = r^{2}, (3, 3) = r^{2} {\sin (θ)}^{2}, (4, 4) = - 2 {k (r, t)}^{2}\}, spaceindices = lowercaselatin_is]$

(8)

To avoid the repeated display of the functionality in $k (r, t)$ in the following examples and have the display of derivatives in jet notation, indexed, use CompactDisplay

>	$CompactDisplay ()$

$k (r, t) will now be displayed as k$

(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)

>	$Weyl [scalarsdefinition]$

$ψ__0 = - C_{}^{μ, ν, α, β} l_{μ} m_{ν} l_{α} m_{β}, ψ__1 = - C_{}^{μ, ν, α, β} l_{μ} n_{ν} l_{α} m_{β}, ψ__2 = - C_{}^{μ, ν, α, β} l_{μ} m_{ν} {\overset{&conjugate0;}{m}}_{α} n_{β}, ψ__3 = - C_{}^{μ, ν, α, β} l_{μ} n_{ν} {\overset{&conjugate0;}{m}}_{α} n_{β}, ψ__4 = - C_{}^{μ, ν, α, β} n_{μ} {\overset{&conjugate0;}{m}}_{ν} n_{α} {\overset{&conjugate0;}{m}}_{β}$

(13)

>	$Weyl [scalars]$

$ψ__0 = 0, ψ__1 = 0, ψ__2 = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, ψ__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

>	$with (Tetrads)$

$_______________________________________________________$

$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_{μ}, {\overset{&conjugate0;}{m}}_{μ}$

$_______________________________________________________$

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

(15)

>	$WeylScalars ()$

$ψ__0 = 0, ψ__1 = 0, ψ__2 = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, ψ__3 = 0, ψ__4 = 0$

(16)

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

$e_[]$

$𝔢_{a, μ} = [\begin{array}{c} \frac{- \frac{I}{2} \sqrt{2}}{k} & 0 & 0 & 0 \\ 0 & −I r & 0 & 0 \\ 0 & 0 & −I r \sin (θ) & 0 \\ \frac{\frac{I}{2} \sqrt{2}}{k} & 0 & 0 & I k \sqrt{2} \end{array}]$

(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 $ψ_{4} \neq 0$ , via

>	$Tetrads :- TransformTetrad (nullrotationwithfixedl_)$

$[\begin{array}{c} 0 & \frac{\sqrt{2} r}{2} & - \frac{I}{2} \sqrt{2} r \sin (θ) & - \overset{&conjugate0;}{E} k \\ 0 & 0 & 0 & - k \\ - \frac{1}{k} & \frac{\sqrt{2} r (E + \overset{&conjugate0;}{E})}{2} & - \frac{I}{2} \sin (θ) \sqrt{2} r (E - \overset{&conjugate0;}{E}) & - k (E \overset{&conjugate0;}{E} + 1) \\ 0 & \frac{\sqrt{2} r}{2} & \frac{I}{2} \sqrt{2} r \sin (θ) & - E k \end{array}]$

(18)

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

>	$e_[definition]$

$𝔢_{a, μ} 𝔢_{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 = - \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{6 r^{2}}, ψ__3 = - \frac{\overset{&conjugate0;}{E} (2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1)}{2 r^{2}}, ψ__4 = - \frac{{\overset{&conjugate0;}{E}}^{}^{2} (2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1)}{r^{2}}$

(21)

We see in this result the expected $ψ_{4} \neq 0$ . 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_{α, β, μ, ν} - \frac{R_{α, κ, μ}^{κ} g_{β, ν}}{2} + \frac{R_{α, κ, ν}^{κ} g_{β, μ}}{2} + \frac{R_{β, κ, μ}^{κ} g_{α, ν}}{2} - \frac{R_{β, κ, ν}^{κ} g_{α, μ}}{2} + \frac{R_{κ, λ}^{λ κ} (g_{α, μ} g_{β, ν} - g_{α, ν} g_{β, μ})}{6}$

(22)

>	$= convert (, Christoffel)$

$C_{α, β, μ, ν} = g_{α, ω} (\partial_{μ} (Γ_{β, ν}^{ω}) - \partial_{ν} (Γ_{β, μ}^{ω}) + Γ_{μ, ζ}^{ω} Γ_{β, ν}^{ζ} - Γ_{ν, ζ}^{ω} Γ_{β, μ}^{ζ}) - \frac{(\partial_{λ} (Γ_{α, μ}^{λ}) - \partial_{μ} (Γ_{α, λ}^{λ}) + Γ_{α, μ}^{σ} Γ_{λ, σ}^{λ} - Γ_{α, λ}^{σ} Γ_{μ, σ}^{λ}) g_{β, ν}}{2} + \frac{(\partial_{λ} (Γ_{α, ν}^{λ}) - \partial_{ν} (Γ_{α, λ}^{λ}) + Γ_{α, ν}^{σ} Γ_{λ, σ}^{λ} - Γ_{α, λ}^{σ} Γ_{ν, σ}^{λ}) g_{β, μ}}{2} + \frac{(\partial_{λ} (Γ_{β, μ}^{λ}) - \partial_{μ} (Γ_{β, λ}^{λ}) + Γ_{β, μ}^{σ} Γ_{λ, σ}^{λ} - Γ_{β, λ}^{σ} Γ_{μ, σ}^{λ}) g_{α, ν}}{2} - \frac{(\partial_{λ} (Γ_{β, ν}^{λ}) - \partial_{ν} (Γ_{β, λ}^{λ}) + Γ_{β, ν}^{σ} Γ_{λ, σ}^{λ} - Γ_{β, λ}^{σ} Γ_{ν, σ}^{λ}) g_{α, μ}}{2} + \frac{g_{}^{χ, κ} (\partial_{σ} (Γ_{χ, κ}^{σ}) - \partial_{χ} (Γ_{κ, σ}^{σ}) + Γ_{χ, κ}^{τ} Γ_{σ, τ}^{σ} - Γ_{κ, σ}^{τ} Γ_{χ, τ}^{σ}) (g_{α, μ} g_{β, ν} - g_{α, ν} g_{β, μ})}{6}$

(23)

>	$= convert (, g_)$

$C_{α, β, μ, ν} = g_{α, ω} (\frac{\partial_{μ} (g_{}^{λ, ω}) (\partial_{ν} (g_{β, λ}) + \partial_{β} (g_{λ, ν}) - \partial_{λ} (g_{β, ν}))}{2} + \frac{g_{}^{λ, ω} (\partial_{μ} (\partial_{ν} (g_{β, λ})) + \partial_{β} (\partial_{μ} (g_{λ, ν})) - \partial_{λ} (\partial_{μ} (g_{β, ν})))}{2} - \frac{\partial_{ν} (g_{}^{κ, ω}) (\partial_{μ} (g_{β, κ}) + \partial_{β} (g_{κ, μ}) - \partial_{κ} (g_{β, μ}))}{2} - \frac{g_{}^{κ, ω} (\partial_{μ} (\partial_{ν} (g_{β, κ})) + \partial_{β} (\partial_{ν} (g_{κ, μ})) - \partial_{κ} (\partial_{ν} (g_{β, μ})))}{2} + \frac{g_{}^{ω, σ} (\partial_{ζ} (g_{μ, σ}) + \partial_{μ} (g_{σ, ζ}) - \partial_{σ} (g_{μ, ζ})) g_{}^{χ, ζ} (\partial_{ν} (g_{β, χ}) + \partial_{β} (g_{χ, ν}) - \partial_{χ} (g_{β, ν}))}{4} - \frac{g_{}^{ω, τ} (\partial_{ζ} (g_{ν, τ}) + \partial_{ν} (g_{τ, ζ}) - \partial_{τ} (g_{ν, ζ})) g_{}^{υ, ζ} (\partial_{μ} (g_{β, υ}) + \partial_{β} (g_{μ, υ}) - \partial_{υ} (g_{β, μ}))}{4}) - \frac{(\frac{\partial_{λ} (g_{}^{λ, τ}) (\partial_{μ} (g_{α, τ}) + \partial_{α} (g_{μ, τ}) - \partial_{τ} (g_{α, μ}))}{2} + \frac{g_{}^{λ, τ} (\partial_{λ} (\partial_{μ} (g_{α, τ})) + \partial_{α} (\partial_{λ} (g_{μ, τ})) - \partial_{λ} (\partial_{τ} (g_{α, μ})))}{2} - \frac{\partial_{μ} (g_{}^{κ, λ}) \partial_{α} (g_{κ, λ})}{2} - \frac{g_{}^{κ, λ} \partial_{α} (\partial_{μ} (g_{κ, λ}))}{2} + \frac{g_{}^{ω, σ} (\partial_{μ} (g_{α, ω}) + \partial_{α} (g_{μ, ω}) - \partial_{ω} (g_{α, μ})) g_{}^{λ, υ} \partial_{σ} (g_{λ, υ})}{4} - \frac{g_{}^{ψ, σ} (\partial_{λ} (g_{α, ψ}) + \partial_{α} (g_{λ, ψ}) - \partial_{ψ} (g_{α, λ})) g_{}^{χ, λ} (\partial_{σ} (g_{χ, μ}) + \partial_{μ} (g_{χ, σ}) - \partial_{χ} (g_{μ, σ}))}{4}) g_{β, ν}}{2} + \frac{(\frac{\partial_{λ} (g_{}^{λ, τ}) (\partial_{ν} (g_{α, τ}) + \partial_{α} (g_{ν, τ}) - \partial_{τ} (g_{α, ν}))}{2} + \frac{g_{}^{λ, τ} (\partial_{λ} (\partial_{ν} (g_{α, τ})) + \partial_{α} (\partial_{λ} (g_{ν, τ})) - \partial_{λ} (\partial_{τ} (g_{α, ν})))}{2} - \frac{\partial_{ν} (g_{}^{κ, λ}) \partial_{α} (g_{κ, λ})}{2} - \frac{g_{}^{κ, λ} \partial_{α} (\partial_{ν} (g_{κ, λ}))}{2} + \frac{g_{}^{ω, σ} (\partial_{ν} (g_{α, ω}) + \partial_{α} (g_{ν, ω}) - \partial_{ω} (g_{α, ν})) g_{}^{λ, υ} \partial_{σ} (g_{λ, υ})}{4} - \frac{g_{}^{ψ, σ} (\partial_{λ} (g_{α, ψ}) + \partial_{α} (g_{λ, ψ}) - \partial_{ψ} (g_{α, λ})) g_{}^{χ, λ} (\partial_{σ} (g_{χ, ν}) + \partial_{ν} (g_{χ, σ}) - \partial_{χ} (g_{ν, σ}))}{4}) g_{β, μ}}{2} + \frac{(\frac{\partial_{λ} (g_{}^{λ, τ}) (\partial_{μ} (g_{β, τ}) + \partial_{β} (g_{μ, τ}) - \partial_{τ} (g_{β, μ}))}{2} + \frac{g_{}^{λ, τ} (\partial_{λ} (\partial_{μ} (g_{β, τ})) + \partial_{β} (\partial_{λ} (g_{μ, τ})) - \partial_{λ} (\partial_{τ} (g_{β, μ})))}{2} - \frac{\partial_{μ} (g_{}^{κ, λ}) \partial_{β} (g_{κ, λ})}{2} - \frac{g_{}^{κ, λ} \partial_{β} (\partial_{μ} (g_{κ, λ}))}{2} + \frac{g_{}^{ω, σ} (\partial_{μ} (g_{β, ω}) + \partial_{β} (g_{μ, ω}) - \partial_{ω} (g_{β, μ})) g_{}^{λ, υ} \partial_{σ} (g_{λ, υ})}{4} - \frac{g_{}^{ψ, σ} (\partial_{λ} (g_{β, ψ}) + \partial_{β} (g_{λ, ψ}) - \partial_{ψ} (g_{β, λ})) g_{}^{χ, λ} (\partial_{σ} (g_{χ, μ}) + \partial_{μ} (g_{χ, σ}) - \partial_{χ} (g_{μ, σ}))}{4}) g_{α, ν}}{2} - \frac{(\frac{\partial_{λ} (g_{}^{λ, τ}) (\partial_{ν} (g_{β, τ}) + \partial_{β} (g_{ν, τ}) - \partial_{τ} (g_{β, ν}))}{2} + \frac{g_{}^{λ, τ} (\partial_{λ} (\partial_{ν} (g_{β, τ})) + \partial_{β} (\partial_{λ} (g_{ν, τ})) - \partial_{λ} (\partial_{τ} (g_{β, ν})))}{2} - \frac{\partial_{ν} (g_{}^{κ, λ}) \partial_{β} (g_{κ, λ})}{2} - \frac{g_{}^{κ, λ} \partial_{β} (\partial_{ν} (g_{κ, λ}))}{2} + \frac{g_{}^{ω, σ} (\partial_{ν} (g_{β, ω}) + \partial_{β} (g_{ν, ω}) - \partial_{ω} (g_{β, ν})) g_{}^{λ, υ} \partial_{σ} (g_{λ, υ})}{4} - \frac{g_{}^{ψ, σ} (\partial_{λ} (g_{β, ψ}) + \partial_{β} (g_{λ, ψ}) - \partial_{ψ} (g_{β, λ})) g_{}^{χ, λ} (\partial_{σ} (g_{χ, ν}) + \partial_{ν} (g_{χ, σ}) - \partial_{χ} (g_{ν, σ}))}{4}) g_{α, μ}}{2} + \frac{g_{}^{χ, κ} (\frac{\partial_{σ} (g_{}^{λ, σ}) (\partial_{κ} (g_{χ, λ}) + \partial_{χ} (g_{κ, λ}) - \partial_{λ} (g_{χ, κ}))}{2} + \frac{g_{}^{λ, σ} (\partial_{κ} (\partial_{σ} (g_{χ, λ})) + \partial_{χ} (\partial_{σ} (g_{κ, λ})) - \partial_{λ} (\partial_{σ} (g_{χ, κ})))}{2} - \frac{\partial_{χ} (g_{}^{ψ, σ}) \partial_{κ} (g_{ψ, σ})}{2} - \frac{g_{}^{ψ, σ} \partial_{χ} (\partial_{κ} (g_{ψ, σ}))}{2} + \frac{g_{}^{τ, ζ} (\partial_{κ} (g_{χ, ζ}) + \partial_{χ} (g_{κ, ζ}) - \partial_{ζ} (g_{χ, κ})) g_{}^{ω, σ} \partial_{τ} (g_{ω, σ})}{4} - \frac{g_{}^{ι, τ} (\partial_{σ} (g_{ι, κ}) + \partial_{κ} (g_{ι, σ}) - \partial_{ι} (g_{κ, σ})) g_{}^{σ, υ} (\partial_{τ} (g_{χ, υ}) + \partial_{χ} (g_{τ, υ}) - \partial_{υ} (g_{χ, τ}))}{4}) (g_{α, μ} g_{β, ν} - g_{α, ν} g_{β, μ})}{6}$

(24)

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

>	$Weyl [1, 2, 2, 4]$

$\frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}$

(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]$

$C_{2, 2, 4}^{1}$

(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 ~)

>	$Weyl [nonzero]$

$C_{α, β, μ, ν} = \{(1, 2, 2, 4) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (1, 2, 4, 2) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (1, 3, 3, 4) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (1, 3, 4, 3) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (1, 4, 1, 4) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}}, (1, 4, 4, 1) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}}, (2, 1, 2, 4) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (2, 1, 4, 2) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (2, 3, 2, 3) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) r^{2} {\sin (θ)}^{2}}{3}, (2, 3, 3, 2) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) r^{2} {\sin (θ)}^{2}}{3}, (2, 4, 1, 2) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (2, 4, 2, 1) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (2, 4, 2, 4) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3}, (2, 4, 4, 2) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3}, (3, 1, 3, 4) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (3, 1, 4, 3) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (3, 2, 2, 3) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) r^{2} {\sin (θ)}^{2}}{3}, (3, 2, 3, 2) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) r^{2} {\sin (θ)}^{2}}{3}, (3, 4, 1, 3) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (3, 4, 3, 1) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (3, 4, 3, 4) = - \frac{2 k^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) {\sin (θ)}^{2}}{3}, (3, 4, 4, 3) = \frac{2 k^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) {\sin (θ)}^{2}}{3}, (4, 1, 1, 4) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}}, (4, 1, 4, 1) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}}, (4, 2, 1, 2) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (4, 2, 2, 1) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (4, 2, 2, 4) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3}, (4, 2, 4, 2) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3}, (4, 3, 1, 3) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (4, 3, 3, 1) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (4, 3, 3, 4) = \frac{2 k^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) {\sin (θ)}^{2}}{3}, (4, 3, 4, 3) = - \frac{2 k^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) {\sin (θ)}^{2}}{3}\}$

(29)

>	$Weyl [~alpha, β, μ, ν, nonzero]$

$C_{β, μ, ν}^{α} = \{(1, 1, 1, 4) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}}, (1, 1, 4, 1) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}}, (1, 2, 1, 2) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (1, 2, 2, 1) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (1, 3, 1, 3) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (1, 3, 3, 1) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (1, 4, 1, 4) = \frac{4 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (1, 4, 4, 1) = - \frac{4 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (2, 1, 2, 4) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, (2, 1, 4, 2) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{6 r^{2}}, (2, 3, 2, 3) = - \frac{2 {\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (2, 3, 3, 2) = \frac{2 {\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (2, 4, 1, 2) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{6 r^{2}}, (2, 4, 2, 1) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, (2, 4, 2, 4) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (2, 4, 4, 2) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (3, 1, 3, 4) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, (3, 1, 4, 3) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{6 r^{2}}, (3, 2, 2, 3) = \frac{2 {k_{r}}^{2} r^{2}}{3} + \frac{2 k_{r, r} k r^{2}}{3} - \frac{4 k r k_{r}}{3} + \frac{2 k^{2}}{3} - \frac{1}{3}, (3, 2, 3, 2) = - \frac{2 {k_{r}}^{2} r^{2}}{3} - \frac{2 k_{r, r} k r^{2}}{3} + \frac{4 k r k_{r}}{3} - \frac{2 k^{2}}{3} + \frac{1}{3}, (3, 4, 1, 3) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{6 r^{2}}, (3, 4, 3, 1) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{6 r^{2}}, (3, 4, 3, 4) = - \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (3, 4, 4, 3) = \frac{2 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}}, (4, 2, 2, 4) = - \frac{{k_{r}}^{2} r^{2}}{3} - \frac{k_{r, r} k r^{2}}{3} + \frac{2 k r k_{r}}{3} - \frac{k^{2}}{3} + \frac{1}{6}, (4, 2, 4, 2) = \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6}, (4, 3, 3, 4) = - \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (4, 3, 4, 3) = \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3}, (4, 4, 1, 4) = \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}}, (4, 4, 4, 1) = \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}}\}$

(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 $C_{1}^{1}$ we have

>	$C [1, 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]$

$C_{β, 1, ν}^{1} = [\begin{array}{c} 0 & 0 & 0 & \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}} \\ 0 & \frac{{k_{r}}^{2} r^{2}}{3} + \frac{k_{r, r} k r^{2}}{3} - \frac{2 k r k_{r}}{3} + \frac{k^{2}}{3} - \frac{1}{6} & 0 & 0 \\ 0 & 0 & \frac{{\sin (θ)}^{2} (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2})}{3} & 0 \\ 0 & 0 & 0 & \frac{4 (- \frac{1}{2} + {k_{r}}^{2} r^{2} + k_{r, r} k r^{2} - 2 k r k_{r} + k^{2}) k^{2}}{3 r^{2}} \end{array}]$

(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]$

$C_{β, 1, ν}^{1} + C_{β, 2, ν}^{2} + C_{β, 3, ν}^{3} + C_{β, 4, ν}^{4} = [\begin{array}{c} 0 & 0 & 0 & \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}} + \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}} \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ \frac{2 {k_{r}}^{2} r^{2} + 2 k_{r, r} k r^{2} - 4 k r k_{r} + 2 k^{2} - 1}{3 r^{2}} + \frac{- 2 {k_{r}}^{2} r^{2} - 2 k_{r, r} k r^{2} + 4 k r k_{r} - 2 k^{2} + 1}{3 r^{2}} & 0 & 0 & 0 \end{array}]$

(34)

>	$simplify ()$

$C_{β, 1, ν}^{1} + C_{β, 2, ν}^{2} + C_{β, 3, ν}^{3} + C_{β, 4, ν}^{4} = [\begin{array}{c} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{array}]$

(35)

In different contexts it is sometimes necessary to classify the spacetime, for instance performing a PetrovType classification, or the SegreType classification of the Ricci tensor. These classification commands are in the Tetrads package. For example, for the spacetime considered, we have

>	$Tetrads :- PetrovType ()$

$D$

(36)

>	$Tetrads :- SegreType ()$

$D, [2,(11)]$

(37)

	See Also
	`.`, Array, ArrayElems, Christoffel, CompactDisplay, Coordinates, DifferentialGeometry[Tensor][WeylTensor], Einstein, g_, l_, m_, mb_, n_, Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Ricci, Riemann, Setup, Tetrads, TransformTetrad, Typesetting, value

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.

Compatibility

•	The Physics[Weyl] command was introduced in Maple 16.

•	For more information on Maple 16 changes, see Updates in Maple 16.

•	The Physics[Weyl] and Physics[Tetrads][WeylScalars] commands were updated in Maple 2017.

•	The Physics[Tetrads][WeylScalars] command was introduced in Maple 2021.

•	For more information on Maple 2021 changes, see Updates in Maple 2021.

Download Help Document

Maple

Maple Add-Ons

MapleSim

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim

Maple

Powerful math software that is easy to use

Maple Add-Ons

MapleSim

Advanced System Level Modeling

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

Online Help

All Products Maple MapleSim