Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Formation of Fuzzy Patterns in Logical Analysis of Data Using a Multi-Criteria Genetic Algorithm
Previous Article in Journal
3D Copyright Protection Based on Binarized Computational Ghost Imaging Encryption and Cellular Automata Transform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalization of Quaternions and Their Applications

1
Department of Electrical Engineering and Computer Science, University of Cincinnati, Cincinnati, OH 45221, USA
2
Independent Researcher, Lockesley Avenue, Conisbrough, South Yorkshire DN12 2AA, UK
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(3), 599; https://doi.org/10.3390/sym14030599
Submission received: 26 November 2021 / Revised: 23 December 2021 / Accepted: 8 March 2022 / Published: 17 March 2022
(This article belongs to the Section Physics)

Abstract

:
There are a total of 64 possible multiplication rules that can be defined starting with the generalized imaginary units first introduced by Hamilton. Of these sixty-four choices, only eight lead to non-commutative division algebras: two are associated to the left- and right-chirality quaternions, and the other six are generalizations of the split-quaternion concept first introduced by Cockle. We show that the 4 × 4 matrix representations of both the left- and right-chirality versions of the generalized split-quaternions are algebraically isomorphic and can be related to each other by 4 × 4 permutation matrices of the C 2 × C 2 group. As examples of applications of the generalized quaternion concept, we first show that the left- and right-chirality quaternions can be used to describe Lorentz transformations with a constant velocity in an arbitrary spatial direction. Then, it is shown how each of the generalized split-quaternion algebras can be used to solve the problem of quantum-mechanical tunneling through an arbitrary one-dimensional (1D) conduction band energy profile. This demonstrates that six different spinors ( 4 × 4 matrices) can be used to represent the amplitudes of the left and right propagating waves in a 1D device.

1. Introduction

Hamilton first introduced the concept of a quaternion in 1843 [1]. Over the years, quaternions have been applied in many areas including control theory, signal processing, robotics, number theory, computer vision and graphics, flight dynamics and navigation systems of aircraft, orbital mechanics, bioinformatics, molecular dynamics, robotics, quantum mechanics, and crystallographic structure analysis, among others [2,3,4,5,6].
Girard described the use of the quaternion group in modern physics [7]. Adler has analyzed in great detail how quaternions can be used to outline a new formulation of quantum mechanics and quantum field theory [5]. Recently, Cahay and Morris have shown how a quaternion representation of the Pauli spinor [8,9] can be used to reformulate a quaternionic approach to quantum computing [10] and the Pauli–Schrödinger equation [11,12].
A few years after the pioneering work on quaternions by Hamilton, James Cockle introduced the concept of split-quaternions, also referred to as coquaternions, which are also elements of a four-dimensional associative algebra [13]. Contrary to quaternions, split-quaternions have nontrivial zero divisors, idempotents, and nilpotent elements [14,15,16]. Split-quaternions have found important applications in many fields such as spatial geometry, physics, quantum mechanics, signal processing, and cryptography, among others [17,18,19,20,21,22,23,24,25,26]. Thakur and Tripathi proposed a new multidimensional public key cryptosystem using split-quaternion algebra [27]. Babaarslan and Yayli derived some relations between split-quaternions and space-like contact slope surfaces in Minkowski three-space [28]. Pop and Cretu used a split-quaternion formalism to formulate a transfer matrix approach to 1D longitudinal elastic wave propagation through multilayered media [29].
In Section 2, it is shown that there are a total of 64 possible multiplication rules that can be defined starting with the imaginary units first introduced by Hamilton to define the fundamental concept of quaternion. Of these sixty-four possibilities, only eight lead to non-commutative division algebras, two of which are associated to the left- and right-chirality quaternions, while the other six are generalizations of the split-quaternion concept first introduced by Cockle [13]. For the latter, as shown by Morris [2], their 4 × 4 real matrix representations form the six four-dimensional representations of the order eight dihedral group. The 4 × 4 of both the left- and right-chirality versions of the generalized split-quaternions are algebraically isomorphic and can be related to each other by 4 × 4 permutation matrices of the C 2 × C 2 group.
In Section 3, both left- and right-chirality quaternions are used to describe Lorentz transformations or boosts in an arbitrary direction. In Section 4, the split-quaternion formalism is used to reformulate the transmission matrix technique, which is a well-known approach to tackle tunneling problems through an arbitrary conduction band energy profile in quantum mechanics. In this approach, the transmission and reflection probabilities through an arbitrary conduction band energy profile can then be easily derived by multiplying real 4 × 4 matrices. The approach is illustrated by revisiting the tunneling problems bt using a one-dimensional (1D) delta-scatterer and by examining the resonant tunneling structure composed of two and three delta-scatterers separated by a distance L. The split-quaternion formalism is also related to the transfer matrix approach to solve the tunneling problem [30]. The split-quaternion associated to a barrier with finite width and height is then derived and used to calculate transmission probability by using realistic resonant tunneling structures. It is shown that six different spinors ( 4 × 4 matrices) can be used to represent the amplitudes of the left and right propagating waves in a device and tackle tunneling problems in quantum mechanics.
In Section 5, an interesting connection between the non-commutative split-quaternion division algebras and the uniqueness of our four-dimensional spacetime is pointed out by introducing a super split-quaternion formed of a superposition of either the left- or right-chirality versions of the generalized split-quaternions. It is shown that the norm of this super split-quaternion has a form identical to the distance function of our four-dimensional spacetime. This paves a way to an interpretation of our four-dimensional spacetime and its physical properties based on a superposition principle of fundamental structures of our physical world, which can be represented by division algebras.
Finally, Section 5 contains our conclusions.

2. Generalization of the Quaternion Concept

We first describe for the sake of completeness the well-known 4 × 4 matrix representations of the quaternion and split-quaternion concepts first introduced by Hamilton [1] and Cockle [13], respectively. We then show that out of the sixty-four possible generalizations of the quaternion concept, only eight can be used to define algebras that satisfy the axioms of proper division algebras.

2.1. Quaternions and Their 4 × 4 Matrix Representations

In the original work of Hamilton [1], quaternions are defined as follows:
Q = a + i b + j c + k d ,
where a, b, c, and d are real numbers; the quaternion units i, j, k obey the following multiplication rules.
i 2 = j 2 = k 2 = 1 ,
i j = k , j i = k , j k = i , k j = i , k i = j , i k = j .
These commutation rules are isomorphic to SU(2). These quaternions are usually referred to as left-chirality quaternions [2,11].
Expressing the quaternions in a row format with the four components a, b, c, and d, i.e., Q = ( a , b , c , d ) , and using the multiplication rules (Equations (2) and (3)), the following quaternion multiplication rule is obtained.
Q 1 Q 2 = ( a 1 , b 1 , c 1 , d 1 ) ( a 2 , b 2 , c 2 , d 2 ) = ( a 1 a 2 b 1 b 2 c 1 c 2 d 1 d 2 , a 1 b 2 + b 1 a 2 + c 1 d 2 d 1 c 2 , a 1 c 2 b 1 d 2 + c 1 a 2 + d 1 b 2 , a 1 d 2 + b 1 c 2 c 1 b 2 + d 1 a 2 ) .
This last equation shows that quaternion multiplication is non-commutative, as can readily be seen by swapping the indices 1 and 2 on the right-hand side of Equation (4).
The quaternion Q = ( a , b , c , d ) defined above is referred to as the left-chirality quaternion [2] and is sometimes represented by using a 4 × 4 matrix [2,11].
Q L < = > a b c d b a d c c d a b d c b a .
The determinant of the 4 × 4 matrix in Equation (5) is equal to d e t ( Q L ) = ( a 2 + b 2 + c 2 + d 2 ) 2 and is related to the norm of the quaternion Q, which is equal to Q L Q L ¯ = a 2 + b 2 + c 2 + d 2 , where Q L ¯ = a b i c j d k is the conjugate of Q L .
If, on the other hand, we assume that the quaternion units i, j, and k satisfy the following multiplication rules:
i 2 = j 2 = k 2 = 1
i j = k , j i = k , j k = i , k j = i , k i = j , i k = j ,
then the signs of the products of the complex units are opposite to the signs used in Equation (3) to define the left-chirality quaternions, and we obtain the so-called right-chirality quaternions, which were already known to Hamilton [1,2].
The right-chirality quaternion has a 4 × 4 matrix representation given by the following.
Q R < = > a b c d b a d c c d a b d c b a .
The elements marked as bold in the expression of Q R have the opposite signs to the same matrix elements in the expression of Q L in Equation (5).

2.2. Split-Quaternions and Their 4 × 4 Matrix Representation

