Calculation of elements of spin groups using method of averaging in Clifford's geometric algebra

TL;DR

This paper introduces a generalized averaging method in Clifford's geometric algebra to compute spin group elements in arbitrary dimensions, extending Hestenes' approach for four dimensions.

Contribution

It develops explicit formulas for spin group elements in any dimension, broadening the applicability of geometric algebra in rotational computations.

Findings

Derived explicit formulas for spin group elements in arbitrary dimensions

Extended Hestenes' method from 4D to higher dimensions

Provided computational tools for rotors connecting different frames

Abstract

We present a method of computing elements of spin groups in the case of arbitrary dimension. This method generalizes Hestenes method for the case of dimension 4. We use the method of averaging in Clifford's geometric algebra previously proposed by the author. We present explicit formulas for elements of spin group that correspond to the elements of orthogonal groups as two-sheeted covering. These formulas allow us to compute rotors, which connect two different frames related by a rotation in geometric algebra of arbitrary dimension.

Equations

$$e_a{e}_b+e_b{e}_a=2{\eta }_{ab}e,$$

$$e_A=e_{a_1\dots a_k},1\le a_1<\cdots <a_k\le n.$$

$$C{\ell }_{p,q}=\underset{k=0}{\overset{n}{⨁}}C{\ell }_{p,q}^k,C{\ell }_{p,q}^{(0)}=\underset{k=0\mathrm{mod}2}{⨁}C{\ell }_{p,q}^k,C{\ell }_{p,q}^{(1)}=\underset{k=1\mathrm{mod}2}{⨁}C{\ell }_{p,q}^k.$$

$$\mathrm{O}(p,q):=\{P\in \mathrm{Mat}(n,\mathbb{R}):P^{\mathrm{T}}\eta P=\eta \}.$$

$$p_B^A=p_{b_1\dots b_k}^{a_1\dots a_k},a_1<\cdots <a_k,b_1<\cdots <b_k,$$

$${\beta }_a{\beta }_b+{\beta }_b{\beta }_a=2{\eta }_{ab}e$$

$$P=||p_a^b||\in \mathrm{O}(p,q)=\{P\in \mathrm{Mat}(n,\mathbb{R}):P^{\mathrm{T}}\eta P=\eta \}.$$

$${\beta }_a{\beta }_b+{\beta }_b{\beta }_a=p_a^c{p}_b^d(e_c{e}_d+e_d{e}_c)=2p_a^c{p}_b^d{\eta }_{cd}.$$

$${\beta }_a:=p_a^b{e}_b,P=||p_a^b||\in \mathrm{O}(p,q),$$

$${\beta }_A={\beta }_{a_1\dots a_k}:={\beta }_{a_1}\cdots {\beta }_{a_k},1\le a_1<\cdots <a_k\le n.$$

$${\beta }_{a_1\dots a_k}=p_{a_1\dots a_k}^{b_1\dots b_k}{e}_{b_1\dots b_k},$$

$${\beta }_A=p_A^B{e}_B,$$

$${\beta }_{a_1{a}_2}$$

$$p_{a_1}^b{p}_{a_2}^b{\eta }_{bb}=0,a_1<a_2,P=||p_a^b||\in \mathrm{O}(p,q).$$

$${\beta }_{a_1\dots a_m}={\beta }_{a_1\dots a_{m-1}}{\beta }_{a_m}=(p_{a_1\dots a_{m-1}}^{c_1\dots c_{m-1}}{e}_{c_1\dots c_{m-1}})(p_{a_m}^c{e}_c).$$

$$p_{a_1\dots a_m}^{b_1\dots b_m}=\sum _{j=1}^m(-1{)}^{m+j}{p}_{a_1\dots a_{m-1}}^{b_1\dots \check{b_j}\dots b_m}{p}_{a_m}^{b_j},$$

$${\beta }_{1\dots n}=\det (P)e_{1\dots n},\det P=\pm 1.$$

$$\det P=\pm 1,|p_{1\dots p}^{1\dots p}|\ge 1,|p_{p+1\dots n}^{p+1\dots ,n}|\ge 1,p_{1\dots p}^{1\dots p}=\frac{p_{p+1\dots n}^{p+1\dots n}}{\det P},$$

