Next Article in Journal
Comparison Methods for Solving Non-Linear Sturm–Liouville Eigenvalues Problems
Next Article in Special Issue
A Spacetime Symmetry Approach to Relativistic Quantum Multi-Particle Entanglement
Previous Article in Journal
Anatomical Considerations and Study of the Fractal Dimension around the Posterior Superior Alveolar Artery
Previous Article in Special Issue
Ordered Gyrovector Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Binary Operations in the Unit Ball: A Differential Geometry Approach

by
Nikita E. Barabanov
and
Abraham A. Ungar
*
Department of Mathematics, North Dakota State University, Fargo, ND 58108, USA
*
Author to whom correspondence should be addressed.
Symmetry 2020, 12(7), 1178; https://doi.org/10.3390/sym12071178
Submission received: 7 May 2020 / Revised: 11 June 2020 / Accepted: 1 July 2020 / Published: 16 July 2020
(This article belongs to the Special Issue Symmetry and Geometry in Physics)

Abstract

:
Within the framework of differential geometry, we study binary operations in the open, unit ball of the Euclidean n-space R n , n N , and discover the properties that qualify these operations to the title addition despite the fact that, in general, these binary operations are neither commutative nor associative. The binary operation of the Beltrami-Klein ball model of hyperbolic geometry, known as Einstein addition, and the binary operation of the Beltrami-Poincaré ball model of hyperbolic geometry, known as Möbius addition, determine corresponding metric tensors in the unit ball. For a variety of metric tensors, including these two, we show how binary operations can be recovered from metric tensors. We define corresponding scalar multiplications, which give rise to gyrovector spaces, and to norms in these spaces. We introduce a large set of binary operations that are algebraically equivalent to Einstein addition and satisfy a number of nice properties of this addition. For such operations we define sets of gyrolines and co-gyrolines. The sets of co-gyrolines are sets of geodesics of Riemannian manifolds with zero Gaussian curvatures. We also obtain a special binary operation in the ball, which is isomorphic to the Euclidean addition in the Euclidean n-space.

1. Introduction

Let B be the unit, open ball in the Euclidean n-space R n ,
B = { x R n : x < 1 } ,
n N , where · is the Euclidean norm.
Einstein addition is a binary operation, E , in the ball B R n that stems from his velocity composition law in the ball B R 3 of relativistically admissible velocities. Seemingly structureless, Einstein addition is neither commutative nor associative. However, Einstein addition turns out to be both gyrocommutative and gyroassociative, thus giving rise to the rich algebraic structures that became known as a (gyrocommutative) gyrogroup and a gyrovector space, the definitions of which are presented in Definitions 1–3, Section 2.
Einstein addition, E , and its isomorphic copy, Möbius addition, M , are studied in the literature algebraically, along with applications to the hyperbolic geometry of Lobachevsky and Bolyai. Naturally, one may expect that the rich algebraic structure of Einstein addition can find home in differential geometry, giving rise to a novel branch called Binary Operations in the Ball.
Accordingly, the aim of this paper is to develop a differential geometry approach to Einstein addition and, hence, to discover the resulting novel branch of differential geometry that involves binary operations in the ball. We thus begin with the study of an arbitrary binary operation in the ball that satisfies some general conditions.
A binary operation in B is a function f: B × B B . We consider functions f of class C 2 , that is, functions f having continuous second derivatives. This operation determines a metric tensor G ( x ) in B given by
g ( x ) = f ( x , y ) y | y = x , G ( x ) = g ( x ) g ( x ) R n × n ,
x B , where denotes transposition, and R n × n is the space of all real n × n -matrices. Then, B is a Riemannian manifold with a metric tensor G.
We pay special attention to the following three binary operations in the ball, along with their associated scalar multiplication:
  • Einstein addition E in the ball, presented in (136), and the scalar multiplication E that it admits, presented in (138), are recovered in Section 5 within the framework of differential geometry. The triple ( B , E , E ) is an Einstein gyrovector space that forms the algebraic setting for the Beltrami-Klein ball model of hyperbolic geometry.
  • Möbius addition M in the ball, presented in (153), and the scalar multiplication M that it admits, presented in (154), are recovered in Section 6 within the framework of differential geometry. The triple ( B , M , M ) is a Möbius gyrovector space that forms the algebraic setting for the Beltrami–Poincaré ball model of hyperbolic geometry.
  • A novel, interesting binary operation ⊕ in the ball, presented in (199), and the scalar multiplication ⊗ that it admits, presented in (205), are discovered in Section 7 within the framework of differential geometry. Remarkably, the triple ( B , , ) is a vector space isomorphic to the Euclidean vector space ( R n , + , · ) . As such, the binary operation ⊕ in B is commutative, associative and distributive. Accordingly, the vector space ( B , , ) forms the algebraic setting for a novel n-dimensional Euclidean geometry ball model.
A procedure that we present in this paper enables binary operations in the ball, like E , M and ⊕, to be obtained from corresponding metric tensors G. Interestingly, a metric tensor G is determined by the behavior of the function f in a neighborhood of the set D = { ( x , x ) : x B } , rather than in the whole of the space B × B . Hence, the global operations turn out to be determined by local properties of the functions f and their second derivatives in the set D .
The procedure is formulated in terms of geodesics and a parallel transport. It is applied to a wide class of metric tensors that satisfy three properties: (i) smoothness, (ii) rotation invariance, and (iii) plane invariance.
For Einstein addition E the operation of scalar multiplication t E a with a B , t R is well defined [1]. It gives rise to a structure called a gyrovector space. For each metric tensor considered in this paper we define an operation t a of scalar multiplication, which leads to corresponding gyrovector spaces.
For every a , b B the set { f ( b , t a ) : t R } is called a gyroline. The set of gyrolines defines a metric tensor G. The sets { f ( t a , b ) : t R } are called co-gyrolines. It is proved in the paper that for operations f isomorphic to Einstein addition, the co-gyrolines are gyrolines for other binary operations, f c o , with metric tensors that we denote by G c o . The Gaussian curvatures of spaces with metric tensors G c o are equal to zero.
Various algebraic and geometric properties of the operations E and M have been intensively studied in recent papers and monographs; see, for instance, References [1,2,3,4,5,6,7,8,9,10]).
Einstein addition and Möbius addition are neither commutative nor associative. As such, they do not form group operations. Yet, the lack of the commutative and associative laws is compensated by the gyrocommutative and gyroassociative laws that these binary operations obey. As such, these binary operations give rise to the algebraic objects known as gyrogroups and gyrovector spaces. Remarkably, gyrovector spaces form the algebraic setting for various models of analytic hyperbolic geometry, just as the standard vector spaces form the algebraic setting for analytic Euclidean geometry; see, for instance, References [4,5,9,11,12,13].
The special interest of our study of both Einstein addition and Möbius addition within the framework of differential geometry stems from the result that they are gyrocommutative gyrogroup operations. Indeed they give rise to (i) Einstein gyrogroups ( B , E ) and Einstein gyrovector spaces ( B , E , E ) ; and to (ii) Möbius gyrogroups ( B , M ) and Möbius gyrovector spaces ( B , M , M ) . The definitions of the abstract (gyrocommutative) gyrogroup and gyrovector space are presented in Section 2.
The organization of the paper is as follows. In Section 2 we present the definitions of the abstract (gyrocommutative) gyrogroup and gyrovector space. In Section 3 we introduce metric tensors satisfying three properties: (i) smoothness; (ii) rotation invariance; and (iii) plane invariance. Then we derive the equations for geodesics and parallel transport. In Section 4 we (i) introduce binary operations defined by metric tensors G; (ii) define corresponding operations of scalar multiplication; and (iii) introduce distances and gyronorms. In Section 5 we prove that the operation f, defined in terms of the metric tensor of Einstein addition, coincides with Einstein addition. In Section 6 and Section 7 similar results are obtained for Möbius addition and for a novel addition. Section 8 is devoted to properties of binary operations similar to those of Einstein addition: (i) left cancellation law; (ii) existence and properties of unitary gyrators; and (iii) a gyrocommutative law.
In the paper, we use the following notation— R is the set of real numbers. I is the identity matrix. A real square matrix U is called unitary [14] if U U = I . The function sign is the signum function in R . Components of an n-vector x are denoted by x j , j = 1 , , n .

2. Gyrogroups and Gyrovector Spaces

Definition 1.
((Gyrogroup), [4] Definition 2.5). A groupoid ( S , ) is a nonempty set S with a binary operation. A groupoid ( S , ) is a gyrogroup if its binary operation satisfies the following axioms. In S there is at least one element, 0, called a left identity, satisfying
(G1)    0 a = a ,
for all a S . There is an element 0 S satisfying axiom ( G 1 ) such that for each a S there is an element a S , called a left inverse of a, satisfying
(G2)    a a = 0 .
Moreover, for any a , b , c S there exists a unique element gyr [ a , b ] c S such that the binary operation obeys the left gyroassociative law
(G3)    a ( b c ) = ( a b ) gyr [ a , b ] c .
The map gyr [ a , b ] : S S given by c gyr [ a , b ] c is an automorphism of the groupoid ( S , ) , that is,
(G4)    gyr [ a , b ] Aut ( S , ) ,
and the automorphism gyr [ a , b ] of S is called the gyroautomorphism, or the gyration, of S generated by a , b S . The operator gyr : S × S Aut ( S , ) is called the gyrator of S. Finally, the gyroautomorphism gyr [ a , b ] generated by any a , b S possesses the left reduction property
(G5)    gyr [ a , b ] = gyr [ a b , b ] ,
called the reduction axiom.
The gyrogroup axioms ( G 1 )–( G 5 ) in Definition 1 split up into three classes:
  • The first pair of axioms, ( G 1 ) and ( G 2 ) , is a reminiscent of the group axioms.
  • The last pair of axioms, ( G 4 ) and ( G 5 ) , presents the gyrator axioms.
  • The middle axiom, ( G 3 ) , is a hybrid axiom linking the two pairs of axioms in (1) and (2).
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.
Definition 2.
((Gyrocommutative Gyrogroup), [4] Definition 2.6). A gyrogroup ( S , ) is gyrocommutative if its binary operation obeys the gyrocommutative law
(G6)    a b = gyr [ a , b ] ( b a ) ,
for all a , b S .
The first concrete example of a gyrogroup was discovered in 1988 [15]. It became known as an Einstein gyrogroup. Einstein gyrogroups are employed, for instance, in References [5,16,17,18,19,20,21]. Möbius gyrogroups are employed as well, for instance, in References [4,22,23,24]. In full analogy with groups, there are topological gyrogroups, studied in References [25,26]. Gyrogroups share remarkable analogies with groups studied, for instance, in References [27,28,29,30,31,32]. Applications of gyrogroups in harmonic analysis are found in References [17,23,33]. For other interesting studies of gyrogroups see References [34,35,36,37,38,39]. Einstein gyrogroups and gyrovector spaces are extended to Einstein bi-gyrogroups and bi-gyrovector spaces in References [10,40], along with an application to relativistic quantum entanglement of multi-particle systems.
The gyrocommutative gyrogroups that we study in this paper admit scalar multiplication, turning themselves into gyrovector spaces, the formal definition of which follows.
Definition 3.
((Gyrovector space), [4] Definition 6.2). A real inner product gyrovector space ( G , , ) (gyrovector space in short) is a gyrocommutative gyrogroup ( G , ) that obeys the following axioms:
(1) G is a subset of a real inner product vector space V called the carrier of G, G V , from which it inherits its inner product, ·, and norm, · , which are invariant under gyroautomorphisms, that is,
(V0)    g y r [ u , v ] a · g y r [ u , v ] b = a · b for all a , b , u , v G .
(2) G admits a scalar multiplication, , possessing the following properties. For all real numbers r , r 1 , r 2 R and all points a G :
(V1)    1 a = a ,
(V2)    ( r 1 + r 2 ) a = r 1 a r 2 a ,
(V3)    ( r 1 r 2 ) a = r 1 ( r 2 a ) ,
(V4)    | r | a r a = a a , a 0 , r 0 ,
(V5)    g y r [ u , v ] ( r a ) = r g y r [ u , v ] a ,
(V6)    g y r [ r 1 v , r 2 v ] = I ,
(V7)    r a = | r | a ,
(V8)    a b a b .
Like gyrogroups, also gyrovector spaces are studied in the literature. The papers [41,42,43,44,45] are devoted to various aspects of Möbius gyrovector spaces. Einstein and Möbius gyrovector spaces in the context of a gyrovector space approach to hyperbolic geometry are the subject of Reference [6]. Generalized gyrovector spaces are studied in References [46,47,48]. Interesting results about the differential geometry of some gyrovector spaces may be found in Reference [49]. Other interesting studies of gyrovector spaces are found in References [50,51,52,53,54,55].
In this paper we introduce, within the framework of differential geometry, a large number of gyrocommutative gyrogroup operations, which enjoy key properties of Einstein addition and Möbius addition. Furthermore, we present corresponding scalar multiplications that turn these gyrocommutative gyrogroups into gyrovector spaces.

3. Metric Tensors

In this section we consider metric tensors that satisfy certain conditions, and derive equations for geodesics and parallel transports for these tensors.

3.1. Parametrization of Metric Tensors

