Physics for Maple 2015 - Maple Help

 Physics

Maple provides a state-of-the-art environment for algebraic computations in Physics, with emphasis on ensuring that the computational experience is as natural as possible. The theme of the Physics project for Maple 2015 has been Vector Analysis, symbolic Tensor manipulations, Quantum Mechanics, and General Relativity. With more than 400 enhancements throughout the entire package to increase robustness and versatility, two new commands, Assume and SubstituteTensor, a new Tetrads subpackage with 13 commands, as well as 26 new Physics:-Library commands to support further explorations and extensions, and an enlargement of the database of solutions to Einstein's equations with more than 100 new metrics, Maple 2015 extends again the range of Physics-related algebraic computations that can be done, in a natural way, using computer algebra software.

As part of its commitment to providing the best possible environment for algebraic computations in Physics, Maplesoft launched a Maple Physics: Research and Development web site with Maple 18, which enabled users to download research versions, ask questions, and provide feedback. The results from this accelerated exchange with people around the world have been incorporated into the Physics package in Maple 2015. Examples illustrating the use of the new capabilities in the context of more general problems are found in the MaplePrimes post Computer Algebra for Theoretical Physics.

Simplification

Simplification is perhaps the most common operation performed in a computer algebra system. In Physics, this typically entails simplifying tensorial expressions, or expressions involving noncommutative operators that satisfy certain commutator/anticommutator rules, or sums and integrals involving quantum operators and Dirac delta functions in the summands and integrands. Relevant enhancements were introduced in Maple 2015 for all these cases, including enhancements in the simplification of:

 • Products of LeviCivita tensors in curved spacetimes when LeviCivita represents the Galilean pseudo-tensor (related to Setup(levicivita = Galilean)), instead of its generalization to curved spaces (related to Setup(levicivita = nongalilean)).
 • Tensorial expressions in general that have spacetime, space, and/or tetrad contracted indices, possibly at the same time.
 • New option tryhard, that resolves zero recognition in an important number of nontrivial situations.
 • Expressions involving the Dirac function.
 • Vectorial expressions involving cylindrical or spherical coordinates and related unit vectors.
 • Expressions simplified with respect to side relations (equations) in the presence of quantum vectorial equations.
 • Expressions involving products of quantum operators entering parameterized algebra rules.
 • Expressions involving vectorial quantum operators simplified with respect to other vectorial equations.
 • Add support for the simplification and integration of spherical harmonics (SphericalY ) relevant in quantum mechanics.

Examples

 • Enhancements in the simplification of tensorial expressions.
 >
 ${\mathrm{* Partial match of \text{'}notation\text{'} against keyword \text{'}mathematicalnotation\text{'}}}$
 ${\mathrm{* Partial match of \text{'}coordinates\text{'} against keyword \text{'}coordinatesystems\text{'}}}$
 ${\mathrm{Default differentiation variables for d_, D_ and dAlembertian are:}}\left\{{X}{=}\left({\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right)\right\}$
 ${\mathrm{Systems of spacetime Coordinates are:}}\left\{{X}{=}\left({\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right)\right\}$
 $\left[{\mathrm{coordinatesystems}}{=}\left\{{X}\right\}{,}{\mathrm{mathematicalnotation}}{=}{\mathrm{true}}\right]$ (1)
 • Define a tensor ${l}_{\mathrm{\mu }}\left(X\right)$ and define the metric in terms of the components of this tensor so that ${g}_{\mathrm{\mu },\mathrm{\nu }}={\mathrm{\eta }}_{\mathrm{\mu },\mathrm{\nu }}+{l}_{\mathrm{\mu }}{l}_{\mathrm{\nu }}$ and ${\mathrm{\eta }}_{\mathrm{\mu },\mathrm{\nu }}$ is a Minkowski metric.
 > ${l}_{\mathrm{\mu }},{\mathrm{\eta }}_{\mathrm{\mu },\mathrm{\nu }}=-\mathrm{rhs}\left({\mathrm{g_}}_{\mathrm{Minkowski}}\right)$
 ${\mathrm{The Minkowski metric in cartesian coordinates}}$
 >
 ${\mathrm{Defined objects with tensor properties}}$
 $\left\{{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{X}}_{{\mathrm{\mu }}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{\mathrm{η}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{l}}_{{\mathrm{\mu }}}{,}{{\mathrm{\delta }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (2)
 • Avoid redundant display:
 > $\mathrm{PDEtools}:-\mathrm{declare}\left(l\left(X\right)\right)$
 ${l}\left({\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right){}{\mathrm{will now be displayed as}}{}{l}$ (3)
 • The metric:
 >
 ${{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{\mathrm{η}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{+}{{l}}_{{\mathrm{\mu }}}{}{{l}}_{{\mathrm{\nu }}}$ (4)
 • New: you can define it directly using a tensorial equation like (4), using either Define or Setup.
 > $\mathrm{Define}\left(\right)$
 ${\mathrm{Defined objects with tensor properties}}$
 $\left\{{{𝒟}}_{{\mathrm{\mu }}}{,}{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{C}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{X}}_{{\mathrm{\mu }}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{\mathrm{η}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{l}}_{{\mathrm{\mu }}}{,}{{\mathrm{\Gamma }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}}{,}{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\delta }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (5)
 • Verify the resulting ${g}_{\mathrm{\mu },\mathrm{\nu }}$:
 > $\mathrm{g_}\left[\right]$
 ${{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}\left[\begin{array}{cccc}{1}{+}{{{l}}_{{1}}}^{{2}}& {{l}}_{{1}}{}{{l}}_{{2}}& {{l}}_{{1}}{}{{l}}_{{3}}& {{l}}_{{1}}{}{{l}}_{{4}}\\ {{l}}_{{1}}{}{{l}}_{{2}}& {1}{+}{{{l}}_{{2}}}^{{2}}& {{l}}_{{2}}{}{{l}}_{{3}}& {{l}}_{{2}}{}{{l}}_{{4}}\\ {{l}}_{{1}}{}{{l}}_{{3}}& {{l}}_{{2}}{}{{l}}_{{3}}& {1}{+}{{{l}}_{{3}}}^{{2}}& {{l}}_{{3}}{}{{l}}_{{4}}\\ {{l}}_{{1}}{}{{l}}_{{4}}& {{l}}_{{2}}{}{{l}}_{{4}}& {{l}}_{{3}}{}{{l}}_{{4}}& {-}{1}{+}{{{l}}_{{4}}}^{{2}}\end{array}\right]$ (6)
 • New: you can query about the definition you gave for any tensor directly from the tensor itself.
 > $\mathrm{g_}\left[\mathrm{definition}\right]$
 ${{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{\mathrm{η}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{+}{{l}}_{{\mathrm{\mu }}}{}{{l}}_{{\mathrm{\nu }}}$ (7)
 • Enhanced simplification capabilities: show that this expression involving derivatives of ${l}_{\mathrm{\mu }}\left(X\right)$ is equal to zero.
 > $-\left({\mathrm{d_}}_{\mathrm{\alpha }}\left({l}_{\mathrm{\lambda }}\left(X\right)\right){\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\beta }}}^{\phantom{}\mathrm{\alpha },\mathrm{\beta }}{\mathrm{d_}}_{\mathrm{\beta }}\left({l}_{\mathrm{\kappa }}\left(X\right)\right){l}_{\phantom{}\phantom{\mathrm{\kappa }}}^{\phantom{}\mathrm{\kappa }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\lambda }}}^{\phantom{}\mathrm{\lambda }}\left(X\right)\right)+\left({\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\beta }}}^{\phantom{}\mathrm{\alpha },\mathrm{\beta }}{\mathrm{d_}}_{\mathrm{\beta }}\left({l}_{\mathrm{\kappa }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\lambda }}\left({l}_{\mathrm{\alpha }}\left(X\right)\right){l}_{\phantom{}\phantom{\mathrm{\kappa }}}^{\phantom{}\mathrm{\kappa }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\lambda }}}^{\phantom{}\mathrm{\lambda }}\left(X\right)\right)+\left({\mathrm{d_}}_{\mathrm{\lambda }}\left({l}_{\mathrm{\kappa }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right)\right){l}_{\phantom{}\phantom{\mathrm{\kappa }}}^{\phantom{}\mathrm{\kappa }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\lambda }}}^{\phantom{}\mathrm{\lambda }}\left(X\right)\right)+\left({\mathrm{d_}}_{\mathrm{\beta }}\left({\mathrm{d_}}_{\mathrm{\lambda }}\left({l}_{\mathrm{\kappa }}\left(X\right)\right)\right){l}_{\phantom{}\phantom{\mathrm{\beta }}}^{\phantom{}\mathrm{\beta }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\kappa }}}^{\phantom{}\mathrm{\kappa }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\lambda }}}^{\phantom{}\mathrm{\lambda }}\left(X\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\mu }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({\mathrm{d_}}_{\mathrm{\mu }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\nu }}\left({l}_{\mathrm{\mu }}\left(X\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({\mathrm{d_}}_{\mathrm{\nu }}\left({l}_{\mathrm{\mu }}\left(X\right)\right)\right)+\left(\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\alpha },\mathrm{\rho }}{\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\mu }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\alpha }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\alpha },\mathrm{\rho }}{\mathrm{d_}}_{\mathrm{\mu }}\left({l}_{\mathrm{\alpha }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)+\left(\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\alpha },\mathrm{\rho }}{\mathrm{d_}}_{\mathrm{\alpha }}\left({l}_{\mathrm{\mu }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){\mathrm{η}}_{\phantom{}\phantom{\mathrm{\alpha }}\phantom{,}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\alpha },\mathrm{\rho }}{\mathrm{d_}}_{\mathrm{\nu }}\left({l}_{\mathrm{\alpha }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\mu }}\left(X\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\alpha }}}^{\phantom{}\mathrm{\alpha }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\nu }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\alpha }}\left({l}_{\mathrm{\mu }}\left(X\right)\right)+\left(\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\alpha }}}^{\phantom{}\mathrm{\alpha }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\mu }}\left({l}_{\mathrm{\alpha }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)\right)-\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\alpha }}}^{\phantom{}\mathrm{\alpha }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\mu }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\alpha }}\left({l}_{\mathrm{\nu }}\left(X\right)\right)+\left(\frac{1}{2}{l}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\nu }}}^{\phantom{}\mathrm{\nu }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\alpha }}}^{\phantom{}\mathrm{\alpha }}\left(X\right){l}_{\phantom{}\phantom{\mathrm{\rho }}}^{\phantom{}\mathrm{\rho }}\left(X\right){\mathrm{d_}}_{\mathrm{\nu }}\left({l}_{\mathrm{\alpha }}\left(X\right)\right){\mathrm{d_}}_{\mathrm{\rho }}\left({l}_{\mathrm{\mu }}\left(X\right)\right)\right)$
 ${-}\left({{\partial }}_{{\mathrm{\alpha }}}\left({{l}}_{{\mathrm{\lambda }}}\right)\right){}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}\left({{\partial }}_{{\mathrm{\beta }}}\left({{l}}_{{\mathrm{\kappa }}}\right)\right){}{{l}}_{\phantom{{}}\phantom{{\mathrm{\kappa }}}}^{\phantom{{}}{\mathrm{\kappa }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\lambda }}}}^{\phantom{{}}{\mathrm{\lambda }}}{+}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}\left({{\partial }}_{{\mathrm{\beta }}}\left({{l}}_{{\mathrm{\kappa }}}\right)\right){}\left({{\partial }}_{{\mathrm{\lambda }}}\left({{l}}_{{\mathrm{\alpha }}}\right)\right){}{{l}}_{\phantom{{}}\phantom{{\mathrm{\kappa }}}}^{\phantom{{}}{\mathrm{\kappa }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\lambda }}}}^{\phantom{{}}{\mathrm{\lambda }}}{+}\left({{\partial }}_{{\mathrm{\lambda }}}\left({{l}}_{{\mathrm{\kappa }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}\right)\right){}{{l}}_{\phantom{{}}\phantom{{\mathrm{\kappa }}}}^{\phantom{{}}{\mathrm{\kappa }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\lambda }}}}^{\phantom{{}}{\mathrm{\lambda }}}{+}\left({{\partial }}_{{\mathrm{\beta }}}\left({{\partial }}_{{\mathrm{\lambda }}}\left({{l}}_{{\mathrm{\kappa }}}\right)\right)\right){}{{l}}_{\phantom{{}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\beta }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\kappa }}}}^{\phantom{{}}{\mathrm{\kappa }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\lambda }}}}^{\phantom{{}}{\mathrm{\lambda }}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}\right)\right){}\left({{\partial }}_{{\mathrm{\mu }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{\partial }}_{{\mathrm{\mu }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}\right)\right){}\left({{\partial }}_{{\mathrm{\nu }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{\partial }}_{{\mathrm{\nu }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right)\right)}{{2}}{+}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right){}\left({{\partial }}_{{\mathrm{\alpha }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\mu }}}\left({{l}}_{{\mathrm{\alpha }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{+}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\alpha }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{\mathrm{η}}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\alpha }}{,}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\nu }}}\left({{l}}_{{\mathrm{\alpha }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}}^{\phantom{{}}{\mathrm{\alpha }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right){}\left({{\partial }}_{{\mathrm{\alpha }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right)}{{2}}{+}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}}^{\phantom{{}}{\mathrm{\alpha }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\mu }}}\left({{l}}_{{\mathrm{\alpha }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{-}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}}^{\phantom{{}}{\mathrm{\alpha }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right){}\left({{\partial }}_{{\mathrm{\alpha }}}\left({{l}}_{{\mathrm{\nu }}}\right)\right)}{{2}}{+}\frac{{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\nu }}}}^{\phantom{{}}{\mathrm{\nu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}}^{\phantom{{}}{\mathrm{\alpha }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\rho }}}}^{\phantom{{}}{\mathrm{\rho }}}{}\left({{\partial }}_{{\mathrm{\nu }}}\left({{l}}_{{\mathrm{\alpha }}}\right)\right){}\left({{\partial }}_{{\mathrm{\rho }}}\left({{l}}_{{\mathrm{\mu }}}\right)\right)}{{2}}$ (8)
 > $\mathrm{Simplify}\left(\right)$
 ${0}$ (9)
 • A new tryhard option significantly improves zero recognition, which also got improved regardless of the new option.