$$P=\left(\begin{array}[]{cc}A_{p\times p}&B_{p\times q}\\ C_{q\times p}&D_{q\times q}\\ \end{array}\right)\in{\rm O}(p,q)$$

$$P^{\mathrm{T}}\eta P=\eta ,B^{\mathrm{T}}B-D^{\mathrm{T}}D=-1,|\det (D)|=|p_{p+1\dots n}^{p+1\dots ,n}|\ge 1;$$

$$P\eta P^{\mathrm{T}}=\eta ,A{A}^{\mathrm{T}}-B{B}^{\mathrm{T}}=1,|\det (A)|=|p_{1\dots p}^{1\dots p}|\ge 1.$$

$$p_{1\dots p}^{1\dots p}=(p^{-1}{)}_{1\dots p}^{1\dots p}=\frac{p_{p+1\dots n}^{p+1\dots n}}{\det P},$$

$${\mathrm{O}}_{+}(p,q)$$

$${\mathrm{O}}_{-}(p,q)$$

$$\mathrm{SO}(p,q)$$

$${\mathrm{SO}}_{+}(p,q)$$

$$\mathrm{O}(n)$$

$$\mathrm{SO}(n)$$

$$\hat{U}:=U{|}_{e_a\to -e_a},U\in C{\ell }_{p,q}$$

$$\tilde{U}:=U{|}_{e_{a_1\dots a_k}\to e_{a_k}\dots e_{a_1}},U\in C{\ell }_{p,q}.$$

$${\Gamma }_{p,q}^{\pm }$$

$${\Gamma }_{p,q}^{+}$$

Full text

1. Introduction

In Clifford’s geometric algebra, it is convenient to describe rotations using elements of spin groups. Spin groups of arbitrary dimension are naturally realized in this algebra. Nowadays Clifford’s geometric algebra is widely used in physics, computer science, engineering, and other sciences.

Professor D. Hestenes presented the method of computing elements of spin group ${\rm Spin}_{+}(1,3)$ in the case of dimension $n=4$ in [7] (pp. 52–53). This method is mentioned in other papers and books (see, for example, [13], p. 130). We generalize this method to the case of arbitrary dimension $n=p+q$ and all spin groups ${\rm Spin}_{+}(p,q)$, ${\rm Spin}(p,q)$, ${\rm Pin}_{+}(p,q)$, ${\rm Pin}_{-}(p,q)$, ${\rm Pin}(p,q)$. We use the method of averaging in Clifford’s geometric algebra previously proposed by the author [17], [18] to do this.

There are other methods of calculation of elements of spin groups using exponentials and exterior exponentials of bivectors (see [7], [8], [4], [13], [14]). But all these methods work only in the cases of fixed dimensions, especially $n=3$ or $n=4$. In this paper, we present explicit formulas for the elements of spin groups, which work in the case of arbitrary $n$.

Note that some years ago we presented another method of computing elements of spin groups using generalized Pauli’s theorem. These results were presented at the conference AGACSE 2012 (La Rochelle, France, July 2012) and published in the Conference Proceedings in AACA [19]. In [19], we presented an algorithm (we had no explicit formulas) for computing elements of spin groups. In the new method presented in the current paper, using the method of averaging, we obtain explicit formulas for the elements of spin groups.

The paper is organized as follows. In Section 2, we discuss a formalism of Clifford’s geometric algebra and use it for the consideration of pseudo-orthogonal groups. We prove some auxiliary lemmas. In Section 3, we present a complete picture of five orthogonal groups and five corresponding spin groups in the case of arbitrary dimension. In Section 4, we discuss Hestenes method of computing elements of the group ${\rm Spin}_{+}(1,3)$. In Section 5, we generalize this method for the case of arbitrary dimension using the method of averaging in Clifford’s geometric algebra. In Section 6, we present an explicit formula for computing rotors, which connect two different frames related by a rotation in geometric algebra of arbitrary dimension.

2. Clifford algebras and pseudo-orthogonal group

Let us consider the real Clifford algebra ${C}\!\ell_{p,q}$, $p+q=n$, with the identity element $e$ and the generators $e_{a}$, $a=1,\ldots,n$, satisfying

[TABLE]