A split-quaternion is defined as follows:
S Q = w + i x + j y + k z ,
where w, x, y, and z are real numbers and the split-quaternion units i, j, and k satisfy the following multiplication rules:
i 2 = 1 , j 2 = k 2 = 1 ,
and the Hamilton multiplication rules for left-chirality quaternions given by Equation (3).
The set (1, i, j, k, −1, −i, −j, −k) forms a group under split-quaternion multiplication, which is isomorphic to dihedral group D 4 and is the symmetry group of the square [2]. If we represent the split-quaternions in a row format with the four components w, x, y, and z, i.e., S Q = ( w , x , y , z ) , and take into account the multiplication rules (Equations (3) and (10)) given above; the following rule for multiplying split-quaternions can easily be derived:
S Q 1 S Q 2 = ( w 1 , x 1 , y 1 , z 1 ) ( w 2 , x 2 , y 2 , z 2 ) = ( w 1 w 2 x 1 x 2 + y 1 y 2 + z 1 z 2 , w 1 x 2 + w 2 x 1 y 1 z 2 + y 2 z 1 , w 1 y 2 x 1 z 2 + y 1 w 2 + z 1 x 2 , w 1 z 2 + x 1 y 2 y 1 x 2 + z 1 w 2 ) ,
showing that the product of two split-quaternions is non-commutative.
The split-quaternion S Q = ( w , x , y , z ) can be represented using the following 4 × 4 real matrix [2].
S Q < = > w x y z x w z y y z w x z y x w .
The determinant of this last 4 × 4 matrix is equal to the following:
d e t ( S Q ) = ( w 2 + x 2 y 2 z 2 ) 2 ,
and is related to the norm of the split-quaternion SQ , which is equal to S Q S Q ¯ = w 2 + x 2 y 2 z 2 , where S Q ¯ = w i x j y k z is the conjugate of SQ .
The multiplication rule of split-quaternions given in Equation (11) can then be easily generated by multiplying 4 × 4 matrices associated to two split-quaternions S Q 1 and S Q 2 using the rule of matrix multiplication.

2.3. The Generalized Quaternion Concept

Hereafter, we address the question whether it is possible to extend the concept of quaternion, as defined by Equation (1), for which the quaternion units would satisfy the only possible multiplication rules listed in Table 1.
There are only 64 possible combinations for the products i 2 , j 2 , k 2 , ij, jk, and ki, corresponding to the fact that i 2 , j 2 , and k 2 can only take the values ± 1 and the products ij, jk, and ki can only be equal to ± k , ± i , and ± j , respectively.
We proceed to find a 4 × 4 matrix representation of all 64 possible quaternion generalizations using the following procedure. The most general form G of the product of two generalized quaternions will be of the following form:
G = G 1 G 2 ,
G = ( a 1 + b 1 i + c 1 j + d 1 k ) ( a 2 + b 2 i + c 2 j + d 2 k ) = a 1 a 2 + a 1 b 2 i + a 1 c 2 j + a 1 d 2 k + b 1 a 2 i + b 1 b 2 i 2 + b 1 c 2 i j + b 1 d 2 i k + c 1 a 2 j + c 1 b 2 j i + c 1 c 2 j 2 + c 1 d 2 j k + d 1 a 2 k + d 1 b 2 k i + d 1 c 2 k j + d 1 d 2 k 2 ,
where the different terms can be written explicitly according to which one of the sixty-four possible multiplication rules we select from the table above for imaginary units.
Proceeding with the procedure for the 64 potential generalized quaternions with the multiplication rules listed in Table 1 results in only eight acceptable division algebras that are closed under multiplication. Their multiplication rules and 4 × 4 matrix representations are provided in Table 2. The last column in the table below gives the values of the determinant of the 4 × 4 matrix associated with each of the generalized quaternions. The latter are equal to the square of the norm associated with each generalized quaternion. These norms are defined as follows: for each of the eight generalized quaternions, its complex conjugate is defined as follows:
Q ¯ = ( a , b , c , d ) ,
and norm N of the generalized quaternion is then given by the following.
N 2 = Q Q ¯ .
Hence, for Q = a + ib + jc + kd, Q ¯ = aibjckd, we obtain the following.
N 2 = a 2 i 2 b 2 j 2 c 2 k 2 d 2 i j ( b c ) j i ( b c ) i k ( b d ) k i ( b d ) j k ( c d ) k j ( c b ) .
Taking into account the fact that, for each of the eight division algebras identified in Table 2, the following anti-commutation rules apply.
i j = j i ; i k = k i ; j k = k j .
Hence, we obtain the following.
N 2 = a 2 i 2 b 2 j 2 c 2 k 2 d 2 .
In Table 2, the two top generalized quaternions have 4 × 4 matrix representations given by Equation (5) and Equation (8), respectively, and are referred to as left- and right-chirality quaternions, respectively. The characterization of the other six algebras as the generalization of the split-quaternion introduced in Section 2.3 and their left- or right-chirality characters will be discussed in Section 4.5. The values of the determinants of the 4 × 4 matrices listed in Table 2 are found to be equal to the fourth power of the norm of the generalized quaternions.
To prevent zero-divisors in the division algebras, the norms of the generalized quaternions must be greater than zero. It turns out that only two out of the eight 4 × 4 matrices always have a norm greater than zero, with both of their determinants being described as follows.
d e t ( G ) = ( a 2 + b 2 + c 2 + d 2 ) 2 , { a , b , c , d 0 }
The associated matrices correspond to the left- and right-chirality quaternions, which are discussed in Section 2.1 and Section 2.2. The other six division algebras in Table 2 have the square of their norms equal to one of the following forms: a 2 b 2 + c 2 d 2 , a 2 b 2 c 2 + d 2 , and a 2 + b 2 c 2 d 2 ; they constitute generalizations of the split-quaternion introduced by Cockle [13].
Since the norms associated with these six algebras contain terms with negative signs, these division algebras have zero divisors. As discussed by Morris [2], the exponential forms of their 4 × 4 matrix representations have a positive determinant and can be, therefore, used as proper division algebras. The 4 × 4 real matrix representations of those six division algebras consist of the six four-dimensional representations of the order eight dihedral group [2].
Next, we show on a specific example why any of the other 56 possible generalized quaternions cannot be used as elements of a proper division algebra. The proof for all other 55 cases proceed in a similar fashion. Particularly, for those 56 generalized quaternions, their multiplication rules are such that their quaternion units (i, j, k) do not obey the anti-commuting principle, i.e., they do not satisfy the following commutation rules, which are satisfied by the eight algebras listed in Table 2.
i j = j i , j k = k j , k i = i k .

2.4. Proof That 56 Other Generalized Quaternions Are Not Proper Division Algebras

