A Spacetime Symmetry Approach to Relativistic Quantum Multi-Particle Entanglement

A Lorentz transformation group SO(m, n) of signature (m, n), m, n ∈ N, in m time and n space dimensions, is the group of pseudo-rotations of a pseudo-Euclidean space of signature (m, n). Accordingly, the Lorentz group SO(1, 3) is the common Lorentz transformation group from which special relativity theory stems. It is widely acknowledged that special relativity and quantum theories are at odds. In particular, it is known that entangled particles involve Lorentz symmetry violation. We, therefore, review studies that led to the discovery that the Lorentz group SO(m, n) forms the symmetry group by which a multi-particle system of m entangled n-dimensional particles can be understood in an extended sense of relativistic settings. Consequently, we enrich special relativity by incorporating the Lorentz transformation groups of signature (m, 3) for all m ≥ 2. The resulting enriched special relativity provides the common symmetry group SO(1, 3) of the (1 + 3)-dimensional spacetime of individual particles, along with the symmetry group SO(m, 3) of the (m + 3)-dimensional spacetime of multi-particle systems of m entangled 3-dimensional particles, for all m ≥ 2. A unified parametrization of the Lorentz groups SO(m, n) for all m, n ∈ N, shakes down the underlying matrix algebra into elegant and transparent results. The special case when (m, n) = (1, 3) is supported experimentally by special relativity. It is hoped that this review article will stimulate the search for experimental support when (m, n) = (m, 3) for all m ≥ 2.


Introduction
Nature organizes itself using the language of symmetries. Thus, in particular, the underlying symmetry group by which Einstein's special relativity theory can be understood is the Lorentz group SO c (1, 3). A physical system obeys the Lorentz symmetry if the relevant laws of physics are invariant under Lorentz transformations. Lorentz symmetry is one of the cornerstones of modern physics. However, it is known that entangled particles in relativistic quantum mechanics involve Lorentz symmetry violation [1][2][3][4][5]. Indeed, several explorers exploit entangled particles to observe Lorentz symmetry violation; see, for instance, [6][7][8][9][10][11].
Quantum entanglement [12] was named by Einstein, "spooky action at a distance." It is a physical phenomenon that occurs when groups of particles interact in ways such that the quantum state of each particle cannot be described independently of the others, even when the particles are separated by a large distance. Instead, a quantum state must be described for a system of particles as a whole.
Understanding entanglement in relativistic settings has been a key question in relativistic quantum mechanics. Some results show that entanglement is observer-dependent [2]. The aim of this review article is, therefore, (i) to present results that demonstrate in extended relativistic settings that the Lorentz transformation groups SO c (m, n), m, n ∈ N, are the missing symmetry groups of multi-particle systems of entangled particles; and (ii) to stimulate the search for experimental support when m ≥ 2 and n = 3, in addition to the available experimental support when (m, n) = (1, 3) that Einstein's

Einstein Velocity Addition
Einstein's special relativity stems from a transformation law, the Lorentz transformation, for the spacetime coordinates of a point particle. Guided by mathematical patterns and by analogies with the intuitively clear Galilei transformation, we extend both Galilei and Einstein boosts of particles to multi-boosts, called bi-boosts, of multi-particle systems. This extension enriches special relativity to a relativistic theory that accommodates multi-particle entanglement. The journey to our enriched special relativity theory begins with the review of Einstein velocity addition law of relativistically admissible velocities and its gyroformalism.
Let c > 0 be any positive constant and let R n = (R n , +, ·) be the Euclidean n-space, n ∈ N, endowed with the common vector addition, +, and inner product, ·. The space of all n-dimensional relativistically admissible velocities is the c-ball R n c = {v ∈ R n : v < c} .
Einstein velocity addition is a binary operation, ⊕, in the c-ball R n c given by [20][21][22][23][24] u⊕v = 1 1 + u·v for all u, v ∈ R n c . In physical applications n = 3, but in geometry n ∈ N. Here, γ u is the Lorentz gamma factor, where u·v and v are the inner product and the norm in the ball, which the ball R n c inherits from its ambient space R n , and v 2 = v·v = v 2 . A nonempty set with a binary operation is a groupoid, so that the pair (R n c , ⊕) is an Einstein groupoid. The constant c > 0 represents the vacuum speed of light. In the Euclidean-Newtonian limit of large c, c → ∞, the ball R n c expands to the whole of its ambient space R n , as we see from (1), and the Einstein addition ⊕ in R n c descends to vector addition + in R n , as we see from (2) and (3). When the nonzero vectors u and v in the ball R n c of R n are parallel in R n , u v, that is, u = λv for some λ ∈ R \ 0, Einstein addition (2) descends to the Einstein addition of parallel velocities, u, v ∈ R n c , which was partially confirmed experimentally by the Fizeau's 1851 experiment [25]. The restricted Einstein addition (4) is both commutative and associative, so that it is a group operation, as Einstein noted in [26]; see [27], p. 142. However, the general Einstein addition is neither commutative nor associative, so that it is not a group operation. Rather, it is a gyrocommutative gyrogroup operation, a rich algebraic structure discovered more than 80 years later, in 1988 [20,[28][29][30], which we review in Section 5.
Einstein addition (2) with n = 3 was introduced by Einstein in his 1905 paper [26] ( [27], p. 141) that founded the special theory of relativity, where the magnitudes of the two sides of Einstein addition (2) are presented.
We use the notation u v = u⊕(−v) for Einstein subtraction, so that, for instance, v v = 0 and Einstein addition and subtraction satisfy the equations (u⊕v) = u v (6) and u⊕(u⊕v) = v (7) for all u, v in the ball R n c . Identity (6) expresses the gyroautomorphic inverse property of Einstein addition, and Identity (7) expresses the left cancellation law of Einstein addition.
Einstein addition does not obey the naive right counterpart of the left cancellation law (7) since, in general, (u⊕v) This seeming lack of a right cancellation law of Einstein addition is remedied in (76) by the incorporation of a second, dual binary operation, .
Einstein addition and the gamma factor are linked by the gamma identity, for all u, v ∈ R n c . The gamma identity signaled the emergence of hyperbolic geometry into special relativity ( [31], p. 79).
The structure of Einstein addition gives rise to our gyrolanguage in which we prefix a gyro to any term that describes a concept in Euclidean geometry and in associative algebra to mean the analogous concept in hyperbolic geometry and in nonassociative algebra. The prefix "gyro" stems from "gyration", which is the mathematical abstraction of the special relativistic effect known as "Thomas precession". In gyrolanguage, thus, Einstein addition is a gyroaddition. The gamma factor (3) associated with Einstein gyroaddition will be extended in (129) to the bi-gamma pair associated with Einstein bi-gyroaddition as a step in our program to extend the Lorentz group SO c (1, 3) to SO c (m, n) for all m, n ∈ N under the same umbrella.
A useful identity that follows immediately from (3) is Einstein addition is noncommutative. Indeed, while Einstein addition is commutative under the norm, u⊕v = v⊕u , (11) in general, it is noncommutative, u⊕v = v⊕u , (12) u, v ∈ R n c . Einstein addition is nonassociative as well since, in general, (u⊕v)⊕w = u⊕(v⊕w) , (13) u, v, w ∈ R n c . As an application of the gamma identity (9), we present the Einstein gyrotriangle inequality in the following theorem.
Proof. By the gamma identity (9) and the Cauchy-Schwarz inequality [32], we have for all u, v in an Einstein groupoid (R n c , ⊕), n ∈ N. The function γ x = γ x is a monotonically increasing function of x , 0 ≤ x < c. Hence (15) implies (14).
The gyrotriangle inequality (14) will be extended from Einstein gyrovector spaces (R n c , ⊕, ⊗) to Einstein bi-gyrovector spaces (R n×m c , ⊕ E , ⊗) of any signature (m, n), m, n ∈ N, in (175). Signature (m, n)). Einstein addition ⊕ in (2) is of signature (1, n). It will be generalized in Section 22 to the Einstein addition ⊕ E = ⊕ E,(m,n),c of signature (m, n) in the ball R n×m c for all m, n ∈ N. Einstein addition ⊕ in (2) will, then, be recognized as the Einstein addition ⊕ E,(1,n),c of signature (1, n) in the ball R n×1 c = R n c , that is, ⊕ = ⊕ E = ⊕ E,(1,n),c .