where $\eta=||\eta_{ab}||$ is the diagonal matrix with its first $p$ entries equal to $1$ and the last $q$ entries equal to $-1$ on the diagonal.

We use notation with ordered multi-indices $A$ for the basis elements of the Clifford algebra ${C}\!\ell_{p,q}$:

[TABLE]

We denote the length of multi-index $A$ by $|A|$. In the case of the identity element $e$, we have empty multi-index $\o$ of length [math]. We call the subspace of ${C}\!\ell_{p,q}$ of Clifford algebra elements, which are linear combinations of basis elements with multi-indices of length $|A|=k$, the subspace of grade $k$ and denote it by ${C}\!\ell^{k}_{p,q}$. We denote the projection operator onto subspace of the grade $k$ by $\pi_{k}$. We denote inverses of generators by $e^{a}:=\eta^{ab}e_{b}=(e_{a})^{-1}$, $a=1,\ldots,n$, and inverses of basis elements by $e^{A}=(e_{A})^{-1}$. The frame $e^{a}$, $a=1,\ldots,n$ is often called reciprocal frame for the frame $e_{a}$, $a=1,\ldots,n$.

Even and odd subspaces we denote by ${C}\!\ell^{(0)}_{p,q}$ and ${C}\!\ell^{(1)}_{p,q}$. We have

[TABLE]

Let us consider the pseudo-orthogonal group ${\rm O}(p,q)$, $p+q=n$:

[TABLE]

Denote by

[TABLE]

the minors of the matrix $P=||p^{a}_{b}||$. They are determinants of the submatrices formed by rows $a_{1}$, …, $a_{k}$ and columns $b_{1}$, …, $b_{k}$. In the particular case of multi-indices of length $1$ ($A=a$, $B=b$), the corresponding minor is just an element $p^{a}_{b}$ of the matrix $P$. In the case of empty multi-indices $A$ and $B$, the corresponding minor equals $1$ by definition. We use Einstein summation convention for ordered multi-indices too.

We have the following simple fact.

Lemma 2.1**.**

The set $\beta_{a}:=p_{a}^{b}e_{b}\in{C}\!\ell_{p,q}$ satisfies the following conditions

[TABLE]

if and only if

[TABLE]

Proof.

We have

[TABLE]

This means that $\beta_{a}\beta_{b}+\beta_{b}\beta_{a}=2\eta_{ab}e$ if and only if $p_{a}^{c}p_{b}^{d}\eta_{cd}=\eta_{ab}$. From the last condition, we obtain the definition of ${\rm O}(p,q)$ (2.2). ∎

For

[TABLE]

let us consider the elements

[TABLE]

Lemma 2.2**.**

For (2.3), we have

[TABLE]

where $p_{a_{1}\ldots a_{k}}^{b_{1}\ldots b_{k}}$ are minors of the matrix $P=||p_{a}^{b}||\in{\rm O}(p,q)$ and we have a sum over all ordered multi-indices $b_{1}\ldots b_{k}$ of length $k$ in (2.5).

Using our notations, we can rewrite (2.5) in the following way

[TABLE]

where we have a sum over all ordered multi-indices $B$ of the same length as the length of the multi-index $A$.

Proof.

For $k=1$, we have $\beta_{a}=p_{a}^{b}e_{b}$ by the definition. For $k=2$, we have

[TABLE]

where we use

[TABLE]

In the general case, the proof is by induction on $k$. Suppose that we have (2.5) for $k=m-1$. Let us prove it for $k=m$. We have

[TABLE]

Multiplying two sums, we obtain the element of grade $m-2$ (in the case $c\in\{c_{1},\ldots,c_{m-1}\}$) and the element of grade $m$ (in the case $c\notin\{c_{1},\ldots,c_{m-1}\}$).

The corresponding element of grade $m-2$ equals zero. Using (2.7) and Laplace expansion, we can prove it again by induction. We omit detailed proof because of its cumbersomeness.

The corresponding element of grade $m$ equals $p_{a_{1}\ldots a_{m}}^{b_{1}\ldots b_{m}}e_{b_{1}\ldots b_{m}}$ because of the Laplace expansion along one column of the corresponding minor:

[TABLE]