Hereafter, we consider the product of two generalized quaternions of the following form:
Q 1 = a 1 + i b 1 + j c 1 + k d 1 , Q 2 = a 2 + i b 2 + j c 2 + k d 2 ,
for which the following multiplication rules are used:
i 2 = 1 , j 2 = 1 , k 2 = 1 , i j = k , j k = i , k i = j ,
from which we derive that result that the following additional relations are satisfied.
j i = i j = k , j k = k j = i , k i = i k = j .
There is actually an inconsistency in the multiplication rules (23), as evidenced by the fact that they result in the following contradiction: j i = ( k i ) ( j k ) = k 2 k = k and j i = j 2 k = k .
Using a similar simple argument for the other 55 possible quaternion generalizations, it can be shown that the only generalized quaternions for which their algberas are closed under multiplication are listed in Table 2.
The simple argument does not explicitly show the impossibility to find a 4 × 4 matrix expression for the generalized quaternions satisfying multiplication rules (23), which render them closed under multiplication. The procedure outlined below was actually used to derive the explicit 4 × 4 matrix representations of the only eight divisions algebras listed in Table 2, which are closed under multiplication.
Using the multiplication rules (23), the product, P = Q 1 Q 2 , is found to be as follows.
P = ( a 1 + i b 1 + j c 1 + k d 1 ) ( a 2 + i b 2 + j c 2 + k d 2 ) = ( a 1 a 2 + b 1 b 2 c 1 c 2 + d 1 d 2 ) + i ( c 1 d 2 d 1 c 2 + b 1 a 2 + a 1 b 2 ) j ( a 2 c 1 + a 1 c 2 + d 1 b 2 + b 1 d 2 ) + k ( a 2 d 1 + a 1 d 2 c 1 b 2 b 1 c 2 ) .
According to this last equation, the real and imaginary components of P can be obtained as the first row of the product of the following two 4 × 4 matrices.
a 1 b 1 c 1 d 1 × × × × × × × × × × × × a 2 b 2 c 2 d 2 b 2 a 2 d 2 c 2 c 2 d 2 a 2 b 2 d 2 c 2 b 2 a 2
Indeed, the elements of the top row of the 4 × 4 matrix resulting from the product above are given by the following:
P ( 1 , 1 ) = a 1 a 2 + b 1 b 2 c 1 c 2 + d 1 d 2 , P ( 1 , 2 ) = a 1 b 2 + b 1 a 2 + c 1 d 2 d 1 c 2 , P ( 1 , 3 ) = a 1 c 2 + b 1 d 2 + c 1 a 2 + d 1 b 2 , P ( 1 , 4 ) = a 1 d 2 b 1 c 2 c 1 b 2 + d 1 a 2 ,
where P(1,1) agrees with the real part of the product in Equation (25), and P(1,2), P(1,3), and P(1,4) agree with the three imaginary components found above.
For the type of generalized quaternions to form a proper division algebra, the rest of the matrix elements represented by crosses in Equation (26) should be of the following form:
a 1 b 1 c 1 d 1 b 1 a 1 d 1 c 1 c 1 d 1 a 1 b 1 d 1 c 1 b 1 a 1
of the same canonical form as the 4 × 4 matrix associated to Q 2 . Since the product must be closed under multiplication, the 4 × 4 matrix elements associated with the product of Q 1 and Q 2 must satisfy the following relations.
P ( 1 , 1 ) = P ( 2 , 2 ) = P ( 3 , 3 ) = P ( 4 , 4 ) , P ( 2 , 1 ) = P ( 1 , 2 ) , P ( 3 , 1 ) = P ( 1 , 3 ) , P ( 4 , 1 ) = P ( 1 , 4 ) , P ( 2 , 3 ) = P ( 3 , 2 ) = P ( 1 , 4 ) , P ( 2 , 4 ) = P ( 1 , 3 ) = P ( 4 , 2 ) , P ( 3 , 4 ) = P ( 1 , 2 ) = P ( 4 , 3 ) .
By performing the matrix multiplication of the 4 × 4 matrices associated with Q 1 and Q 2 :
a 1 b 1 c 1 d 1 b 1 a 1 d 1 c 1 c 1 d 1 a 1 b 1 d 1 c 1 b 1 a 1 a 2 b 2 c 2 d 2 b 2 a 2 d 2 c 2 c 2 d 2 a 2 b 2 d 2 c 2 b 2 a 2 ,
we obtain the following expressions for the diagonal elements associated with product P:
P ( 1 , 1 ) = a 1 a 2 + b 1 b 2 c 1 c 2 + d 1 d 2 , P ( 2 , 2 ) = a 1 a 2 + b 1 b 2 + c 1 c 2 + d 1 d 2 , P ( 3 , 3 ) = a 1 a 2 b 1 b 2 c 1 c 2 + d 1 d 2 , P ( 4 , 4 ) = a 1 a 2 b 1 b 2 + c 1 c 2 + d 1 d 2 .
They are not equal to the following.
P ( 1 , 1 ) = a 1 a 2 + b 1 b 2 c 1 c 2 + d 1 d 2 .
Similarly, we find the following.
P ( 2 , 1 ) = a 2 b 1 + a 1 b 2 c 2 d 1 c 1 d 2 , P ( 1 , 2 ) = a 1 b 2 + a 2 b 1 + c 1 d 2 d 1 c 2 , P ( 3 , 1 ) = a 2 c 1 + b 2 d 1 a 1 c 2 b 1 d 2 , P ( 1 , 3 ) = a 1 c 2 + c 1 a 2 b 1 d 2 + d 1 b 2 .
Thus, P ( 2 , 1 ) P ( 1 , 2 ) and P ( 3 , 1 ) P ( 1 , 3 ) .
Therefore, the 4 × 4 matrices of form (28) are not closed under multiplication.
The fact that the generalized quaternions of the form (28) with the multiplication rules (23) and (24) do not form a proper division algebra is related to the fact that the following is the case.
i j = j i , k i = i k , j k = k j ,
In other words, the (i, j) and (i, k) complex units commute but the (j, k) units anti-commute.
Using a similar proof for the other 55 possible quaternion generalizations, it can be shown that no 4 × 4 matrix representations can be found for those generalized quaternions, which would render them closed under multiplication. The only division algebras for which its 4 × 4 matrix representations are closed under multiplication are those listed in Table 2.

2.5. Isomorphism between the Non-Commutative Division Algebra

In this section, we establish the isomorphism between the split-quaternions listed in Table 2. Introducing the permutation matrix p 23 :
p 23 = 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1 ,
we have the following.
a c b d c a d b b d a c d b c a = 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1 a b c d b a d c c d a b d c b a 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1
Since p 23 is equal to its own inverse, Equation (36) is a similarity transformation connecting Q + + k , i , j ( a , b , c , d ) and Q + + k , i , j ( a , c , b , d ) split-quaternions. The latter has the 4 × 4 matrix representation listed in Table 2, with elements b and c interchanged.
The Q + + k , i , j ( a , b , c , d ) split-quaternion can be transformed into all five others using the following similarity transformation:
Q + + k , i , j ( a , c , b , d ) = p 23 Q + + k , i , j ( a , b , c , d ) p 23 , Q + + k , i , j ( a , c , d , b ) = p 34 p 23 Q + + k , i , j ( a , b , c , d ) p 23 p 34 , Q + + k , i , j ( a , d , c , b ) = p 24 Q + + k , i , j ( a , b , c , d ) p 24 , Q + + k , i , j ( a , b , d , c ) = p 34 Q + + k , i , j ( a , b , c , d ) p 34 , Q + + k , i , j ( a , d , b , c ) = p 24 p 23 Q + + k , i , j ( a , b , c , d ) p 23 p 24 ,
with the permutation matrices are given by the following.
p 34 = 1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0 ; p 24 = 1 0 0 0 0 0 0 1 0 0 1 0 0 1 0 0 .
A similar isomorphism can be established between the left-chirality and right-chirality quaternions in Table 2. In fact, using the permutation matrices defined above, we have the following.
H r i g h t ( a , c , b , d ) = p 23 H l e f t ( a , b , c , d ) p 23 , H r i g h t ( a , d , c , b ) = p 24 H l e f t ( a , b , c , d ) p 24 , H r i g h t ( a , b , d , c ) = p 34 H l e f t ( a , b , c , d ) p 34 .
Notice that the following is the case.
p 23 p 23 = p 24 p 24 = p 34 p 34 = I .
Notice that all the transformations listed in Equations (37) and (39) are similar transformations since each permutation matrix listed above is equal to its own inverse.
Finally, as shown in Table 2, the determinant of each of 4 × 4 matrix representation of the eight proper division algebras is equal to the fourth power of its norm of the corresponding generalized split-quaternion. This last property can be used to show that the left- and right-chirality quaternions cannot be mapped by using a unitary transformation to any of the six generalized quaternions listed at the end of Table 2.
Indeed, if such a transformation existed, there would be a 4 × 4 unitary matrix U such that the following is the case:
U 1 Q U = H l e f t o r H r i g h t
for each of the six Q’s associated with the split-quaternions listed in Table 2. Hence, this would imply the following:
d e t ( U 1 Q U ) = d e t ( Q ) = d e t ( H l e f t ) = d e t ( H r i g h t )
but d e t ( Q ) is equal to the following:
( a 2 + b 2 c 2 d 2 ) 2 , ( a 2 b 2 + c 2 d 2 ) 2 , o r ( a 2 b 2 c 2 + d 2 ) 2 ,
for the six generalized split-quaternions, whereas the following is the case.
d e t ( H l e f t ) = d e t ( H r i g h t ) = ( a 2 + b 2 + c 2 + d 2 ) 2 .
The unitary transformation U only exists for the trivial case where b = c = d = 0 .

3. Examples of Generalized Quaternion Applications

Recently, Cahay et al. [11] used a quaternion formalism to solve the time dependent Pauli–Schrödinger equation. They showed that both the unitary operator describing the time evolution of the Pauli spinor and the Pauli spinor itself can be written either as left- or right-chirality quaternions [11]. They also show a method for reformulating different quantum gates and the evolution of a qubit on a Bloch sphere using quaternions. Quaternions were also used to derive the Larmor precession of an electron spin in an uniform magnetic field and revisit the issue of measurements in quantum mechanics.
In [12], it was shown that a quaternion algebraic approach can be used to solve the spatially invariant time-dependent Pauli–Schrödinger equation. It was applied to derivation of the Rabi formula and to the calculation of the transmission probabilities by using a one-dimensional delta scatterer and by using a resonant tunneling structure composed of two one-dimensional quaternionic delta-scatterers in series.
Hereafter, we illustrate how both the left- and right-chirality quaternion formalisms can be used to describe Lorentz boosts in special relativity [31]. The Lorentz transformation or boost is a linear transformation from one coordinate frame in spacetime to another moving at a constant velocity relative to the first. The connections between the four dimensional coordinates of two Lorentzian frames where the second (with prime coordinates) moves with constant velocity confined in the x-direction compared to the first are given by the following relations [31]:
c t = γ ( c t β x ) ,
x = γ ( x β c t ) ,
y = y ; z = z ,
where β = v / c , γ = 1 / 1 β 2 , and c is the speed of light in vacuum.
The more general relations between the coordinates in the two frames when the second moves at constant velocity in a direction characterized by unit vector n in the first frame are given by the following [31].
c t = γ c t β γ ( n ^ · r )
r = r β γ c t n ^ + ( γ 1 ) ( n ^ · r ) n ^ .
Hereafter, we show that these can be obtained fairly easily by using a quaternion formalism. First, we describe the approach developed by Sweetser to derive the equations describing the Lorentz boost along the x-direction using left-chirality quaternions [32]. The derivation can be readily extended to describe boosts along y and z directions. Next, we show how the use of the multiplication rule for quaternions can be used to generate relations (48) and (49). Finally, we show how the formulation of a general Lorentz boost can also be described using right-chirality quaternions.

3.1. Lorentz Boost along the x, y, or z Directions Using Left-Chirality Quaternions