Set spacetimeindices to be represented with lowercase Latin letters to make the input simpler.

 >
 $\left[{\mathrm{spacetimeindices}}{=}{\mathrm{lowercaselatin}}\right]$ (10)
 >
 ${\mathrm{Defined objects with tensor properties}}$
 $\left\{{{A}}_{{a}}{,}{{B}}_{{a}{,}{b}}{,}{{𝒟}}_{{\mathrm{\mu }}}{,}{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{F}}_{{a}{,}{b}{,}{c}}{,}{{H}}_{{a}{,}{b}{,}{c}{,}{d}}{,}{{J}}_{{a}{,}{b}{,}{c}{,}{d}{,}{e}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{C}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{X}}_{{\mathrm{\mu }}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{\mathrm{η}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{l}}_{{\mathrm{\mu }}}{,}{{\mathrm{\Gamma }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}}{,}{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\delta }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (11)

Consider the following tensorial expressions.

 > $\mathrm{e__1}≔{A}_{c}{H}_{d,h,c,b}{F}_{h,d,f}{H}_{a,e,e,f}-{A}_{h}{H}_{a,c,c,f}{F}_{e,d,f}{H}_{d,e,h,b}$
 $\mathrm{e__1}{≔}{{H}}_{{a}{,}{e}\phantom{{e}}{f}}^{\phantom{{a}}\phantom{{,}}\phantom{{e}}{e}\phantom{{f}}}{}{{H}}_{{d}{,}{h}{,}{c}{,}{b}}{}{{F}}_{\phantom{{}}\phantom{{h}}\phantom{{,}}\phantom{{d}}\phantom{{,}}\phantom{{f}}}^{\phantom{{}}{h}{,}{d}{,}{f}}{}{{A}}_{\phantom{{}}\phantom{{c}}}^{\phantom{{}}{c}}{-}{{H}}_{{a}{,}{c}\phantom{{c}}{f}}^{\phantom{{a}}\phantom{{,}}\phantom{{c}}{c}\phantom{{f}}}{}{{H}}_{{d}{,}{e}{,}{h}{,}{b}}{}{{F}}_{\phantom{{}}\phantom{{e}}\phantom{{,}}\phantom{{d}}\phantom{{,}}\phantom{{f}}}^{\phantom{{}}{e}{,}{d}{,}{f}}{}{{A}}_{\phantom{{}}\phantom{{h}}}^{\phantom{{}}{h}}$ (12)
 > $\mathrm{e__2}≔{J}_{f,a,b,g,c}{F}_{e,f,d}{F}_{e,d,g}-{J}_{e,a,b,f,c}{F}_{d,e,g}{F}_{d,g,f}$
 $\mathrm{e__2}{≔}{-}{{J}}_{{e}{,}{a}{,}{b}{,}{f}{,}{c}}{}{{F}}_{{d}\phantom{{e}}{g}}^{\phantom{{d}}{e}\phantom{{g}}}{}{{F}}_{\phantom{{}}\phantom{{d}}\phantom{{,}}\phantom{{g}}\phantom{{,}}\phantom{{f}}}^{\phantom{{}}{d}{,}{g}{,}{f}}{+}{{J}}_{{f}{,}{a}{,}{b}{,}{g}{,}{c}}{}{{F}}_{{e}{,}{d}\phantom{{g}}}^{\phantom{{e}}\phantom{{,}}\phantom{{d}}{g}}{}{{F}}_{\phantom{{}}\phantom{{e}}\phantom{{,}}\phantom{{f}}\phantom{{,}}\phantom{{d}}}^{\phantom{{}}{e}{,}{f}{,}{d}}$ (13)
 > $\mathrm{e__3}≔{J}_{g,c,k,j,e}{H}_{i,h,a,h}{H}_{e,b,i,g}{F}_{d,k,j}-{H}_{h,b,g,j}{F}_{d,k,e}{J}_{j,c,k,e,h}{H}_{g,i,a,i}$
 $\mathrm{e__3}{≔}{-}{{J}}_{{j}{,}{c}{,}{k}{,}{e}{,}{h}}{}{{H}}_{{g}{,}{i}{,}{a}\phantom{{i}}}^{\phantom{{g}}\phantom{{,}}\phantom{{i}}\phantom{{,}}\phantom{{a}}{i}}{}{{H}}_{\phantom{{}}\phantom{{h}}{b}\phantom{{g}}\phantom{{,}}\phantom{{j}}}^{\phantom{{}}{h}\phantom{{b}}{g}{,}{j}}{}{{F}}_{{d}\phantom{{k}}\phantom{{,}}\phantom{{e}}}^{\phantom{{d}}{k}{,}{e}}{+}{{J}}_{{g}{,}{c}{,}{k}{,}{j}{,}{e}}{}{{H}}_{\phantom{{}}\phantom{{e}}{b}{,}{i}\phantom{{g}}}^{\phantom{{}}{e}\phantom{{b}}\phantom{{,}}\phantom{{i}}{g}}{}{{H}}_{\phantom{{}}\phantom{{i}}{h}{,}{a}\phantom{{h}}}^{\phantom{{}}{i}\phantom{{h}}\phantom{{,}}\phantom{{a}}{h}}{}{{F}}_{{d}\phantom{{k}}\phantom{{,}}\phantom{{j}}}^{\phantom{{d}}{k}{,}{j}}$ (14)

There are two, three and four free indices, and five, four and six repeated indices respectively in each of these three expressions, $\mathrm{e__1},\mathrm{e__2},\mathrm{e__3}$ (see Check).

 >
 $\left[\left[\left[\left\{{c}{,}{d}{,}{e}{,}{f}{,}{h}\right\}{,}\left\{{c}{,}{d}{,}{e}{,}{f}{,}{h}\right\}\right]{,}\left\{{a}{,}{b}\right\}\right]{,}\left[\left[\left\{{d}{,}{e}{,}{f}{,}{g}\right\}{,}\left\{{d}{,}{e}{,}{f}{,}{g}\right\}\right]{,}\left\{{a}{,}{b}{,}{c}\right\}\right]{,}\left[\left[\left\{{e}{,}{g}{,}{h}{,}{i}{,}{j}{,}{k}\right\}{,}\left\{{e}{,}{g}{,}{h}{,}{i}{,}{j}{,}{k}\right\}\right]{,}\left\{{a}{,}{b}{,}{c}{,}{d}\right\}\right]\right]$ (15)
 > $\mathrm{Simplify}\left(\mathrm{e__1},\mathrm{tryhard}\right)$
 ${0}$ (16)
 > $\mathrm{Simplify}\left(\mathrm{e__2},\mathrm{tryhard}\right)$
 ${0}$ (17)
 > $\mathrm{Simplify}\left(\mathrm{e__3}\right)$
 ${0}$ (18)
 • Improve the simplification with respect to side relations (equations) in the presence of quantum vectorial equations.
 > $\mathrm{with}\left(\mathrm{Vectors}\right):$
 >
 $\left[{\mathrm{hermitianoperators}}{=}\left\{{L}{,}{p}{,}\stackrel{{\to }}{{p}}{,}\stackrel{{\to }}{{r}}{,}{x}{,}{y}{,}{z}\right\}{,}{\mathrm{quantumoperators}}{=}\left\{{A}{,}{B}{,}{L}{,}{p}{,}\stackrel{{\to }}{{p}}{,}\stackrel{{\to }}{{r}}{,}{x}{,}{y}{,}{z}\right\}\right]$ (19)
 • The identification of the vectorial character of expressions got improved. For example, this commutator is actually a non-projected vector.
 > $\mathrm{Commutator}\left(\mathrm{p_}·\mathrm{p_},\mathrm{r_}\right)$
 ${\left[{‖\stackrel{{\to }}{{p}}‖}^{{2}}{,}\stackrel{{\to }}{{r}}\right]}_{{-}}$ (20)

 > $\mathrm{Identify}\left(\right)$
 ${5}$ (21)
 • Consider an expression involving vectorial quantum operators.
 >
 ${I}$ (22)
 > $\mathrm{_i}\cdot \left(i{x}{{p}}_{{z}}-i{z}{{p}}_{{x}}\right)+\mathrm{_j}\cdot \left(i{y}{{p}}_{{z}}-i{z}{{p}}_{{y}}\right)$
 $\stackrel{{\wedge }}{{i}}{}\left({i}{}{x}{}{{p}}_{{z}}{-}{i}{}{z}{}{{p}}_{{x}}\right){+}\stackrel{{\wedge }}{{j}}{}\left({i}{}{y}{}{{p}}_{{z}}{-}{i}{}{z}{}{{p}}_{{y}}\right)$ (23)
 • Simplify (23) taking into account the definition of angular momentum.
 > $\left\{{y}{{p}}_{{z}}-{z}{{p}}_{{y}}={{L}}_{{x}},{z}{{p}}_{{x}}-{x}{{p}}_{{z}}={{L}}_{{y}}\right\}$
 $\left\{{y}{}{{p}}_{{z}}{-}{z}{}{{p}}_{{y}}{=}{{L}}_{{x}}{,}{z}{}{{p}}_{{x}}{-}{x}{}{{p}}_{{z}}{=}{{L}}_{{y}}\right\}$ (24)
 > $\mathrm{simplify}\left(,\right)$
 ${-}{i}{}\stackrel{{\wedge }}{{i}}{}{{L}}_{{y}}{+}{i}{}\stackrel{{\wedge }}{{j}}{}{{L}}_{{x}}$ (25)
 • Enhancements in the simplification and integration of spherical harmonics (SphericalY) relevant in quantum mechanics.
 > $\frac{1}{8648640}\frac{\sqrt{7293}{ⅇ}^{6I\mathrm{theta}}{\left(1-{\mathrm{cos}\left(\mathrm{phi}\right)}^{2}\right)}^{3}\left(\frac{2027025}{2}{\mathrm{cos}\left(\mathrm{phi}\right)}^{2}-\frac{135135}{2}\right)}{\sqrt{\mathrm{\pi }}}+\mathrm{SphericalY}\left(8,6,\mathrm{phi},\mathrm{theta}\right)$
 $\frac{\sqrt{{7293}}{}{{ⅇ}}^{{6}{}{i}{}{\mathrm{\theta }}}{}{\left({1}{-}{{\mathrm{cos}}\left({\mathrm{\phi }}\right)}^{{2}}\right)}^{{3}}{}\left(\frac{{2027025}{}{{\mathrm{cos}}\left({\mathrm{\phi }}\right)}^{{2}}}{{2}}{-}\frac{{135135}}{{2}}\right)}{{8648640}{}\sqrt{{\mathrm{\pi }}}}{+}{{Y}}_{{8}}^{{6}}\left({\mathrm{\phi }}{,}{\mathrm{\theta }}\right)$ (26)
 > $\mathrm{simplify}\left(\right)$
 ${0}$ (27)
 > $\mathrm{conjugate}\left(\mathrm{SphericalY}\left(8,6,\mathrm{phi},\mathrm{theta}\right)\right)*\left(1+1/5*\mathrm{sin}\left(6*\mathrm{theta}\right)*\mathrm{sin}\left(5*\mathrm{phi}\right)\right)*\mathrm{sin}\left(\mathrm{phi}\right)$
 $\stackrel{{&conjugate0;}}{{{Y}}_{{8}}^{{6}}\left({\mathrm{\phi }}{,}{\mathrm{\theta }}\right)}{}\left({1}{+}\frac{{\mathrm{sin}}\left({6}{}{\mathrm{\theta }}\right){}{\mathrm{sin}}\left({5}{}{\mathrm{\phi }}\right)}{{5}}\right){}{\mathrm{sin}}\left({\mathrm{\phi }}\right)$ (28)
 >
 ${{\int }}_{{0}}^{{\mathrm{\pi }}}{{\int }}_{{0}}^{{2}{}{\mathrm{\pi }}}\stackrel{{&conjugate0;}}{{{Y}}_{{8}}^{{6}}\left({\mathrm{\phi }}{,}{\mathrm{\theta }}\right)}{}\left({1}{+}\frac{{\mathrm{sin}}\left({6}{}{\mathrm{\theta }}\right){}{\mathrm{sin}}\left({5}{}{\mathrm{\phi }}\right)}{{5}}\right){}{\mathrm{sin}}\left({\mathrm{\phi }}\right)\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{\theta }}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}{ⅆ}{\mathrm{\phi }}$ (29)
 > $\mathrm{value}\left(\right)$
 ${-}\frac{{37}{}{i}}{{81920}}{}{{\mathrm{\pi }}}^{{3}}{{2}}}{}\sqrt{{7293}}$ (30)