For every smooth binary operation f C 1 : B × B B ,
f ( a , b ) = a b ,
we determine a matrix function G R n × n by the equations
g ( x ) = y [ ( x ) y ] | y = x G ( x ) = g ( x ) g ( x ) ,
for all x B .
We consider three conditions on the tensors G. These are—(i) rotation invariance; (ii) two-dimensional plane invariance; and (iii) smoothness.
Rotation invariance. For every unitary n × n -matrix U and x B we have
G ( U x ) = U G ( x ) U .
Plane invariance. For every a , b B and every two-dimensional plane that contains a and b, the whole geodesic x ( · ) with initial conditions x ( 0 ) = a and x ˙ ( 0 ) = b belongs to this plane.
Smoothness. The function G is differentiable in B .
The first property means invariance with respect to rotations. The second property means that every geodesic belongs to a plane that contains its initial data. In particular, if a and b belong to the same line l, then f ( a , b ) l . The third property is a standard assumption that allows to derive the differential equations for geodesics.
These three properties are valid if there exist differentiable functions m 0 , m 1 : [ 0 , 1 ) ( 0 , ) such that m 0 ( 0 ) = m 1 ( 0 ) , m 0 ( s ) > 0 , m 1 ( s ) > 0 for all s [ 0 , ) , and
G ( x ) = { g i j ( x ) } i , j = 1 n = m 0 ( x 2 ) ( I x x x 2 ) + m 1 ( x 2 ) x x x 2 .
The inverse matrix G 1 ( x ) of G ( x ) in (6) exists, being given by
G 1 ( x ) = { g j s ( x ) } j , s = 1 n = 1 m 0 ( x 2 ) ( I x x x 2 ) + 1 m 1 ( x 2 ) x x x 2 .
We say that the metric tensor G ( x ) in (6) is parametrized by the two functions m 0 and m 1 . Special attention will be paid to two parametrizations of G ( x ) .
The first parametrization of G ( x ) is when m 0 and m 1 are given by
m 0 ( r 2 ) = 1 1 r 2 , m 1 ( r 2 ) = 1 ( 1 r 2 ) 2 , r [ 0 , 1 ) ,
which will be associated with Einstein addition. Einstein addition, in turn, is the binary operation that results from Einstein’s addition law of relativistically admissible velocities in special relativity theory [1].
The second parametrization of G ( x ) is when m 0 and m 1 are given by
m 0 ( r 2 ) = 1 ( 1 r 2 ) 2 , m 1 ( r 2 ) = 1 ( 1 r 2 ) 2 , r [ 0 , 1 ) ,
which will be associated with Möbius addition. Möbius addition, in turn, is a natural generalization to the ball of a well known Möbius transformation of the disc [4].
In order to construct binary operations that result from metric tensors G we will find the differential equations for geodesics and parallel transports. To this end we have to calculate the Christoffel symbols in Section 3.3.

3.2. Geodesics

Assume x: [ 0 , 1 ] B is a curve in B . The length of this curve is
l = 0 1 x ˙ ( s ) G ( x ( s ) ) x ˙ ( s ) d s .
Denote G ( x ) = { g i j ( x ) } i , j = 1 n , a = x ( 0 ) , b = x ( 1 ) , and
F ( x ( s ) , x ˙ ( s ) ) = x ˙ ( s ) G ( x ( s ) ) x ˙ ( s ) .
Assume x ( · ) minimizes l over all smooth curves x such that x ( 0 ) = a , x ( 1 ) = b . Use a new parametrization
d t = F ( x ( s ) , d x ( s ) d s ) d s 0 1 F ( x ( u ) , x ˙ ( u ) ) d u .
Then, F ( x ( t ) , d x d t ( t ) ) c o n s t .
According to Euler-Lagrange equations for all k = 1 , , n we have
F x k d d t F x ˙ k = 0 .
This equation takes the form
i = 1 n j = 1 n g i j ( x ) x k d x i ( t ) d t d x j ( t ) d t d d t j = 1 n g k j ( x ) d x j ( t ) d t = 0 .
Denote the entries of the matrix G 1 ( x ) by g k p ( x ) . Then (14) is equivalent to the following. For all p = 1 , , n
x ¨ p + 1 2 k = 1 n i = 1 n j = 1 n { g k j ( x ) x i + g k i ( x ) x j g i j ( x ) x k } g k p ( x ) x ˙ i x ˙ j = 0 .
Denote
Γ i j p ( x ) = 1 2 k = 1 n { g k j ( x ) x i + g k i ( x ) x j g i j ( x ) x k } g k p ( x ) .
Then (15) and (16) yield the well known equations of geodesics [56],
x ¨ p + i = 1 n j = 1 n Γ i j p ( x ) x ˙ i x ˙ j = 0 .
Note that along every geodesic the function
F ( x ( t ) , x ˙ ( t ) ) = i = 1 n j = 1 n g i j ( x ( t ) ) x ˙ i ( t ) x ˙ j ( t )
is constant.

3.3. Christoffel Symbols

For the sake of brevity we omit the arguments of the functions m j : m j = m j ( x 2 ) , and of the derivatives of these functions,
m j ( r 2 ) = d m j ( z ) d z | z = r 2 .
We also use the Kronecker delta δ i j which is equal to one if i = j , and to zero otherwise. The elements g i j ( x ) of the metric tensor G ( x ) in (6) have the form
g i j ( x ) = m 0 ( δ i j x i x j x 2 ) + m 1 x i x j x 2 .
Then
m j x k = m j 2 x k x k x i x j x 2 = 2 x i x j x k x 4 + δ i k x j + δ j k x i x 2 .
The partial derivatives of g i j ( x ) are given by
g i j ( x ) x k = m 0 2 x k ( δ i j x i x j x 2 ) + m 1 2 x k x i x j x 2 + ( m 1 m 0 ) ( δ i k x j + δ j k x i x 2 2 x i x j x k x 4 ) .
Now we introduce the functions Γ i j k ( x ) ,
Γ i j k ( x ) = 1 2 ( g i k x j + g j k x i g i j x k ) = m 0 ( δ i k x j + δ j k x i δ i j x k x i x j x k x 2 ) + m 1 x i x j x k x 2 + ( m 1 m 0 ) ( δ i j x k x 2 x i x j x k x 2 ) .
Therefore, Christoffel symbols Γ i j s ( x ) are given by
Γ i j s ( x ) = k = 1 n Γ i j k ( x ) g k s ( x ) = m 1 m 1 x i x j x s x 2 + m 0 m 0 ( δ i s x j + δ j s x i 2 x i x j x s x 2 ) + m 1 m 0 x 2 m 0 x 2 m 1 ( δ i j x s x i x j x s x 2 ) .
Having the Christoffel symbols Γ i j s ( x ) , we are in the position to calculate the equation of geodesics in Section 3.4.

3.4. Geodesics

The well known differential equation of geodesics is (17). Substituting the values of Christoffel symbols Γ i j s ( x ) from (24) into (17) yields
x ¨ + m 0 m 0 2 ( x x ˙ ) x ˙ + [ ( m 1 m 1 2 m 0 m 0 ) ( x x ˙ ) 2 x 2 + m 1 m 0 x 2 m 0 x 2 m 1 ( x ˙ 2 ( x x ˙ ) 2 x 2 ) ] x = 0 .
In the special case when the initial values x ( 0 ) and x ˙ ( 0 ) are parallel, that is, x ( 0 ) = λ 1 a and x ˙ ( 0 ) = λ 2 a for some vector a and numbers λ 1 , λ 2 , the solution of (25) has the form
x ( t ) = r ( t ) a .
Here, the scalar function r, to be determined, satisfies the initial conditions r ( 0 ) = λ 1 and r ˙ ( 0 ) = λ 2 , along with the differential equation
r ¨ + m 0 m 0 2 a 2 r r ˙ 2 + ( m 1 m 1 2 m 0 m 0 ) a 2 r r ˙ 2 = 0
or, equivalently,
r ¨ + m 1 m 1 a 2 r r ˙ 2 = 0 .
Noting the notation m 1 = m 1 ( x 2 ) = m 1 ( a 2 r ( t ) 2 ) that we use, (28) can be integrated, obtaining
ln | r ˙ ( t ) | + 1 2 ln | m 1 ( a 2 r ( t ) 2 ) | = C 1 ,
where C 1 is an arbitrary constant.
The resulting first order Equation (29) can be written as
r ˙ ( t ) = C m 1 ( a 2 r ( t ) 2 ) ,
where
C = r ˙ ( 0 ) m 1 ( x ( 0 ) 2 ) .
Equation (30) is separable, giving rise to the equation
m 1 ( a 2 r 2 ) d r = C d t .
Integrating (32) yields
r ( 0 ) r ( 1 ) m 1 ( a 2 s 2 ) d s = C ,
so that, finally, by (31) and (33),
r ˙ ( 0 ) = r ( 0 ) r ( 1 ) m 1 ( a 2 s 2 ) m 1 ( a 2 r ( 0 ) 2 ) d s .
Let us now consider the important special case when r ( 0 ) = 0 and r ( 1 ) = 1 (this means x ( 0 ) = 0 and x ( 1 ) = a ). In this special case (34) reduces to
r ˙ ( 0 ) = 0 1 m 1 ( a 2 s 2 ) m 1 ( 0 ) d s = 1 a 0 a m 1 ( s 2 ) m 1 ( 0 ) d s .
Formula (35) will prove useful in the definition of special binary operations and in the introduction of norms associated with these operations.

3.5. Parallel Transport

The well known differential equation for a parallel transport of a vector X = { X s } s = 1 n along a curve x is [56]
X ˙ s + i = 1 n j = 1 n Γ i j s x ˙ i X j = 0 , s = 1 , , n .
Substituting the values of Christoffel symbols Γ i j s ( x ) from (24) into (36) yields
X ˙ + m 1 m 1 ( x x ˙ ) ( x X ) x 2 x + { m 0 m 0 [ ( x X ) ( I x x x 2 ) x ˙ + ( x x ˙ ) ( I x x x 2 ) X ] + m 1 m 0 x 2 m 0 x 2 m 1 x ˙ ( I x x x 2 ) X } x = 0 .
If the curve x is linear, joining the origin and a point a R n , then x ( t ) = a t , and (37) takes the form
X ˙ + m 1 m 1 t a X a + m 0 m 0 [ t a X a + t a 2 X ] = 0 .
With the initial condition X ( 0 ) = X 0 , (38) has the solution
X ( t ) = m 0 ( x ( 0 ) 2 ) m 0 ( x ( t ) 2 ) ( I a a a 2 ) X 0 + m 1 ( x ( 0 ) 2 ) m 1 ( x ( t ) 2 ) a a a 2 X 0 .
In particular, at t = 1 , we have
X ( 1 ) = m 0 ( 0 ) m 0 ( a 2 ) ( I a a a 2 ) X 0 + m 1 ( 0 ) m 1 ( a 2 ) a a a 2 X 0 .
If the vectors X 0 and a are not parallel, then the first term in (39) describes the evolution of X ( t ) in a direction orthogonal to a in the plane determined by the vectors X 0 and a, and the second term in (39) describes the evolution of X ( t ) in a direction parallel to a.
If we multiply the functions m 0 and m 1 by the same positive number, then equation (25) of geodesics and (37) of parallel transport remain unchanged. Hence, without loss of generality we may assume that m 0 ( 0 ) = m 1 ( 0 ) = 1 .
Assumption 1.
The functions m 0 and m 1 satisfy m 0 ( 0 ) = m 1 ( 0 ) = 1 .

4. Binary Operations

We are now in the position to define a binary operation ⊕ in B that results from the metric tensor G in (6).

4.1. Vector Addition

In this subsection we define a binary operation ⊕ in the ball B . Let vectors a , b B be given. If b = 0 , then a b = a . If a = 0 , then a b = b . For the case b 0 , a 0 we perform the following four steps that lead to the definition of a b .
Step 1. We calculate the value r ˙ ( 0 ) using Formula (35) with a replaced by b,
r ˙ ( 0 ) = 1 b 0 b m 1 ( s 2 ) d s .
Step 2. We calculate the value X ( 1 ) using Formula (40) with initial condition X 0 = r ˙ ( 0 ) b ,
X ( 1 ) = 1 m 0 ( a 2 ) ( I a a a 2 ) X 0 + 1 m 1 ( a 2 ) a a a 2 X 0 .
Step 3. We solve the second order differential Equation (25),
x ¨ + m 0 m 0 2 ( x x ˙ ) x ˙ + [ ( m 1 m 1 2 m 0 m 0 ) ( x x ˙ ) 2 x 2 + m 1 m 0 x 2 m 0 x 2 m 1 ( x ˙ 2 ( x x ˙ ) 2 x 2 ) ] x = 0
with initial conditions x ( 0 ) = a , x ˙ ( 0 ) = X ( 1 ) .
Step 4. A binary operation a b is defined by the equation
a b = x ( 1 ) .
We say that the binary operation ⊕ in B is generated by the metric tensor G given by (2).

4.2. Elementary Properties of the Binary Operations

Lemma 1.
For every vector v B
0 v = v ,
( v ) v = 0 .
Proof. 
The first equation stems from the definition of ⊕. Assume a = v , b = v , v 0 . Applying steps 1–4, we have
X 1 = 1 m 1 ( v 2 ) v .
Consider a solution x of Equation (43) with initial conditions x ( 0 ) = v , x ˙ ( 0 ) = X 1 . Since the vectors x ( 0 ) , x ˙ ( 0 ) are parallel to v, the solution x has the form x ( t ) = g ( t ) v , where g is a scalar function to be determined. Equation (43) shows that g satisfies the conditions
g ¨ + m 1 m 1 ( g ˙ ) 2 g v 2 = 0 , g ( 0 ) = 1 , g ˙ ( 0 ) = 1 m 1 ( v 2 ) 0 1 m 1 ( s 2 v 2 ) d s .
The solution g of this initial value problem satisfies the equation
1 g ( t ) m 1 ( s 2 v 2 ) d s = t 0 1 m 1 ( s 2 v 2 ) d s .
Obviously, g ( 1 ) = 0 . Hence, ( a ) a = x ( 1 ) = 0 . □

4.3. Metric Tensors Associated with Binary Operations