The Gyroformalism of Einstein Addition
Vector addition, +, in R n is both commutative and associative; that is, for all u, v, w ∈ R n . In contrast, Einstein addition, ⊕, in R n c is neither commutative nor associative, where the deviation from both commutativity and associativity is controlled by gyrations, as evidenced from (18).
Gyrations gyr[u, v] ∈ Aut(R n c , ⊕), u, v ∈ R n c , are given in terms of Einstein addition by the gyrator equation for all u, v, w ∈ R n c . Equation (17) presents the application to w of the gyration gyr[u, v] generated by u and v. Gyrations are automorphisms of the Einstein groupoid (R n c , ⊕). An automorphism of a groupoid (S, ⊕) is a bijective map f of S onto itself that respects the binary operation, that is, The set of all automorphisms of a groupoid (S, ⊕) forms a group, denoted by Aut(S, ⊕), where the group operation is given by automorphism composition. To emphasize that the gyrations of an Einstein groupoid (R n c , ⊕) are automorphisms of the groupoid, gyrations are also called gyroautomorphisms.
Possessing their own rich structure, gyrations measure the extent to which Einstein addition deviates from commutativity and associativity as we see from the following list of identities [20,21,31]: Gyrocommutative Law u⊕(v⊕w) = (u⊕v)⊕gyr [u, v] Gyration Inversion Law (18) for all u, v, w ∈ R n c . It is clear from the gyrocommutative and gyroassociative laws in (18) that the departure of Einstein addition, ⊕, from commutativity and associativity is controlled by gyrations.
The properties of Einstein gyroaddition and resulting gyrations in (18) give rise to a mathematical formalism and its associated gyrolanguage. These properties will be extended in Section 22 to Einstein bi-gyroaddition and bi-gyrations of signature (m, n) for all m, n ∈ N, where the special signature (m, n) = (1, 3) corresponds to Einstein's special relativity in m = 1 time dimension and n = 3 space dimensions.
In Definition 3, the left reduction property is elevated to the reduction axiom. As noted by F. Chatelin in [41], the reduction axiom triggers remarkable reduction in complexity as, for instance, in (72).
Gyrations are expressed in (17) in terms of Einstein addition. Explicitly, they are given by the equation where for all u, v, w ∈ R n c , as shown in [19]. In the three special cases when (i) u = 0, or (ii) v = 0, or (iii) u and v are parallel in R n , u v, we have Au + Bv = 0, so that gyr[u, v] is trivial; that is, for all u, v ∈ R n c , where u v in the third equation, and for all w ∈ R n . Following (20) for all u, v ∈ R n c . Gyrations keep the inner product invariant; that is, for all a, b, u, v ∈ R n c . Hence, gyr [u, v] is an isometry of R n c , keeping the norm of elements of the ball R n c invariant.
Accordingly, gyr[u, v] represents a rotation of the ball R n c about its origin for any u, v ∈ R n c .  (28) for all a, b, u, v ∈ R n c . Hence, gyr [u, v] is an automorphism of the Einstein groupoid (R n c , ⊕).
Example 1. As an example that illustrates the use of the invariance of the norm under gyrations, we note that Indeed, we have the following chain of equations, which are numbered for subsequent derivation.

From Einstein Velocity Addition to Gyrogroups
Guided by analogies with groups, the key features of Einstein groupoids (R n c , ⊕), n = 1, 2, 3, . . ., suggest the formal gyrogroup definition in which gyrogroups form a most natural generalization of groups.

Definition 1. (Binary Operations).
A binary operation + in a set S is a function + : S×S → S. We use the notation a + b to denote +(a, b) for any a, b ∈ S.

Definition 2. (Groupoids, Automorphisms).
A groupoid (S, +) is a nonempty set, S, with a binary operation, +. An automorphism φ of a groupoid (S, +) is a bijective self-map of S which respects its groupoid operation; that is, φ(a + b) = φ(a) + φ(b) for all a, b ∈ S. The automorphisms of a groupoid (S, +) form a group denoted by Aut(S, +).

Definition 3. (Gyrogroup ([31] Definition 2.5)).
A groupoid (G, ⊕) is a gyrogroup if its binary operation satisfies the following axioms. In G there is at least one element, 0, called a left identity, satisfying (G1) 0⊕a = a for all a ∈ G. There is an element 0 ∈ G satisfying axiom (G1) such that for each a ∈ G there is an element a ∈ G, called a left inverse of a, satisfying (G2) a⊕a = 0 . Moreover, for any a, b, c ∈ G there exists an automorphism gyr[a, b] ∈ Aut(G, ⊕) such that the binary operation obeys the left gyroassociative law (G3) a⊕(b⊕c) = (a⊕b)⊕gyr[a, b]c .
The automorphism gyr[a, b] of G is called the gyroautomorphism, or the gyration, of G generated by a, b ∈ G. The operator gyr : G × G → Aut(G, ⊕) is called the gyrator of G. Finally, the gyroautomorphism gyr[a, b] generated by any a, b ∈ G obeys the left reduction axiom As in group theory, we use the notation a b = a⊕( b) in gyrogroup theory as well.
In full analogy with groups, gyrogroups split up into gyrocommutative and non-gyrocommutative ones.

Gyrogroup Cooperation (Coaddition)
In order to capture analogies with groups, we introduce into the abstract gyrogroup (G, ⊕) a second binary operation, , called the gyrogroup cooperation, or coaddition. The two binary operations ⊕ and in G are dual to each other in some sense [20].
Naturally, we use the notation The gyrogroup cooperation is commutative if and only if the gyrogroup operation is gyrocommutative, as stated in ([31] Theorem 3.4, p. 50). In particular, Einstein coaddition is commutative since Einstein addition ⊕ is gyrocommutative.
With the emergence of Einstein coaddition, the dream of the famous mathematician Émile Borel to create a commutative variant of Einstein addition came true [55].

