IsTetrad - Maple Help

All Products Maple MapleSim

Home : Support : Online Help : Physics : Tetrads : IsTetrad

Physics[Tetrads][NullTetrad] - compute different forms of a null tetrad for the current spacetime metric

Physics[Tetrads][OrthonormalTetrad] - compute different forms of an orthonormal tetrad for the current spacetime metric

Physics[Tetrads][IsTetrad] - return true or false according to whether a given matrix is a null or orthonormal tetrad and indicate what kind of the tetrad is

	Calling Sequence
	NullTetrad(method = ..., firstvector = ...) OrthonormalTetrad(method = ..., firstvector = ...) IsTetrad(T, orthonormal, null, quiet)

Parameters

method = ..	-	optional, indicates the method to be used, that could be GramSchmidt or Eigenvectors
firstvector = ..	-	optional, a Vector or a list with the vector components, as many as the dimension of spacetime; indicates the first departing vector to be used with the GramSchmidt method
T	-	a Matrix representing a tetrad, it is also optionally the first argument in calls to NullTetrad and OrthonormalTetrad
orthonormal	-	optional, to indicate to IsTetrad to check only for an orthonormal tetrad
null	-	optional, to indicate to IsTetrad to check only for an null tetrad
quiet	-	optional, to indicate to IsTetrad to avoid displaying the type of tetrad on the screen

Description

NullTetrad and OrthonormalTetrad compute tetrads for the spacetime metric set. Recalling, given the metric of a local (tetrad) system of references, there are infinitely many tetrads (represented in Physics by the Tetrads:-e_ command) relating the components of a tensor in the global (spacetime) and local (tetrad) systems of references. These two commands NullTetrad and OrthonormalTetrad compute different tetrads, that may result more convenient in different contexts, for example when computing the Weyl and Ricci scalars or performing a Petrov classification.

•	When the tetrad you get from NullTetrad or OrthonormalTetrad is the one you prefer, say E, you can set it, to be used by the system as the current tetrad, by entering $Setup (tetrad = E)$ . You can always see the components of the current tetrad by entering e_[].

•

When called with no arguments, both NullTetrad and OrthonormalTetrad will depart from the current tetrad and compute a corresponding - null or orthonormal - tetrad. The formula used depends on the signature. For the signatures (---+) and (+++-), given an orthonormal tetrad E, NullTetrad returns a null tetrad constructed as M . E, where

$M = [\begin{array}{c} - \frac{1}{2} \sqrt{2} & 0 & 0 & \frac{1}{2} \sqrt{2} \\ 0 & \frac{1}{2} \sqrt{2} & \frac{1}{2} I \sqrt{2} & 0 \\ 0 & \frac{1}{2} \sqrt{2} & - \frac{1}{2} I \sqrt{2} & 0 \\ \frac{1}{2} \sqrt{2} & 0 & 0 & \frac{1}{2} \sqrt{2} \end{array}]$

and, for the signatures (+---) and (-+++),

$M = [\begin{array}{c} \frac{1}{2} \sqrt{2} & \frac{1}{2} \sqrt{2} & 0 & 0 \\ \frac{1}{2} \sqrt{2} & - \frac{1}{2} \sqrt{2} & 0 & 0 \\ 0 & 0 & \frac{1}{2} \sqrt{2} & \frac{1}{2} I \sqrt{2} \\ 0 & 0 & \frac{1}{2} \sqrt{2} & - \frac{1}{2} I \sqrt{2} \end{array}]$

Given a null tetrad E, OrthonormalTetrad returns an orthonormal tetrad constructed as $M^{- 1} . E$ .

•	Both NullTetrad and OrthonormalTetrad can also be called with a first argument, T, a matrix representing a tetrad, in which case instead of departing from the current components of Tetrads:-e_, they apply the formulas just described to T.

•	Additionally, instead of constructing tetrads multiplying by M or its inverse, NullTetrad and OrthonormalTetrad can construct tetrads from scratch using two different methods: GramSchmidt or Eigenvectors - see the Examples section.

•