Tensors

A number of relevant changes happened in the tensor routines of the Physics package, towards making the routines pack more functionality, the simplification more powerful, and the handling of symmetries, substitutions, and other operations more flexible and natural.

 • Physics now works with four kinds of Minkowski spaces (different signatures) to accommodate the typical situations seen in textbooks; to these, correspond the signatures +---, ---+, -+++ and +++-.
 • Allow setting the metric by specifying the signature directly, as in g_[-] or g_[+---], or $\mathrm{Setup}\left(\mathrm{metric}=\mathrm{---+}\right)$ or $\mathrm{Setup}\left({g}_{\mathrm{\mu },\mathrm{\nu }}=\mathrm{---+}\right)$.
 • The signature keyword of the Physics Setup is now in use, to set the metric and to indicate the form of the orthonormal tetrad, in turn used to derive the form of a null tetrad.
 • Automatic detection of the position of t as the time variable when you set the coordinates automatically sets the signature of the default Minkowski spacetime metric accordingly to ---+ or +---.
 • New keywords with special meaning when indexing the Physics (also the user defined) tensors: · ~; for example g_[~] returns the all-contravariant matrix form of the metric. · definition; for example Ricci[definition] returns the definition of the Ricci tensor; works also with user-defined tensors. · scalars; for example Weyl[scalars] and Ricci[scalars] return the five Weyl and seven Ricci scalars used to perform a Petrov classification and in the Newman-Penrose formalism. · scalarsdefinition, and invariantsdefinition; for example Weyl[scalarsdefinition] or Riemann[invariantsdefinition] return the corresponding definitions for the scalars and invariants. · nullvectors; for example, when the new Tetrads subpackage is loaded, e_[nullvectors] returns a sequence of null vectors with their products normalized according to the Newman-Penrose formalism. · matrix; this keyword was introduced in previous releases, and in Maple 2015 it can appear after a space index (not spacetime), in which case a matrix with only the space components is returned.
 • Tensorial expressions can now have spacetime indices (related to a global system of references) and tetrad indices (related to a local system of references) at the same time, or they be rewritten in one (spacetime) or the other (tetrad) frames.
 • The matrix keyword can be used with spacetime, space, or tetrad indices, resulting in the corresponding matrix
 • Implement automatic determination of symmetry under permutation of tensor indices when the tensor is defined as a matrix.
 • New conversions from the Weyl to the Ricci tensors, and from Weyl to the Christoffel symbols.
 • New option evaluatetrace = true or false within convert/Ricci, to avoid automatically evaluating the Ricci trace when performing conversions that involve this trace.
 • New option 'evaluate' to convert/g_, convert/Christoffel and convert/Ricci. With this option set to false, it is possible to see the algebraic form of the result (that is, of the tensors involved) before evaluating it.
 • The Maple 18 Library:-SubstituteTensor command, got enhanced and transformed into one of the main Physics commands, that substitutes tensorial equation(s) Eqs into an expression, taking care of the free and repeated indices, such that: 1) equations in Eqs are interpreted as mappings having the free indices as parameters, 2) repeated indices in Eqs do not clash with repeated indices in the expression and 3) spacetime, space, and tetrad indices are handled independently, so they can all be present in Eqs and in the expression at the same time. This new command can also substitute algebraic sub-expressions of type product or sum within the expression, generalizing and unifying the functionality of the subs and algsubs commands for algebraic tensor expressions.