First Gyrogroup Properties
It is clear how to define a right identity and a right inverse in a gyrogroup. The existence of such elements is not presumed since the existence of a unique identity and a unique inverse, both left and right, is a consequence of the gyrogroup axioms. This result is stated in the following theorem, along with other important basic results in gyrogroup theory.
Hence, by Item (1) above we have x = gyr[0, a]x for all x ∈ G so that gyr[0, a] = I, I being the trivial (identity) map. 3. By the left reduction property and by Item (2) above we have 4. That follows from an application of the left reduction property and Item (2) above. Thus, 5. Let x be a left inverse of a corresponding to a left identity, 0, of G. Then by left gyroassociativity and Item (3) above, Hence, by (1), a⊕0 = a for all a ∈ G so that 0 is a right identity. 6. Suppose 0 and 0 * are two left identities, one of which, say 0, is also a right identity. Then 0 = 0 * ⊕0 = 0 * . 7. Let x be a left inverse of a. Then by left gyroassociativity, (G2) of Definition 3 and Items (3), (5), (6) above. By Item (1) we have a⊕x = 0 so that x is a right inverse of a. 8. Suppose x and y are left inverses of a. By Item (7) above, they are also right inverses, so a⊕x = 0 = a⊕y.
By Item (1), x = y. Let a be the resulting unique inverse of a. Then a⊕a = 0 so that the inverse ( a) of a is a. 9. By left gyroassociativity and by Item (3) we have 10. By an application of the left cancellation law in Item (9) to the left gyroassociative law (G3) in Definition 3 we obtain the result in Item (10) . 11. We obtain Item (11) from Item (10) with x = 0. 12. Since gyr[a, b] is an automorphism of (G, ⊕), we have from Item (11) gyr and hence the result. 13. We obtain Item (13) from Item (10) with b = 0, and a left cancellation, Item (9).

Elements of Gyrogroup Theory
Einstein gyrogroups (G, ⊕) possess the gyroautomorphic inverse property, according to which (a⊕b) = a b for all a, b ∈ G. In general, however, (a⊕b) = a b in some gyrogroups. Hence, the following theorem is interesting.
Proof. The first identity in (44) follows from Theorem 2 (10) with x = b, and Theorem 2 (12), and the second part of Theorem 2 (8). The second identity in (44) follows from the first one by replacing b by b, noting Theorem 2 (12).
A nested gyroautomorphism is a gyration generated by points that depend on another gyration. Thus, for instance, some gyrations in (45) Proof. By two successive applications of the left gyroassociative law in two different ways, we obtain the following two chains of equations for all a, b, c, x ∈ G, and By comparing the extreme right-hand sides of these two chains of equations, and by employing the left cancellation law, Theorem 2 (1), we obtain the identity for all x ∈ G, thereby verifying (45).
In the special case when c = b, (45) descends to (46); note that the left-hand side of (45) becomes trivial owing to Items (2) and (3) of Theorem 2.
Identity (47) results from the following chain of equations, which are numbered for subsequent derivation: Derivation of the numbered equalities in (53): 1. Follows from (46). 2. Follows from (1) by the left reduction property. 3. Follows from (2) by the left gyroassociative law. Indeed, an application of the left gyroassociative law to the first entry of the left gyration in (3) gives the first entry of the left gyration in (2); that is, To verify (48) we consider the special case of (45) where the second identity in (54) follows from Items (2) and (3) of Theorem 2. Through replacing a by a and c by b in (54), we obtain (48). Finally, (49) is derived from (48) by an application of the left reduction property to the first gyroautomorphism in (48) followed by a left cancellation, Theorem 2 (9). Accordingly, as desired.
The nested gyroautomorphism identity (47) in Theorem 5 allows the equation that defines the coaddition to be dualized with its corresponding equation in which the roles of the binary operations and ⊕ are interchanged, as shown in the following theorem. Then for all a, b ∈ G.
Proof. Let a and b be any two elements of G. By (56) and (47) we have thus verifying (57).
Hence, a a = a a = 0 .
Identity (60) implies the equality between the inverses of a ∈ G with respect to ⊕ and , for all a, b, c ∈ G.
Proof. By the left gyroassociative law and the left cancellation law, and using the notation d = b⊕c, we have (63)
Proof. Identity (64) is a rearrangement of (62) obtained by left-gyroadding the term ( a⊕b) to both sides of (62) followed by a left cancellation.
The significance of Identity (64) stems from the analogy it shares with its group counterpart, −(−a + b) + (−a + c) = −b + c in any group with group operation +.

Basic Gyrogroup Equations
The basic equations of gyrogroup theory are and a, b, x ∈ G, each equation for the unknown x in a gyrogroup (G, ⊕). Theorem 9 below asserts that each of the two basic equations possesses a unique solution. The unique solutions give rise in Section 10 to the basic gyrogroup cancellation laws. in G for the unknown x is x = a⊕b (68) and the unique solution of the equation Proof. Let x be a solution of the first basic equation, (65). Then we have by (65) and the left cancellation law, Theorem 2 (9), p. 9, Hence, if a solution x of (65) exists then it must be given by x = a⊕b, as we see from (71).
Conversely, x = a⊕b is, indeed, a solution of (65) as wee see by substituting x = a⊕b into (65) and applying the left cancellation law in Theorem 2 (9). Hence, the gyrogroup equation (65) possesses the unique solution x = a⊕b.
The solution of the second basic gyrogroup equation, (66), is quiet different from that of the first, (65), owing to the noncommutativity of the gyrogroup operation. Let x be a solution of (66). Then we have the following chain of equations, which are numbered for subsequent derivation: Derivation of the numbered equalities in (72): 1. Follows from the existence of a unique identity element, 0, in the gyrogroup (G, ⊕) by Theorem 2, p. 9. 2. Follows from the existence of a unique inverse element a of a in the gyrogroup (G, ⊕) by Theorem 2. 3. Follows from (2) by the left gyroassociative law in Axiom (G3) of gyrogroups in Definition 3, p. 7. 4. Follows from (3) by Theorem 2 (12). 5. Follows from (4) by the left reduction property (G4) of gyrogroups in Definition 3. 6. Follows from (5) by the assumption that x is a solution of (66). 7. Follows from (6) by (59).
Hence, if a solution x of (66) exists then it must be given by x = b a, as we see from (72).
Conversely, x = b a is, indeed, a solution of (66), as we see from the following chain of equations: Derivation of the numbered equalities in (73): 1. Follows from the assumption that x = b a. (1) by (59). 3. Follows from (2) by Identity (47) (5) since 0 is the unique identity element of the gyrogroup (G, ⊕).

Basic Gyrogroup Cancellation Laws
Basic cancellation laws of gyrogroup theory are obtained from the basic equations of gyrogroups solved in Section 9. By substituting the solution (68) into its equation (Equation (67)), we obtain the left cancellation law for all a, b ∈ G, stated in Theorem 2 (9). Similarly, by substituting the solution (70) into its equation (Equation (69)), we get the first right cancellation law for all a, b ∈ G. The latter can be dualized, obtaining the second right cancellation law for all a, b ∈ G. Indeed, (76) follows from the chain of equations where one employs the left gyroassociative law, the left reduction property, and the gyrogroup cooperation definition. Identities (74)-(76) form the three basic cancellation laws of gyrogroup theory ([19] Section 2.11).