Recently, Sweetser has been able to formulate a derivation of a Lorentz boost along the x direction using left-chirality quaternions with real components [32]. The derivation is reproduced here for the sake of completeness and then generalized to a boost along the y and z directions; finally, it is generalized to a boost with constant velocity v in an arbitrary direction characterized by unit vector n ^ .
Sweetser first performed the following transformation of the left-chirality quaternion Q L , L x l Q L L x l ¯ , where Q L and L x l are the left-chirality quaternions given by the following:
Q L = ( c t , x , y , z ) ,
and the following is the case:
L x l = ( cosh α , sinh α , 0 , 0 ) ,
respectively, and parameter α must be determined so that the components of Q L = ( c t , x , y , z ) in the second Lorentz frame moving with a constant velocity v along the x axis of the initial Lorentz frame agree with Equations (45)–(47).
The explicit 4 × 4 matrix representation of the quaternion L x l is given by the following.
L x l cosh α sinh α 0 0 sinh α cosh α 0 0 0 0 cosh α sinh α 0 0 sinh α cosh α .
Using the multiplication rule for left-chirality quaternions given in Equation (4), we obtain the following.
Q L L x l ¯ = ( c t cosh α + x sinh α , c t sinh α + x cosh α , y cosh α z sinh α , y sinh α + z cosh α ) .
Using the multiplication rule (4) one more time and taking into account of the following hyperbolic trigonometric identities:
cosh 2 α + sinh 2 α = cosh 2 α ,
and the following as well:
sinh 2 α = 2 cosh α sinh α ,
we obtain the following.
L x l Q L L x l ¯ = ( c t cosh 2 α , x cosh 2 α , y 2 z cosh α sinh α , z + 2 y sinh α cosh α ) .
If we select α , the following is the case:
cosh 2 α = γ .
since the relations between γ and β and the relation between cosh 2 α and sinh 2 α are as follows.
γ = 1 1 β 2
cosh 2 2 α sinh 2 2 α = 1
We can derive a hyperbolic expression for γ β .
2 cosh α sinh α = sinh 2 α = γ β .
Equation (56) can then be rewritten.
L x l Q L L x l ¯ = ( c t γ , x γ , y β γ z , z + γ β y ) .
In order for this quantity to be equal ( c t , x , y , z ) , where c t , x , y , and z are given by Equations (45)–(47), the following left-chirality quaternion must be added to Equation (61).
( γ β x , γ β c t , γ β z , γ β y ) .
This last quaternion can be generated by first calculating the following quaternion L x l L x l Q L . We find the following.
L x l L x l = ( 1 , sinh 2 α , 0 , 0 ) .
Hence, the following is the case.
L x l L x l Q L = ( c t γ β x , x + γ β c t , y γ β z , z + γ β y ) .
The conjugate of the latter quaternion is equal to the following.
L x l L x l Q L ¯ = ( c t γ β x , x γ β c t , y + γ β z , z γ β y ) .
Substracting Q L ¯ = ( c t , x , y , z ) from the latter, we obtain the following.
L x l L x l Q L ¯ Q L ¯ = ( γ β x , γ β c t , γ β z , γ β y ) .
Adding the latter to Equation (61), we obtain the desired result.
Q L = L x Q L = ( c t , x , y , z ) = L x l Q L L x l ¯ + L x l L x l Q L ¯ Q L ¯ .
This last expression can be rewritten using the property for the conjugate of the product of two quaternions.
Q L = L x Q L = L x l Q L L x l ¯ + Q L ¯ L x l L x l ¯ 1 .
Following a similar derivation, it can be shown that a Lorentz boost along the y-axis can be generated by using the following transformation:
Q L = L y Q L = L y l Q L L y l ¯ + Q L ¯ L y l L y l ¯ 1 .
where the following is the case:
L y l = ( cosh α , 0 , sinh α , 0 ) ,
and a Lorentz boost along the z-direction can be generated by using the following transformation:
Q L = L z Q L = L z l Q L L z l ¯ + Q L ¯ L z l L z l ¯ 1 .
where
L z l = ( cosh α , 0 , 0 , sinh α ) .

3.2. Lorentz Boost along an Arbitrary Direction n ^ Using Left-Chirality Quaternions

In this section, we show that the results of the previous section can be extended to describe a Lorentz boost with constant velocity v along an arbitrary direction n ^ using the following relation link between the left-chirality quaternions Q L = ( c t , x , y , z ) and Q L = ( c t , x , y , z ) associated with the spatio-temporal coordinates of the two Lorentz frames:
Q L = L n ^ l Q L = L n ^ l Q L L n ^ l ¯ + Q L ¯ L n ^ l L n ^ l ¯ 1 ,
where quaternion L n ^ l reduces to L x l , L y l , and L z l for a Lorentz boost along the x-, y-, and z-axis, respectively, as found in the previous section.
To prove that the quaternionic transformation in Equation (73) leads to general Equations (48) and (49), we first rewrite the product of two left-chirality quaternions as follows.
A left-chirality quaternion can be represented as follows:
Q 1 , L = ( a 1 , b 1 , c 1 , d 1 ) = ( s , v ) ,
where s = a 1 and v = ( b 1 , c 1 , d 1 ) are the real and vector parts of the quaternion Q 1 , L , respectively.
Taking the product of this quaternion with another quaternion as follows:
Q 2 , L = ( a 2 , b 2 , c 2 , d 2 ) = ( t , w ) ,
and using the multiplication rule of quaternion Equation (4), product Q 1 , L Q 2 , L can be written as follows:
( s , v ) ( t , w ) = ( s t v · w ) + ( s w + t v + v × w ) ,
where v · w , and v × w is the dot and vector product of the two three-dimensional vectors v and w , respectively.
Using this last equation, we first calculate product P = L n ^ l Q L L n ^ l ¯ :
P = ( cosh α , sinh α n ^ ) ( c t , r ) ( cosh α , sinh α n ^ ) = ( cosh α , sinh α n ^ ) ( c t cosh α + sinh α r · n ^ , c t sinh α n ^ + cosh α r sinh α r × n ^ ) = ( cosh α ( c t cosh α + sinh α r · n ^ ) ( sinh α n ^ ) · ( c t sinh α n ^ + cosh α r sinh α r × n ^ ) , cosh α ( c t sinh α n ^ + cosh α r sinh α r × n ^ ) + sinh α n ^ ( c t cosh α + sinh α r · n ^ ) + ( sinh α n ^ ) × ( c t sinh α n ^ + cosh α r sinh α r × n ^ ) ) ,
where r = ( x , y , z ) .
This last expression can be simplified using the following identities.
n ^ × n ^ = 0 , n ^ · r × n ^ = 0 , n ^ · n ^ = 1 , n ^ × r × n ^ = ( n ^ · n ^ ) r ( n ^ · r ) n ^ .
Then, using hyperbolic function identities, Equation (77) becomes the following.
P = c t ( cosh 2 α + sinh 2 α ) , cosh 2 α r + 2 sinh α cosh α ( n ^ × r ) + sinh 2 α n ^ ( r · n ^ ) sinh 2 α ( ( n ^ · n ^ ) r ( n ^ · r ) n ^ ) = ( c t cosh 2 α , cosh 2 α r + sinh 2 α ( n ^ × r ) + 2 sinh 2 α ( n ^ · r ) n ^ sinh 2 α r ) .
On the other hand, we have the following.
Q L ¯ ( L n ^ l L n ^ l 1 ) = ( sinh 2 α ( n ^ · r ) , c t sinh 2 α n ^ + sinh 2 α ( r × n ^ ) )
Combining Equations (79) and (80), we obtain the following.
Q L = L n ^ l Q L = ( c t cosh 2 α sinh 2 α ( n ^ · r ) , cosh 2 α r + sinh 2 α ( n ^ × r ) + 2 sinh 2 α ( n ^ · r ) n ^ sinh 2 α r c t sinh 2 α n ^ + sinh 2 α ( r × n ^ ) ) = ( c t cosh 2 α sinh 2 α ( n ^ · r ) , r c t sinh 2 α n ^ + 2 sinh 2 α ( n ^ · r ) n ^ ) .
Taking into account the following relations:
cosh 2 α = γ , sinh 2 α = β γ , 2 sinh 2 α = ( sinh 2 α + cosh 2 α ) ( cosh 2 α sinh 2 α ) = cosh 2 α 1 = γ 1 ,
we arrive at the desired quaternionic expression for the Lorentz boost with velocity v in an arbitrary direction n .
Q L = L n ^ l Q L = ( γ c t β γ ( n ^ · r ) , r β γ c t n ^ + ( γ 1 ) ( n ^ · r ) n ^ ) .
Equating the real and vectorial parts of the quaternions on both sides of this last equation results in the following:
c t = γ c t β γ ( n ^ · r ) , r = r β γ c t n ^ + ( γ 1 ) ( n ^ · r ) n ^ ,
which are the well-known relations connecting two Lorentz frames where the second one moves with a constant velocity in direction n ^ in the first frame [31].

3.3. Three Lorentz Boosts along an Arbitrary Direction n ^ Using Right-Chirality Quaternions