where $b_{1}\ldots\check{b_{j}}\ldots b_{m}$ is the ordered multi-index of length $m-1$, which is obtained from $b_{1}\ldots b_{m}$ by discarding $b_{j}$. ∎

Note that as particular case of (2.5) we get

[TABLE]

The conditions $\beta_{1\ldots n}=\pm e_{1\ldots n}$ mean that $2^{n}$ elements $\beta_{A}$ (2.4) are linearity independent and constitute a new basis of ${C}\!\ell_{p,q}$ (see pp. 289–290 in [21] or pp. 127–128 in [16]).

3. Complete picture of orthogonal and spin groups

For the convenience of the reader, we present a complete picture of five orthogonal groups and the corresponding five spin groups in the case of arbitrary dimension. For more details, see [20] and [3]. Often, only the group ${\rm Spin}_{+}(p,q)$ is considered, but sometimes other spin groups are also required for different applications.

Lemma 3.1**.**

For an arbitrary matrix $P\in{\rm O}(p,q)$, we have

[TABLE]

where $p^{1\ldots p}_{1\ldots p}$ and $p^{p+1\ldots,n}_{p+1\ldots n}$ are the corresponding minors of the matrix $P$.

Proof.

The first statement is trivial. For the matrix

[TABLE]

with the blocks $A,B,C,D$ of corresponding sizes, we have

[TABLE]

[TABLE]

From $\eta P^{\rm T}\eta=P^{-1}$ and the well-known formula for the minor of the inverse of a matrix (see, for example, [11], pp. 25 - 27), we obtain

[TABLE]

where $(p^{-1})^{1\ldots p}_{1\ldots p}$ is the corresponding minor of the matrix $P^{-1}$. ∎

The group ${\rm O}(p,q)$ has the following subgroups:

[TABLE]

For example, in the particular case of Minkowski space, we have Lorentz group ${\rm O}(1,3)$, special (or proper) Lorentz group ${\rm SO}(1,3)$, orthochronous (or time preserving) Lorentz group ${\rm O}_{+}(1,3)$, orthochorous (or parity preserving) Lorentz group ${\rm O}_{-}(1,3)$, special orthochronous Lorentz group ${\rm SO}_{+}(1,3)$.

In Euclidean cases ($p=0$ or $q=0$), we have only two orthogonal groups instead of five groups:

[TABLE]

We denote grade involution (main involution) in ${C}\!\ell_{p,q}$ by

[TABLE]

and reversion (anti-involution) by

[TABLE]

Denote by $M^{\times}$ the subset of invertible elements of any set $M$. Let us consider the Lipschitz group

[TABLE]

and its subgroup

[TABLE]

The following groups are called spin groups:

[TABLE]

Let us consider the twisted adjoint representation

[TABLE]

The following homomorphisms are surjective with the kernel $\{\pm 1\}$:

[TABLE]

It means that, for all $P=||p^{a}_{b}||\in{\rm O}(p,q)$, there exist $\pm S\in{\rm Pin}(p,q)$ such that

[TABLE]

and for the other groups similarly. The spin groups (3.2) are two-sheeted coverings of the corresponding orthogonal groups.

Our goal is to find out elements $\pm S\in{\rm Pin}(p,q)$ for each $P\in{\rm O}(p,q)$ in the case of arbitrary $p$ and $q$ using the relation (3.3).

4. Hestenes method

Let us consider the method proposed by D. Hestenes [7] for the case of dimension $n=4$, ${C}\!\ell_{1,3}$.

For each element $P=||p^{a}_{b}||\in{\rm SO}_{+}(1,3)$, there exist two elements $\pm S\in{\rm Spin}_{+}(1,3)$ such that

[TABLE]

The elements $\pm S$ can be found in the following way

[TABLE]

Let us discuss the plan of the proof of the formula (4.2). Multiplying both sides of the first equation (4.1) on the right by $e^{a}=(e_{a})^{-1}$, we obtain

[TABLE]

We have the following well-known formula (see [13], [15])

[TABLE]

We have $S\in{\rm Spin}_{+}(1,3)$, so $S=\pi_{0}(S)+\pi_{2}(S)+\pi_{4}(S)$. Using (4.4), we get from (4.3)

[TABLE]