Automorphisms and Gyroautomorphisms
In this section we find a commuting relation between automorphisms of a gyrogroup and its gyroautomorphisms.
Proof. For any three elements a, b, x ∈ (G, ⊕) and any automorphism A ∈ Aut(G, ⊕) we have by the left gyroassociative law, Hence, by a left cancellation, for all x ∈ G, implying (78).
Theorem 12 enhances the duality symmetry that Einstein addition and coaddition share in Theorem 6 and in (75) and (76).

On the Rich Gyrostructure that Stems from Einstein Addition
The introduction of a coaddition, , into the abstract gyrogroup (G, ⊕) is dictated by our wish to capture analogies with classical results.
Einstein addition, ⊕, admits scalar multiplication, ⊗, which turns Einstein gyrocommutative gyrogroups (R n c , ⊕) to Einstein gyrovector spaces (R n c , ⊕, ⊗). The latter are extended by abstraction to abstract gyrovector spaces, which are fully analogous to vector spaces. Thus, in particular, (i) Einstein gyrovector spaces form the algebraic setting for the Beltrami-Klein ball model of hyperbolic geometry; and (ii) Möbius gyrovector spaces form the algebraic setting for the Poincaré ball model of hyperbolic geometry, just as (iii) vector spaces form the algebraic setting for the standard model of Euclidean geometry, as shown, for instance, in [20,21,31].
Within the frame of gyrovector spaces we study concepts such as gyrolines, cogyrolines, gyroangles and cogyroangles, in full analogy with the study of their classical counterparts. Remarkably, along with Einstein addition, it is Einstein coaddition that enables us to capture analogies between parallelograms in the Euclidean geometry of R n and gyroparallelograms in the hyperbolic geometry of the ball R n c . For details, see [19][20][21]31,54,[56][57][58]), The importance of understanding the rich gyrostructure that stems from Einstein addition to our study of Lorentz symmetry groups rests on the fact that (i) relativistically admissible velocities in the ball R n c parametrize the Lorentz groups; and (ii) compositions of Lorentz transformations involve both Einstein addition and gyroautomorphisms, as we will see in (103).

Galilei and Lorentz Boosts and Multi-boosts
In order to pave the road to extending the Galilei and Lorentz boosts to multi-boosts called bi-boosts of signature (m, n), m, n ∈ N, we present in this section the common Galilei and Lorentz boosts along with an additive decomposition (104) that associates Galilei and Lorentz boosts to each other.
A Galilei transformation without rotations is called a Galilei boost. The Galilei boost is a linear transformation of time-space coordinates in (1 + 3)-dimensional time-space, represented by the matrix The action of a Galilei boost G(v) on the time-space coordinates (t, x) t of a particle with position x at time t gives where and where exponent t denotes transposition.
Here v represents velocity and (t, x) t represents a particle with position x at time t. Accordingly, the Galilei boost B(v) in (96) boosts a particle with position x at time t to a boosted particle with position x + vt at time t.
Extending from 3 to n space dimensions, the Galilei boost is represented by the (n + 1) × (n + 1) matrix and its action on time-space coordinates gives where v = (v 1 , v 2 , . . . , v n ) t ∈ R n and x = (x 1 , x 2 , . . . , x n ) t ∈ R n are column vectors. Here, I n is the n × n identity matrix and 0 m,n is the m × n zero matrix. Galilei boosts are parametrized by the velocity parameter v, and the Galilei boost composition is given by velocity addition, that is, for all v 1 , v 2 ∈ R n . Obviously, the inverse of Anticipating the extension of the Galilei boost to m time and n space dimensions for all m, n ∈ N, we call G(v) in (97)-(99) a Galilei boost of signature (1, n). The Galilei boosts of signature (1, n) form a group, and the group of Galilei boosts of signature (1, 3) forms the symmetry group of classical mechanics.
A Lorentz boost is a Lorentz transformation without rotations. The Lorentz boosts B c (v), v ∈ R n c , that correspond to the Galilei boosts of signature (1, n) are represented by the (n + 1) × (n + 1) matrix parametrized by the velocity vector v ∈ R n c . The inverse of B c (v) is given by We clearly have the limit in which a Lorentz boost descends to its Galilean counterpart, for any v ∈ R n c . The composition of two Galilei boosts of signature (1, n) is, again, a Galilei boost of signature (1, n), as we see from (99). In contrast, the composition of two Lorentz boosts of signature (1, n) is not a Lorentz boost. Rather, it is a Lorentz boost of signature (1, n) preceded (or, followed) by a space gyration (a gyration is a rotation known in special relativity as Thomas precession), for all v 1 , v 2 ∈ R n c , where Einstein addition ⊕ in the ball R n c is given by (2), and the gyration gyr[v 1 , v 2 ] is given in terms of Einstein addition in Theorem 2 (10). The study of the Lorentz boost composition law (103) is presented in ([21] Theorem 4.10).
The Lorentz boost composition law (103) involves Einstein addition along with a gyration. As such, it demonstrates the important role that the gyroalgebra of Einstein addition plays in our understanding of the Lorentz transformation group of signature (1, n). The gyration gyr[v 1 , v 2 ] in (103) represents a space rotation. In the extension of our gyroformalism to Lorentz transformations of signature (m, n), m, n ≥ 2, a second gyration, which represents a time rotation, emerges. Naturally, the pair of a space and a time gyration will be called a bi-gyration, as we will see in the sequel.
The Lorentz transformation group of signature (1, n) is the symmetry group of Einstein's special relativity, where n = 3 in physical applications. In order to understand the Lorentz transformation B c (v), v ∈ R n c , of signature (1, n) in terms of its associated Galilei transformation of signature (1, n), we consider the following additive decomposition of B c (v) ([19] Section 6.3.1), which one can readily prove by means of (10), The additive decomposition (104) of B c (v) has two components. It can be written as Accordingly, it is intuitively clear how to extend the Galilei boost from the common boost that acts on particles individually to a multi-boost that acts on several particles collectively. In contrast, it is not yet clear how to extend the Lorentz boost from the common boost that acts on particles individually to a multi-boost that acts on several particles collectively.
Fortunately, the intuitively clear Galilei multi-boosts come to the rescue. They will guide us in the search for the Lorentz multi-boosts that act on several particles collectively, rather than individually. Remarkably, the Lorentz multi-boosts that the Galilei multi-boosts suggest, will turn out in the sequel to be nothing else but the well-known Lorentz transformations of signature (m, n), m, n > 1.
We, therefore, proceed with the study of Lorentz transformations of any signature (m, n), m, n ∈ N, before turning to study Galilei multi-boosts. Our goal is, accordingly, to discover how to extend the additive decomposition (104) of the Lorentz boost from signature (1, n) to any signature (m, n), m, n ∈ N.