A Lorentz boost in an arbitrary direction n ^ can also be described in terms of right-chirality quaternions. The multiplication rule (11) for right-chirality quaternions is first rewritten as follows:
Q 1 , R Q 2 , R = ( s , v ) ( t , w ) = ( s t v · w ) + ( s w + t v v × w ) ,
which differs from Equation (76) by the sign in front of the vector product, v × w .
A proof similar to the one provided in Section 3.2 can be used to show that the general formula describing the connection between two Lorentz frame for which its 4-dimensional variables are given by the right-chirality quaternions Q R = ( c t , x , y , z ) and Q R = ( c t , x , y , z ) is provided by the following:
Q R = L n ^ l Q L = L n ^ r ¯ Q R L n ^ r + L n ^ r L n ^ r ¯ 1 Q R ¯ ,
where the 4 × 4 expression of the right-chirality quaternion Q R is given by the following:
Q R < = > c t x y z x c t z y y z c t x z y x c t ,
and where the 4 × 4 representation of the quaternion L n ^ r is given by the following.
L n ^ r = cosh α sinh α n x sinh α n y sinh α n z sinh α n x cosh α sinh α n z sinh α n y sinh α n y sinh α n z cosh α sinh α n x sinh α n z sinh α n y sinh α n x cosh α .

3.4. Lorentz Invariant

In this section, starting with Equation (83), we show that the real part of Q L 2 is a Lorentz invariant, i.e., it is identical in two frames related by a Lorentz boost.
We first rewrite Q L = ( s , v ) where the real part is given by the following:
s = c t β γ ( n ^ · r ) ,
and the vectorial part is given by the following.
r = r c t sinh 2 α n ^ + 2 sinh 2 α ( n ^ · r ) n ^ .
Now, using Equation (76), the real part of Q L Q L is given by the following.
r e a l ( Q L Q L ) = s 2 v · v .
It can easily be shown that the following is the case.
s 2 = c 2 t 2 cosh 2 2 α + sinh 2 2 α + ( n ^ · r ) 2 2 c t cosh 2 α sinh 2 α ( n ^ · r ) .
Moreover, we have the following.
v · v = r 2 2 c t ( sinh 2 α + 2 sinh 2 α sinh 2 α ) ( n ^ · r ) + 4 ( sinh 2 α + sinh 4 α ) ( n ^ · r ) 2 .
Using the last two equations and some hyperbolic trigonometric identities, it is then easy to show that the following is the case.
r e a l ( Q L Q L ) = s 2 v · v = c 2 t 2 r 2 = r e a l ( Q L Q L ) ,
In other words, it is indeed invariant.

4. A Generalized Split-Quaternion Approach to Tunneling Problems

The transmission matrix technique is a well-known approach for tackling tunneling problems by using an arbitrary conduction band energy profile in quantum mechanics. Hereafter, the transmission matrix approach is reformulated using split-quaternion formalism. The transmission and reflection probabilities through an arbitrary 1D conduction band energy profile can then be easily derived by multiplying real 4 × 4 matrices. The approach is illustrated by revisiting the tunneling problems by using a 1D delta-scatterer and by using a resonant tunneling structure composed of two and three delta-scatterers separated by distance L. We show that the tunneling problems can be treated using any of the six split-quaternions introduced in Section 2.3 and Table 2.

4.1. Split Quaternion Formulation of Transmission Matrix

Consider a device with an arbitrary 1D conduction band energy profile sandwiched between two contacts, as shown in Figure 1. Hereafter, we assume that the two contacts are identical and that there is no bias across the device (the analysis described hereafter can easily be extended to the case of a finite bias). Assuming a constant effective mass m throughout, the general solution of the 1D Schrödinger equation in the left contact is given as the superposition of left and right propagating plane waves:
ϕ I ( x ) = A e i k x + B e i k x ,
where k = 1 2 m E , and E is the kinetic energy of the electron in the contact ( E c is assumed to be zero in both contacts).
Similarly, in the right contact, we have the following.
ϕ I I ( x ) = C e i k x + D e i k x .
By definition, the transmission matrix T associated with the device between the two contacts is the 2 × 2 complex matrix such that the following is the case.
C D = T A B .
Starting with the fundamental properties of the 1D Schrödinger equation, i.e., current density conservation and time reversal symmetry, it can be shown that the matrix T is of the general form:
T = m n n m ,
where m and n are complex functions of energy and d e t ( T ) = | m | 2 | n | 2 = 1 .
The matrix elements m and n are related to the transmission and reflection amplitudes of an electron incident from the contacts. Considering an electron incident from the left with A = 1 and B = r , the reflection amplitude, and D = 0 and C = t , which are the transmission amplitudes, using Equations (95) and (96), we obtain the following.
t 0 = m n n m 1 r .
It can be easily shown that the following is the case:
t = 1 / m , r = n / m ,
and the following is the case:
| t | 2 + | r | 2 = 1 ,
where | t | 2 and | r | 2 are the transmission and reflection probabilities of the device, respectively. Equation (101) expresses the conservation of the current density in the scattering problem.
Starting with Equation (97), we obtain the following.
C = m A m A + n B n B ,
C = m A + m A + n B + n B ,
D = n A + n A + m B + m B ,
D = n A n A + m B m B .
Replacing the complex elements of the transmission matrix (Equation (98)) by their 2 × 2 matrix representations [11], we obtain the following 4 × 4 real matrix representation of the transmission matrix.
T = m n n m < = > m m n n n m n n n n m m n n m m .
The latter is of the generic form of the 4 × 4 matrix representation of a split-quaternion given by Equation (12) if we select w = Rem, x = Imm, y = Ren, and z = Imn.
If we replace the column vector associated with the amplitudes of the incident and reflected amplitudes in the left contact by the following split-quaternion:
A B < = > A A B B A A B B B B A A B B A A ,
and multiply it on the left by the split-quaternion (106) associated with the transmission matrix; we obtain a split-quaternion P for which its 4 × 4 matrix has the following elements in its top row.
P ( 1 , 1 ) = m A m A + n B n B ,
P ( 1 , 2 ) = m A + m A + n B + n B ,
P ( 1 , 3 ) = n A + n A + m B + m B ,
P ( 1 , 4 ) = n A + n A m B + m B .
Comparing these results with Equations (102)–(105), we obtain the following:
P ( 1 , 1 ) = C , P ( 1 , 2 ) = C , P ( 1 , 3 ) = D , I m ( 1 , 4 ) = D ,
which is the split-quaternion associated with the outgoing and incoming wave amplitudes in the right contact.
C D < = > C C D D C C D D D D C C D D C C .
The rules outlined above establish a one-to-one correspondence between the 2 × 2 complex transmission formalism and the 4 × 4 real split-quaternion approach to the general tunneling problem in 1D.
If a general conduction band energy profile is divided into subsections, the overall split-quaternion associated with the scattering problem describing the tunneling of an electron incident from the left contact can be obtained by multiplying, from right to left, the individual split-quaternions associated with the subsections of the device moving from the left to right through the device. The overall split-quaternion associated with the device will then be of a general form in Equation (106) from which the transmission probability can be obtained as the sum of square magnitudes of the first two components in the top row of the overall split-quaternion.
Starting with Equation (106) in its quaternion form and using Equation (99) results in the following:
| t | 4 = ( | m | 2 | n | 2 ) 2 ( 1 | r | 2 ) 2 = ( 1 | r | 2 ) 2 ,
which results in relation | t | 2 + | r | 2 = 1 , i.e., the conservation of the probability current density. Equation (114) also has solution | r | 2 | t | 2 = 1 , which is not physical.
Next, we illustrate the split-quaternion formalism by revisiting several tunneling problems, including transmission through a one-dimensional delta-scatterer and resonant tunneling structures composed of two and three delta-scatterers separated by free propagating regions of length L.

4.2. Split-Quaternion Associated to a 1D Delta Scatterer

The 1D Schrödinger equation for an electron with effective mass m in the presence of a 1D delta-scatterer is given by the following:
2 2 m d 2 ψ d x 2 + V ( x ) ψ ( x ) = E ψ ( x )
with potential energy profile V ( x ) = Γ δ ( x ) and Γ being the strength of the delta-scatterer.
To calculate the transmission matrix associated with the delta-scatterer, we seek solutions to the 1D Schrödinger equation to the left and right of the scatterer of the following form.
ψ L = A e i k x + B e i k x ( x < 0 ) ψ R = C e i k x + D e i k x ( x > 0 ) .
Imposing continuity of the wavefunction at x = 0 results in the following.
C + D = A + B .
Integrating Equation (115) on both sides from ϵ to + ϵ and letting ϵ tends to zero results in the following.
d ψ d x ( x = 0 + ) d ψ d x ( x = 0 ) = 2 m 2 ψ ( 0 ) .
Making use of Equation (116), this last equation results in another relation between coefficients A, B, C, and D.
C D = A B 2 i m Γ 2 ( A + B ) .
Equations (117) and (119) can be rewritten as follows.
C D = 1 i m Γ k 2 i m Γ k 2 i m Γ k 2 1 + i m Γ k 2 A B .
Hence, the transmission matrix T through a one-dimensional delta scatterer is of the generic form (98) with the following being the case:
m = 1 i m Γ k 2 ,
and the following.
n = i m Γ k 2 .
Using the expressions for m and n, we easily check the following.
d e t ( T ) = | m | 2 | n | 2 = ( 1 + m Γ k 2 2 ) m Γ k 2 2 = 1 .
Using Equations (106) and (120), the 4 × 4 matrix associated with the split-quaternion of a delta-scatterer is given by the following.
T δ = 1 m Γ k 2 0 m Γ k 2 m Γ k 2 1 m Γ k 2 0 0 m Γ k 2 1 m Γ k 2 m Γ k 2 0 m Γ k 2 1 .
In order to calculate the transmission probability through an array of identical delta-scatterers with the same separation L between them, we determine next the split-quaternion associated to a free propagation region. The overall split-quaternion associated with the array of delta-scatterers is then obtained by cascading (i.e., multiplying), from right to left, the 4 × 4 split-quaternion matrices associated with the individual scatterers and free propagation regions between them.