Examples

 • Minkowski spacetime: four different conventions.
 >
 ${\mathrm{Default differentiation variables for d_, D_ and dAlembertian are:}}\left\{{X}{=}\left({\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right)\right\}$
 ${\mathrm{Systems of spacetime Coordinates are:}}\left\{{X}{=}\left({\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right)\right\}$
 $\left\{{X}\right\}$ (31)

By default, the metric has a signature - - - +, with time in the last position.

 > $\mathrm{g_}\left[\right]$
 > $\mathrm{x0}$
 ${\mathrm{x4}}$ (32)

You can now set the signature of a Minkowski spacetime directly from the metric, or using the Setup keyword signature to use any of the four conventions frequently found in textbooks (- - - +), (+ - - -), (+++ -), and (- +++), and to indicate an Euclidean spacetime you can indicate (++++) or just +.

 > $\mathrm{g_}\left[\mathrm{+}\right]$
 ${\mathrm{Changing the signature of the tensor spacetime to: + + + +}}$
 ${\mathrm{The Euclidean metric in coordinates}}\left[{\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right]$
 > $\mathrm{g_}\left[\mathrm{+++-}\right]$
 ${\mathrm{The Minkowski metric, with signature + + + -, in coordinates}}\left[{\mathrm{x1}}{,}{\mathrm{x2}}{,}{\mathrm{x3}}{,}{\mathrm{x4}}\right]$

Using Setup's signature keyword:

 >
 $\left[{\mathrm{signature}}{=}{\mathrm{- + + +}}\right]$ (33)
 > $\mathrm{g_}\left[\right]$
 • The position of the different sign in a Minkowski spacetime always, either the first or the last symbol in the signature, refers to the time variable. So if you indicate $t$ as one of the coordinates, say in the first position, and the signature is not in agreement with time in the first position, it is now automatically corrected and the metric set. For example, in this moment, due to (33), time is implicitly assumed to be in the first position; set coordinates with $t$ in the last (fourth) position.
 > $\mathrm{Coordinates}\left(X=\left[x,y,z,t\right]\right)$
 ${\mathrm{Detected t, the time variable, in position 4. Changing the signature of the spacetime metric accordingly, to: + + + -}}$
 ${\mathrm{Default differentiation variables for d_, D_ and dAlembertian are:}}\left\{{X}{=}\left({x}{,}{y}{,}{z}{,}{t}\right)\right\}$
 ${\mathrm{Systems of spacetime Coordinates are:}}\left\{{X}{=}\left({x}{,}{y}{,}{z}{,}{t}\right)\right\}$
 $\left\{{X}\right\}$ (34)
 > $\mathrm{g_}\left[\right]$
 • New shortcut for the matrix form of the all contravariant components of a tensor.

Since previous releases, you can always request the matrix (or array in the case of more than 2 indices) form of a tensor by adding the keyword matrix as a last index. For example, ${\mathrm{g_}}_{\mathrm{\mu },\mathrm{\nu },\mathrm{matrix}}$ . A shortcut notation for the all covariant components is to omit the indices, as in . A new shortcut notation for the all contravariant components is to pass only they tilde ~. To see the difference, consider a non-Minkowski spacetime, for instance set it in one go via:

 > $\mathrm{g_}\left[\mathrm{sc}\right]$
 ${\mathrm{Systems of spacetime Coordinates are:}}\left\{{X}{=}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}{,}{t}\right)\right\}$
 ${\mathrm{Default differentiation variables for d_, D_ and dAlembertian are:}}\left\{{X}{=}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}{,}{t}\right)\right\}$
 ${\mathrm{The Schwarzschild metric in coordinates}}\left[{r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}{,}{t}\right]$
 ${\mathrm{Parameters:}}\left[{m}\right]$