Pseudo-Euclidean Spaces and Lorentz Transformations of any Signature
We now face the task to discover that the rich gyroalgebra of Einstein addition and Lorentz transformations of signature (1, n) that we have studied in Sections 2-13 within the framework of special relativity, survive unimpaired in their transition from signature (1, n) to any signature (m, n), m, n ∈ N ([19] Section 4.2), [59].
The study of pseudo-Euclidean spaces and pseudo-rotations is well-known in linear algebra [13] and in physics [14]. A pseudo-Euclidean space R m,n , m, n ∈ N, of signature (m, n) is an (m + n)-dimensional space with an orthogonal basis e i , i = 1, . . . , m + n, where δ ij , j = 1, . . . , m + n, is the Kronecker delta, where and where c is an arbitrarily fixed positive constant.
Without loss of generality one may select c = 1. However, we prefer to view c as a free positive parameter, 0 < c < ∞, in order to allow limits as c → ∞ to be available, where the modern and new descends to the classical and familiar. An illustrative point in case is the Lorentz transformation of special relativity theory, which descends to the Galilei transformation of classical mechanics when the vacuum speed of light c tends to infinity, as shown in (149)-(155) with any signature (m, n). The usefulness of viewing c as a positive parameter rather than selecting c = 1 is clearly demonstrated in the additive decomposition (155) of the Lorentz multi-boost, which plays an important role in our study.
A pseudo-Euclidean space R m,n of signature (m, n) is equipped with an inner product of signature (m, n). The inner product x·y of two vectors x, y ∈ R m,n , is given by Accordingly, the squared norm of x ∈ R m,n is Let η be the (m + n) × (m + n) diagonal matrix where I m is the m × m identity matrix and 0 m,n is the m × n zero matrix. Then, the matrix representation of the inner product (110) is given by where x, y ∈ R m,n are the column vectors and exponent t denotes transposition. Let Λ be an (m + n) × (m + n) matrix that leaves the inner product (113) invariant. Then, for all x, y ∈ R m,n , (Λx) t ηΛy = x t ηy , implying x t Λ t ηΛy = x t ηy, so that ( [60] p. 193) The determinant of the matrix equation (116) yields noting that det(Λ t ηΛ) = (detΛ t )(detη)(detΛ) and detΛ t = detΛ. Hence, detΛ = ±1. The special transformations Λ that can be reached continuously from the identity transformation of R m,n constitute the special pseudo-orthogonal group SO c (m, n), called the group of pseudo-rotations of signature (m, n). A pseudo-rotation of signature (m, n) is also called a Lorentz transformation of signature (m, n), or in short, an (m, n)-Lorentz transformation. The Lorentz transformation of signature (1, 3) turns out to be the common homogeneous, proper, orthochronous Lorentz transformation of Einstein's special theory of relativity [21]. The formal definition of SO c (m, n) follows. If Λ ∈ SO c (m, n) then its determinant is equal to 1, and the determinant of its first m rows and columns is positive ( [61] p. 478).

Matrix Balls of Radius c
The concept of relativistically admissible velocities is extended in this section from velocity vectors in a ball R n c ⊂ R n of vectors to velocity matrices in a ball R n×m c ⊂ R n×m of matrices.

Definition 7. (Spectrum, Spectral Radius ([62] p. 35)).
Let R n×m be the set of all real n × m matrices, m, n ∈ N. The set of all complex numbers that are eigenvalues of a square matrix A ∈ R n×n is called the spectrum of A and is denoted by σ(A). The spectral radius of A is the nonnegative number The spectral radius ρ(A) is thus the radius of the smallest disc centered at the origin of the complex plane that includes all the eigenvalues of A.
For Definition 8 below, we note that for any V ∈ R n×m the set of nonzero eigenvalues of VV t ∈ R n×n equals the set of nonzero eigenvalues of V t V ∈ R m×m .

Definition 8. (Matrix Ball, Matrix Spectral Norm ([19] Definition 5.7)
). For any m, n ∈ N and c > 0, the c-ball R n×m c of the ambient space R n×m of all n × m real matrices is given by The matrix spectral norm V of V ∈ R n×m , or norm in short, is defined by It is clear from Definition 8 that so that the matrix ball R n×m c ⊂ R n×m is a natural generalization of the ball R n c = R n×1 c ⊂ R n×1 = R n in (1) from signature (1, n) to any signature (m, n), m, n ∈ N.
Well-known properties of the matrix spectral norm are [62]: for any r ∈ R and square matrices A, B ∈ R k×k , k ∈ N. exists.
The obvious limits will prove useful in (152). We use the abbreviated notation Γ L V and Γ R V in (123) and (124), rather than the full notation Γ L n,V,c and Γ R m,V,c , since the values of m, n and c are always known from the context. We call Γ L V and Γ R V , respectively, the left and right gamma factors of signature (m, n). Collectively, the pair (Γ L V , Γ R V ) is called the bi-gamma factor.
In view of (125), it is interesting to note that for any fixed m, n ∈ N and c > 0, the left and right gamma factors of signature (m, n) satisfy the matrix identities ( [19] Lemma 5.82) In (126) we use the convenient matrix division notation A/B to denote either AB −1 or B −1 A when no confusion may arise, that is, when the matrices A and B satisfy AB −1 = B −1 A. The identities in (126) prove useful in establishing the additive decomposition (155), which plays an important role in the interpretation of the Lorentz transformation of signature (m, n), m, n ∈ N.
In the special case when m = 1, the right gamma factor, Γ R V in (124), descends to the Lorentz gamma factor, γ V in (3), of special relativity theory in a single time dimension and n space dimensions, It can be shown that, for m = 1, the left counterpart of (127) is given by Equation (5.172) in [19].
V ∈ R n×1 c = R n c , which holds for m = 1 and all n ∈ N.

Bi-Gamma Factor
Analogies with the gamma factor of special relativity theory in (3) and (127) suggest the definition of a left gamma factor Γ L V and a right gamma factor Γ R V by the equations V ∈ R n×m c . It clearly follows from (129) that for all V ∈ R n×m c . Following (130) we have interesting matrix identities as, for instance, the relations for any V ∈ R n×m c . The matrices Γ L V and Γ R V are symmetric. Furthermore, they are even in the parameter V ∈ R n×m c , that is, Naturally, the pair (Γ L V , Γ R V ) of a left and a right gamma factor is called a bi-gamma factor.

Theorem 14. (Bi-gamma Commuting Relation).
For any two positive integers m, n ∈ N, the left and right gamma factors are related by the commuting relation Proof. For simplicity, we assume c = 1. For c > 0 the proof is similar. By (129) and (133), we have to prove the matrix identity or, equivalently, The diagonal matrix Hence, so that, by (136) and (138), and By means of (139) we have the chain of equations implying Similarly, Hence, by (142) and (136) we have and by (143) and (136) we have The extreme right-hand sides of (144) and (145) are equal, thereby implying (135), and the proof is complete.
Following (133), the left and right gamma factors are related by the first commuting relation in (146) below. The remaining commuting relations in (146) follow immediately from the first one; note that left and right gamma factors are symmetric matrices. (146)