4.3. Split-Quaternion Associated to a Free Propagation Region of Length L

Consider a free propagation region of length L (see Figure 2) between two contacts and assume that the conduction band energy profile is constant and the same throughout. The following relations are found between the wave amplitudes of the left- and right-propagating waves in the two contacts.
C = A e i k L ,
D = B e i k L .
This leads to the following expression for the transmission matrix associated with a free propagation of length L:
T p r o p = e i k L 0 0 e i k L ,
which is of generic form (98) with m = e i k L and n = 0 . Using Equation (106), the split-quaternion associated with a free propagation region of length L is, therefore, given by the following.
T p r o p < = > cos k L sin k L 0 0 sin k L cos k L 0 0 0 0 cos k L sin k L 0 0 sin k L cos k L .

4.4. Split-Quaternion Associated to a Finite Square Barrier

To derive the split-quaternion associated with a square barrier, we need to establish the connection between the transmission matrix and the transfer matrix, which is another well-known technique for solving tunneling problems in quantum mechanics [30].
For a region of length d sandwiched between two contacts, the transfer matrix is defined as the matrix which relates solution ϕ of the 1D Schrödinger equation and its derivatives on both sides of the device as follows.
ϕ ˙ ( d ) ϕ ( d ) = W 11 W 12 W 21 W 22 ϕ ˙ ( 0 ) ϕ ( 0 ) .
If ϕ 1 ( x ) and ϕ 2 ( x ) are two linearly independent solutions of the Schrödinger equation satisfying the boundary conditions at the left contact, then the following is the case.
ϕ ˙ 1 ( 0 ) ϕ ˙ 2 ( 0 ) ϕ 1 ( 0 ) ϕ 2 ( 0 ) = 1 0 0 1 .
The transfer matrix can be written as follows [30].
W = ϕ ˙ 1 ( d ) ϕ ˙ 2 ( d ) ϕ 1 ( d ) ϕ 2 ( d ) .
For a square barrier with height V 0 and width d, the explicit form of the linearly independent solutions ϕ 1 ( x ) and ϕ 2 ( x ) are well-known, leading to the following transfer matrix [30]:
W = cos ( k d ) k sin ( k d ) 1 k sin ( k d ) cos ( k d ) ,
when E > V 0 . In Equation (132), the following is the case.
k = 1 2 m ( E V 0 ) .
When V 0 > E , the transfer matrix is given by the following [30]:
W = cosh ( κ d ) κ sinh ( κ d ) 1 κ sinh ( κ d ) cosh ( κ d ) ,
with the following being the case.
κ = 1 2 m ( V 0 E ) .
For the 1D scattering problem, the solution to the Schrödinger equation can be written as follows.
ϕ L ( x ) = A + e i k x + A e i k x ( x < 0 ) ϕ ( x ) = c 1 ϕ 1 ( x ) + c 2 ϕ 2 ( x ) ( 0 < x < d ) ϕ R ( x ) = B + e i k x + B e i k x ( x > d ) .
At x = 0, we obtain the following:
c 1 ϕ ˙ 1 ( 0 ) + c 2 ϕ ˙ 2 ( 0 ) = i k ( A + A ) ,
and the following is also obtained.
c 1 ϕ 1 ( 0 ) + c 2 ϕ 2 ( 0 ) = A + + A .
Similarly, at x = d, the following relations must be satisfied:
c 1 ϕ ˙ 1 ( d ) + c 2 ϕ ˙ 2 ( d ) = i k ( B + B ) ,
and the following must be satisfied as well.
c 1 ϕ 1 ( d ) + c 2 ϕ 2 ( d ) = B + + B .
Starting with Equations (137) and (138), we obtain the following relations between cofficients ( A + , A ) and ( B + , B ) .
i k i k 1 1 A + A = W 1 i k i k 1 1 B + B .
Using Equation (97), the connection between the transmission and transfer matrix of the device between the two contacts is then given by the following:
T d e v i c e = i k i k 1 1 1 W i k i k 1 1 = 1 2 i k 1 i k 1 i k W i k i k 1 1 ,
where k is the wavevector in the contact regions. Carrying out matrix multiplication on the right-hand side of Equation (142) leads to the following.
T d e v i c e = 1 2 ( W 11 + W 22 ) i 2 k ( W 12 k 2 W 21 ) 1 2 ( W 11 W 22 ) i 2 k ( W 12 + k 2 W 21 ) 1 2 ( W 11 W 22 ) + i 2 k ( W 12 + k 2 W 21 ) 1 2 ( W 11 + W 22 ) + i 2 k ( W 12 k 2 W 21 ) .
For the particular case of a square barrier, Equations (132) and (134) show that W 11 = W 22 . Therefore, the transmission matrix has the following simplified form:
T b a r r i e r = W 11 i 2 k ( W 12 k 2 W 21 ) i 2 k ( W 12 + k 2 W 21 ) + i 2 k ( W 12 + k 2 W 21 ) W 11 + i 2 k ( W 12 k 2 W 21 ) ,
which, in its 4 × 4 split-quaternion form, is given by the following.
T b a r r i e r = W 11 1 2 k ( W 12 k 2 W 21 ) 0 1 2 k ( W 12 + k 2 W 21 ) 1 2 k ( W 12 k 2 W 21 ) W 11 1 2 k ( W 12 + k 2 W 21 ) 0 0 1 2 k ( W 12 + k 2 W 21 ) W 11 1 2 k ( W 12 k 2 W 21 ) 1 2 k ( W 12 + k 2 W 21 ) 0 1 2 k ( W 12 k 2 W 21 ) W 11 .
The tunneling probability through an RTD composed of two identical square barriers separated by a free propagation region can then be obtained from the split-quaternion formalism by first calculating the overall split-quaternion associated with this structure:
T R T D = T b a r r i e r T p r o p T b a r r i e r ,
where T p r o p and T b a r r i e r are given by Equations (124) and (145), respectively.
Once the overall 4 × 4 matrix T R T D is found, the tunneling probability through RTD can then be found from Equation (100).
Hereafter, we illustrate how an RTD consisting of two 1D delta scatterers separated by distance L can be replaced by a realistic resonant tunneling device composed of two square barriers of height V e f f and width d e f f with a free propagation of length L d e f f between them by increasing the effective barrier height V e f f and reducing its width d e f f while keeping product V e f f d e f f constant and equal to Γ , as illustrated in Figure 3.
Figure 4 shows the energy dependence of the tunneling probability through an RTD composed of two delta scatterers of strength Γ = 15 eV separated by distance L = 100. The effective mass was assumed to be m = 0.067 m 0 . Moreover, the transmission probabilities through an RTD composed of two square barriers of height V e f f and separated by a distance L d e f f for four different V e f f are shown as red dashed lines: (a) 0.3, (b) 1, (c) 5, and (d) 10 eV. In each plot, the effective width of the barrier d e f f was adjusted to keep the product V e f f d e f f = Γ = 15 eV. Figure 4 clearly shows that the transmission probability through the two delta scatterers can be well approximated by tunneling through two square barriers of increasing height and decreasing width.
This analysis can be extended to the case of multiple barriers. Figure 5 shows a comparison of the transmission probability through an array of three delta scatterers of strength Γ = 15 eV and with a separation of 100 between adjacent scatterers. The frames (a), (b), (c), and (d) in Figure 5 correspond to a value of V e f f equal to 0.3, 1, 5, and 10 eV, respectively.

4.5. Solving Tunneling Problems with Alternative Split-Quaternions