The all contravariant components of the metric are:

 > $\mathrm{g_}\left[\mathrm{~}\right]$
 • It is now possible to also use the keyword matrix passing only space (not spacetime) indices, resulting in a matrix with only the space components.
 > $\mathrm{Setup}\left(\mathrm{spaceindices}=\mathrm{lowercase_is}\right)$
 $\left[{\mathrm{spaceindices}}{=}{\mathrm{lowercaselatin_is}}\right]$ (35)
 >
 ${{g}}_{\phantom{{}}\phantom{{i}}\phantom{{,}}\phantom{{j}}}^{\phantom{{}}{i}{,}{j}}{=}\left[\begin{array}{ccc}\frac{{-}{r}{+}{2}{}{m}}{{r}}& {0}& {0}\\ {0}& {-}\frac{{1}}{{{r}}^{{2}}}& {0}\\ {0}& {0}& {-}\frac{{1}}{{{r}}^{{2}}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}}\end{array}\right]$ (36)

 • All the tensor predefined in the Physics package, or defined using defining equations and the Define command, accept a new keyword, definition, that returns their definition. For example:
 > $\mathrm{Ricci}\left[\mathrm{definition}\right]$
 ${{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{\partial }}_{{\mathrm{\alpha }}}\left({{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}{\mathrm{\mu }}{,}{\mathrm{\nu }}}^{\phantom{{}}{\mathrm{\alpha }}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}}\right){-}\left({{\partial }}_{{\mathrm{\nu }}}\left({{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}{\mathrm{\mu }}{,}{\mathrm{\alpha }}}^{\phantom{{}}{\mathrm{\alpha }}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}}\right)\right){+}{{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\beta }}}{\mathrm{\mu }}{,}{\mathrm{\nu }}}^{\phantom{{}}{\mathrm{\beta }}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}}{}{{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}{\mathrm{\beta }}{,}{\mathrm{\alpha }}}^{\phantom{{}}{\mathrm{\alpha }}\phantom{{\mathrm{\beta }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}}{-}{{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\beta }}}{\mathrm{\mu }}{,}{\mathrm{\alpha }}}^{\phantom{{}}{\mathrm{\beta }}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}}{}{{\mathrm{\Gamma }}}_{\phantom{{}}\phantom{{\mathrm{\alpha }}}{\mathrm{\nu }}{,}{\mathrm{\beta }}}^{\phantom{{}}{\mathrm{\alpha }}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}$ (37)
 > $\mathrm{Einstein}\left[\mathrm{definition}\right]$
 ${{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{-}\frac{{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{}{{R}}_{{\mathrm{\alpha }}\phantom{{\mathrm{\alpha }}}}^{\phantom{{\mathrm{\alpha }}}{\mathrm{\alpha }}}}{{2}}$ (38)

Define a tensor using a tensorial equation:

 > $\mathrm{PDEtools}:-\mathrm{declare}\left(A\left(X\right)\right)$
 ${A}\left({r}{,}{\mathrm{\theta }}{,}{\mathrm{\phi }}{,}{t}\right){}{\mathrm{will now be displayed as}}{}{A}$ (39)
 >
 ${{F}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{\partial }}_{{\mathrm{\mu }}}\left({{A}}_{{\mathrm{\nu }}}\right){-}\left({{\partial }}_{{\mathrm{\nu }}}\left({{A}}_{{\mathrm{\mu }}}\right)\right)$ (40)

Note that for consistency of free indices, if you now define the left-hand side ${F}_{\mathrm{\mu },\mathrm{\nu }}$ as a tensor, in the right-hand side ${A}_{\mathrm{\mu }}$ is also a tensor and so is automatically defined as well.

 > $\mathrm{Define}\left(\right)$
 ${\mathrm{Defined objects with tensor properties}}$
 $\left\{{{A}}_{{\mathrm{\mu }}}{,}{{𝒟}}_{{\mathrm{\mu }}}{,}{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{F}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{C}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{X}}_{{\mathrm{\mu }}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\Gamma }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}}{,}{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\delta }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (41)

Check the definition of ${F}_{\mathrm{\mu },\mathrm{\nu }}$:

 > $F\left[\mathrm{definition}\right]$
 ${{F}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{\partial }}_{{\mathrm{\mu }}}\left({{A}}_{{\mathrm{\nu }}}\right){-}\left({{\partial }}_{{\mathrm{\nu }}}\left({{A}}_{{\mathrm{\mu }}}\right)\right)$ (42)
 • A new keyword scalars for the Weyl and Ricci tensors, generates the Weyl and Ricci scalars of the Newman-Penrose formalism. Correspondingly, another new keyword scalarsdefinition return the definition of these scalars.
 > $\mathrm{Weyl}\left[\mathrm{scalars}\right]$
 $\mathrm{ψ__0}{=}{0}{,}\mathrm{ψ__1}{=}{0}{,}\mathrm{ψ__2}{=}\frac{{m}}{{{r}}^{{3}}}{,}\mathrm{ψ__3}{=}{0}{,}\mathrm{ψ__4}{=}{0}$ (43)
 > $\mathrm{Weyl}\left[\mathrm{scalarsdefinition}\right]$
 $\mathrm{ψ__0}{=}{{C}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}{{l}}_{{\mathrm{\mu }}}{}{{m}}_{{\mathrm{\nu }}}{}{{l}}_{{\mathrm{\alpha }}}{}{{m}}_{{\mathrm{\beta }}}{,}\mathrm{ψ__1}{=}{{C}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{{\mathrm{\nu }}}{}{{l}}_{{\mathrm{\alpha }}}{}{{m}}_{{\mathrm{\beta }}}{,}\mathrm{ψ__2}{=}{{C}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}{{l}}_{{\mathrm{\mu }}}{}{{m}}_{{\mathrm{\nu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\alpha }}}{}{{n}}_{{\mathrm{\beta }}}{,}\mathrm{ψ__3}{=}{{C}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{{\mathrm{\nu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\alpha }}}{}{{n}}_{{\mathrm{\beta }}}{,}\mathrm{ψ__4}{=}{{C}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}\phantom{{,}}\phantom{{\mathrm{\nu }}}\phantom{{,}}\phantom{{\mathrm{\alpha }}}\phantom{{,}}\phantom{{\mathrm{\beta }}}}^{\phantom{{}}{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{}{{n}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\nu }}}{}{{n}}_{{\mathrm{\alpha }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\beta }}}$ (44)

In the definition earlier, the tensors ${l}_{\mathrm{\mu }},{n}_{\mathrm{\mu }},{m}_{\mathrm{\mu }},{\stackrel{&conjugate0;}{m}}_{\mathrm{\mu }}$ are the null tensors of the Newman-Penrose formalism, implemented in Maple 2015 within the new Tetrads package.

 > $\mathrm{with}\left(\mathrm{Tetrads}\right)$
 ${\mathrm{Setting lowercaselatin_ah letters to represent tetrad indices}}$
 ${\mathrm{Defined as tetrad tensors}}\left({\mathrm{see ?Physics,tetrads}}\right){,}{{𝔢}}_{{a}{,}{\mathrm{\mu }}}{,}{{\mathrm{\eta }}}_{{a}{,}{b}}{,}{{\mathrm{\gamma }}}_{{a}{,}{b}{,}{c}}{,}{{\mathrm{\lambda }}}_{{a}{,}{b}{,}{c}}$
 ${\mathrm{Defined as spacetime tensors representing the NP null vectors of the tetrad formalism}}\left({\mathrm{see ?Physics,tetrads}}\right){,}{{l}}_{{\mathrm{\mu }}}{,}{{n}}_{{\mathrm{\mu }}}{,}{{m}}_{{\mathrm{\mu }}}{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}$
 $\left[{\mathrm{IsTetrad}}{,}{\mathrm{NullTetrad}}{,}{\mathrm{OrthonormalTetrad}}{,}{\mathrm{SimplifyTetrad}}{,}{\mathrm{TransformTetrad}}{,}{\mathrm{e_}}{,}{\mathrm{eta_}}{,}{\mathrm{gamma_}}{,}{\mathrm{l_}}{,}{\mathrm{lambda_}}{,}{\mathrm{m_}}{,}{\mathrm{mb_}}{,}{\mathrm{n_}}\right]$ (45)

In this new Tetrads package, ${𝔢}_{a,\mathrm{\mu }}$ is the tetrad (vierbein), by default an orthonormal tetrad that can be set to be a null tetrad using $\mathrm{Setup}\left(\mathrm{tetrad}=\mathrm{null}\right)$, ${\mathrm{\gamma }}_{a,b,c},{\mathrm{\lambda }}_{a,b,c}$ are respectively the Ricci rotation coefficients and the lambda tensor defined in the "Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2" (definitions (98.9) and (98.10))

For instance, the form of the orthonormal tetrad for the $\mathrm{Schwarzschild}$ metric set in  is:

 > $\mathrm{e_}\left[\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}\frac{\sqrt{{r}}}{\sqrt{{-}{r}{+}{2}{}{m}}}& {0}& {0}& {0}\\ {0}& {-I}{}{r}& {0}& {0}\\ {0}& {0}& {-I}{}{r}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)& {0}\\ {0}& {0}& {0}& \frac{\sqrt{{-}{r}{+}{2}{}{m}}}{\sqrt{{r}}}\end{array}\right]$ (46)

A new keyword for the tetrad ${𝔢}_{a,\mathrm{\mu }}$ is nullvectors:

 > $\mathrm{e_}\left[\mathrm{nullvectors}\right]$
 ${{l}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}\frac{\sqrt{{2}}{}\sqrt{{r}}}{{2}{}\sqrt{{-}{r}{+}{2}{}{m}}}& {0}& {0}& \frac{\sqrt{{2}}{}\sqrt{{-}{r}{+}{2}{}{m}}}{{2}{}\sqrt{{r}}}\end{array}\right]{,}{{n}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{-}\frac{\sqrt{{2}}{}\sqrt{{r}}}{{2}{}\sqrt{{-}{r}{+}{2}{}{m}}}& {0}& {0}& \frac{\sqrt{{2}}{}\sqrt{{-}{r}{+}{2}{}{m}}}{{2}{}\sqrt{{r}}}\end{array}\right]{,}{{m}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{0}& {-}\frac{{I}}{{2}}{}\sqrt{{2}}{}{r}& \frac{\sqrt{{2}}{}{r}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}{{2}}& {0}\end{array}\right]{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{0}& {-}\frac{{I}}{{2}}{}\sqrt{{2}}{}{r}& {-}\frac{\sqrt{{2}}{}{r}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}{{2}}& {0}\end{array}\right]$ (47)