V-Parametric Realization of Lorentz Transformations of Signature (m,n)
The following theorem realizes the (m, n)-Lorentz group SO c (m, n) parametrically. Each element Λ ∈ SO c (m, n) is a Lorentz transformation of signature (m, n). It is parametrized by a velocity matrix parameter V ∈ R n×m c in the ball, and two orientation parameters O m ∈ SO(m) and O n ∈ SO(n), according to the following Theorem 15. The term velocity matrix will be justified in the sequel, where we will see by analogies with Galilei transformations that the columns of a velocity matrix represent relative velocities of the constituents of a multi-particle system.
parametrized by the velocity matrix V ∈ R n×m c , and the two orientation parameters O m ∈ SO(m) and O n ∈ SO(n).
A slightly different decomposition, called the polar decomposition of Lorentz groups of any signature, is provided by the following theorem.
for any V ∈ R n×m c , O m ∈ SO(m) and O n ∈ SO(n).
It follows from (147), or (148) that a Lorentz transformation Λ of signature (m, n) without rotations is what we call a Lorentz bi-boost, B c (V), of signature (m, n), parametrized by the velocity matrix V ∈ R n×m c . Like (101), the inverse of B c (V) is given by The matrix B c (V) is a bi-boost (as opposed to a boost) in the sense that it admits in (147) and in (148) the bi-rotation (O m , O n ) ∈ SO(m) × SO(n).
Here B c (V) ∈ R (m+n)×(m+n) is an (m + n) × (m + n) block matrix consisting of the four blocks Γ R V ∈ R m×m , Γ L V ∈ R n×n , Γ L V V ∈ R n×m , and c −2 Γ R V V t ∈ R m×n . An elegant, straightforward demonstration that B c (V) ∈ SO c (m, n), as stated in (149), based on the commuting relations in (146), is presented in ([19] Section 5.8).
In the special case when the signature is (m, n) = (1, n), the Lorentz bi-boost B c (V) in (149) descends to the common Lorentz transformation of Einstein's special theory of relativity, as shown in (151), where the parameter V ∈ R n×1 c = R n c represents relativistically admissible n-dimensional velocities.
Indeed, when m = 1 the bi-boost B c (V) ∈ SO(m, n) in (149) can be manipulated by means of (146) and by means of (127) and (128), obtaining the following chain of equations.
where V ∈ R n×1 c ⊂ R n×1 = R n is a column vector in the ball R n×1 c = R n c of R n and where γ V is given by (127).
The extreme right-hand side of (151) turns out to be the standard special relativistic (n + 1) × (n + 1) matrix representation of the Lorentz group in one time dimension and n space dimensions [28] ( [20] p. 254) ( [21] p. 447), presented in (100). Accordingly, it follows from (151) that in the special case when m = 1 the Lorentz group of signature (m, n) specializes to the Lorentz group of special relativity theory, where (m, n) = (1, n).
The bi-gyration decomposition (147) presents the generic Lorentz transformation of signature (m, n) as a bi-boost along with a left rotation O n ∈ SO(n) and a right rotation O m ∈ SO(m) acting on V ∈ R n×m c . Collectively, the pair (O n , O m ) of a left and a right rotation, taking V into O n VO m , is called a bi-rotation.
In the limit of large c, the Lorentz bi-boost B c (V) of signature (m, n) tends to its Galilean counterpart B ∞ (V), called the Galilei bi-boost of signature (m, n). Indeed, by means of (125) we have is the group of all Galilei transformations of signature (m, n). The Galilei bi-boost of signature (1, 3) is the common Galilei boost of classical mechanics. The physical interpretation of the Galilei bi-boost of signature (1, 3) is intuitively clear. We will see in Section 20 that the physical interpretation of the Galilei bi-boost of signature (m, 3) for any m > 1 is intuitively clear as well. The physical interpretation of the Galilei bi-boost of signature (m, 3), in turn, leads to a physical interpretation of the Lorentz bi-boost of signature (m, 3), by means of the additive decomposition (155) of the Lorentz bi-boost, which we present in Section 18.

Theorem 17. (Additive Decomposition of the Lorentz Bi-boost ([19] Theorem 5.83)). Let
V ∈ R n×m c ⊂ R n×m , be the Lorentz bi-boost of signature (m, n), m, n ∈ N, and let V ∈ R n×m , be the Galilei bi-boost of signature (m, n). Then, B c (V) and B ∞ (V) are related to each other by the additive decomposition of the Lorentz bi-boost, for all V ∈ R n×m c .
In (155) we continue using the convenient matrix division notation, as in (126). The additive decomposition (155) of the Lorentz bi-boost B c (V) forms the natural generalization of the additive decomposition (104) of the Lorentz boost to the Lorentz bi-boost, that is, a generalization from signature (1, n) to any signature (m, n), m, n ∈ N.
The additive decomposition (155) follows immediately from (153) and (126), and it provides a correspondence between Lorentz and Galilei bi-boosts. Thus, the set of all Lorentz bi-boosts B c (V), V ∈ R n×m c , of signature (m, n) is associated in (155) with the set of all Galilei bi-boosts B ∞ (V), V ∈ R n×m , of same signature (m, n), and conversely, the set of all Galilei bi-boosts of signature (m, n) is associated in (155) with the set of all Lorentz bi-boosts of same signature (m, n), m, n ∈ N. This association is important, enabling us to interpret Lorentz transformations of signature (m, n) in terms of the intuitively clear interpretation of corresponding Galilei transformations of same signature (m, n) that we will study in Section 20. It suggests the following formal definition.

Definition 9. (Additive Decomposition of the Lorentz Bi-boost ([19] Definition 5.84)). Let
V ∈ R n×m c , and V ∈ R n×m , be the Lorentz bi-boost of signature (m, n) and its corresponding Galilei bi-boost of same signature (m, n), m, n ∈ N. Furthermore, let V ∈ R n×m c , so that, by Theorem 17, we have the Lorentz bi-boost additive decomposition

Application of the Galilei Bi-Boost of Signature (1,3)
In this section about the Galilei bi-boost of signature (1, 3), we make way on the road to the presentation of the Galilei bi-boost of signature (2,3), and hence, of signature (m, 3), m > 1, in Section 20.
A Galilei bi-boost of signature (1, 3) is the common Galilei boost of classical mechanics. Let B ∞ (V) = B ∞ (v) be the Galilei bi-boost of signature (m, n) = (1, 3), parametrized by the velocity parameter V = (v), In order to conform to the formalism of bi-boosts of signature (m, n), m, n > 1, in Section 20, we view V in (160) as a 3 × 1 matrix in R 3×1 , and v as a vector in R 3 , while noting that R 3×1 = R 3 . Accordingly, V = (v) is a matrix V the single column of which is v, so that v is a velocity vector and V is a velocity matrix. Let be a 4 × 1 matrix that represents the time-space coordinates of a particle with position x = (x 1 , x 2 , x 3 ) t ∈ R 3 at time t ∈ R. The point (t, x) is said to be a particle of signature (m, n)=(1, 3) with position x ∈ R 3 at time t ∈ R. A particle of signature (m, n) is also called an (m, n)-particle, in short.
The application of the Galilei bi-boost B ∞ (V) of signature (m, n) = (1, 3) to a (1, 3)-particle (t, x) in m + n = 1 + 3 time-space dimensions is described in the following chain of equations.
Here, B ∞ (V) in (162) is given by (157) with m = 1 and n = 3. Accordingly, the Galilei bi-boost B ∞ (V) of signature (1, 3) is a Galilei boost that keeps the time invariant, t = t, and boosts the position x ∈ R 3 of the particle (t, x) t into the position x = x + vt ∈ R 3 , v ∈ R 3 , of the boosted particle (t, x + vt) t , at time t.