Independently of NullTetrad and OrthonormalTetrad, you can always set the tetrad to have any preferred matrix form, say E, by entering $Setup (tetrad = E)$ , or transform a tetrad into another one using TransformTetrad. Of particular interest is the the ability of TransformTetrad to compute a tetrad in canonicalform, that is a form that results in a canonical form of the Weyl scalars.

•	IsTetrad is a complementary command: it returns true or false depending on whether the given matrix T is a tetrad, optionally indicating the type of tetrad.

Examples

>	$with (Physics) &colon;$ $Setup (mathematicalnotation = true) &semi;$ $with (Tetrads)$

$[mathematicalnotation = true]$

$_______________________________________________________$

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

(1)

Set the spacetime metric to something non-flat, for instance take the metric [13,7,5] of the Exact solutions book referenced at the end of this page:

>	$g_[[13, 7, 5]]$

$_______________________________________________________$

$Systems of spacetime coordinates are: \{X = (x, y, z, t)\}$

$Default differentiation variables for d_, D_ and dAlembertian are: \{X = (x, y, z, t)\}$

$The metric in coordinates [x, y, z, t]$

$Parameters: [ε, A (t), B (t), A1]$

$Comments: &ExponentialE;psⅈlon=1 or &ExponentialE;psⅈlon=-1$

$Domain: [0 < t]$

$Assumptions: [A (t) :: real, B (t) :: real, 0 < A1, 0 < x, 0 < A (t), 0 < B (t)]$

$Resetting the signature of spacetime from (- - - +) to (+ + + -) in order to match the signature in the database of metrics$

$_______________________________________________________$

$Setting lowercaselatin_is letters to represent space indices$

$g_{μ, ν} = [\begin{array}{c} ε {A (t)}^{2} & 0 & 0 & 0 \\ 0 & {B (t)}^{2} {&ExponentialE;}^{2 A1 x} (2 {\cosh (x)}^{2} - 1) & - 2 {B (t)}^{2} {&ExponentialE;}^{2 A1 x} \cosh (x) \sinh (x) & 0 \\ 0 & - 2 {B (t)}^{2} {&ExponentialE;}^{2 A1 x} \cosh (x) \sinh (x) & {B (t)}^{2} {&ExponentialE;}^{2 A1 x} (2 {\cosh (x)}^{2} - 1) & 0 \\ 0 & 0 & 0 & - ε \end{array}]$

(2)

Compact the display of the functions $A (t), B (t)$ :

>	$CompactDisplay (A (t), B (t))$

$A (t) will now be displayed as A$

$B (t) will now be displayed as B$

(3)

The default orthonormal tetrad is:

$e_[]$