In Section 2.3, we found that there are five division algebras that are generalizations of the split-quaternion first introduced by Cockle [13]. In fact, the formulation of tunneling problems can be performed using any of the split-quaternion algebras listed in Table 3 since Equation (37) establishes isomorphism between the six generalized split-quaternions algebras listed in Table 2.
Next, we illustrate how to derive the left- and right-chirality nature of the six split-quaternions listed in Table 2.
In Section 4.1, we formulated tunneling problems using split-quaternion Q + + k , i , j ( a , b , c , d ) for which its 4 × 4 matrix representation is given at the bottom of Table 2. In that case, the split-quaternion associated with the outgoing wave amplitudes (C,D) in Equation (113) was obtained as a product of the split-quaternion associated with the transmission matrix in Equation (106) and the split-quaternion associated with the incoming wave amplitudes (A,B) in Equation (107). Because the 4 × 4 matrix associated with the split-quaternion of the transmission matrix is applied to the left of the split-quaternion associated with the incoming wave amplitudes, we refer to the 4 × 4 matrix at the bottom of Table 2 as the left-chirality split-quaternion. Next, we show that tunneling problems can also be solved using any of the other five split-quaternion algebras listed in Table 2.
Starting with the split-quaternion labelled as Q + + k , i , j ( a , b , c , d ) for which its 4 × 4 matrix representation is given by the following:
Q + + k , i , j ( a , b , c , d ) = a b c d b a d c c d a b d c b a ,
we introduce the following split-quaternion associated with the incoming wave amplitudes (A, B):
A B < = > A A B B A A B B B B A A B B A A ,
with components (a, b, c, d) = ( A , A , B , B ) , and the split-quaternion associated with the transmission matrix.
T < = > m m n n m m n n n n m m n n m m .
Multiplying the split-quaternion in Equation (148) on the right by the split-quaternion in Equation (149), the resulting split-quaternion associated with outgoing amplitudes (C, D) is found to be of the following form:
C B < = > C C D D C C D D D D C C D D C C ,
where ReC, ImC, ReD, and ImD have the same expressions as in Equations (102)–(105), respectively.
Because this result was obtained by applying the split-quaternion associated with the transmission matrix T to the right of the split-quaternion matrix associated with the incoming amplitudes (A, B), we refer to the split-quaternion of the form Q + + k , i , j ( a , b , c , d ) as a right-chirality split-quaternion.
Following similar derivations, it can be shown that the four other split-quaternions division algebras listed in Table 2 have chirality properties listed in Table 3.
Similarly to quaternions, split-quaternions are elements of a 4-dimensional associative algebra [13]. They have found important applications in many fields such as spatial geometry, physics, quantum mechanics, signal processing, and cryptography, among others [16,17,18,19,20,21,22,23,24,25].

5. Conclusions and Suggestions for Future Work

This paper outlines a formalism to study various physical phenomena using proper non-commutative division algebras based on a generalization of the concept of quaternions first introduced by Hamilton. It was shown that there are a total of 64 possible multiplication rules, which can be defined starting with the generalized imaginary units first introduced by Hamilton. Of these sixty-four choices, only eight lead to non-commutative division algebras, two of which are associated with the left- and right-chirality quaternions while the other six are generalizations of the split-quaternion concept first introduced by Cockle. For the latter, their 4 × 4 real matrix representations consist of the six four-dimensional representations of the order eight dihedral group. Next, we showed that the 4 × 4 of both the left- and right-chirality versions of the generalized split-quaternions are algebraically isomorphic and can be related to each other by 4 × 4 permutation matrices of the C 2 × C 2 group.
It was also shown that the left- and right-chirality quaternions can be used to describe Lorentz boosts with constant velocity in an arbitrary direction. Next, the generalized split-quaternion algebras were used to solve the problem of quantum-mechanical tunneling by using an arbitrary one-dimensional conduction band energy profile. This demonstrates that six different spinors ( 4 × 4 matrices) can be used to represent the amplitudes of the left and right propagating waves in a device.
The generalized quaternion formalism introduced in this paper could be used to reformulate a derivation of the Maxwell equations [33,34,35,36,37,38] and a quaternionic form of the Dirac equations [39,40,41,42,43]. These topics will be the subject of future investigations.

Author Contributions

Conceptualization, M.C. and D.M.; methodology, H.-Y.L., M.C. and D.M.; software, H.-Y.L.; formal analysis, H.-Y.L., M.C., B.N.V. and D.M.; writing—original draft preparation, M.C. and H.-Y.L.; visualization, H.-Y.L. and M.C.; supervision, M.C.; project administration, M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Computer codes used to generate the simulation results are available from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hamilton, W.R. Computer codes used to generate the simulation results are available from the corresponding author. Philos. Mag. 1844, 25, 489. [Google Scholar]
  2. Morris, D. Quaternions; Abane and Right: Port Mulgrave, UK, 2015. [Google Scholar]
  3. Ward, J.P. Quaternions and Cayley Numbers: Algebra and Applications; Kluwer Academic, Springer-Science+Business Media, B.V.: Dordrecht, The Netherlands, 1997. [Google Scholar]
  4. Jafari, M. Quaternions Algebra and Its Applications: An Overview. Int. J. Theor. Appl. Math. 2016, 2, 79. [Google Scholar]
  5. Adler, S.L. Quaternionic Quantum Mechanics and Quantum Fields; Oxford University Press: New York, NY, USA, 1995. [Google Scholar]
  6. Arbab, A.I. Derivation of Dirac, Klein-Gordon, Schrödinger, Diffusion and Quantum Heat Transport Equations from a Universal Quantum Wave Equation. Europhys. Lett. 2010, 92, 40001. [Google Scholar] [CrossRef]
  7. Girard, P.R. The Quaternion Group and Modern Physics. Eur. J. Phys. 1984, 5, 25. [Google Scholar] [CrossRef]
  8. Cartan, E. The Theory of Spinors; MIT Press: Cambridge, MA, USA, 1966. [Google Scholar]
  9. Bandyopadhyay, S.; Cahay, M. Introduction to Spintronics, 2nd ed.; CRC Book Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  10. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information; Cambridge University Press: New York, NY, USA, 2000. [Google Scholar]
  11. Cahay, M.; Purdy, G.; Morris, D. On the Quaternion Representation of the Pauli Spinor of An Electron. Phys. Scr. 2019, 94, 085205. [Google Scholar] [CrossRef]
  12. Cahay, M.; Morris, D. On the Quaternionic Form of the Pauli- Schrödinger Equation. Phys. Scr. 2020, 95, 015204. [Google Scholar] [CrossRef]
  13. Cockle, J. On Systems of Algebra Involving More Than One Imaginary; And On Equations of The Fifth Degree. Philos. Mag. 1849, 35, 434. [Google Scholar] [CrossRef] [Green Version]
  14. Kula, L.; Yayli, Y. Split Quaternions and Rotations in Semi Euclidean Space. J. Korean Math. Soc. 2007, 44, 1313. [Google Scholar] [CrossRef] [Green Version]
  15. Machen, C. The Exponential of a Quaternion Matrix. Rose-Hulman Undergr. Math. J. 2011, 12, 29. [Google Scholar]
  16. Özdemir, M. The Roots of a Split Quaternion. Appl. Math. Lett. 2009, 22, 258. [Google Scholar] [CrossRef] [Green Version]
  17. Özdemir, M.; Ergin, A.A. Rotations with Unit Timelike Quaternions in Minkowski 3-space. J. Geom. Phys. 2006, 56, 322. [Google Scholar] [CrossRef]
  18. Ujang, B.C.; Took, C.C.; Mandic, D.P. Split Quaternion Nonlinear Adaptive Filtering. Neural Netw. 2010, 23, 426. [Google Scholar] [CrossRef] [PubMed]
  19. Ata, E.; Yayli, Y. Split Quaternions and Semi-Euclidean Projective Spaces. Chaos Solitons Fractals 2009, 41, 1910. [Google Scholar] [CrossRef]
  20. Jiang, T.; Jiang, Z.; Zhang, Z. Algebraic Techniques for Diagonalization of a Split Quaternion Matrix In Split Quaternion Mechanics. J. Math. Phys. 2015, 56, 083509. [Google Scholar] [CrossRef]
  21. Brody, D.C.; Graefe, E.M. On Complexified Mechanics and Coquaternions. J. Phys. A Math. Theor. 2011, 44, 072001. [Google Scholar] [CrossRef] [Green Version]
  22. Brody, D.C.; Graefe, E.M. Coquaternions Quantum Dynamics for Two-Level Systems. Acta Polytech. 2011, 51, 14. [Google Scholar] [CrossRef]
  23. Bender, C.M. Coquaternions Quantum Dynamics for Two-Level Systems. Rep. Prog. Phys. 2007, 70, 947. [Google Scholar] [CrossRef] [Green Version]
  24. Gogberashvili, M. Split Quaternions and Particles in (2+1)-Space. Eur. Phys. J. C 2014, 74, 3200. [Google Scholar] [CrossRef] [Green Version]
  25. Yormaz, C.; Elamz, S.N.; Simsek, S. Hamiltonian Mechanical System with Split Quaternions. Univers. J. Appl. Math. 2018, 6, 17. [Google Scholar] [CrossRef]
  26. Antonuccio, F. Split-Quaternions and the Dirac Equation. Adv. Appl. Clifford Algebra 2015, 25, 13. [Google Scholar] [CrossRef] [Green Version]
  27. Thakur, K.; Tripathi, B.P. A Variant of NTRU with Split Quaternions Algebra. Palest. J. Math. 2017, 6, 598610. [Google Scholar]
  28. Babaarslan, M.; Yayli, Y. Split Quaternions and Spacelike Constant Slope Surfaces in Minkowski 3-space. Int. J. Geom. 2013, 2, 23. [Google Scholar]
  29. Pop, M.-I.; Cretu, N. Intrinsic Transfer Matrix Method and Split Quaternion Formalism for Multilayer Media. Wave Motion 2016, 65, 105. [Google Scholar] [CrossRef]
  30. Cahay, M.; Bandyopadhyay, S. Problems in Quantum Mechanics, for Material Scientists, Applied Physicists, and Device Engineers; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2017. [Google Scholar]
  31. Jackson, J.D. Classical Electrodynamics, 2nd ed.; Wiley: New York, NY, USA, 1975. [Google Scholar]
  32. Is This a Quaternion Lorentz Boost? Available online: https://physics.stackexchange.com/questions/28797/is-this-a-quaternion-lorentz-boost (accessed on 11 October 2021).
  33. Maxwell, J.C. A Dynamic Theory of The Electromagnetic Field, Philosophical Transactions of the Royal Society of London. R. Soc. Trans. 1865, 155, 459. [Google Scholar]
  34. Maxwell, J.C. A Treatise on Electricity and Magnetism; Clarendon Press: Oxford, UK, 1973. [Google Scholar]
  35. Morris, D. The Physics of Empty Space; Abana and Right: Port Mulgrave, UK, 2015. [Google Scholar]
  36. Kravchenko, V. Quaternionic equation for electromagnetic fields in homogeneous media. arXiv 2002, arXiv:math-phys/0202010. [Google Scholar]
  37. Jack, P.M. Physical space as a quaternion structure, I: Maxwell equations. A Brief Note. arXiv 2003, arXiv:math-phys/0307038v1. [Google Scholar]
  38. Edmonds, J.D. Quaternionic Quantum Theory: New Physics or Number Misticism? Am. J. Phys. 1978, 46, 430. [Google Scholar] [CrossRef]
  39. Davies, A.J. Quaternionic Dirac Equation. Phys. Rev. D 1990, 41, 2628. [Google Scholar] [CrossRef] [PubMed]
  40. Schwartz, C. Relativistic Quaternionic Wave Equation. J. Math. Phys. 2006, 47, 122301. [Google Scholar] [CrossRef] [Green Version]
  41. Sobhani, H.; Hassanabadi, H. Quaternionic Dirac Equation. Can. J. Phys. 2016, 94, 262. [Google Scholar] [CrossRef]
  42. Schuricht, D.; Greiter, M. Contemplations on Dirac’s Equation in Quaternionic Coordinates. Eur. J. Phys. 2004, 25, 755. [Google Scholar] [CrossRef] [Green Version]
  43. Peres, N.M.R. Scattering in One-Dimensional Heterostructures Described by the Dirac Equation. J. Phys. Condens. Matter 2009, 21, 095501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Illustration of a device with an arbitrary conduction band energy profile between two contacts. The transmission matrix of the device relates the amplitudes of the left- and right-propagating plane waves (A,B) and (C,D) in both contacts.