Application of the Galilei Bi-boost of any Signature
We are now in the position to consider the Galilei bi-boost of signature (m, n) for all m, n ∈ N, while paying special attention to the case when (m, n) = (2, 3) as an illustrative example. Let of two velocity vectors v k = (v 1k , v 2k , v 3k ) t ∈ R 3 , k = 1, 2. These two velocity vectors form the two columns of the velocity matrix V, in analogy with (160). Furthermore, let be a 5 × 2 matrix that represents a (2, 3)-particle consisting of the time-space coordinates of two subparticles, (t k , x k ), k = 1, 2, with positions x k = (x 1k , x 2k , x 3k ) t ∈ R 3 , at time t k ∈ R, respectively. Here t 1 , t 2 > 0, is a 2 × 2 diagonal matrix that represents the times t 1 and t 2 when two subparticles are observed at positions x 1 and x 2 in R 3 , respectively; and is a 3 × 2 matrix the columns of which represent the positions x 1 , x 2 ∈ R 3 of two subparticles at times t 1 , t 2 ∈ R, respectively.
Accordingly, the point T X = (T, X) ∈ R 5×2 represents a (2, 3)-particle, which is a multi-particle system consisting of two subparticles (t 1 , x 1 ) and (t 2 , x 2 ) with positions x 1 and x 2 in R 3 at times t 1 and t 2 , respectively. Here we use the displayed notation, T X , and the inline notation, (T, X), interchangeably.
The collective application of the Galilei bi-boost B ∞ (V) of signature (2,3) to the pair of subparticles (T, X) in m + n = 2 + 3 time-space dimensions yields which is described in the following chain of equations, Here, B ∞ (V) in (167) and (168) is given by (157) with m = 2 and n = 3. The chain of equations (168) describes the application of a Galilei bi-boost B ∞ (V) of signature (2,3) to collectively bi-boost two subparticles, (t 1 , x 1 ) and (t 2 , x 2 ), into the two bi-boosted subparticles, (t 1 , x 1 + v 1 t 1 ) and (t 2 , x 2 + v 2 t 2 ), by 3-dimensional velocity vectors v 1 = (v 11 , v 21 , v 31 ) t and v 2 = (v 12 , v 22 , v 32 ) t in R 3 . It is important to note that the two collectively bi-boosted subparticles are not entangled in the sense that the boost of each boosted subparticle is independent of the boost of the other boosted subparticle. Interestingly, this observation fails when we replace Galilei bi-boosts of signature (m, 3), m ≥ 2, by corresponding Lorentz bi-boosts of same signature (m, 3), as we will see in the sequel.
Each of the two particles (t 1 , x 1 ) and (t 2 , x 2 ) possesses a one-dimensional time, t 1 ∈ R and t 2 ∈ R, respectively. Accordingly, the system consisting of the two particles possesses the two-dimensional time, (t 1 , t 2 ) ∈ R 2 . Each of the two particles possesses its own clock, so that the two-dimensional time of the system is measured by two clocks. In general, a multi-particle system consisting of m particles possesses an m-dimensional time, measured by m clocks, m ∈ N.
The extension of (163)-(168) from signature (2, 3) to signature (m, n), m, n ∈ N, is now obvious. The Galilei bi-boost B ∞ (V) of signature (m, n) is parametrized by a velocity matrix V ∈ R n×m of order n × m that consists of m columns, V = (v 1 v 2 . . . v m ), which represent the m velocity vectors v 1 , v 2 , . . . , v m ∈ R n relative to some inertial frame. Furthermore, when B ∞ (V) is applied to collectively bi-boost m particles in R n (i) it keeps invariant each of the times t k , k = 1, . . . , m of the m particles (t k , x k ), that is, t k = t k , and (ii) it bi-boosts their positions x k ∈ R n into the bi-boosted positions x k + v k t k ∈ R n at times t k , respectively. The m collectively bi-boosted particles are not entangled in the sense that (i) the boost of each boosted particle is independent of the boosts and times of the other boosted particles and (ii) the time of each boosted particle is independent of the times and boosts of the other boosted particles.
A Galilei bi-boost of signature (m, 3), applied collectively to the m subparticles of an (m, 3)-particle is thus equivalent to m Galilei boosts applied individually to each subparticle. Hence, a Galilei bi-boost of signature (m, n), m, n ≥ 2, can be viewed as a Galilei multi-boost acting on multi-particle systems. While Galilei multi-boosts involve no entanglement, we will see that Lorentz multi-boosts accommodate entanglement of space and time coordinates of multi-particle systems.
The chain of equations (168) for the action of Galilei bi-boosts of signature (2, 3) and its obvious extension to the action of Galilei bi-boosts of any signature (m, n), m, n ∈ N, demonstrate that the extension of the common Galilei boost of signature (1, 3) to Galilei bi-boosts of any signature (m, n) is quite natural and intuitively clear. The additive decomposition (155) provides a correspondence between Galilei bi-boosts, B ∞ (V), V ∈ R n×m , of signature (m, n), and Lorentz bi-boosts, B c (V), V ∈ R n×m c , of same signature (m, n). This correspondence indicates that the extension of the common Lorentz boost of signature (1, 3) to Lorentz bi-boosts of any signature (m, n) is quite natural as well. Yet, unlike Lorentz bi-boosts, Galilei bi-boosts of signature (m, n) are intuitively clear. We, therefore, explore the Lorentz bi-boosts of any signature (m, n) in Section 21.