$𝔢_{a, μ} = [\begin{array}{c} A \sqrt{ε} & 0 & 0 & 0 \\ 0 & B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x} & - \frac{B \sinh (2 x) {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \\ 0 & 0 & \frac{B {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \\ 0 & 0 & 0 & - \sqrt{ε} \end{array}]$

(4)

It is possible to test whether these tetrads satisfy the tetrad definition

>	$e_[definition]$

$𝔢_{a, μ} 𝔢_{b}^{μ} = η_{a, b}$

(5)

using IsTetrad

>	$IsTetrad ()$

$Type of tetrad: orthonormal$

$true$

(6)

To compute with a null tetrad instead of an orthonormal tetrad, set the tetrad or tetradmetric to null using Setup.

>	$Setup (tetradmetric = null)$

$[tetradmetric = \{(1, 4) = −1, (2, 3) = 1\}]$

(7)

Note: this form of the null tetrad metric is consistent with time in position 4 (default) - and not with time in position 1. To see the more frequent form of the tetradmetric

>	$Matrix (4, \{(1, 2) = - 1, (3, 4) = 1\}, shape = symmetric)$

$[\begin{array}{c} 0 & −1 & 0 & 0 \\ −1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 1 & 0 \end{array}]$

(8)

you can either set the signature of the underlying inertial system of references (Minkowski) indicating time (different sign) in the first place, or set the tetradmetric directly using Setup(tetradmetric = the_matrix_above).

So now the metric of the local (tetrad) system of references is:

$eta_[]$

$η_{a, b} = [\begin{array}{c} 0 & 0 & 0 & −1 \\ 0 & 0 & 1 & 0 \\ 0 & 1 & 0 & 0 \\ −1 & 0 & 0 & 0 \end{array}]$

(9)

and the null tetrad is:

$e_[]$

$𝔢_{a, μ} = [\begin{array}{c} - \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & - \frac{\sqrt{2} \sqrt{ε}}{2} \\ 0 & \frac{\sqrt{2} B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x}}{2} & \frac{\sqrt{2} B {&ExponentialE;}^{A1 x} (- \sinh (2 x) + I)}{2 \sqrt{\cosh (2 x)}} & 0 \\ 0 & \frac{\sqrt{2} B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x}}{2} & - \frac{\sqrt{2} B {&ExponentialE;}^{A1 x} (\sinh (2 x) + I)}{2 \sqrt{\cosh (2 x)}} & 0 \\ \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & - \frac{\sqrt{2} \sqrt{ε}}{2} \end{array}]$

(10)

>	$IsTetrad ()$

$Type of tetrad: null$

$true$

(11)

You can compute different forms of an orthonormal tetrad or a null tetrad using the OrthonormalTetrad and NullTetrad commands, requesting the use of different methods, or passing a starting vector. For instance, this is the default form of an orthonormal tetrad for the current spacetime metric:

>	$OrthonormalTetrad ()$

$[\begin{array}{c} A \sqrt{ε} & 0 & 0 & 0 \\ 0 & B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x} & - \frac{B \sinh (2 x) {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \\ 0 & 0 & \frac{B {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \\ 0 & 0 & 0 & - \sqrt{ε} \end{array}]$

(12)

There are two methods available for computing these tetrads: the GramSchmidt (generalization to curved spaces) and Eigenvectors (based on computing the Eigenvectors of a related Matrix). By default, the Physics package routines decide on what method is more convenient to use for the spacetime metric set. In this case the routines use the GramSchmidt method; you can check the other method via:

>	$OrthonormalTetrad (method = Eigenvectors)$

$[\begin{array}{c} A \sqrt{ε} & 0 & 0 & 0 \\ 0 & - \frac{B \sqrt{2} {&ExponentialE;}^{A1 x}}{2 \sqrt{\cosh (2 x) - \sinh (2 x)}} & \frac{B \sqrt{2} {&ExponentialE;}^{A1 x}}{2 \sqrt{\cosh (2 x) - \sinh (2 x)}} & 0 \\ 0 & \frac{B {&ExponentialE;}^{A1 x} \sqrt{2}}{2 \sqrt{\cosh (2 x) + \sinh (2 x)}} & \frac{B {&ExponentialE;}^{A1 x} \sqrt{2}}{2 \sqrt{\cosh (2 x) + \sinh (2 x)}} & 0 \\ 0 & 0 & 0 & \sqrt{ε} \end{array}]$

(13)

If this is an orthonormal tetrad you prefer, you can set it via Setup(e_ = ...), where ... is the matrix you prefer (see the previous output). To see what the null tetrads corresponding to these two orthonormal tetrads would be:

>	$NullTetrad ()$

$[\begin{array}{c} - \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & - \frac{\sqrt{2} \sqrt{ε}}{2} \\ 0 & \frac{\sqrt{2} B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x}}{2} & \frac{\sqrt{2} B {&ExponentialE;}^{A1 x} (- \sinh (2 x) + I)}{2 \sqrt{\cosh (2 x)}} & 0 \\ 0 & \frac{\sqrt{2} B \sqrt{\cosh (2 x)} {&ExponentialE;}^{A1 x}}{2} & - \frac{\sqrt{2} B {&ExponentialE;}^{A1 x} (\sinh (2 x) + I)}{2 \sqrt{\cosh (2 x)}} & 0 \\ \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & - \frac{\sqrt{2} \sqrt{ε}}{2} \end{array}]$

(14)

>	$NullTetrad (method = Eigenvectors)$

$[\begin{array}{c} - \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & \frac{\sqrt{2} \sqrt{ε}}{2} \\ 0 & \frac{B {&ExponentialE;}^{A1 x} (I \sqrt{\cosh (2 x) - \sinh (2 x)} - \sqrt{\cosh (2 x) + \sinh (2 x)})}{2 \sqrt{\cosh (2 x) - \sinh (2 x)} \sqrt{\cosh (2 x) + \sinh (2 x)}} & \frac{B {&ExponentialE;}^{A1 x} (I \sqrt{\cosh (2 x) - \sinh (2 x)} + \sqrt{\cosh (2 x) + \sinh (2 x)})}{2 \sqrt{\cosh (2 x) - \sinh (2 x)} \sqrt{\cosh (2 x) + \sinh (2 x)}} & 0 \\ 0 & - \frac{B {&ExponentialE;}^{A1 x} (I \sqrt{\cosh (2 x) - \sinh (2 x)} + \sqrt{\cosh (2 x) + \sinh (2 x)})}{2 \sqrt{\cosh (2 x) - \sinh (2 x)} \sqrt{\cosh (2 x) + \sinh (2 x)}} & - \frac{B {&ExponentialE;}^{A1 x} (I \sqrt{\cosh (2 x) - \sinh (2 x)} - \sqrt{\cosh (2 x) + \sinh (2 x)})}{2 \sqrt{\cosh (2 x) - \sinh (2 x)} \sqrt{\cosh (2 x) + \sinh (2 x)}} & 0 \\ \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & 0 & \frac{\sqrt{2} \sqrt{ε}}{2} \end{array}]$

(15)

When using the GramSchmidt method you can also specify the first vector. Recalling, the method works iterating the computation of the next vector (line of the tetrad) starting from a first vector. For example, instead of starting from $[\sqrt{ε} A, 0, 0, 0]$ , start from $[0, 0, 0, 1]$ compare with default result of OrthonormalTetrad lines above.

>	$OrthonormalTetrad (firstvector = [0, 0, 0, 1])$

$[\begin{array}{c} A \sqrt{ε} & 0 & 0 & 0 \\ 0 & 0 & 0 & \sqrt{- ε} \\ 0 & - \frac{B {&ExponentialE;}^{A1 x}}{\sqrt{sech (2 x)}} & \frac{B {&ExponentialE;}^{A1 x} (- csch (4 x) + \coth (4 x))}{\sqrt{sech (2 x)}} & 0 \\ 0 & 0 & \frac{−I B {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \end{array}]$

(16)

>	$IsTetrad ()$

$Type of tetrad: orthonormal$

$true$

(17)

>	$NullTetrad (firstvector = [0, 0, 0, 1])$

$[\begin{array}{c} - \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & \frac{- \frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \\ 0 & \frac{- \frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x}}{\sqrt{sech (2 x)}} & \frac{- \frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x} (csch (4 x) - \coth (4 x))}{\sqrt{sech (2 x)}} & \frac{I}{2} \sqrt{2} \sqrt{ε} \\ 0 & \frac{\frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x}}{\sqrt{sech (2 x)}} & \frac{\frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x} (csch (4 x) - \coth (4 x))}{\sqrt{sech (2 x)}} & \frac{I}{2} \sqrt{2} \sqrt{ε} \\ \frac{\sqrt{2} A \sqrt{ε}}{2} & 0 & \frac{- \frac{I}{2} \sqrt{2} B {&ExponentialE;}^{A1 x}}{\sqrt{\cosh (2 x)}} & 0 \end{array}]$

(18)

>	$IsTetrad ()$

$Type of tetrad: null$

$true$

(19)

	See Also
	d_, D_, e_, eta_, g_, gamma_, l_, lambda_, m_, mb_, n_, Physics, Physics conventions, Physics examples, Physics Updates, Tensors - a complete guide, Mini-Course Computer Algebra for Physicists, Ricci, Setup, Tetrads, TransformTetrad, Weyl

References

Textbooks

Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2, fourth revised English edition. Elsevier, 1975.

Stephani, H., Kramer, D., MacCallum, M., Hoenselaers, C. Herlt, E. Exact Solutions of Einstein's Field Equations, Cambridge Monographs on Mathematical Physics, second edition. Cambridge University Press, 2003.

Compatibility

•	The Physics[Tetrads][NullTetrad], Physics[Tetrads][OrthonormalTetrad] and Physics[Tetrads][IsTetrad] commands were introduced in Maple 2015.

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

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