In this section we show that the metric tensor G that generates the operation ⊕ can be recovered from the operation ⊕.
Theorem 1.
Let G be given by (6) and let be the binary operation (44) generated by G. Furthermore, for every x 0 B let g ( x 0 ) be the n × n -matrix such that
( x 0 ) ( x 0 + Δ x ) = g ( x 0 ) Δ x + o ( Δ x ) .
Then
g ( x 0 ) g ( x 0 ) = G ( x 0 ) = m 0 ( x 0 2 ) I x 0 x 0 x 0 2 + m 1 ( x 0 2 ) x 0 x 0 x 0 2 .
Proof. 
The matrix
g ( x 0 ) = [ ( x 0 ) y ] y | y = x 0
exists since the functions m 0 and m 1 are smooth. It is sufficient to prove that
g ( x 0 ) Δ x = m 0 ( x 0 2 ) Δ x
when x 0 Δ x = 0 , and
g ( x 0 ) Δ x = m 1 ( x 0 2 ) Δ x
when x 0 Δ x .
Let’s assume that the vectors x 0 and Δ x are parallel, that is, for some number δ we have Δ x = δ x 0 . We use the procedure described in Section 4.1 with a = x 0 and b = ( 1 + δ ) x 0 . Then, following (40),
X ( 1 ) = 0 ( 1 + δ ) x 0 m 1 ( s 2 ) m 1 ( x 0 2 ) d s x 0 x 0 .
The vectors x ( t ) and x ˙ ( t ) are parallel. Hence (25) takes the form
x ¨ + m 1 m 1 x ˙ 2 x = 0
with initial conditions x ( 0 ) = x 0 , x ˙ ( 0 ) = X ( 1 ) . The unique solution of this initial value problem has the form x ( t ) = q ( t ) x 0 , where the scalar function q satisfies the second order differential equation
q ¨ + q ˙ 1 2 d d t ln m 1 ( x 0 2 q 2 ) = 0
and the initial conditions
q ( 0 ) = 1 q ˙ ( 0 ) = 0 1 + δ m 1 ( x 0 2 s 2 ) m 1 ( x 0 2 ) d s .
Integrating Equation (57) yields
q ˙ ( t ) = q ˙ ( 0 ) m 1 ( x 0 2 ) m 1 ( x 0 2 q ( t ) 2 ) .
Equation (59) is separable. Integrating it over the interval t [ 0 , 1 ] yields
1 q ( 1 ) m 1 ( x 0 2 q 2 ) m 1 ( x 0 2 ) d q = q ˙ ( 0 ) = 0 1 + δ m 1 ( x 0 2 s 2 ) m 1 ( x 0 2 ) d s .
Therefore
q ( 1 ) 1 m 1 ( x 0 2 q 2 ) d q = 0 1 + δ m 1 ( x 0 2 s 2 ) d s
and
q ( 1 ) 0 m 1 ( x 0 2 q 2 ) d q = 1 1 + δ m 1 ( x 0 2 s 2 ) d s .
Since m 1 ( 0 ) = 1 , we get
q ( 1 ) = δ m 1 ( x 0 2 ) + o ( δ ) .
Thus, for Δ x = δ x 0 , and x ( 1 ) = ( x 0 ) ( x 0 + Δ x ) = q ( 1 ) x 0 we have
x ( 1 ) = m 1 ( x 0 2 ) Δ x + o ( Δ x ) ,
so that (53) is proved.
We now assume that x 0 Δ x = 0 . Then x 0 + Δ x x 0 = o ( Δ x ) , and
X ( 1 ) = Δ x x 0 1 m 0 ( x 0 2 ) 0 x 0 m 1 ( s 2 ) d s + x 0 x 0 0 x 0 m 1 ( s 2 ) m 1 ( x 0 2 ) d s + o ( Δ x ) .
Let us consider the solution x of (25) with initial conditions x ( 0 ) = x 0 and x ˙ ( 0 ) = X ( 1 ) . For every t 0 the vector x ( t ) belongs to a two-dimensional plane that contains x 0 and X ( 1 ) . We introduce an orthonormal basis with the first unit vector e 1 = x 0 x 0 , and the second unit vector e 2 = Δ x Δ x , and use the notation x ( t ) = e 1 x ( t ) and x ( t ) = e 2 x ( t ) .
Then x ( t ) = x ( t ) e 1 + x ( t ) e 2 for all t 0 . Owing to (25), the functions x and x satisfy the same second order differential equation
y ¨ + h 1 ( x , x ˙ ) y ˙ + h 2 ( x , x ˙ ) y = 0
with initial conditions
x ( 0 ) = x 0 , x ˙ ( 0 ) = 0 x 0 m 1 ( s 2 ) m 1 ( x 0 2 ) d s + o ( Δ x ) ,
x ( 0 ) = 0 , x ˙ ( 0 ) = Δ x x 0 1 m 0 ( x 0 2 ) 0 x 0 m 1 ( s 2 ) d s + o ( Δ x ) ,
where
h 1 ( x , x ˙ ) = m 0 ( z ) m 0 ( z ) | z = x 2 2 ( x x ˙ ) , h 2 ( x , x ˙ ) = [ ( m 1 ( z ) m 1 ( z ) 2 m 0 ( z ) m 0 ( z ) ) ( x x ˙ ) 2 x 2 + m 1 ( z ) m 0 ( z ) z m 0 ( z ) z m 1 ( z ) ( x ˙ 2 ( x x ˙ ) 2 x 2 ) ] | z = x 2 .
Since the functions h 1 ( x ( · ) , x ˙ ( · ) ) and h 2 ( x ( · ) , x ˙ ( · ) ) are bounded on [ 0 , 1 ] and x ( 0 ) = 0 , x ˙ ( 0 ) = O ( Δ x ) , we have
x ( t ) = O ( Δ x ) , x ˙ ( t ) = O ( Δ x )
uniformly for t [ 0 , 1 ] . Denote by x 0 a solution of Equation (25) with initial data
x 0 ( 0 ) = x 0 , x ˙ 0 ( 0 ) = x 0 x 0 0 x 0 m 1 ( s 2 ) m 1 ( x 0 2 ) d s .
Then x 0 ( t ) = 0 , x 0 ( t ) = x ( t ) for all t [ 0 , 1 ] , x 0 ( · ) is increasing, and
x 0 ( 1 ) = 0 , x ˙ 0 ( 1 ) = x 0 x 0 0 x 0 m 1 ( s 2 ) d s .
Denote by t ^ the value t [ 0 , 1 ] such that x 0 ( t ) = Δ x 2 / 3 . Then 1 t ^ = O ( Δ x 2 / 3 ) . On the interval [ 0 , t ^ ] we have
| ( x ( t ) x ˙ ( t ) ) 2 x ( t ) 2 x ˙ 2 ( t ) | = | ( x ( t ) x ˙ ( t ) + x ( t ) x ˙ ( t ) ) 2 x 2 ( t ) + x 2 ( t ) x ˙ 2 ( t ) | = O ( Δ x 4 / 3 ) .
Since x 2 = x 2 + x 2 , x x ˙ = x x ˙ + x x ˙ , and x ˙ 2 = x ˙ 2 + x ˙ 2 , Equations (69) and (70) imply
| h j ( x ( t ) , x ˙ ( t ) ) h j ( x ( t ) e 1 , x ˙ ( t ) e 1 ) | = O ( Δ x 4 / 3 ) , j = 1 , 2 .
On the other hand for every differentiable scalar function w we have
h 1 ( w ( t ) e 1 , w ˙ ( t ) e 1 ) = m 0 ( z ) m 0 ( z ) | z = w ( t ) 2 2 ( w ( t ) w ˙ ( t ) ) , h 2 ( w ( t ) e 1 , w ˙ ( t ) e 1 ) = [ ( m 1 ( z ) m 1 ( z ) 2 m 0 ( z ) m 0 ( z ) ) | z = w ( t ) 2 w ˙ ( t ) 2 .
Therefore, there exist constants C 1 and C 2 that are independent of bounded Δ x , and such that
| h j ( x ( t ) e 1 , x ˙ ( y ) e 1 ) h j ( x 0 ( t ) , x ˙ 0 ( t ) ) | C 1 x ( t ) e 1 x 0 ( t ) + C 2 x ˙ ( t ) e 1 x ˙ 0 ( t ) , j = 1 , 2 .
Since x and x 0 are solutions of Equation (66), there exist constants C 3 , C 4 and C 5 independent on bounded Δ x and such that
x ¨ e 1 x ¨ 0 C 3 x ˙ e 1 x ˙ 0 + C 4 x e 1 x 0 + C 5 Δ x 4 / 3 .
Taking into account that x ( 0 ) e 1 = x 0 = x 0 , x ˙ ( 0 ) e 1 x ˙ 0 ( 0 ) = o ( Δ x ) , we get
x ( t ^ ) e 1 x 0 ( t ^ ) = o ( Δ x ) , x ˙ ( t ^ ) e 1 x ˙ 0 ( t ^ ) = o ( Δ x ) .
Since 1 t ^ = O ( Δ x 2 / 3 ) , for t [ t ^ , 1 ] we have
x ˙ ( t ) x ˙ 0 ( t ) = o ( Δ x ) + O ( t t ^ ) = O ( Δ x 2 / 3 ) ,
x ( t ) x 0 ( t ) = o ( Δ x ) + O ( t t ^ ) 2 = o ( Δ x ) .
In particular,
x ( 1 ) = o ( Δ x ) , x ˙ ( 1 ) = 0 x 0 m 1 ( s 2 ) d s + O ( Δ x 2 / 3 ) .
In order to prove (54) we consider the Wronskian
W ( t ) = d e t x ( t ) x ( t ) x ˙ ( t ) x ˙ ( t ) .
It is well known that
W ( 1 ) = W ( 0 ) e 0 1 m 0 m 0 2 ( x x ˙ ) d t .
First, we calculate the exponent. Since m 1 ( 0 ) = 1 and x ( 1 ) = O ( Δ x ) , we have
0 1 m 0 m 0 2 ( x x ˙ ) d t = 0 1 d d t ( ln ( m 0 ( x ( t ) 2 ) ) d t = ln m 0 ( x ( 1 ) 2 ) m 0 ( x ( 0 ) 2 ) = ln 1 m 0 ( x 0 2 ) + O ( Δ x ) .
Hence,
e 0 1 m 0 m 0 2 ( x x ˙ ) d t = m 0 ( x 0 2 ) + O ( Δ x ) .
Then, we calculate the values of the Wronskian:
W ( 0 ) = x ( 0 ) x ˙ ( 0 ) x ˙ ( 0 ) x ( 0 ) = Δ x 1 m 0 ( x 0 2 ) 0 x 0 m 1 ( s 2 ) d s ,
W ( 1 ) = x ( 1 ) x ˙ ( 1 ) x ˙ ( 1 ) x ( 1 ) = x ( 1 ) 0 x 0 m 1 ( s 2 ) d s + o ( Δ x ) .
Thus
x ( 1 ) = m 0 ( x 0 2 ) Δ x + o ( Δ x ) .
Taking into account that x ( 1 ) = o ( Δ x ) and Δ x | | e 2 , we get
( x 0 ) ( x 0 + Δ x ) = x ( 1 ) = x ( 1 ) e 1 + x ( 1 ) e 2 = m 0 ( x 0 2 ) Δ x + o ( Δ x ) ,
and
g ( x 0 ) Δ x = m 0 ( x 0 2 ) Δ x .
This completes the proof of the theorem. □
Theorem 1 shows that there is a one-to-one correspondence between metric tensors G in the form (6), and binary operations ⊕ that G generates, as defined in Section 4.1.

4.4. Multiplication of Vectors by Numbers

In this subsection we define a function ⊗: R × B B (multiplication of vectors by numbers) that is compatible with the binary operation ⊕ in the sense that
( t 1 a ) ( t 2 a ) = ( t 1 + t 2 ) a ( t 1 ( t 2 a ) ) = ( t 1 t 2 ) a
for all t 1 , t 2 R , a B . Notice that t 1 + t 2 and t 1 t 2 in (91) are the common addition and multiplication of t 1 and t 2 in R . The first identity in (91) is called the scalar distributive law, and the second identity in (91) is called the scalar associative law.
Let a , b B be two nonzero parallel vectors (that is, belong to the same line passing through the origin). We calculate a b using the four steps described in Section 4.1.
Step 1. We define
r ˙ ( 0 ) = 1 b 0 b m 1 ( s 2 ) d s .
Step 2. We compute X ( 1 ) using (40) with X 0 = r ˙ ( 0 ) b , noting that
( I a a a 2 ) b = 0 and a b = sign ( a b ) a b .
Then
X ( 1 ) = sign ( a b ) a a 0 b m 1 ( s 2 ) m 1 ( a 2 ) d s .
Step 3. We integrate Equation (25) with initial data x ( 0 ) = a and x ˙ ( 0 ) = X ( 1 ) . We parametrize the solution x as follows: x ( t ) = p ( t ) a with scalar function p. Then p ( 0 ) = 1 , p ˙ ( 0 ) a = X ( 1 ) , x ( 1 ) = p ( 1 ) a = a b . Equation (34) implies
p ˙ ( 0 ) = p ( 0 ) p ( 1 ) m 1 ( a 2 s 2 ) m 1 ( a 2 p ( 0 ) 2 ) d s .
Therefore
sign ( a b ) a 0 b m 1 ( s 2 ) m 1 ( a 2 ) d s = p ( 0 ) p ( 1 ) m 1 ( a 2 s 2 ) m 1 ( a 2 p ( 0 ) 2 ) d s .
Following obvious simplifications we get
0 p ( 1 ) a m 1 ( s 2 ) d s = 0 a m 1 ( s 2 ) d s + sign ( a b ) 0 b m 1 ( s 2 ) d s .
Step 4. We have a b = p ( 1 ) a .
Therefore p ( 1 ) a = sign ( ( a b ) a ) a b . We define a function h as follows. For every number p we set
h ( p ) = 0 p m 1 ( s 2 ) d s .
Then, h ( p ) = sign ( p ) h ( | p | ) and
sign ( ( a b ) a ) h ( a b ) = h ( a ) + sign ( a b ) h ( b ) .
The function h is monotonically increasing. Therefore, h is invertible, the inverse of which is denoted by h 1 .
Assumption 2.
We assume that
0 1 m 1 ( s 2 ) d s = .
Under this assumption the function h is a bijection ( 1 , 1 ) ( , ) . In particular, h ( 0 ) = 0 , and for every t R and p ( 1 , 1 ) there exists h 1 ( t h ( p ) ) .
Now for every number t we define
t a = a a h 1 ( t h ( a ) ) .
Then h ( t a ) = | t | h ( a ) and ( t ) a = ( t a ) . We assume t 1 > 0 , t 2 > 0 . Then
sign ( a ( t 1 a ) ) = sign ( a ( t 2 a ) ) = sign ( a ( ( t 1 + t 2 ) a ) = 1 .
Now we have
( t 1 a ) ( t 2 a ) = a a h 1 ( h ( ( t 1 a ) ( t 2 a ) ) ) = a a h 1 ( h ( t 1 a ) + h ( t 2 a ) ) = a a h 1 ( t 1 h ( a ) + t 2 h ( a ) ) = a a h 1 ( ( t 1 + t 2 ) h ( a ) ) = ( t 1 + t 2 ) a
and
( t 1 ( t 2 a ) ) = a a h 1 ( t 1 h ( t 2 a ) ) = a a h 1 ( t 1 t 2 h ( a ) ) = ( t 1 t 2 ) a .
Owing to the property ( t ) a = ( t a ) , the cases with t 1 and t 2 of arbitrary signs are considered similarly. Accordingly, the operation of multiplication by a number is well defined.
From the definition of scalar product (101) it is easy to see that
1 a = a
and for all a B \ { 0 } , r 0
| r | a r a = a a .

4.5. Distances and Norms

We introduce the standard definition of distance between points a and b of the unit ball B .
Definition 4.
The distance between points a , b B is the minimal length of a curve connecting a and b,
d i s t ( a , b ) = min { 0 1 j = 1 n i = 1 n g i j ( x ( t ) ) x ˙ i ( t ) x ˙ j ( t ) d t } ,
where the minimum is taken over all smooth functions x: [ 0 , 1 ] B with boundary conditions x ( 0 ) = a and x ( 1 ) = b .
Obviously, d i s t ( a , b ) 0 for all a , b B , and d i s t ( a , b ) = 0 iff a = b since G ( x ) > 0 for all x B . Besides, we have the triangle inequality: for all a , b , c B
d i s t ( a , b ) + d i s t ( b , c ) d i s t ( a , c ) .
The value d i s t ( 0 , a ) is called the norm of a, denoted by a .
Lemma 2.
For every a , b B there exists z B such that a z = b .
Proof. 
Fix arbitrary points a , b B . Since B is convex, there exist curves in B connecting a and b. The minimum of the lengths of such curves does exist since G ( x ) as x 1 . Let this minimum be attained at a curve q: [ 0 , 1 ] B . Then q is a geodesic which connects a and b: q ( 0 ) = a , q ( 1 ) = b . Consider the vector q ˙ ( 0 ) . Make a parallel transport of this vector along the interval connecting a and the origin. Denote the vector at the origin by y. Consider a geodesic w: [ 0 , 1 ] B with initial conditions w ( 0 ) = 0 , w ˙ ( 0 ) = y . Denote z = w ( 1 ) . Then according to the steps 1,2,3 described in the Section 4.1 we have b = a z . The lemma is proved. □
Consider a geodesic x: [ 0 , 1 ] B with boundary conditions x ( 0 ) = a and x ( 1 ) = b . It is known that d i s t ( a , b ) is equal to the length of this geodesic, and that along the geodesic the integrand on the right-hand side of (107) is constant [56]. Therefore
d i s t ( a , b ) = j = 1 n i = 1 n g i j ( x ( 0 ) ) x ˙ i ( 0 ) x ˙ j ( 0 ) .
Let us consider a geodesic y: [ 0 , 1 ] B with boundary points y ( 0 ) = 0 and y ( 1 ) = z . The value of j = 1 n i = 1 n g i j ( y ( t ) ) y ˙ i ( t ) y ˙ j ( t ) is constant over t [ 0 , 1 ] . Hence,
z = j = 1 n i = 1 n g i j ( y ( 0 ) ) y ˙ i ( 0 ) y ˙ j ( 0 ) .
But the vector x ˙ ( 0 ) is a parallel transport of the vector y ˙ ( 0 ) along the curve { a s : s [ 0 , 1 ] } . Therefore,
j = 1 n i = 1 n g i j ( x ( 0 ) ) x ˙ i ( 0 ) x ˙ j ( 0 ) = j = 1 n i = 1 n g i j ( y ( 0 ) ) y ˙ i ( 0 ) y ˙ j ( 0 ) .
Finally,
d i s t ( a , b ) = j = 1 n i = 1 n g i j ( y ( 0 ) ) y ˙ i ( 0 ) y ˙ j ( 0 ) = z .
In Section 3.4 we found that y ( t ) = r ( t ) z , and
r ˙ ( 0 ) = 1 z 0 z m 1 ( s 2 ) d s .
In view of the equation
j = 1 n i = 1 n g i j ( y ( 0 ) ) y ˙ i ( 0 ) y ˙ j ( 0 ) = r ˙ ( 0 ) z
we have
z = d i s t ( a , b ) = 0 z m 1 ( s 2 ) d s = h ( z ) .
Interestingly, the norm · does not depend on the function m 0 . In particular, if we have two spaces with metric tensors that have the same function m 1 , then the distances between points (and hence the norms) in these two spaces coincide, and the operations of multiplication by a number also coincide. We’ll see this result in several examples, including the spaces with Einstein and Möbius additions.
Now let us consider a binary operation ⊕ applied to numbers. That is, given functions m 0 , m 1 , we consider the tensor G given by (6), and its resulting binary operation ⊕ introduced in Section 4.1. Also, let’s consider a tensor with the same functions m 0 and m 1 in the one dimensional space. The corresponding operation between numbers is denoted by the same symbol, ⊕.
For arbitrary numbers p, q the value of p q is defined by the four-steps procedure presented in Section 4.1. If q = 0 , then p q = p . Assume q 0 .
Step 1. We calculate
r ˙ ( 0 ) = 1 | q | 0 | q | m 1 ( s 2 ) d s .
Step 2. Since 1 q q | q | 2 = 0 , we have
X 0 = 0 q m 1 ( s 2 ) d s , X 1 = 1 m 1 ( q 2 ) 0 q m 1 ( s 2 ) d s .
Step 3. If x is a scalar then Equation (25) for geodesics takes the form
x ¨ + m 1 m 1 ( x ˙ ) 2 x = 0 .
A solution x of (118) with the initial values x ( 0 ) = p , x ˙ ( 0 ) = X 1 satisfies the condition
p x ( t ) m 1 ( s 2 ) d s = t 0 q m 1 ( s 2 ) d s .
for all t > 0 .
Step 4. Since p q = x ( 1 ) , we have
0 p q m 1 ( s 2 ) d s = 0 p m 1 ( s 2 ) d s + 0 q m 1 ( s 2 ) d s .
Thus, taking into account the definition of the function h, we get
h ( p q ) = h ( p ) + h ( q ) .
Similarly, the operation ⊗ for numbers is defined by (101). In particular, for arbitrary numbers r and p we have
r p = h 1 ( r h ( p ) ) .
Thus, for every number r and vector a B
r a = | h 1 ( r h ( a ) ) | = | r | a .
Let a , b B \ { 0 } . The triangle inequality implies
d i s t ( a , 0 ) + d i s t ( 0 , b ) d i s t ( a , b )
and equality is attained only if the vectors a and b belong to the same ray, that is, if there exists a positive number λ such that a = λ b .
In terms of the function h, inequality (124) yields
h ( a ) + h ( b ) h ( a b ) .
Employing (121), we get
a b = h 1 [ h ( a ) + h ( b ) ] h 1 [ h ( a b ) ] = a b .
We have thus proved the following theorem.
Theorem 2.
For all vectors a , b B we have the gyrotriangle inequality
a b a b ,
and equality is attained iff the vectors a , b lie on the same ray, that is, for some nonnegative number λ we have a = λ b or b = λ a .
Inequality (127) may be considered as a triangle inequality in the spaces of vectors and numbers with the same addition operation ⊕.
It is shown below that for Einstein addition and all additions isomorphic to Einstein addition (for example, Möbius addition) a solution of the equation a x = b is given by x = ( a ) b . Hence, we have the following result.
Theorem 3.
If a solution x of the equation a x = b is given by x = ( a ) b , then
d i s t ( a , b ) = h ( ( a ) b ) .

5. Spaces with Einstein Addition

In this section, we consider a special binary operation in the Beltrami-Klein ball model of hyperbolic geometry that turns out to be Einstein addition [1,5].
We consider the metric tensor G in (6)–(8) parametrized by the functions
m 0 ( r 2 ) = 1 1 r 2 and m 1 ( r 2 ) = 1 ( 1 r 2 ) 2 ,
so that
G ( x ) = 1 1 x 2 I x x x 2 + 1 ( 1 x 2 ) 2 x x x 2 .

5.1. Einstein Addition

We follow the four steps described in Section 4:
Step 1. We evaluate
r ˙ ( 0 ) = 1 b 0 b d s 1 s 2 = atanh ( b ) b .
Step 2. We find a solution of Equation (40) with X 0 = atanh ( b ) b b ,
X ( 1 ) = atanh ( b ) b [ 1 a 2 ( I a a a 2 ) + ( 1 a 2 ) a a a 2 ] b .
Step 3. We solve Equation (25) with initial conditions x ( 0 ) = a , x ˙ ( 0 ) = X ( 1 ) . Noticing that
m 1 ( r 2 ) m 1 ( r 2 ) 2 m 0 ( r 2 ) m 0 ( r 2 ) = 0 m 1 ( r 2 ) m 0 ( r 2 ) r 2 m 0 ( r 2 ) = 0 ,
Equation (25) takes the form
x ¨ ( t ) + d ( ln ( m 0 ( x ( t ) 2 ) ) d t x ˙ ( t ) = 0 .
The solution to the initial value problem is given by
x ( t ) = a + tanh ( t atanh ( b ) ) atanh ( b ) 1 + a b b tanh ( t atanh ( b ) ) X ( 1 )
Step 4. Finally, we obtain the addition
a E b = a + 1 a 2 1 + a b [ b a b 1 a 2 + 1 a ] = 1 1 + a b { [ 1 + a b a 2 ( 1 1 a 2 ) ] a + 1 a 2 b } ,
for all a , b B , which is recognized as Einstein addition.
The binary operation described in Section 4 enables Einstein addition (136) to be recovered from the metric tensor (130). Thus, the operation a E b is determined by its local properties in a neighborhood of the set { ( a , b ) : a , b B , a = b } .
The binary operation E in (136) turns out to be the well known Einstein addition in the ball, studied in special relativity theory, where the speed of light in empty space is normalized to c = 1 [1,4,5].

5.2. Einstein Multiplication by a Number

To define an operation of multiplication by a number we evaluate the function h in (98) with the function m 1 in (129), obtaining,
h ( p ) = 0 p 1 1 s 2 d s = atanh ( p ) .
Therefore, by (101),
t E a = atanh ( t tanh ( a ) ) a a .

5.3. Derivation of the Metric Tensor Associated with Einstein Addition

We have derived Einstein addition from the metric tensor (130). According to Theorem 1 we can make also the reverse procedure, that is, we can recover the metric tensor (130) from Einstein addition. In this subsection we demonstrate it directly. Let us consider the Einstein difference in the following chain of equations.
( x ) E ( x + Δ x ) = x + 1 1 x ( x + Δ x ) [ 1 x 2 ( x + Δ x x ( x + Δ x ) x 2 x ) + x ( x + Δ x ) x 2 ( 1 x 2 ) x ] = 1 1 x ( x + Δ x ) [ 1 x 2 Δ x + ( 1 1 x 2 ) x Δ x x 2 x ] = g ( x ) Δ x + o ( Δ x ) ,
where the n × n -matrix g ( x ) is given by
g ( x ) = 1 1 x 2 [ 1 x 2 ( I x x x 2 ) + x x x 2 ] .
Then, the metric tensor G ( x ) is given by the equation
G ( x ) = g ( x ) g ( x ) = 1 1 x 2 ( I x x x 2 ) + 1 ( 1 x 2 ) 2 x x x 2 ,
thus recovering (130).

6. Spaces with Möbius Addition

Contrasting (129), let’s assume that
m 0 ( r 2 ) = m 1 ( r 2 ) = 1 ( 1 r 2 ) 2 ,
so that the metric tensor G ( x ) in (6) and (9) takes the form
G ( x ) = 1 ( 1 x 2 ) 2 I .

6.1. Möbius Addition

Again, we follow the four steps described in Section 4.
Step 1. Similar to the case of Einstein addition, we evaluate
r ˙ ( 0 ) = 1 b 0 b 1 1 s 2 d s = atanh ( b ) b .
Step 2. With
X 0 = atanh ( b ) b ,
we find a solution X ( 1 ) to Equation (40),
X ( 1 ) = ( 1 a 2 ) atanh ( b ) b b .
Step 3. We solve the initial value problem
x ¨ + 2 1 x 2 [ 2 ( x x ˙ ) x ˙ x ˙ 2 x ] = 0 x ( 0 ) = a x ˙ ( 0 ) = X ( 1 ) ,
seeking a solution of the form
x ( t ) = P + R 1 cos φ ( t ) + R 2 sin φ ( t ) ,
with constant vectors P, R 1 , R 2 .
The initial value problem (147) has a solution (148) with the following parameters and a function φ , as may be checked directly,
R 1 = ( 1 a 2 ) b 2 2 ( a 2 b 2 ( a b ) 2 ) a + ( 1 a 2 ) a b 2 ( a 2 b 2 ( a b ) 2 ) b R 2 = 1 a 2 2 a 2 b 2 ( a b ) 2 b P = a R 1 φ ( t ) = 2 arctan a 2 b 2 ( a b ) 2 tanh ( atanh ( b ) t ) b + a b tanh ( atanh ( b ) t ) .
Therefore,
x ( 1 ) = P + R 1 cos φ ( 1 ) + R 2 sin φ ( 1 ) = ( 1 + 2 a b + b 2 ) a + ( 1 a 2 ) b 1 + 2 a b + a 2 b 2 .
We notice that
R 1 = R 2 = P x ( t ) , t 0 ,
and
P 2 = R 1 2 + 1 .
Therefore, each geodesic is a circular arc centered at P with radius R 1 , which intersects the unit circle orthogonally, and it tends to one of these two intersections as t and as t .
Step 4. Finally, we get
a M b = ( 1 + 2 a b + b 2 ) a + ( 1 a 2 ) b 1 + 2 a b + a 2 b 2 ,
for all a , b B , which is the well known formula for Möbius addition in the ball [1,4].
Möbius addition is recovered here from the metric tensor (143). Thus, the operation a M b is determined by its local properties in a neighborhood of the set { ( a , b ) : a , b B , a = b } .

6.2. Möbius Multiplication by a Number

Since the function m 1 is the same for both Einstein and Möbius metric tensors, the multiplication by a scalar is the same as in (138) in Section 5.2,
t M a = atanh ( t tanh ( a ) ) a a .

6.3. Derivation of the Metric Tensor Associated with Möbius Addition

As for Einstein addition, to recover the underlying metric tensor from its associated Möbius addition it is possible to use Theorem 1. In this subsection we calculate directly the metric tensor G that determines the binary operation M . Let us consider the Möbius difference
( x ) M ( x + Δ x ) = 1 2 x ( x + Δ x ) + x 2 x + Δ x 2 ) ( x ) + ( 1 x 2 ) ( x + Δ x ) 1 2 x ( x + Δ x ) + x 2 x + Δ x 2 = ( 1 x 2 ) Δ x Δ x 2 x ( 1 x 2 ) 2 2 x ( 1 x 2 ) + x 2 Δ x 2 = g M ( x ) Δ x + o ( Δ x ) ,
where
g M ( x ) = 1 1 x 2 I .
Then, the metric tensor G ( x ) is given by the equation
G ( x ) = g M ( x ) g M ( x ) = 1 ( 1 x 2 ) 2 I ,
thus recovering (143).

7. A Space with an Operation Isomorphic to Euclidean Addition

Assume
m 0 ( r 2 ) = 1 1 r 2 , m 1 ( r 2 ) = 1 ( 1 r 2 ) 3 .
Then the metric tensor (6) takes the form
G ( x ) = 1 1 x 2 I x x x 2 + 1 ( 1 x 2 ) 3 x x x 2 .

7.1. Binary Operation

We follow the four steps described in Section 4.
Step 1. We calculate
r ˙ ( 0 ) = 1 b 0 b 1 ( 1 s 2 ) 3 d s = 1 1 b 2 .
Step 2. We find a solution of Equation (40) with X 0 = b 1 b 2 ,
X ( 1 ) = 1 a 2 1 b 2 b ( a b ) a .
Step 3. We solve Equation (25) with initial conditions x ( 0 ) = a , x ˙ ( 0 ) = X ( 1 ) . Noticing that
m 1 ( r 2 ) m 1 ( r 2 ) 2 m 0 ( r 2 ) m 0 ( r 2 ) = 1 1 r 2 , m 1 ( r 2 ) m 0 ( r 2 ) r 2 m 0 ( r 2 ) r 2 m 1 ( r 2 ) = 1 ,
we see that Equation (25) takes the form
x ¨ + 2 x x ˙ 1 x 2 x ˙ + x ˙ 2 + ( x x ˙ ) 2 1 x 2 x = 0 .
Assume vectors a and b are not parallel. The solution of the initial value problem lies entirely in the plane containing the vectors a and b. Denote by
b = b a b a 2 a
the component of vector b orthogonal to a. Denote r ( t ) = x ( t ) . Then the vector function x may be presented in the polar coordinates as follows:
x = r cos φ a a + r sin φ b b ,
where φ is a scalar function, the angle of x in the plane with unit vectors a a and b b . Then r ( 0 ) = a , φ ( 0 ) = 0 , x ˙ ( 0 ) = r ˙ ( 0 ) a a + r ( 0 ) φ ˙ ( 0 ) b b , and
b r ( 0 ) φ ˙ ( 0 ) = b x ˙ ( 0 ) = b X ( 1 ) = 1 a 2 1 b 2 b b > 0 .
Therefore φ ˙ ( 0 ) > 0 . Moreover,
x x ˙ = r r ˙ , x ˙ 2 = r ˙ 2 + r 2 φ ˙ 2 , x ˙ 2 ( x x ˙ ) 2 x 2 = r 2 φ ˙ 2 ,
and Equation (163) is equivalent to the following two scalar equations:
r ¨ + 3 1 r 2 r ˙ 2 r + r ( r 2 1 ) φ ˙ 2 = 0 ,
and
r φ ¨ + 2 φ ˙ r ˙ 1 r 2 = 0 .
Equation (169) can be written as
d d t r 2 φ ˙ 1 r 2 = 0 .
Hence, there exists a number C such that
r 2 φ ˙ 1 r 2 = C .
Since φ ˙ ( 0 ) > 0 , we have C > 0 , and φ ˙ ( t ) > 0 for all t [ 0 , 1 ] .
Substituting (171) into (168) yields
r ¨ + 3 r r ˙ 1 r 2 ( 1 r 2 ) 3 C 2 r 3 = 0 .
It follows that if r ˙ ( 0 ) 0 , then r ˙ ( t ) > 0 for all t > 0 . If r ˙ ( 0 ) < 0 , then there exists a number t ^ > 0 such that r ˙ ( t ) < 0 for all t [ 0 , t ^ ) , r ˙ ( t ^ ) = 0 , and r ˙ ( t ) > 0 for all t > t ^ . Notice that
r ( 0 ) r ˙ ( 0 ) = x ( 0 ) x ˙ ( 0 ) = a X ( 1 ) = a b ( 1 a 2 ) 3 / 2 ( 1 b 2 ) 1 / 2 .
Therefore, sign r ˙ ( 0 ) = sign a b . Below we consider the case r ˙ ( 0 ) > 0 . The case r ˙ ( 0 ) 0 may be studied in a very similar manner and leads to the same formula for a b .
Let z = r 2 . Then
z ˙ = 2 r r ˙ , z ¨ = 2 r r ¨ + 2 r ˙ 2 ,
and Equation (172) acquires the form
z ¨ 2 z ˙ 2 4 z + 3 z ˙ 2 4 ( 1 z ) ( 1 z ) 3 C 2 z = 0 .
Since
d d t z ˙ 2 z ( 1 z ) 3 + 4 C 2 ( 1 z ) z = 4 z ˙ z ( 1 z ) 3 z ¨ 2 z ˙ 2 4 z + 3 z ˙ 2 4 ( 1 z ) ( 1 z ) 3 C 2 z = 0 ,
there exists a constant C 1 such that
z ˙ 2 z ( 1 z ) 3 + 4 C 2 ( 1 z ) z = C 1 .
Equation (177) is autonomous, and therefore separable,
d z ( 1 z ) 3 ( C 1 z 4 C 2 ( 1 z ) ) = d t .
Let’s calculate C and C 1 in terms of a and b. Let
p = 0 b d s ( 1 s 2 ) 3 = b 1 b 2 ,
and let α ( 0 , π / 2 ) be the angle between a and b:
cos α = a b a b .
Then
x ( 0 ) = a x ˙ ( 0 ) = X 1 = p 1 a 2 b b a b b a .
We have r ( 0 ) = x ( 0 ) = a ,
2 x ( 0 ) x ˙ ( 0 ) = d ( z ) 2 d t = 2 r ( 0 ) r ˙ ( 0 ) .
Therefore,
r ˙ ( 0 ) = p ( 1 a 2 ) 3 a b a b = p ( 1 a 2 ) 3 cos α .
For every t we have
φ ˙ 2 = 1 x 2 x ˙ 2 ( x x ˙ ) 2 x 2 .
Therefore
φ ˙ ( 0 ) = p 1 a 2 a sin α .
From (171) we get
C = r ( 0 ) 2 φ ˙ ( 0 ) 1 r ( 0 ) 2 = p a sin α 1 a 2 .
Now we calculate the constant C 1 using (177),
C 1 = 4 r ( 0 ) 2 r ˙ ( 0 ) 2 r ( 0 ) 2 ( 1 r ( 0 ) 2 ) 3 + 4 C 2 ( 1 r ( 0 ) 2 ) r ( 0 ) 2 = 4 p 2 .
We integrate Equation (178) over t [ 0 , 1 ] using these values of C, C 1 and the equation z = r 2 , obtaining
0 1 r ˙ d t ( 1 r 2 ) 3 1 a 2 ( 1 r 2 ) r 2 ( 1 a 2 ) sin 2 α = p .
Changing the variable of integration on the left-hand side, and using the initial condition r ( 0 ) = a , we have
a r ( 1 ) d s ( 1 s 2 ) 3 1 a 2 ( 1 s 2 ) s 2 ( 1 a 2 ) sin 2 α = 0 b d s ( 1 s 2 ) 3 .
Both integrals can be evaluated in terms of elementary functions,
r ( 1 ) 2 1 r ( 1 ) 2 a 2 sin 2 α 1 a 2 a 2 cos 2 α 1 a 2 = b 1 b 2 .
In order to calculate φ ( 1 ) , we solve Equation (171) for φ ˙ ,
φ ˙ = C 1 r 2 r 2 .
Now we rewrite Equation (178) in terms of r,
r ˙ = p ( 1 r 2 ) 3 ( 1 a 2 ( 1 r 2 ) r 2 ( 1 a 2 ) sin 2 α .
Then, we multiply and divide the right-hand side of (190) by d r d t , use Equation (192) and integrate (190) over [ 0 , 1 ] taking into account that φ ( 0 ) = 0 , obtaining
φ ( 1 ) = 0 1 sin α r ˙ d t r 2 ( 1 r 2 ) 2 r 2 1 r 2 a 2 1 a 2 sin 2 α a 2 1 a 2 = a r ( 1 ) sin α d r r 2 ( 1 r 2 ) 2 r 2 1 r 2 a 2 1 a 2 sin 2 α a 2 1 a 2 = α arctan a 2 sin 2 α ( 1 r ( 1 ) 2 ) ( 1 a 2 ) r ( 1 ) 2 ( 1 r ( 1 ) 2 ) a 2 1 a 2 sin 2 α .
Noticing that
r ( 1 ) 2 1 r ( 1 ) 2 = a 1 a 2 + b 1 b 2 2
and
r ( 1 ) 2 1 r ( 1 ) 2 a 2 1 a 2 sin 2 α = a 2 cos 2 α 1 a 2 + b 2 1 b 2 2 ,
we have
cos φ ( 1 ) = r ( 1 ) 2 1 r ( 1 ) 2 a 2 1 a 2 + b 2 cos 2 α 1 b 2 sin φ ( 1 ) = r ( 1 ) 2 1 r ( 1 ) 2 b 2 sin 2 α 1 b 2 .
Taking into account (165), (195) and (196), we finally get
x ( 1 ) = a 1 a 2 + b 1 b 2 1 + a 1 a 2 + b 1 b 2 2 = γ a a + γ b b 1 + γ a a + γ b b 2 ,
where
γ a = 1 1 a 2 .
Step 4. We set
a b = x ( 1 ) = γ a a + γ b b 1 + γ a a + γ b b 2 .
To see that the novel operation ⊕ in the ball of R n is isomorphic to the Euclidean vector addition in R n , let us consider the mapping f: B R n given by
f ( x ) = γ x x .
Then
1 a b 2 = 1 1 + γ a a + γ b b 2
and
f ( a b ) = γ a a + γ b b = f ( a ) + f ( b ) .
Equation (202) implies that the novel operation ⊕ in the ball is isomorphic to the Euclidean operation, +, in R n and hence it is commutative and associative.

7.2. Multiplication by Numbers

To define the operation ⊗ of multiplication by numbers satisfying properties (91) in the Section 4.4 we have to introduce the function h: [ 0 , 1 ) R , (98),
h ( p ) = 0 p 1 ( 1 s 2 ) 3 d s = p 1 p 2 .
Obviously, h is an increasing bijection. Therefore the inverse function h 1 exists. It is easy to see that
h 1 ( q ) = q 1 + q 2 .
The operation ⊗ is given by Formula (101), which in our case has the form
t a = h 1 ( t h ( a ) a a = t 1 + ( t 2 1 ) a 2 a a .
Combining (202) and (203), (203) and (204) we get
f ( t a ) = t a 1 t a 2 = h ( t a ) t a a = t h ( a ) a a = t f ( a ) .
Equations (202) and (206) imply the distributive property of the operations ⊕ and ⊗: for any numbers t 1 , t 2 , and vectors a , b B we have
f ( ( t 1 a ) ( t 2 b ) ) = t 1 f ( a ) + t 2 f ( b ) .
Thus, the operation ⊕ is isomorphic to Euclidean addition.

7.3. Derivation of the Metric Tensor Associated with the Binary Operation

To recover the underlying metric tensor from the operation ⊕ we can use Theorem 1. Noticing that
( x ) ( x + Δ x ) = x 1 x 2 + x + Δ x 1 x + Δ x 2 1 + x 1 x 2 + x + Δ x 1 x + Δ x 2 2 = Δ x 1 x 2 + x 1 1 x + Δ x 2 1 1 x 2 + o ( Δ x ) = 1 1 x 2 I + x x ( 1 x 2 ) 3 Δ x + o ( Δ x ) ,
we have
g ( x ) = 1 1 x 2 I + x x ( 1 x 2 ) 3 = 1 1 x 2 I x x x 2 + 1 ( 1 x 2 ) 3 x x x 2 . .
The metric tensor G associated with the operation ⊕ is given in (51),
G ( x ) = g ( x ) g ( x ) = 1 1 x 2 I x x x 2 + 1 ( 1 x 2 ) 3 x x x 2 .
It has the form (6) parametrized by the functions
m 0 ( r 2 ) = 1 1 r 2 m 1 ( r 2 ) = 1 ( 1 r 2 ) 3 ,
thus recovering (158)–(159).

8. Properties of Einstein Addition

Seemingly structureless, Einstein addition is neither commutative nor associative. The Einstein groupoid ( B , E ) formed by the ball B equipped with Einstein addition E is a gyrocommutative gyrogroup. A gyrogroup is a group-like object that possesses a rich structure [1,5], formally defined in Definitions 1 and 2, Section 2. Einstein addition is a gyrocommutative gyrogroup operation, several properties of which are presented below.
  • Left Cancellation Law:
    a E ( ( a ) E b ) = b , a , b B .
  • Existence of Gyrations: for every a , b B there exists a unitary matrix denoted by g y r [ a , b ] such that for all c B we have the following gyroassociative law:
    a E ( b E c ) = ( a E b ) E ( g y r [ a , b ] c ) , a , b , c B .
  • Gyrocommutative Law:
    a E b = g y r [ a , b ] ( b E a ) , a , b B ,
    such that
    a E b = b E a .
  • Left Reduction Property:
    g y r [ a E b , b ] = g y r [ a , b ] , a , b B .
The validity of these properties was proved by using computer algebra (specifically, the software MATHEMATICA for symbolic manipulation). Below we prove these properties by hand using the definition of Einstein addition (136), and provide the matrix representation of gyrations g y r [ a , b ] .

8.1. Left Cancellation Law

By Einstein addition (136), we have
( a ) E b = a + 1 1 a b [ 1 a 2 ( b a b a 2 a ) + a b a 2 ( 1 a 2 ) a ] .
We notice that
a [ ( a ) E b ] = a 2 + 1 a 2 1 a b a b = a b a 2 1 a b
and
1 + a [ ( a ) E b ] = 1 a 2 1 a b .
Therefore,
a E ( ( a ) E b ) = a + 1 a 2 1 a 2 ( 1 a b ) { a + 1 1 a b [ 1 a 2 ( b a b a 2 a ) + a b a 2 ( 1 a 2 ) a + ( a b a 2 ) ( 1 a 2 1 ) a 2 a ] } = a + 1 1 a 2 [ a + a b a + 1 a 2 b 1 a 2 a b a 2 a + a b a 2 ( 1 a 2 ) a + ( a b a 2 ) ( 1 a 2 1 ) a a 2 ] = b ,
thus proving the left cancellation law.
The property established in this subsection may be reformulated as follows. For every vectors a , b B there exists a unique vector x B such that
a E x = b
and, moreover,
x = ( a ) E b .
Eliminating x between (221) and (222) yields the left cancellation law (212).

8.2. Existence of Gyrations

Evaluating the squared norm a E b 2 yields an elegant result,
a E b 2 = 1 a 2 1 + a b ( I a a a 2 ) b + a 2 + a b a 2 ( 1 + a b ) a 2 = 1 a 2 ( 1 + a b ) 2 ( b 2 ( a b ) 2 a 2 ) + ( a 2 + a b ) 2 ( 1 + a b ) 2 a 2 = 1 ( 1 a 2 ) ( 1 b 2 ) ( 1 + a b ) 2 .
The function on the extreme right-hand side of (223) is symmetric with respect to a , b . Hence, a E b = b E a . Owing to this symmetry, there exists a unitary matrix which is the identity matrix in the subspace orthogonal to both a and b, and which maps b E a to a E b . This matrix is said to be the gyration generated by a and b, and is denoted by g y r [ a , b ] . Gyrations play an important role in the geometric theory of spaces with the binary operation a E b . Hence, we introduce them along with their properties in Section 8.3.

8.3. Definition of Gyrations

For every three vectors a , b , c B we consider a vector d B given by
d = ( a E b ) E ( a E ( b E c ) ) .
First we check that d depends on c linearly, that is, d is equal to some matrix depending on a, b only, times c. This matrix is destined to be called the gyration matrix generated by a and b.
Owing to the left cancellation law, Equation (224) is equivalent to
( a E b ) E d = a E ( b E c ) .
Using (136), we get,
( a E b ) E d = 1 1 + ( a E b ) d { ( 1 + ( a E b ) d a E b 2 ) a E b + 1 a E b 2 [ d ( a E b ) d a E b 2 ( a E b ) ] }
and
a E ( b E c ) = 1 1 + a ( b E c ) { [ 1 + a ( b E c ) a 2 ( 1 1 a 2 ) ] a + 1 a 2 ( b E c ) } .
Furthermore, we need identity (228) of the following Lemma.
Lemma 3.
For every vectors a , b , c B we have
( 1 + b c ) ( 1 + a ( b E c ) ) = ( 1 + a b ) ( 1 + ( b E a ) c ) .
Proof. 
By means of (136) we have,
( 1 + b c ) ( 1 + a ( b E c ) ) = 1 + b c + ( 1 + b c b 2 ) a b + 1 b 2 ( a c b c b 2 a b )
and
( 1 + a b ) ( 1 + ( b E a ) c ) = 1 + a b + ( 1 + b a b 2 ) b c + 1 b 2 ( a c b a b 2 b c ) .
The right-hand sides of (229) and (230) are identical as desired, and the proof is complete. □
We solve Equation (225) for d under the assumption that
( a E b ) d = ( b E a ) c .
This assumption will be verified later, below (237).
Then, taking into account Lemma 3, we see that (225) is equivalent to
[ 1 + ( b E a ) c a E b 2 ( 1 1 a E b 2 ) ] ( a E b ) + 1 a E b 2 d = 1 + b c 1 + a b { [ 1 + a ( b E c ) a 2 ( 1 1 a 2 ) ] a + 1 a 2 ( b E c ) } .
Equation (232) allows us to determine d,
d = 1 1 a E b 2 [ 1 + b c 1 + a b { [ 1 + a ( b E c ) a 2 ( 1 1 a 2 ) ] a + 1 a 2 ( b E c ) } [ 1 + ( b E a ) c a E b 2 ( 1 1 a E b 2 ) ] ( a E b ) ] .
To prove that (233) is a solution of (225) for the unknown d, we need only to show that (231) holds.
By means of (136), we have
d = 1 1 a E b 2 ( 1 + a b ) [ a + b c a + ( 1 1 a 2 ) [ ( 1 + b c b 2 ) a b + 1 b 2 ( a c b c b 2 a b ) ] a + 1 a 2 [ ( 1 + b c b 2 ) b + 1 b 2 ( c b c b 2 b ) ] ( 1 + a b ) [ ( a E b ) 1 a E b 2 ( b E a ) c a E b 2 ( a E b ) + ( b E a ) c a E b 2 ( a E b ) ] ] .
The terms without c in the big square brackets cancel each other, and we get the following expression for d, which is linear in c,
d = ( b E a ) c a E b 2 ( a E b ) + 1 1 a E b 2 ( 1 + a b ) [ b c a + ( 1 1 a 2 ) [ b c b 2 a b + 1 b 2 ( a c b c b 2 a b ) ] a + 1 a 2 [ b c b 2 b + 1 b 2 ( c b c b 2 b ) ] ( 1 + a b ) ( a E b ) a E b 2 [ ( 1 + b a b 2 ) b c + 1 b 2 ( a c b c b 2 a b ) ] ] .
Now we check condition (231). We use the following four formulas,
a ( a E b ) = | a 2 + a b 1 + a b b ( a E b ) = 1 a b [ ( 1 + a b a 2 ) a b + 1 a 2 ( b 2 ( a b ) 2 a 2 ) ] c ( a E b ) = 1 1 + a b [ ( 1 + a b a 2 ) c b + 1 a 2 ( c b a b a 2 b c ) ] 1 a E b 2 ( 1 + a b ) = ( 1 a 2 ) ( 1 b 2 ) ,
in the representation
( a E b ) d ( b E a ) c = 1 ( 1 a 2 ) ( 1 b 2 ) [ b c a ( a E b ) + ( 1 1 a 2 ) [ b c b 2 a b + 1 b 2 ( a c b c b 2 a b ) ] a ( a E b ) + 1 a 2 [ b c b 2 b ( a E b ) + 1 b 2 ( c ( a E b ) b c b 2 b ( a E b ) ) ] ( 1 + a b ) [ ( 1 + b a b 2 ) b c + 1 b 2 ( a c b c b 2 a b ) ] ]
to show by direct calculations that the right-hand side of (237) vanishes. Hence, Equation (231) is true. Therefore, the vector d satisfies (225). Equation (235) shows that d depends linearly on c,
d = M c ,
where M is a matrix depending only on a and b. This matrix is called the gyration matrix generated by a and b, and is denoted by g y r [ a , b ] .
In particular, (225) implies the Gyroassociative law (213):
a E ( b E c ) = ( a E b ) E ( M c ) a , b , c B .

8.4. Properties of the Gyration

We have already shown that the gyration operator g y r [ a , b ] , defined by the identity
( a E b ) E g y r [ a , b ] c = a E ( b E c ) ,
is linear, that is, g y r [ a , b ] c = M c for some matrix M. Below we determine the matrix M and prove that M is unitary.
From (235) and (238) we obtain the following representation of the matrix M.
M = g y r [ a , b ] = I + C a a a a + C a b a b + C b a b a + C b b b b ,
where the scalar functions C a a , C a b , C b a , C b b are given by
C a a = 1 a 2 ( 1 1 a 2 ) ( 1 1 b 2 ) 1 + a b + ( 1 a 2 ) ( 1 b 2 ) C a b = 1 + 2 a b ( 1 + 1 a 2 ) ( 1 + 1 b 2 ) 1 + a b + ( 1 a 2 ) ( 1 b 2 ) C b a = 1 1 + a b + ( 1 a 2 ) ( 1 b 2 ) C b b = 1 b 2 ( 1 1 a 2 ) ( 1 1 b 2 ) 1 + a b + ( 1 a 2 ) ( 1 b 2 ) .
Lemma 4.
The matrix M is unitary,
M M = I .
Proof. 
We have
M M = I + α a a a a + α a b a b + α b a b a + α b b b b ,
where the numbers α a a , α a b , α b a , α b b are given by
α a a = 2 C a a + a 2 C a a 2 + 2 C a b C a a a b + C a b 2 b 2 α a b = α b a = C a b + C b a + a b ( C a a C b b + C a b C b a ) + C a a C b a a 2 + C a b C b b b 2 α b b = 2 C b b + b 2 C b b 2 + 2 C b a C b b a b + C b a 2 b 2 .
Straightforward calculations show that all these numbers vanish. To show these results we use the notation x = 1 a 2 and y = 1 b 2 , so that a 2 = 1 x 2 and b 2 = 1 y 2 . Indeed, we have
α a a ( 1 + a b + x y ) a 2 = 2 ( 1 x ) ( 1 y ) ( 1 + a b + x y ) + ( 1 x ) 2 ( 1 y ) 2 2 a b ( 1 x ) ( 1 y ) ( 1 + 2 a b ( 1 x ) ( 1 y ) a 2 b 2 + ( 1 x 2 ) ( 1 y 2 ) + 4 a b ( 1 x ) ( 1 y ) + 4 ( a b ) 2 ( 1 x ) 2 ( 1 y ) 2 a 2 b 2 = 0 ,
α b b ( 1 + a b + x y ) b 2 = 2 ( 1 x ) ( 1 y ) ( 1 + a b + x y ) + ( 1 x ) 2 ( 1 y ) 2 + 2 a b ( 1 x ) ( 1 y ) + ( 1 x 2 ) ( 1 y 2 ) = 0
and
α a b ( 1 + a b + x y ) a 2 b 2 = ( ( 1 x 2 ) ( 1 y 2 ) + 2 a b ( 1 x ) ( 1 y ) ) ( 1 + a b + x y ) ( 1 x 2 ) ( 1 y 2 ) ( 1 + a b + x y ) + a b [ ( 1 x ) 2 ( 1 y ) 2 ( ( 1 x 2 ) ( 1 y 2 ) + 2 a b ( 1 x ) ( 1 y ) ) ] + ( 1 x ) ( 1 y ) ( 1 x 2 ) ( 1 y 2 ) [ ( 1 x 2 ) ( 1 y 2 ) + 2 a b ( 1 x ) ( 1 y ) ] ( 1 x ) ( 1 y ) = 0 .
The proof of the Lemma is, thus, complete. □
Lemma 5.
If vectors a and b are parallel, that is, there exists a non zero pair of numbers ( λ 1 , λ 2 ) such that
λ 1 a + λ 2 b = 0 ,
then
g y r [ a , b ] = I .
Proof. 
If λ 2 = 0 , then a = 0 , and the statement follows from (241). Assume λ 2 0 . Set λ = λ 1 λ 2 . Then b = λ a , and the representation (241) amounts to
g y r [ a , b ] = I + d a a ,
where
d = C a a + λ ( C a ( λ a ) + C ( λ a ) a ) + λ 2 C ( λ a ) ( λ a ) = 0
 □
Corollary 1.
For every numbers r 1 , r 2 and a vector v B we have
g y r [ r 1 v , r 2 v ] = I .
Proof. 
By the definition of the scalar multiplication in (101), the vectors r 1 v and r 2 v are parallel. Hence, (253) follows from Lemma 5. □

8.5. Reduction Property

In this subsection we establish property (216) of Einstein addition. Let
γ v = 1 1 v 2
be the well known Lorentz gamma factor of special relativity theory [1,5]. Then, in view of (242) we have
C a a = γ a 2 γ a + 1 γ b 1 γ a γ b ( 1 + a b ) + 1 C a b = γ a γ b + 2 a b γ a 2 γ b 2 ( γ a + 1 ) ( γ b + 1 ) γ a γ b ( 1 + a b ) + 1 C b a = γ a γ b γ a γ b ( 1 + a b ) + 1 C b b = γ b 2 γ b + 1 γ a 1 γ a γ b ( 1 + a b ) + 1 .
Lemma 6.
For every vectors a , b B we have
g y r [ a E b , b ] = g y r [ a , b ] .
Proof. 
We have by (241),
g y r [ a E b , b ] = I + C ˜ ( a E b ) ( a E b ) ( a E b ) ( a E b ) + C ˜ ( a E b ) b ( a E b ) b + C ˜ b ( a E b ) b ( a E b ) + C ˜ b b b b
where
C ˜ ( a E b ) ( a E b ) = γ a E b 2 γ a E b + 1 γ b 1 γ a E b γ b ( 1 + ( a E b ) b ) + 1 C ˜ ( a E b ) b = γ a E b γ b + 2 ( a E b ) b γ a E b 2 γ b 2 ( γ a E b + 1 ) ( γ b + 1 ) γ a E b γ b ( 1 + ( a E b ) b ) + 1 C ˜ b ( a E b ) = γ a E b γ b γ a E b γ b ( 1 + ( a E b ) b ) + 1 C ˜ b b = γ b 2 γ b + 1 γ a E b 1 γ a E b γ b ( 1 + ( a E b ) b ) + 1 .
By straightforward calculations we have the following identities,
( a E b ) ( a E b ) = ( 1 + a b γ a γ a + 1 ) 2 ( 1 + a b ) 2 a a + 1 + a b γ a γ a + 1 1 + a b ( a b + b a ) + 1 γ a 2 ( 1 + a b ) 2 b b
( a E b ) b = 1 + a b γ a γ a + 1 1 + a b a b + 1 γ a ( 1 + a b ) b b
γ a E b = γ a γ b ( 1 + a b )
and
γ a E b γ b ( 1 + ( a E b ) b ) + 1 = γ b 2 γ a + 1 ( 1 + γ a ( 1 + a b ) ) 2 .
Hence, we have
C ˜ ( a E b ) ( a E b ) = γ a 2 ( 1 + a b ) 2 γ a γ b ( 1 + a b ) + 1 ( γ b 1 ) ( γ a 1 ) ( 1 + γ a ( a + a b ) ) 2
C ˜ ( a E b ) b = γ a ( 1 + a b ) [ ( γ a + 1 ) ( γ b 1 ) ( 1 + γ a 2 γ b ( 1 + a b ) ) + 2 γ b 2 ( 1 + γ a ( 1 + a b ) ) 2 ] ( 1 + γ a γ b 2 ( 1 + a b ) ) ( γ b + 1 ) ( 1 + γ a ( 1 + a b ) ) 2
C ˜ b ( a E b ) = γ a ( γ a + 1 ) ( 1 + a b ) ( 1 + γ a ( 1 + a b ) ) 2
and
C ˜ b b = ( γ a + 1 ) ( γ a γ b ( 1 + a b ) 1 ) ( γ b + 1 ) ( 1 + γ a ( 1 + a b ) ) 2 .
Therefore,
g y r [ a E b , b ] = I + D a a a a + D a b a b + D b a b a + D b b b b ,
where the coefficients are evaluated straightforwardly,
D a a = C ˜ ( a E b ) ( a E b ) ( 1 + a b γ a γ a + 1 ) 2 ( 1 + a b ) 2 = C a a D a b = C ˜ ( a E b ) ( a E b ) 1 + a b γ a γ a + 1 1 + a b + C ˜ ( a E b ) b 1 + a b γ a γ a + 1 1 + a b = C a b D b a = C ˜ ( a E b ) ( a E b ) 1 + a b γ a γ a + 1 1 + a b + C ˜ b ( a E b ) 1 + a b γ a γ a + 1 1 + a b = C b D b b = C ˜ ( a E b ) ( a E b ) 1 γ a 2 ( 1 + a b ) 2 + ( C ˜ ( a E b ) ( a E b ) + C ˜ b ( a E b ) ) 1 γ a ( 1 + a b ) + C ˜ b b = C b b .
Thus,
g y r [ a E b , b ] = I + C a a a a + C a b a b + C b a b a + C b b b b = g y r [ a , b ] .
The proof of the Lemma is, thus, complete. □

8.6. Gyrocommutative Law

We are now in the position to prove the gyrocommutative law (214) of Einstein addition.
Lemma 7.
For every vectors a , b B we have
a E b = g y r [ a , b ] ( b E a ) .
Proof. 
Equation (225) with c = b E a yields the equation
( a E b ) E d = a E ( b E ( b E a ) ,
where d = M c .
The value of d exists and is unique owing to the left cancellation law. Furthermore, we show that d = a E b satisfies (271). Indeed, for every vector v B we have
v E v = 1 1 + v v [ ( 1 + v v v 2 ) v ] = 2 1 + v 2 v .
In particular,
( a E b ) E ( a E b ) = 2 1 + a E b 2 ( a E b ) .
We now evaluate the right-hand side of (271), using the identities in (236) and in (223),
b ( b E a ) = b 2 + a b 1 + a b a E b 2 = 1 ( 1 a 2 ) ( 1 b 2 ) ( 1 a b ) 2 .
Then,
b E ( b E a ) = 1 1 + b ( b E a ) [ ( 1 + b ( b E a ) b 2 ) b + 1 b 2 ( b E a b ( b E a ) b 2 b ) ] = 1 + a b 1 + 2 a b | b 2 [ 1 1 + a b ( 1 + a b + b 2 + a b b 2 ) b + 1 b 2 1 + a b [ ( 1 + a b b 2 ) b + 1 b 2 ( a a b b 2 b ) b 2 + a b b 2 ] ] = 1 1 + 2 a b + b 2 ( 1 b 2 ) a + 2 ( 1 + a b ) b .
Using the formula
a ( b E ( b E a ) ) = ( 1 b 2 ) a 2 + 2 a b ( 1 + a b ) 1 + 2 a b + b 2 a ( b E ( b E a ) ) ,
we get
a E ( b E ( b E a ) ) = 1 1 + a ( b E ( b E a ) ) [ ( 1 + a ( b E ( b E a ) ) a 2 ) a + 1 a 2 [ b E ( b E a ) a ( b E ( b E a ) ) a 2 a ] ] = 1 1 + 2 a b + b 2 + ( 1 b 2 ) a 2 + 2 a b ( 1 + a b ) [ ( 1 + 2 a b + b 2 + [ ( 1 b 2 ) a 2 + 2 a b ( 1 + a b ) ] a + 1 a 2 [ ( 1 b 2 ) a + 2 ( 1 + a b ) b ( ( 1 b 2 ) a 2 + 2 a b ( 1 + a b ) ) a a 2 ] ] = 2 ( 1 + a b ) 2 ( 1 + a b ) 2 ( 1 a 2 ) ( 1 b 2 ) ( 1 + a b b 2 ) a + 1 a 2 ( b a b b 2 ) = 2 1 a E b 2 a E b = ( a E b ) E ( a E b ) .
Therefore, Equation (271) can be written as
( a b ) d = ( a b ) ( a b ) .
Now we add ( a b ) from the left to both sides of (278), obtaining
d = a b .
Recalling that d = M c = g y r [ a , b ] ( ( b a ) , we get the assertion (270) of the lemma. The proof of the Lemma is, thus, complete. □

8.7. Gyrations Preserve Einstein Addition and Multiplication

For any a , b B , the operator c g y r [ a , b ] c is linear, being given by a multiplication by a unitary matrix. Interestingly, it is also gyrolinear in the sense that
g y r [ a , b ] ( t 1 E c 1 ) E ( t 2 E c 2 ) = ( t 1 E ( g y r [ a , b ] c 1 ) ) E ( t 2 E ( g y r [ a , b ] c 2 ) )
for all c 1 , c 2 B and t 1 , t 2 [ , ) .
A result that generalizes (280) is stated in the following Lemma.
Lemma 8.
For every unitary matrix U, for every vectors c 1 , c 2 B , and every numbers t 1 , t 2 R we have
U ( t 1 E c 1 ) E ( t 2 E c 2 ) = ( t 1 E ( U c 1 ) ) E ( t 2 E ( U c 2 ) ) .
Proof. 
The proof follows from properties of Einstein addition and multiplication. For every c 1 , c 2 B
( U c 1 ) E ( U c 2 ) = 1 1 + ( U c 1 ) ( U c 2 ) { [ 1 + ( U c 1 ) ( U c 2 ) U c 1 2 ( 1 1 U c 1 2 ) ] U c 1 + 1 U c 1 2 U c 2 } = U 1 1 + c 1 c 2 { [ 1 + c 1 c 2 c 1 2 ( 1 1 c 1 2 ) ] c 1 + 1 c 1 2 c 2 } = U ( c 1 E c 2 ) ,
and for every vector c B and positive number t,
t E ( U c ) = atanh ( t tanh ( U c ) ) U c U c = U atanh ( t tanh ( c ) ) c c = U [ t E c ] .
The proof of the Lemma is, thus, complete. □
The following two corollaries are immediate consequences of Lemmas 8 and 4.
Corollary 2.
For every unitary matrix U, for every vectors a , b , c 1 , c 2 B , and every numbers t 1 , t 2 R we have
g y r [ a , b ] ( t 1 E c 1 ) E ( t 2 E c 2 ) = ( t 1 E ( g y r [ a , b ] c 1 ) ) E ( t 2 E ( g y r [ a , b ] c 2 ) ) .
Corollary 3.
For every vectors a , b , u , v B we have
( g y r [ u , v ] a ) ( g y r [ u , v ] b ) = a b .

8.8. Gyrogroups and Gyrovector Spaces

Assume E , E are the Einstein addition and scalar multiplication respectively.
Theorem 4.
The groupoid ( B , E ) is a gyrocommutative gyrogroup.
Proof. 
We need to check properties (G1)–(G6) in Definitions 1 and 2. Property (G1) follows from (45), property (G2) follows from (46), property (G3) follows from (240), property (G4) follows from (241), property (G5) follows from (256), and property (G6) follows from (270). □
Theorem 5.
The triple ( B , , ) is a gyrovector space.
Proof. 
We need to check properties (V0)–(V8) in Definition 3. Property (V0) follows from (285), property (V1) follows from (105), property (V2) follows from (103), property (V3) follows from (104), property (V4) follows from (106), property (V5) follows from (280), property (V6) follows from (253), property (V7) follows from (123), and property (V8) follows from (127). □

8.9. Einstein Coaddition

In Section 8.1 we have found a solution to the equation
a E x = b
in B for the unknown x.
We now consider the similar equation
x E a = b
in B for the unknown x.
It is known ([1] Chap. 1, Section 6) that the unique solution of (287) is given by
x = b E g y r [ b , a ] a .
Motivated by (288) we introduce the Einstein coaddition
a E b = a E g y r [ a , b ] b
and use the notation
a E b = a E ( b ) ,
so that the solution of equation (287) is given by
x = b E a .
Accordingly, in this section we calculate
a E b = a E g y r [ a , b ] b .
By the representation of g y r [ a , b ] in (241), we have
g y r [ a , b ] = I + C a a a a + C a b a b + C b a b a + C b b b b ,
where the numbers C a a , C a b , C b a and C b b are given by (242).
Let
Δ = 1 + a b + ( 1 a 2 ) ( 1 b 2 ) .
Then,
g y r [ a , b ] b = D a a + D b b ,
where
D a = C a a a b + C a b b 2 = 1 Δ b 2 + a b a 2 ( 1 1 a 2 ) ( 1 1 b 2 )
and
D b = 1 + C b a a b + C b b b 2 = 1 a 2 + 1 b 2 Δ .
We have
a c = 1 1 a c ( 1 a c a 2 ) a + 1 a 2 ( a c a 2 a c )
for any a , c B .
Calculating (298) with c = g y r [ a , b ] b , noticing that
1 a g y r [ a , b ] b = 1 Δ 1 a 2 b 2 + ( 1 a b ) 1 a 2 1 b 2 a g y r [ a , b ] b a 2 g y r [ a , b ] b = 1 a 2 + 1 b 2 Δ a b a 2 a b ,
we have
a E g y r [ a , b ] b = ( 1 a 2 + 1 b 2 ) ( a 1 b 2 b 1 a 2 ) 1 a 2 b 2 + ( 1 a b ) 1 a 2 1 b 2 .
Using the standard notation for the Lorentz factor in (254), we have the equation
a b = γ a a + γ b b γ a + γ b γ a b + 1 γ a + γ b ,
which is well known ([5] p. 92). An extension of (301) to more than two summands is presented in ([5] p. 425). It is clear from (301) that Einstein coaddition in B is commutative,
a b = b a .
In order to calculate the metric tensor for the binary operation ⊞ in B we obtain the equation
( x ) ( x + Δ x ) = I + x x 1 x 2 Δ x + o ( Δ x .
Hence, the metric tensor for the binary operation ⊞ has the form
G ( x ) = I + x x 1 x 2 2 = I x x x 2 + 1 ( 1 x 2 ) 2 x x x 2 .
The metric tensor G ( x ) possesses the form of (6) with
m 0 ( x 2 ) = 1 m 1 ( x 2 ) = 1 ( 1 x 2 ) 2 .
We should note that the function m 1 is the same as for Einstein addition and Möbius addition. Hence, the distances for all these three binary operations are the same (115):
b E ( a ) E = b M ( a ) M = b ( a ) = atanh ( b a ) ,
where · is the Euclidean norm. Moreover, the multiplication by a number is the same for all these three binary operations,
t E a = t M a = t a .
We can use also the standard Euclidean norm · . Notice that for numbers (that is, for n = 1 ) all binary operations coincide,
x E y = x M y = x y = x + y 1 + x y = tanh ( atanh ( x ) + atanh ( y ) ) .
Therefore, if ⊕ is one of the operations E , M , ⊞, then for all a , b B ,
a b a b ,
and the equality occurs if and only if for some nonnegative number t we have a = t b or b = t a .

Author Contributions

All authors contributed equally to this article. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ungar, A.A. Beyond the Einstein Addition Law and Its Gyroscopic Thomas Precession: The Theory of Gyrogroups and Gyrovector Spaces; Fundamental Theories of Physics, Kluwer Academic Publishers Group: Dordrecht, The Netherlands, 2001; Volume 117. [Google Scholar]
  2. Ungar, A.A. From Möbius to gyrogroups. Am. Math. Mon. 2008, 115, 138–144. [Google Scholar] [CrossRef]
  3. Ungar, A.A. The intrinsic beauty, harmony and interdisciplinarity in Einstein velocity addition law: Gyrogroups and gyrovector spaces. Math. Interdisc. Res. 2016, 1, 5–51. [Google Scholar]
  4. Ungar, A.A. Analytic Hyperbolic Geometry: Mathematical Foundations and Applications; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2005. [Google Scholar]
  5. Ungar, A.A. Analytic Hyperbolic Geometry and Albert Einstein’s Special Theory of Relativity; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2008. [Google Scholar]
  6. Ungar, A.A. A Gyrovector Space Approach to Hyperbolic Geometry; Morgan & Claypool Pub.: San Rafael, CA, USA, 2009. [Google Scholar]
  7. Ungar, A.A. Barycentric Calculus in Euclidean and Hyperbolic Geometry: A Comparative Introduction; World Scientific Publishing Co. Pte. Ltd.: Hackensack, NJ, USA, 2010. [Google Scholar]
  8. Ungar, A.A. Hyperbolic Triangle Centers: The Special Relativistic Approach; Springer: New York, NY, USA, 2010. [Google Scholar]
  9. Ungar, A.A. Analytic Hyperbolic Geometry in n Dimensions: An Introduction; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  10. Ungar, A.A. Beyond Pseudo-Rotations in Pseudo-Euclidean Spaces: An Introduction to the Theory of Bi-Gyrogroups and Bi-Gyrovector Spaces: Mathematical Analysis and its Applications; Elsevier: Amsterdam, The Netherlands; Academic Press: Cambridge, MA, USA, 2018. [Google Scholar]
  11. Rassias, T.M. Book Review: Analytic Hyperbolic Geometry and Albert Einstein’s Special Theory of Relativity, by Abraham A. Ungar. Nonlinear Funct. Anal. Appl. 2008, 13, 167–177. [Google Scholar]
  12. Rassias, T.M. Book Review: A gyrovector space approach to hyperbolic geometry, by Abraham A. Ungar. J. Geom. Symmetry Phys. 2010, 18, 93–106. [Google Scholar]
  13. Walter, S. Book Review: Beyond the Einstein Addition Law and its Gyroscopic Thomas Precession: The Theory of Gyrogroups and Gyrovector Spaces, by Abraham A. Ungar. Found. Phys. 2002, 32, 327–330. [Google Scholar] [CrossRef]
  14. Lang, S. Linear Algebra; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  15. Ungar, A.A. Thomas rotation and the parametrization of the Lorentz transformation group. Found. Phys. Lett. 1988, 1, 57–89. [Google Scholar] [CrossRef]
  16. Abe, T. Gyrometric preserving maps on Einstein gyrogroups, Möbius gyrogroups and proper velocity gyrogroups. Nonlinear Funct. Anal. Appl. 2014, 19, 1–17. [Google Scholar]
  17. Ferreira, M. Harmonic analysis on the Einstein gyrogroup. J. Geom. Symmetry Phys. 2014, 35, 21–60. [Google Scholar]
  18. Frenkel, P.E. On endomorphisms of the Einstein gyrogroup in arbitrary dimension. J. Math. Phys. 2016, 57, 032301. [Google Scholar] [CrossRef]
  19. Molnár, L.; Virosztek, D. On algebraic endomorphisms of the Einstein gyrogroup. J. Math. Phys. 2015, 56, 082302. [Google Scholar] [CrossRef]
  20. Suksumran, T.; Wiboonton, K. Einstein gyrogroup as a B-loop. Rep. Math. Phys. 2015, 76, 63–74. [Google Scholar] [CrossRef]
  21. Urribarri, D.K.; Castro, S.M.; Martig, S.R. Gyrolayout: A hyperbolic level-of-detail tree layout. J. Univers. Comput. Sci. 2013, 19, 132–156. [Google Scholar]
  22. Ferreira, M. Factorizations of Möbius gyrogroups. Adv. Appl. Clifford Algebr. 2009, 19, 303–323. [Google Scholar] [CrossRef] [Green Version]
  23. Ferreira, M. Harmonic Analysis on the Möbius Gyrogroup. J. Fourier Anal. Appl. 2015, 21, 281–317. [Google Scholar] [CrossRef]
  24. Ferreira, M.; Ren, G. Möbius gyrogroups: A Clifford algebra approach. J. Algebra 2011, 328, 230–253. [Google Scholar] [CrossRef] [Green Version]
  25. Atiponrat, W. Topological gyrogroups: Generalization of topological groups. Topol. Appl. 2017, 224, 73–82. [Google Scholar] [CrossRef]
  26. Lal, R.; Yadav, A.C. Topological Right Gyrogroups and Gyrotransversals. Commun. Algebra 2013, 41, 3559–3575. [Google Scholar] [CrossRef]
  27. Suksumran, T. The algebra of gyrogroups: Cayley’s theorem, Lagrange’s theorem and isomorphism theorems. In Essays in Mathematics and Its Applications: In Honor of Vladimir Arnold; Springer: New York, NY, USA, 2016; pp. 369–437. [Google Scholar]
  28. Suksumran, T. Gyrogroup actions: A generalization of group actions. J. Algebra 2016, 454, 70–91. [Google Scholar] [CrossRef]
  29. Suksumran, T. Involutive groups, unique 2-divisibility, and related gyrogroup structures. J. Algebra Appl. 2017, 16, 1750114. [Google Scholar] [CrossRef]
  30. Suksumran, T.; Ungar, A.A. Gyrogroups and the Cauchy property. Quasigroups Relat. Syst. 2016, 24, 277–286. [Google Scholar]
  31. Suksumran, T.; Wiboonton, K. Isomorphism theorems for gyrogroups and L-subgyrogroups. J. Geom. Symmetry Phys. 2015, 37, 67–83. [Google Scholar]
  32. Suksumran, T.; Wiboonton, K. Lagrange’s theorem for gyrogroups and the cauchy property. Quasigroups Relat. Syst. 2015, 22, 283–294. [Google Scholar]
  33. Ferreira, M. Harmonic analysis on the proper velocity gyrogroup. Banach J. Math. Anal. 2017, 11, 21–49. [Google Scholar] [CrossRef]
  34. Abe, T.; Hatori, O. On a characterization of commutativity for C*-algebras via gyrogroup operations. Period. Math. Hungar. 2016, 72, 248–251. [Google Scholar] [CrossRef]
  35. Ferreira, M. Gyrogroups in projective hyperbolic Clifford analysis. In Hypercomplex Analysis and Applications, Trends Math.; Birkhäuser/Springer Basel AG: Basel, Switzerland, 2011; pp. 61–80. [Google Scholar]
  36. Jaiyéọlá, T.G.; Sòlárìn, A.R.T.; Adéníran, J.O. Some Bol-Moufang characterization of the Thomas precession of a gyrogroup. Algebr. Groups Geom. 2014, 31, 341–362. [Google Scholar]
  37. Lal, R.; Yadav, A.C. Twisted automorphisms and twisted right gyrogroups. Commun. Algebra 2015, 43, 3442–3458. [Google Scholar] [CrossRef]
  38. Park, J.; Kim, S. Hilbert projective metric on a gyrogroup of qubit density matrices. Rep. Math. Phys. 2015, 76, 389–400. [Google Scholar] [CrossRef]
  39. Tărnxaxuceanu, M. A note on a class of gyrogroups. Quasigroups Relat. Syst. 2017, 25, 151–154. [Google Scholar]
  40. Ungar, A.A. Symmetry groups of systems of entangled particles. J. Geom. Symmetry Phys. 2018, 48, 47–77. [Google Scholar] [CrossRef]
  41. Demirel, O. The first sharp gyrotriangle inequality in Möbius gyrovector space (𝔻,⊕,⊗). Forum Geom. 2017, 17, 439–447. [Google Scholar]
  42. Demirel, O.; Seyrantepe, E.S. The cogyrolines of Möbius gyrovector spaces are metric but not periodic. Aequationes Math. 2013, 85, 185–200. [Google Scholar] [CrossRef]
  43. Ungar, A.A. Möbius gyrovector spaces in quantum information and computation. Comment. Math. Univ. Carolin. 2008, 49, 341–356. [Google Scholar]
  44. Watanabe, K. A confirmation by hand calculation that the Möbius ball is a gyrovector space. Nihonkai Math. J. 2016, 27, 99–115. [Google Scholar]
  45. Watanabe, K. Orthogonal gyroexpansion in Möbius gyrovector spaces. J. Funct. Spaces 2017, 13, 1518254. [Google Scholar] [CrossRef] [Green Version]
  46. Abe, T.; Hatori, O. Generalized gyrovector spaces and a Mazur-Ulam theorem. Publ. Math. Debr. 2015, 87, 393–413. [Google Scholar] [CrossRef]
  47. Hatori, O. Examples and applications of generalized gyrovector spaces. Results Math. 2017, 71, 295–317. [Google Scholar] [CrossRef]
  48. Hatori, O. Extension of isometries in generalized gyrovector spaces of the positive cones. Contemp. Math. Am. Math. Soc. 2017, 687, 145–156. [Google Scholar]
  49. Ungar, A.A. Gyrovector spaces and their differential geometry. Nonlinear Funct. Anal. Appl. 2005, 10, 791–834. [Google Scholar]
  50. Abe, T.; Watanabe, K. Finitely generated gyrovector subspaces and orthogonal gyrodecomposition in the Möbius gyrovector space. J. Math. Anal. Appl. 2017, 449, 77–90. [Google Scholar] [CrossRef]
  51. Chen, J.-L.; Ungar, A.A. From the group sl(2,C) to gyrogroups and gyrovector spaces and hyperbolic geometry. Found. Phys. 2001, 31, 1611–1639. [Google Scholar] [CrossRef]
  52. Chen, J.-L.; Ungar, A.A. The Bloch gyrovector. Found. Phys. 2001, 32, 531–565. [Google Scholar] [CrossRef]
  53. Kasparian, A.K.; Ungar, A.A. Lie gyrovector spaces. J. Geom. Symmetry Phys. 2004, 1, 3–53. [Google Scholar]
  54. Kim, S. Distributivity on the Gyrovector spaces. Kyungpook Math. J. 2015, 55, 13–20. [Google Scholar] [CrossRef] [Green Version]
  55. Ungar, A.A. Möbius transformations of the ball, Ahlfors’ rotation and gyrovector spaces. In Nonlinear Analysis in Geometry and Topology; Rassias, T.M., Ed.; Hadronic Press: Palm Harbor, FL, USA, 2000; pp. 241–287. [Google Scholar]
  56. Novikov, S.P.; Fomenko, A.T. Elements of Differential Geometry and Topology; Kluwer Academic Publishers: Boston, MA, USA, 1990. [Google Scholar]

Share and Cite

MDPI and ACS Style

Barabanov, N.E.; Ungar, A.A. Binary Operations in the Unit Ball: A Differential Geometry Approach. Symmetry 2020, 12, 1178. https://doi.org/10.3390/sym12071178

AMA Style

Barabanov NE, Ungar AA. Binary Operations in the Unit Ball: A Differential Geometry Approach. Symmetry. 2020; 12(7):1178. https://doi.org/10.3390/sym12071178

Chicago/Turabian Style

Barabanov, Nikita E., and Abraham A. Ungar. 2020. "Binary Operations in the Unit Ball: A Differential Geometry Approach" Symmetry 12, no. 7: 1178. https://doi.org/10.3390/sym12071178

APA Style

Barabanov, N. E., & Ungar, A. A. (2020). Binary Operations in the Unit Ball: A Differential Geometry Approach. Symmetry, 12(7), 1178. https://doi.org/10.3390/sym12071178

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