Application of the Lorentz Bi-Boost of any Signature
We are now in the position to explore the interpretation of the Lorentz bi-boost of signature (m, n) for all m, n ∈ N, paying special attention to the case when (m, n) = (2, 3) as an illustrative example.
The collective application of the Lorentz bi-boost B c (V) of signature (2,3) to the pair of particles (T, X) in 2 + 3 time-space dimensions is obtained by replacing B ∞ (V) by B c (V) in (168), where B ∞ (V) and B c (V) are related by the Lorentz bi-boost additive decomposition (159). Accordingly, in full analogy with (168) we obtain the chain of equations where T and X are given by (165) and (166), V ∈ R 3×2 c ⊂ R 3×2 is given by (163) and E(V) is given by (158) with (m, n) = (2, 3). Owing to the presence of the factor 1/c 2 in (169), the interaction terms that appear in the entangled time components t 21 and t 12 in (169) are directly noticeable only at very high speeds.
The relativistic squared bi-norm (111) of each subparticle of the (2, 3)-particle (T, X) remains invariant under the application in (169) of the bi-boost B c (V) ∈ SO c (2, 3), that is, where we use the notation x 2 = x·x for vectors x ∈ R n .
Moreover, the relativistic bi-inner product (110) of the two subparticles (t k , x k ), k = 1, 2, of the (2, 3)-particle (T, X) remains invariant under a bi-boost application as well. (T, X) ∈ R (m+n)×m , m, n ∈ N, in the special case when (m, n) = (2, 3). Hence, it is important to note the following immediate generalization of the theorem. The two theorems in [19] (Theorem 6.14) and [19] (Theorem 6.15) are stated for (T, X) ∈ R (m+n)×m for all m, n ∈ N, where T ∈ R m×m and X ∈ R n×m . However, these two theorems and their proofs can be unified under a generalization of the theorem in [19] (Theorem 6.15) to the case when (T, X) ∈ R (m+n)×k for all m, n, k ∈ N, so that T ∈ R m×k and X ∈ R n×k . Then, each invariance in (170) corresponds to (m, n, k) = (2, 3, 1), and the invariance in (171) corresponds to (m, n, k) = (2, 3, 2).
According to (169), the collective application of a Lorentz bi-boost B c (V) of signature (m, n) = (2, 3) to the constituents of a system of two 3-dimensional subparticles generates classical effects and bi-relativistic effects.
The classical effects boost the two subparticles of the system individually by velocity vectors v 1 and v 2 , as we see from the first term on the extreme right-hand side of (169), which is the application of the Galilean component of the Lorentz bi-boost of signature (2,3).
The bi-relativistic effects involve collective entanglement of the times and positions of the system constituents, as we see from the second term on the extreme right-hand side of (169), which is the application of the relativistic entanglement component of the Lorentz bi-boost of signature (2,3).
The extension of (169) to the application of Lorentz bi-boosts of any signature (m, n), m, n ∈ N, to multi-particle systems of m n-dimensional subparticles is now clear. It suggests that Lorentz bi-boosts of signature (m, n) form the symmetry group by which multi-particle systems of m n-dimensional moving entangled particles can be understood, just as the Lorentz group of signature (1, 3) forms the symmetry group by which Einstein's special relativity can be understood.

Lorentz Bi-boost Composition Law
The special relativistic Lorentz bi-boost composition law for signature (1, n) is presented in (103). In full analogy, the Lorentz bi-boost composition law for any signature (m, n), n, m ∈ N, is given by ([19] Theorem 6.7) for all V 1 , V 2 ∈ R n×m c .
Here ⊕ E = ⊕ E,(1,n),c is the Einstein addition of signature (m, n) in the ball R n×m c , given by each of the following mutually equivalent equations (Equations 6.54 [19]): V 1 , V 2 ∈ R n×m c . In bi-gyrolanguage, the Einstein addition of signature (m, n) is said to be a bi-gyroaddition. In the special case when m = 1 Einstein bi-gyroaddition (174) descends to the Einstein gyroaddition (2) of special relativity, as shown in ([19] Section 5.17). Not unexpectedly, Einstein bi-gyroaddition (174) in R n×m c obeys the so called bi-gyrotriangle inequality ( [19] Theorem 5.108), V 1 , V 2 ∈ R n×m c , which descends to the gyrotriangle inequality (14) when m = 1. The Lorentz bi-boost composition law (172) involves two kinds of gyrations generated by We say that the multi-particle system (182) is integrated (as opposed to disintegrated) if its constituent particles have equal positions at equal times, that is, if We may consider the integrated system (T, X) in (185) as a single particle capable of being disintegrated into three subparticles. To disintegrate the particle (T, X) we replace V in (183) by V 123 = (v 1 , v 2 , v 3 ) where v 1 , v 2 , v 3 ∈ R 2 are distinct velocity vectors. The bi-boost B c (V 123 ) takes the integrated (3, 2)-particle (T, X) in (185) into a disintegrated multi-particle system (T , X ) = B c (V 123 )(T, X) of three entangled subparticles that move with distinct relative velocities. In this case we say that the bi-boost B c (V 123 ) disintegrates the integrated particle (T, X), resulting in the disintegrated multi-particle system (T , X ) of entangled subparticles.
A multi-particle system is partially integrated if several of its subsystems are integrated. Clearly, if an inertial observer observes that a given multi-particle system is partially integrated in some sense, then all inertial observers observe that the system is partially integrated in the same sense. Hence, we naturally raise the following property (Property 1) of Lorentz bi-boosts that are parametrized by equi-column velocity matrices. Property 1. The property of a multi-particle system to be partially integrated is observer-independent.
Formally, Property 1 asserts that any bi-boost B c (V) parametrized by an equi-column velocity matrix V ∈ R n×m c preserves partially integrated structures of multi-particle systems. Clearly, Property 1 must hold true if the Lorentz bi-boost B c (V) of signature (m, n) possesses physical significance as the symmetry group of multi-particle systems.
Property 1 implies, for instance, that the system (T, X) in (184) is integrated if and only if the bi-boosted system (T , X ) in (184) is integrated.
Property 1 is corroborated numerically in ( [19] Chap. 6) by several examples. However, we do not have a proof of Property 1. Property 1, if it holds true, supports the claim that the Lorentz group SO c (m, n) is the symmetry group of the time and space coordinates of any multi-particle system of m n-dimensional particles, for any m, n ∈ N.
It is, therefore, hoped that this review article will stimulate the search for experimental evidence that establishes the physical significance of Lorentz groups SO c (m, 3) for all m ≥ 2, just as Einstein's special relativity establishes the physical significance of the common Lorentz group SO c (1, 3).

Conclusions
Einstein addition of relativistically admissible velocities, presented in Section 2, is neither commutative nor associative, and as such, it is not a group operation. However, it turns out in Section 3 that, being a gyrocommutative gyrogroup operation, Einstein addition possesses a rich algebraic structure called a gyrocommutative gyrogroup. In particular, it turns out that the departure of Einstein addition from commutativity and associativity is strictly controlled by gyrations, which are automorphisms studied in Section 4. In Sections 5-12, we therefore pave the road from Einstein addition to gyrogroups, and review elements of gyrogroup theory in order to demonstrate the power and elegance that gyrogroup theory shares with group theory.
In particular, we pay attention to the elegant gyro-identities (18) that Einstein addition and its gyrations obey. These remarkable gyro-identities survive unimpaired in (181) in their transition from Einstein addition (2) of relativistically admissible velocity vectors in the ball R n c to Einstein addition (174) of bi-relativistically admissible velocity matrices in the ball R n×m c . Our plan to extend the Lorentz boost from signature (1,3) in special relativity to signature (m, n) for all m, n ∈ N, is guided by analogies with the Galilei boost, the analogous extension of which is clear and intuitive. Therefore, in Section 13 we present both Galilei and Lorentz boosts, paying special attention to the additive decomposition (105) that expresses the Lorentz boost as the sum of two components: (i) a Galilei boost and (ii) a relativistic component, which is directly noticeable only in high speeds. We indicate in Section 13 that Galilei and Lorentz boosts, which act on particles, are to be extended to Galilei and Lorentz multi-boosts, which act on multi-particle systems.