The matrix keyword can be used with different kinds of indices, representing different objects. For example, the spacetime components of the electromagnetic tensor defined in equation (40) are:

 > $F\left[\right]$
 ${{F}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}\left[\begin{array}{cccc}{0}& {\left({{A}}_{{2}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\theta }}}& {\left({{A}}_{{3}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\phi }}}& {\left({{A}}_{{4}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{t}}\\ {\left({{A}}_{{1}}\right)}_{{\mathrm{\theta }}}{-}{\left({{A}}_{{2}}\right)}_{{r}}& {0}& {\left({{A}}_{{3}}\right)}_{{\mathrm{\theta }}}{-}{\left({{A}}_{{2}}\right)}_{{\mathrm{\phi }}}& {\left({{A}}_{{4}}\right)}_{{\mathrm{\theta }}}{-}{\left({{A}}_{{2}}\right)}_{{t}}\\ {\left({{A}}_{{1}}\right)}_{{\mathrm{\phi }}}{-}{\left({{A}}_{{3}}\right)}_{{r}}& {\left({{A}}_{{2}}\right)}_{{\mathrm{\phi }}}{-}{\left({{A}}_{{3}}\right)}_{{\mathrm{\theta }}}& {0}& {\left({{A}}_{{4}}\right)}_{{\mathrm{\phi }}}{-}{\left({{A}}_{{3}}\right)}_{{t}}\\ {\left({{A}}_{{1}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{r}}& {\left({{A}}_{{2}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{\mathrm{\theta }}}& {\left({{A}}_{{3}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{\mathrm{\phi }}}& {0}\end{array}\right]$ (48)

The components of this tensor in the local inertial (tetrad) system of references are:

 > $F\left[a,b,\mathrm{matrix}\right]$
 ${{F}}_{{a}{,}{b}}{=}\left[\begin{array}{cccc}{0}& \frac{{I}{}\sqrt{{-}{r}{+}{2}{}{m}}{}\left({\left({{A}}_{{2}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\theta }}}\right)}{{{r}}^{{3}}{{2}}}}& \frac{{I}{}\sqrt{{-}{r}{+}{2}{}{m}}{}\left({\left({{A}}_{{3}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\phi }}}\right)}{{{r}}^{{3}}{{2}}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}& {\left({{A}}_{{1}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{r}}\\ \frac{{-I}{}\sqrt{{-}{r}{+}{2}{}{m}}{}\left({\left({{A}}_{{2}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\theta }}}\right)}{{{r}}^{{3}}{{2}}}}& {0}& \frac{{\left({{A}}_{{2}}\right)}_{{\mathrm{\phi }}}{-}{\left({{A}}_{{3}}\right)}_{{\mathrm{\theta }}}}{{{r}}^{{2}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}& \frac{{I}{}\left({\left({{A}}_{{2}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{\mathrm{\theta }}}\right)}{\sqrt{{-}{r}{+}{2}{}{m}}{}\sqrt{{r}}}\\ \frac{{-I}{}\sqrt{{-}{r}{+}{2}{}{m}}{}\left({\left({{A}}_{{3}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{\mathrm{\phi }}}\right)}{{{r}}^{{3}}{{2}}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}& \frac{{\left({{A}}_{{3}}\right)}_{{\mathrm{\theta }}}{-}{\left({{A}}_{{2}}\right)}_{{\mathrm{\phi }}}}{{{r}}^{{2}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}& {0}& \frac{{I}{}\left({\left({{A}}_{{3}}\right)}_{{t}}{-}{\left({{A}}_{{4}}\right)}_{{\mathrm{\phi }}}\right)}{\sqrt{{-}{r}{+}{2}{}{m}}{}\sqrt{{r}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right)}\\ {\left({{A}}_{{4}}\right)}_{{r}}{-}{\left({{A}}_{{1}}\right)}_{{t}}& \frac{{I}{}\left({\left({{A}}_{{4}}\right)}_{{\mathrm{\theta }}}{-}{\left({{A}}_{{2}}\right)}_{{t}}\right)}{\sqrt{{r}}{}\sqrt{{-}{r}{+}{2}{}{m}}}& \frac{{I}{}\left({\left({{A}}_{{4}}\right)}_{{\mathrm{\phi }}}{-}{\left({{A}}_{{3}}\right)}_{{t}}\right)}{\sqrt{{r}}{}{\mathrm{sin}}\left({\mathrm{\theta }}\right){}\sqrt{{-}{r}{+}{2}{}{m}}}& {0}\end{array}\right]$ (49)

Likewise, the nonzero components of the Riemann tensor in the global (spacetime, Greek indices) and local (tetrad, lowercase roman indices from a to h) system of references are:

 > $\mathrm{Riemann}\left[\mathrm{nonzero}\right]$
 ${{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{=}\left\{\left({1}{,}{2}{,}{1}{,}{2}\right){=}{-}\frac{{m}}{{-}{r}{+}{2}{}{m}}{,}\left({1}{,}{2}{,}{2}{,}{1}\right){=}\frac{{m}}{{-}{r}{+}{2}{}{m}}{,}\left({1}{,}{3}{,}{1}{,}{3}\right){=}{-}\frac{{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}}{{-}{r}{+}{2}{}{m}}{,}\left({1}{,}{3}{,}{3}{,}{1}\right){=}\frac{{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}}{{-}{r}{+}{2}{}{m}}{,}\left({1}{,}{4}{,}{1}{,}{4}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({1}{,}{4}{,}{4}{,}{1}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({2}{,}{1}{,}{1}{,}{2}\right){=}\frac{{m}}{{-}{r}{+}{2}{}{m}}{,}\left({2}{,}{1}{,}{2}{,}{1}\right){=}{-}\frac{{m}}{{-}{r}{+}{2}{}{m}}{,}\left({2}{,}{3}{,}{2}{,}{3}\right){=}{-}{2}{}{r}{}{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{,}\left({2}{,}{3}{,}{3}{,}{2}\right){=}{2}{}{r}{}{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{,}\left({2}{,}{4}{,}{2}{,}{4}\right){=}\frac{\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}{,}\left({2}{,}{4}{,}{4}{,}{2}\right){=}\frac{{-}{2}{}{{m}}^{{2}}{+}{m}{}{r}}{{{r}}^{{2}}}{,}\left({3}{,}{1}{,}{1}{,}{3}\right){=}\frac{{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}}{{-}{r}{+}{2}{}{m}}{,}\left({3}{,}{1}{,}{3}{,}{1}\right){=}{-}\frac{{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}}{{-}{r}{+}{2}{}{m}}{,}\left({3}{,}{2}{,}{2}{,}{3}\right){=}{2}{}{r}{}{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{,}\left({3}{,}{2}{,}{3}{,}{2}\right){=}{-}{2}{}{r}{}{m}{}{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{,}\left({3}{,}{4}{,}{3}{,}{4}\right){=}\frac{{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{}\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}{,}\left({3}{,}{4}{,}{4}{,}{3}\right){=}{-}\frac{{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{}\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}{,}\left({4}{,}{1}{,}{1}{,}{4}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({4}{,}{1}{,}{4}{,}{1}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({4}{,}{2}{,}{2}{,}{4}\right){=}\frac{{-}{2}{}{{m}}^{{2}}{+}{m}{}{r}}{{{r}}^{{2}}}{,}\left({4}{,}{2}{,}{4}{,}{2}\right){=}\frac{\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}{,}\left({4}{,}{3}{,}{3}{,}{4}\right){=}{-}\frac{{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{}\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}{,}\left({4}{,}{3}{,}{4}{,}{3}\right){=}\frac{{{\mathrm{sin}}\left({\mathrm{\theta }}\right)}^{{2}}{}\left({-}{r}{+}{2}{}{m}\right){}{m}}{{{r}}^{{2}}}\right\}$ (50)
 > $\mathrm{Riemann}\left[a,b,c,d,\mathrm{nonzero}\right]$
 ${{R}}_{{a}{,}{b}{,}{c}{,}{d}}{=}\left\{\left({1}{,}{2}{,}{1}{,}{2}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({1}{,}{2}{,}{2}{,}{1}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({1}{,}{3}{,}{1}{,}{3}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({1}{,}{3}{,}{3}{,}{1}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({1}{,}{4}{,}{1}{,}{4}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({1}{,}{4}{,}{4}{,}{1}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({2}{,}{1}{,}{1}{,}{2}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({2}{,}{1}{,}{2}{,}{1}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({2}{,}{3}{,}{2}{,}{3}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({2}{,}{3}{,}{3}{,}{2}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({2}{,}{4}{,}{2}{,}{4}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({2}{,}{4}{,}{4}{,}{2}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({3}{,}{1}{,}{1}{,}{3}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({3}{,}{1}{,}{3}{,}{1}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({3}{,}{2}{,}{2}{,}{3}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({3}{,}{2}{,}{3}{,}{2}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({3}{,}{4}{,}{3}{,}{4}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({3}{,}{4}{,}{4}{,}{3}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({4}{,}{1}{,}{1}{,}{4}\right){=}{-}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({4}{,}{1}{,}{4}{,}{1}\right){=}\frac{{2}{}{m}}{{{r}}^{{3}}}{,}\left({4}{,}{2}{,}{2}{,}{4}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({4}{,}{2}{,}{4}{,}{2}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}{,}\left({4}{,}{3}{,}{3}{,}{4}\right){=}\frac{{m}}{{{r}}^{{3}}}{,}\left({4}{,}{3}{,}{4}{,}{3}\right){=}{-}\frac{{m}}{{{r}}^{{3}}}\right\}$ (51)
 • Implement automatic determination of symmetry under permutation of tensor indices when the tensor is defined as a matrix.



Define a tensor using a symmetric matrix on the right-hand-side:

 >
 >
 ${\mathrm{Defined objects with tensor properties}}$
 $\left\{{{A}}_{{\mathrm{\mu }}}{,}{{𝒟}}_{{\mathrm{\mu }}}{,}{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{F}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{M}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{R}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{C}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}{,}{\mathrm{\beta }}}{,}{{X}}_{{\mathrm{\mu }}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{𝔢}}_{{a}{,}{\mathrm{\mu }}}{,}{{\mathrm{\eta }}}_{{a}{,}{b}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\gamma }}}_{{a}{,}{b}{,}{c}}{,}{{l}}_{{\mathrm{\mu }}}{,}{{\mathrm{\lambda }}}_{{a}{,}{b}{,}{c}}{,}{{m}}_{{\mathrm{\mu }}}{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}{,}{{n}}_{{\mathrm{\mu }}}{,}{{\mathrm{\Gamma }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}{,}{\mathrm{\alpha }}}{,}{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\delta }}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (52)

Check whether the system detected that ${\mathrm{\eta }}_{\mathrm{\alpha },\mathrm{\beta }}$ is symmetric under permutation of $\mathrm{\alpha },\mathrm{\beta }$.

 > $\mathrm{Library}:-\mathrm{IsTensorialSymmetric}\left({M}_{\mathrm{\alpha },\mathrm{\beta }}\right)$
 ${\mathrm{true}}$ (53)

Hence, the indices of this tensor are automatically normalized taking this symmetry into account, so that:

 > ${M}_{\mathrm{\alpha },\mathrm{\beta }}$
 ${{M}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}}$ (54)
 > ${M}_{\mathrm{\beta },\mathrm{\alpha }}$
 ${{M}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}}$ (55)
 > $-$
 ${0}$ (56)

The formalism of tetrads in general relativity got implemented within Physics as a new package, Physics:-Tetrads, with 13 commands, mainly the null vectors of the Newman-Penrose formalism, the tetrad tensors ${𝔢}_{a,\mathrm{\mu }},{\mathrm{\eta }}_{a,b},{\mathrm{\gamma }}_{a,b,c},{\mathrm{\lambda }}_{a,b,c}$, respectively: the tetrad, the tetrad metric, the Ricci rotation coefficients, and the lambda tensor, plus five algebraic manipulation commands: IsTetrad, NullTetrad, OrthonormalTetrad, SimplifyTetrad, and TransformTetrad to construct orthonormal and null tetrads of different forms and using different methods.

Examples

The new Tetrads package contains 13 commands for computing in a local (tetrad) frame. Tensor components in the local frame are represented with tetrad indices using a type of letter different than the one representing global spacetime indices. You can set the type of letter using Setup, or just load the new package and the type of letter will be set automatically.

 > $\mathrm{restart};\mathrm{with}\left(\mathrm{Physics}\right):\mathrm{with}\left(\mathrm{Tetrads}\right);$
 ${\mathrm{Setting lowercaselatin letters to represent tetrad indices}}$
 ${\mathrm{Defined as tetrad tensors}}\left({\mathrm{see ?Physics,tetrads}}\right){,}{{𝔢}}_{{a}{,}{\mathrm{\mu }}}{,}{{\mathrm{\eta }}}_{{a}{,}{b}}{,}{{\mathrm{\gamma }}}_{{a}{,}{b}{,}{c}}{,}{{\mathrm{\lambda }}}_{{a}{,}{b}{,}{c}}$
 ${\mathrm{Defined as spacetime tensors representing the NP null vectors of the tetrad formalism}}\left({\mathrm{see ?Physics,tetrads}}\right){,}{{l}}_{{\mathrm{\mu }}}{,}{{n}}_{{\mathrm{\mu }}}{,}{{m}}_{{\mathrm{\mu }}}{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}$
 $\left[{\mathrm{IsTetrad}}{,}{\mathrm{NullTetrad}}{,}{\mathrm{OrthonormalTetrad}}{,}{\mathrm{SimplifyTetrad}}{,}{\mathrm{TransformTetrad}}{,}{\mathrm{e_}}{,}{\mathrm{eta_}}{,}{\mathrm{gamma_}}{,}{\mathrm{l_}}{,}{\mathrm{lambda_}}{,}{\mathrm{m_}}{,}{\mathrm{mb_}}{,}{\mathrm{n_}}\right]$ (57)

The most relevant commands are ${𝔢}_{a,\mathrm{\mu }}$ and ${\mathrm{\eta }}_{a,b}$, respectively representing the tetrad (also vierbein; by default, this is an orthonormal tetrad) and the tetrad metric (that is, the metric of the local - by default inertial - frame)

 > $\mathrm{e_}\left[\mathrm{definition}\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{}{{𝔢}}_{{b}\phantom{{\mathrm{\mu }}}}^{\phantom{{b}}{\mathrm{\mu }}}{=}{{\mathrm{\eta }}}_{{a}{,}{b}}$ (58)

Then ${\mathrm{\gamma }}_{a,b,c}$ are the Ricci rotation coefficients and ${\mathrm{\lambda }}_{a,b,c}$ is a linear combination of them, according to the definitions in the "Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2" (definitions (98.9) and (98.10)).

 > $\mathrm{gamma_}\left[\mathrm{definition}\right]$
 ${{\mathrm{\gamma }}}_{{a}{,}{b}{,}{c}}{=}\left({{𝒟}}_{{\mathrm{\nu }}}\left({{𝔢}}_{{a}{,}{\mathrm{\mu }}}\right)\right){}{{𝔢}}_{{b}\phantom{{\mathrm{\mu }}}}^{\phantom{{b}}{\mathrm{\mu }}}{}{{𝔢}}_{{c}\phantom{{\mathrm{\nu }}}}^{\phantom{{c}}{\mathrm{\nu }}}$ (59)
 > $\mathrm{lambda_}\left[\mathrm{definition}\right]$
 ${{\mathrm{\lambda }}}_{{a}{,}{b}{,}{c}}{=}\left({{𝒟}}_{{\mathrm{\nu }}}\left({{𝔢}}_{{a}{,}{\mathrm{\mu }}}\right){-}\left({{𝒟}}_{{\mathrm{\mu }}}\left({{𝔢}}_{{a}{,}{\mathrm{\nu }}}\right)\right)\right){}{{𝔢}}_{{b}\phantom{{\mathrm{\mu }}}}^{\phantom{{b}}{\mathrm{\mu }}}{}{{𝔢}}_{{c}\phantom{{\mathrm{\nu }}}}^{\phantom{{c}}{\mathrm{\nu }}}$ (60)

The ${l}_{\mathrm{\mu }},{n}_{\mathrm{\mu }},{m}_{\mathrm{\mu }},{\stackrel{&conjugate0;}{m}}_{\mathrm{\mu }}$ tensors are the null tensors of the Newman-Penrose formalism and the commands  $\mathrm{IsTetrad},\mathrm{NullTetrad},\mathrm{OrthonormalTetrad},\mathrm{SimplifyTetrad},\mathrm{TransformTetrad}$ are for manipulating and exploring different forms of tetrads.

In a flat space, the spacetime and tetrad metrics are the same, so the orthonormal tetrad is just the identity.

 > $\mathrm{g_}\left[\right]$
 > $\mathrm{eta_}\left[\right]$
 > $\mathrm{e_}\left[\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{1}& {0}& {0}& {0}\\ {0}& {1}& {0}& {0}\\ {0}& {0}& {1}& {0}\\ {0}& {0}& {0}& {1}\end{array}\right]$ (61)

In a curved spacetime, for instance, set a $"Local Rotational Symmetry metric"$ metric:

 > $\mathrm{g_}\left[\left[13,7,5\right]\right]$
 ${\mathrm{Systems of spacetime Coordinates are:}}\left\{{X}{=}\left({x}{,}{y}{,}{z}{,}{t}\right)\right\}$
 ${\mathrm{Default differentiation variables for d_, D_ and dAlembertian are:}}\left\{{X}{=}\left({x}{,}{y}{,}{z}{,}{t}\right)\right\}$
 ${\mathrm{The metric in coordinates}}\left[{x}{,}{y}{,}{z}{,}{t}\right]$
 ${\mathrm{Parameters:}}\left[{\mathrm{\epsilon }}{,}{A}\left({t}\right){,}{B}\left({t}\right)\right]$
 ${\mathrm{Comments:}}{_ⅇpsilon=1 or _ⅇpsilon=-1}$
 >
 ${A}\left({t}\right){}{\mathrm{will now be displayed as}}{}{A}$
 ${B}\left({t}\right){}{\mathrm{will now be displayed as}}{}{B}$ (62)

 > $\mathrm{e_}\left[\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{-I}{}\sqrt{{\mathrm{\epsilon }}}{}{A}& {0}& {0}& {0}\\ {0}& {-I}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}{}{B}{}{{ⅇ}}^{{A}{}{x}}& \frac{{I}{}{\mathrm{sinh}}\left({2}{}{x}\right){}{B}{}{{ⅇ}}^{{A}{}{x}}}{\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\\ {0}& {0}& \frac{{-I}{}{B}{}{{ⅇ}}^{{A}{}{x}}}{\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\\ {0}& {0}& {0}& {I}{}\sqrt{{\mathrm{\epsilon }}}\end{array}\right]$ (63)

The following null vectors correspond to this tetrad:

 > $\mathrm{e_}\left[\mathrm{nullvectors}\right]$
 ${{l}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{-}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}{}{A}& {0}& {0}& \frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}\end{array}\right]{,}{{n}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}{}{A}& {0}& {0}& \frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}\end{array}\right]{,}{{m}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{0}& {-}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}{}{B}{}{{ⅇ}}^{{A}{}{x}}& \frac{\sqrt{{2}}{}{B}{}{{ⅇ}}^{{A}{}{x}}{}\left({I}{}{\mathrm{sinh}}\left({2}{}{x}\right){+}{1}\right)}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\end{array}\right]{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{0}& {-}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}{}{B}{}{{ⅇ}}^{{A}{}{x}}& \frac{\sqrt{{2}}{}{B}{}{{ⅇ}}^{{A}{}{x}}{}\left({I}{}{\mathrm{sinh}}\left({2}{}{x}\right){-}{1}\right)}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\end{array}\right]$ (64)

You can compute these null vectors directly since these are also part of the Tetrads package:

 >
 ${{l}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{,}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{,}{{l}}_{{\mathrm{\mu }}}{}{{m}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}$ (65)
 >
 $\left[{{l}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{1}{,}{{l}}_{{\mathrm{\mu }}}{}{{m}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}\right]$ (66)

You can query about their definition in the same way you can now query any other tensor:

 > $\mathrm{m_}\left[\mathrm{definition}\right]$
 ${{m}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{{n}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{{m}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{-1}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{{\mathrm{\nu }}}{+}{{l}}_{{\mathrm{\nu }}}{}{{n}}_{{\mathrm{\mu }}}{-}{{m}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\nu }}}{-}{{m}}_{{\mathrm{\nu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}$ (67)

Now that this definition depends on the signature:

 > $\mathrm{Setup}\left(\mathrm{signature}\right)$
 $\left[{\mathrm{signature}}{=}{\mathrm{- - - +}}\right]$ (68)

Change the signature, both in sign and placing time in position 1:

 > $\mathrm{Setup}\left(\mathrm{signature}=\mathrm{- +++}\right)$
 $\left[{\mathrm{signature}}{=}{\mathrm{- + + +}}\right]$ (69)

So, now ${m}_{\mathrm{\mu }}{\stackrel{&conjugate0;}{m}}_{\phantom{}\phantom{\mathrm{\mu }}}^{\phantom{}\mathrm{\mu }}=1$ instead of $-1$ and there is a change in the sign of the definition in terms of the metric ${g}_{\mathrm{\mu },\mathrm{\nu }}$ if compared with (67).

 > $\mathrm{m_}\left[\mathrm{definition}\right]$
 ${{m}}_{{\mathrm{\mu }}}{}{{l}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{{n}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{{m}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{0}{,}{{m}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{\phantom{{}}\phantom{{\mathrm{\mu }}}}^{\phantom{{}}{\mathrm{\mu }}}{=}{1}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{-}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{{\mathrm{\nu }}}{-}{{l}}_{{\mathrm{\nu }}}{}{{n}}_{{\mathrm{\mu }}}{+}{{m}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\nu }}}{+}{{m}}_{{\mathrm{\nu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}$ (70)

You can verify these tensorial identities using TensorArray; for example, for the last equation:

 > $\left[-1\right]$
 ${{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}{-}{{l}}_{{\mathrm{\mu }}}{}{{n}}_{{\mathrm{\nu }}}{-}{{l}}_{{\mathrm{\nu }}}{}{{n}}_{{\mathrm{\mu }}}{+}{{m}}_{{\mathrm{\mu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\nu }}}{+}{{m}}_{{\mathrm{\nu }}}{}{\stackrel{{&conjugate0;}}{{m}}}_{{\mathrm{\mu }}}$ (71)
 > $\mathrm{TensorArray}\left(,\mathrm{simplifier}=\mathrm{simplify}@\mathrm{expand}\right)$
 $\left[\begin{array}{cccc}{\mathrm{\epsilon }}{}{{A}}^{{2}}{=}{\mathrm{\epsilon }}{}{{A}}^{{2}}& {0}{=}{0}& {0}{=}{0}& {0}{=}{0}\\ {0}{=}{0}& {{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}\left({2}{}{{\mathrm{cosh}}\left({x}\right)}^{{2}}{-}{1}\right){=}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}\left({2}{}{{\mathrm{cosh}}\left({x}\right)}^{{2}}{-}{1}\right)& {-}{2}{}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}{\mathrm{cosh}}\left({x}\right){}{\mathrm{sinh}}\left({x}\right){=}{-}{2}{}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}{\mathrm{cosh}}\left({x}\right){}{\mathrm{sinh}}\left({x}\right)& {0}{=}{0}\\ {0}{=}{0}& {-}{2}{}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}{\mathrm{cosh}}\left({x}\right){}{\mathrm{sinh}}\left({x}\right){=}{-}{2}{}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}{\mathrm{cosh}}\left({x}\right){}{\mathrm{sinh}}\left({x}\right)& {{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}\left({2}{}{{\mathrm{cosh}}\left({x}\right)}^{{2}}{-}{1}\right){=}{{B}}^{{2}}{}{{ⅇ}}^{{2}{}{A}{}{x}}{}\left({2}{}{{\mathrm{cosh}}\left({x}\right)}^{{2}}{-}{1}\right)& {0}{=}{0}\\ {0}{=}{0}& {0}{=}{0}& {0}{=}{0}& {-}{\mathrm{\epsilon }}{=}{-}{\mathrm{\epsilon }}\end{array}\right]$ (72)

 >
 $\left[{\mathrm{tetradmetric}}{=}\left\{\left({1}{,}{2}\right){=}{-1}{,}\left({3}{,}{4}\right){=}{1}\right\}\right]$ (73)

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

 > $\mathrm{eta_}\left[\right]$

Note: this form of the null tetrad metric is consistent with time in position 1, a change done in (69)- and not with time in position 4 (default). You can in any case redefine the tetrad metric in any particular way also using .

 > $\mathrm{e_}\left[\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{=}\left[\begin{array}{cccc}{-}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}{}{A}& \frac{\sqrt{{2}}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}{}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}}& {-}\frac{\sqrt{{2}}{}{\mathrm{sinh}}\left({2}{}{x}\right){}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\\ {-}\frac{{I}}{{2}}{}\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}{}{A}& {-}\frac{\sqrt{{2}}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}{}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}}& \frac{\sqrt{{2}}{}{\mathrm{sinh}}\left({2}{}{x}\right){}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {0}\\ {0}& {0}& \frac{\sqrt{{2}}{}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& {-}\frac{\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}}{{2}}\\ {0}& {0}& \frac{\sqrt{{2}}{}{B}{}{{ⅇ}}^{{A}{}{x}}}{{2}{}\sqrt{{\mathrm{cosh}}\left({2}{}{x}\right)}}& \frac{\sqrt{{2}}{}\sqrt{{\mathrm{\epsilon }}}}{{2}}\end{array}\right]$ (74)

Compare this result with the orthonormal tetrad (63).

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

 > $\mathrm{e_}\left[\mathrm{definition}\right]$
 ${{𝔢}}_{{a}{,}{\mathrm{\mu }}}{}{{𝔢}}_{{b}\phantom{{\mathrm{\mu }}}}^{\phantom{{b}}{\mathrm{\mu }}}{=}{{\mathrm{\eta }}}_{{a}{,}{b}}$ (75)

Using the manipulation commands of the package:

 > $\mathrm{IsTetrad}\left(\right)$
 ${\mathrm{Type of tetrad: orthonormal}}$