Let us take reversion of both sides of (4.5). We get

[TABLE]

Multiplying both sides of (4.6) by both sides of (4.5) and using $\widetilde{S}S=e$, we obtain

[TABLE]

Both sides of this equation belong to ${C}\!\ell^{0}_{1,3}\oplus{C}\!\ell^{4}_{1,3}\cong{\mathbb{C}}$. Taking square root of both sides of (4.7), we get

[TABLE]

Substituting (4.8) into (4.5), we obtain (4.2).

Note that this method works only in the case of dimension $n=4$ for the matrices $P=||p^{a}_{b}||\in{\rm SO}_{+}(1,3)$ with additional condition

[TABLE]

The condition (4.9) is equivalent to the condition

[TABLE]

for the corresponding element $S\in{\rm Spin}_{+}(1,3)$.

In the next section, we will generalize this method to the case of arbitrary $n=p+q$. We will use some other operators instead of (4.4) to do this.

5. Generalization of Hestenes method

We have the following new theorems.

Theorem 5.1**.**

Let us consider the real Clifford algebra ${C}\!\ell_{p,q}$ with even $n=p+q$. Let $P\in{\rm SO}(p,q)$ be an orthogonal matrix such that

[TABLE]

Then we can find the elements $\pm S\in{\rm Spin}(p,q)$ that correspond to $P=||p^{b}_{a}||\in{\rm SO}(p,q)$ as two-sheeted covering $Se_{a}S^{-1}=p_{a}^{b}e_{b}$ in the following way:

[TABLE]

where

[TABLE]

and the sign

[TABLE]

depends on the component of the orthogonal group ${\rm SO}(p,q)$ (or the corresponding component of the spin group ${\rm Spin}(p,q)$).

Theorem 5.2**.**

Let us consider the real Clifford algebra ${C}\!\ell_{p,q}$ with odd $n=p+q$. Let $P\in{\rm O}(p,q)$ be an orthogonal matrix such that

[TABLE]

Then we can find the elements $\pm S\in{\rm Pin}(p,q)$ that correspond to $P=||p^{b}_{a}||\in{\rm O}(p,q)$ as two-sheeted covering $\widehat{S}e_{a}S^{-1}=p_{a}^{b}e_{b}$ in the following way:

[TABLE]

where

[TABLE]

and the sign

[TABLE]

depends on the component of the orthogonal group ${\rm O}(p,q)$ (or the corresponding component of the group ${\rm Pin}(p,q)$).

The conditions (5.1) and (5.3) are equivalent to the condition $\pi_{{\rm Cen}}(S)\neq 0$ for the corresponding element $S\in{\rm Pin}(p,q)$, where $\pi_{{\rm Cen}}$ is the projection onto the center of the Clifford algebra

[TABLE]

Proof.

In the proof of Theorems 5.1 and 5.2, we use Reynolds operators of Salingaros vee group (see [17])

[TABLE]

We have the following relation (3.3) between orthogonal matrix $P=||p^{b}_{a}||\in{\rm O}(p,q)$ and the corresponding two elements of spin group $\pm S\in{\rm Pin}(p,q)$. We can rewrite this relation in the following way

[TABLE]

because of the relation between parity of the element of spin group and the determinant of the corresponding orthogonal matrix (see Section 3). Multiplying (5.9) by itself several times and using Lemma 2.2, we get

[TABLE]

Multiplying both sides of (5.10) on the right by $e^{A}=(e_{A})^{-1}$, we obtain

[TABLE]

Here we have a sum over multi-indices $A$, $B$ of the same length $|A|=|B|$. We denote the right side of (5.11) by

[TABLE]

and get

[TABLE]

Using $M\in{C}\!\ell^{(0)}_{p,q}$ and $S\in{C}\!\ell^{(0)}_{p,q}\cup{C}\!\ell^{(1)}_{p,q}$, we get $\pi_{{\rm Cen}}(S^{-1})\in{C}\!\ell^{(0)}_{p,q}\cup{C}\!\ell^{(1)}_{p,q}$. The condition $M\neq 0$ is equivalent to the condition $\pi_{{\rm Cen}}(S)\neq 0$ for the element $S\in{\rm Pin}(p,q)$ because of (5.13) and these facts.

We have the following well-known fact (see, for example [13]) about the center of the Clifford algebra ${C}\!\ell_{p,q}$:

[TABLE]

We have

[TABLE]

and

[TABLE]

Note that ${\mathbb{R}}$ and ${\mathbb{C}}$ are fields, but the set of double numbers (split-complex numbers) ${\mathbb{R}}\oplus{\mathbb{R}}$ is not a field. It is an associative commutative algebra of dimension 2 over real numbers and it has zero divisors.

Now let us consider the cases of even and odd $n$ separately.

In the case of even $n$, we have $\pi_{{\rm Cen}}(S^{-1})=\pi_{0}(S^{-1})$. Let us take reversion of both sides of (5.13). We obtain

[TABLE]

Multiplying both sides of (5.17) by both sides of (5.13) we obtain

[TABLE]

We have $\alpha:=\widetilde{S}S=\pm e\in{C}\!\ell^{0}_{p,q}$ and get

[TABLE]

Taking square root and substituting this expression into (5.13), we get (5.2).

In the case of odd $n$, we have $\pi_{{\rm Cen}}(S^{-1})=\pi_{0}(S^{-1})+\pi_{n}(S^{-1})$ and

[TABLE]

In the case $n=1\mod 4$, taking reversion of both sides of (5.13), we get

[TABLE]

Multiplying both sides of (5.26) by both sides of (5.13), we obtain

[TABLE]

We have $\alpha:=\widetilde{S}S=\pm e\in{C}\!\ell^{0}_{p,q}$ and get

[TABLE]

Taking square root and substituting this expression into (5.13), we get (5.4) for the first case (5.7).

In the case $n=3\mod 4$, taking superposition of reversion and grade involution (it is called Clifford conjugation) of both sides of (5.13), we get

[TABLE]

Multiplying both sides of (5.29) by both sides of (5.13), we obtain

[TABLE]

We have $\alpha:=\widehat{\widetilde{S}}S=\pm e\in{C}\!\ell^{0}_{p,q}$ and get

[TABLE]

Taking square root and substituting this expression into (5.13), we get (5.4) for the second case (5.7).

The theorems are proved. ∎

6. Calculation of rotors in geometric algebra

Let us consider the particular case of Theorems 5.1 and 5.2 for the elements of the group ${\rm Spin}_{+}(p,q)$ and the corresponding group ${\rm SO}_{+}(p,q)$. Elements of ${\rm Spin}_{+}(p,q)$ are often called rotors and have wide application in geometric algebra (see, for example, [7], [8], [4], [2], [6], [5], [12]). We obtain the following statement.

Corollary 6.1**.**

Let we have

[TABLE]

where two frames $e_{a}$ and $\beta_{a}$, $a=1,\ldots,n$, are related by a rotation.

If

[TABLE]

then

[TABLE]

Using our previous notation, we can write

[TABLE]

[TABLE]

In the particular case of $n=3$, we have

[TABLE]

because $\beta_{123}=e_{123}\in{\rm Cen}({C}\!\ell_{p,q})$. We can remove scalar “$2$” in (6.2) because of normalization and finally obtain the following well-known formula for the rotor $S$ (see, for example, p. 103 in [4] or p. 72 in [2])

[TABLE]

This formula is widely used in different applications of geometric algebra. Corollary 6.1 generalizes (6.3) to the case of arbitrary $n$.

The results of this paper were reported at the conference AGACSE 2018 (Campinas, Brazil, July 2018). We hope that these results will be useful for different applications in computer science, robotics, and engineering. There are well-known methods of calculating of rotors in dimensions $n=3$ and $4$, but we often need geometric algebra of higher dimensions. For example, the conformal geometric algebra ${C}\!\ell_{4,1}$ of dimension $n=5$ is widely used in different applications (see, for example, [1], [6], [9], [10]). The results of this paper allow us to calculate rotors in arbitrary dimension $n=p+q\geq 1$.

Acknowledgment

The author is grateful to N. Marchuk and participants of the AGACSE 2018 Conference for fruitful discussions. The author is grateful to two anonymous reviewers for their careful reading of the paper and helpful comments. This work is supported by the Russian Science Foundation (project 18-71-00010).