Figure 1. Illustration of a device with an arbitrary conduction band energy profile between two contacts. The transmission matrix of the device relates the amplitudes of the left- and right-propagating plane waves (A,B) and (C,D) in both contacts.
Symmetry 14 00599 g001
Figure 2. Illustration of the left- and right-propagating waves in the two contacts sandwiching a free propagation region of length L. The conduction band energy profile is assumed to be constant throughout.
Figure 2. Illustration of the left- and right-propagating waves in the two contacts sandwiching a free propagation region of length L. The conduction band energy profile is assumed to be constant throughout.
Symmetry 14 00599 g002
Figure 3. A resonant tunneling device consisting of two delta scatterers separated by a distance L can be replaced by a realistic resonant tunneling device composed of two square barriers of height V e f f and width d e f f with a free propagation of length L d e f f between them, by increasing the effective barrier height V e f f and reducing its width d e f f while keeping product V e f f d e f f constant and equal to strength Γ of the delta scatterers.
Figure 3. A resonant tunneling device consisting of two delta scatterers separated by a distance L can be replaced by a realistic resonant tunneling device composed of two square barriers of height V e f f and width d e f f with a free propagation of length L d e f f between them, by increasing the effective barrier height V e f f and reducing its width d e f f while keeping product V e f f d e f f constant and equal to strength Γ of the delta scatterers.
Symmetry 14 00599 g003
Figure 4. Tunneling probability through an RTD composed of two delta scatterers of strength Γ = 15 eV separated by distance L = 100. The effective mass was assumed to be set equal to m = 0.067 m 0 . Moreover, shown as red dashed lines are the the transmission probabilities through a RTD composed of two square barriers of height V e f f and separated by distance L d e f f for V e f f equal to (a) 0.3, (b) 1, (c) 5, and (d) 10 eV, respectively. In each plot, the effective width of the barrier d e f f was adjusted to keep product V e f f d e f f = Γ = 15 eV.
Figure 4. Tunneling probability through an RTD composed of two delta scatterers of strength Γ = 15 eV separated by distance L = 100. The effective mass was assumed to be set equal to m = 0.067 m 0 . Moreover, shown as red dashed lines are the the transmission probabilities through a RTD composed of two square barriers of height V e f f and separated by distance L d e f f for V e f f equal to (a) 0.3, (b) 1, (c) 5, and (d) 10 eV, respectively. In each plot, the effective width of the barrier d e f f was adjusted to keep product V e f f d e f f = Γ = 15 eV.
Symmetry 14 00599 g004
Figure 5. Transmission probability through an array of three delta scatterers of strength Γ = 15 eV and with a separation of 100 between adjacent scatterers. Moreover, shown as red dashed lines are the the transmission probabilities through a RTD composed of three square barriers of height V e f f and separated by a distance L d e f f for a value of V e f f equal to (a) 0.3, (b) 1, (c) 5, and (d) 10 eV, respectively.
Figure 5. Transmission probability through an array of three delta scatterers of strength Γ = 15 eV and with a separation of 100 between adjacent scatterers. Moreover, shown as red dashed lines are the the transmission probabilities through a RTD composed of three square barriers of height V e f f and separated by a distance L d e f f for a value of V e f f equal to (a) 0.3, (b) 1, (c) 5, and (d) 10 eV, respectively.
Symmetry 14 00599 g005
Table 1. Multiplication rules of generalized quaternion units.
Table 1. Multiplication rules of generalized quaternion units.
MultiplicationPossible Products
i 2 ± 1
j 2 ± 1
k 2 ± 1
i j ± k
j k ± i
i k ± j
Table 2. List of properties of 8 non-commutative division algebras.
Table 2. List of properties of 8 non-commutative division algebras.
Algebra i 2 j 2 k 2 ij jk ki ji kj ik Equivalent Matrix Determinant
H l e f t −1−1−1kijkij a b c d b a d c c d a b d c b a ( a 2 + b 2 + c 2 + d 2 ) 2
H r i g h t −1−1−1kijkij a b c d b a d c c d a b d c b a
Q + + k , i , j 11−1kijkij a b c d b a d c c d a b d c b a ( a 2 b 2 c 2 + d 2 ) 2
Q + + k , i , j 11−1kijkij a b c d b a d c c d a b d c b a
Q + + k , i , j 1−11kijkij a b c d b a d c c d a b d c b a ( a 2 b 2 + c 2 d 2 ) 2
Q + + k , i , j 1−11kijkij a b c d b a d c c d a b d c b a
Q + + k , i , j −111kijkij a b c d b a d c c d a b d c b a ( a 2 + b 2 c 2 d 2 ) 2
Q + + k , i , j −111kijkij a b c d b a d c c d a b d c b a
Table 3. Chilarity properties of the six split-quaternion algebras.
Table 3. Chilarity properties of the six split-quaternion algebras.
Algebra Chirality Equivalent Matrix Determinant
Q + + k , i , j ( a , b , c , d ) Left a b c d b a d c c d a b d c b a ( a 2 + b 2 c 2 d 2 ) 2
Q + + k , i , j ( a , b , c , d ) Right a b c d b a d c c d a b d c b a
Q + + k , i , j ( a , b , c , d ) Left a b c d b a d c c d a b d c b a ( a 2 b 2 + c 2 d 2 ) 2
Q + + k , i , j ( a , b , c , d ) Right a b c d b a d c c d a b d c b a
Q + + k , i , j ( a , b , c , d ) Left a b c d b a d c c d a b d c b a ( a 2 b 2 c 2 + d 2 ) 2
Q + + k , i , j ( a , b , c , d ) Right a b c d b a d c c d a b d c b a
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, H.-Y.; Cahay, M.; Vellambi, B.N.; Morris, D. A Generalization of Quaternions and Their Applications. Symmetry 2022, 14, 599. https://doi.org/10.3390/sym14030599

AMA Style

Lin H-Y, Cahay M, Vellambi BN, Morris D. A Generalization of Quaternions and Their Applications. Symmetry. 2022; 14(3):599. https://doi.org/10.3390/sym14030599

Chicago/Turabian Style

Lin, Hong-Yang, Marc Cahay, Badri N. Vellambi, and Dennis Morris. 2022. "A Generalization of Quaternions and Their Applications" Symmetry 14, no. 3: 599. https://doi.org/10.3390/sym14030599

APA Style

Lin, H. -Y., Cahay, M., Vellambi, B. N., & Morris, D. (2022). A Generalization of Quaternions and Their Applications. Symmetry, 14(3), 599. https://doi.org/10.3390/sym14030599

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop