# How Nonassociative Geometry Describes a Discrete Spacetime

- CUCEI, Department of Physics, University of Guadalajara, Guadalajara, Mexico

Nonassociative geometry, providing a unified description of discrete and continuum spaces, is a valuable candidate for the study of discrete models of spacetime. Within the framework of nonassociative geometry we propose a model of emergent spacetime. In our approach, the evolution of spacetime geometry is governed by a random/stochastic process. This leads to a natural appearance of causal structure and arrow of time. We apply our approach to study a toy model of discrete (2+1)-D spacetime and Friedmann-Robertson-Walker cosmological model. We show that in a continuous limit the evolution of the discrete spacetime corresponds to a radiation epoch of the standard cosmological model.

## 1. Introduction

Any attempt to create a quantum gravitation theory faces the challenge of comprehension the concept of quantizing of spacetime and of describing the quantum nature of space and time [1–3]. As part of the general problem of time, one of the most challenging and open problems is the origin of irreversibility in our universe, also known as the problem of the arrow of time [1, 4–9].

One can distinguish two general strategies to achieve this goals [10, 11]. The first strategy consists in quantizing a classical structure, that later is recovered as a limit of the quantum theory. The second strategy assumes that the classical structures are emergent from the other theory, which is more fundamental theory from the very beginning. In the second approach, a formulation of the quantum theory may require omitting use of continuum concepts *a priori*. This means that at the Planck scale the standard concept of spacetime must be replaced by some discrete structure [10, 12].

The proposal that spacetime may be fundamentally discrete has been adopted by numerous strategies [13, 14]. This discreteness is motivated by several heuristic arguments, but so far has not been deduced rigorously. In Nesterov and Sabinin [15–18], we have proposed a new unified approach, based on *nonassociative geometry*, for describing both continuum and discrete spacetime. Among advanced models that propose discreteness, two that will be relevant to our work are Causal Sets (CS)[19–21] and Causal Dynamical Triangulations (CDT)[14, 22–28]. These approaches are based on the hypothesis that spacetime is discrete and causality is a fundamental principle.

Nonassociative geometry, being based on the theory of quasigroups and loops, provides a unified algebraic description of manifolds. In the case of smooth manifolds it is equivalent to conventional differential geometry [29, 30]. The main algebraic structures emerging in nonassociative geometry may be described as follows. Let 𝔐 be a manifold with an affine connection. Then, in a neighborhood of an arbitrary point *a* ∈ 𝔐 on one can introduce the geodesic local loop ${Q}$_{a} with a smooth partial ternary operation, ${L}_{x}^{a}y$, *a* being the neutral element of ${Q}$_{a}, and *a, x, y* ∈ 𝔐. The operation ${L}_{x}^{a}y$ is uniquely defined through the parallel transport of the geodesic 〈*ay*〉, connecting the point *a* with the point *y*, along the geodesic 〈*ax*〉. The family of local loops {${Q}$_{a}} constructed in this way satisfies some natural algebraic identities and uniquely defines the space with affine connection. Next, one can introduce the left-invariant metric, ${g}^{b}({L}_{b}^{a}x,{L}_{b}^{a}y)={g}^{a}(x,y)$, enriching the manifold 𝔐 with a *metric structure*. In nonassociative geometry, the (nonlocal) equivalent of curvature is elementary holonomy, ${h}_{(x,y)}^{a}={L}_{a}^{y}\circ {L}_{y}^{x}\circ {L}_{x}^{a}$.

In the present paper, we introduce a discrete model of emergent spacetime based on nonassociative geometry. In our approach events, being “atoms” of spacetime, are elements of a local loop with a partial operation. We assume that time is quantized and the distance between two neighboring atoms can not be less than the fundamental length, ℓ_{p}. The evolution of spacetime geometry is governed by a random/stochastic process conserving spacetime geometry. The result is a partially ordered set of events, and spacetime geometry is encoded in the nonassociative structure of spacetime. We apply our approach to study a toy model of (2 + 1)-D discrete spacetime and a discrete Friedmann-Robertson-Walker (FRW) cosmological model. We show that the evolution of the discrete spacetime corresponds to a radiation era in the continuous limit with the scale factor growing as $~\sqrt{t}$. The irreversibility of evolution is governed by a random process and by information loss as a consequence of the fundamental Planck length.

The paper is organized as follows. In section 2 we survey the foundations of the theory of smooth quasigroups and loops. In section 3, we give a brief introduction to nonassociative geometry, with detailed examples and applications to the de Sitter and Friedmann-Robertson-Walker (FRW) models of spacetime. In section 4 we introduce nonassociative discrete geometry. In section 5 we present a discrete model of spacetime based on nonassociative geometry with applications to (2+1)-D spacetime and the FRW model. In section 6 we discuss the results of numerical simulations and theoretical models for the implementation of our framework. In the Supplemental Material, we present technical details of our method and algorithms.

## 2. Smooth Quasigroups and Loops

During the last five decades, nonassociative structures have appeared in various fields of physics, and nowadays they play an important role in modern physics. Among others, one may mention the employment of octonions, Lie groupoids, and algebroids in the context of Yang-Mills theories and gauge theories on commutative but nonassociative fuzzy spaces [31–37], and attempts to create nonassociative quantum mechanics [38–42].

The nonassociative objects, such as 3-cocycles, emerge in quantum mechanics with the Dirac monopole and quantum field theory with chiral anomalies [43–55]. The chiral anomalies appear as the Schwinger terms in current algebra [53–56]. These terms, having a cohomological nature, are related to the failure of the Jacobi identity, *J*(*A, B, C*) = 0, where for a given triple of operators (*A, B, C*) we define

In particular, it has been shown that electric fields in chiral gauge theories satisfy the fourth-order identity, *J*(*E*_{i}, *E*_{j}, [*E*_{k}, *E*_{i}]) = [*J*(*E*_{i}, *E*_{j}, *E*_{k}), *E*_{i}], instead of the Jacobi identity, and form the so-called Mal'cev algebra [53–57]. The Mal'cev algebra, being a tangent algebra of the analytic Moufang loop, generalizes the concept of Lie algebra, and for a given tangent Mal'cev algebra there exists a unique simply connected analytical Moufang loop.

Failure of the Jacobi identity implies violation of associativity and forces one to make use of nonassociative generalization of the Lie groups, such as smooth quasigroups and loops. The nonassociative algebraic structures, such as smooth quasigroups and loops, are of considerable potential interest for modern physics. Recently, the theory of smooth loops has been employed for description of Thomas precession, coherent states and geometric phases in quantum mechanics [58–61]. In general relativity, unlike special relativity, the Poincaré group does not exist in general spacetime. However, this problem can be overcome by introducing the Poincaré quasigroup [62, 63].

The theory of fiber bundles, being one of the important applications of Lie groups not only in mathematics but in physics as well, leads to gauge field theory. The nonassociative generalization of Lie groups allows us to construct nonassociative gauge theories with higher “nonassociative gauge symmetries" by using the theory of smooth loops [64–66]. In general terms, a consequence of nonassociativity is that the structure constants of the gauge algebra have to be changed by structure functions.

In what follows we outline the theory of smooth quasigroups and loops. While there is not a consensus about the birth date of the theory of smooth loops, one can consider the seminal work by Mal'cev [57] as a starting point for the development of this theory. General results on the subject, including reviews, may be found in [30, 67–69], Kikkawa [70], Belousov [71], Pflugfelder [72], and Batalin [73] (see also references therein). The point of departure is an identity of quasiassociativity leading to the infinitesimal theory of smooth loops.

**Definition 2.1**. Let 〈*Q*, ·〉 be a groupoid, i.e., a set with a binary operation (*a, b*) ↦ *a* · *b*. A groupoid 〈*Q*, ·〉 is called a *quasigroup* if each of the equations *a*·*x* = *b*, and *y*·*a* = *b* has a unique solution, *x* = *a*\*b* and *y* = *b*/*a*, defining left (“\”) and right (“/”) division. A *loop* is a quasigroup with a two-sided identity, *a*·*e* = *e*·*a* = *a*, where *e* is a neutral element. A loop that is also a differentiable manifold and in which the operation *a*·*b* is a smooth map is called a *smooth loop*.

One can show that the following identities hold:

where *l*_{(a, b)} is a *left associator*, *r*_{(b, c)} is a *right associator*, *L*_{a} is a *left translation* and *R*_{b} is a *right translation*. The left and right translations are defined as follows: *L*_{a}*b* = *R*_{b}*a* = *a*·*b*. As one can see, the loop becomes a group, if the associator is the identity operator. Thus, the associator can be considered as a measure of nonassociativity of the operation in the loop.

Let *T*_{e}(*Q*) be the tangent space at the neutral element *e* ∈ *Q*. Then, for each *X*_{e} ∈ *T*_{e}(*Q*), one can define a smooth vector field on *Q* as

where *L*_{b*} : *T*_{e}(*Q*) ↦ *T*_{b}(*Q*) is the differential of the left translation.

**Definition 2.2**. A vector field *X* on *Q* which satisfies the relation ${L}_{a*}{X}_{b}={\widehat{l}}_{\left(a,b\right)*}{X}_{a\xb7b}$ for any *a, b* ∈ *Q* is called a *left fundamental* or *left quasi-invariant* vector field.

Let ${\text{\Gamma}}_{i}={R}_{i}^{j}\partial /\partial {a}^{j},i=1,2,\dots ,r$, be a basis of left quasi-invariant vector fields. Then, for the commutator, [Γ_{i}, Γ_{j}], we obtain

where ${C}_{ij}^{p}(a)$ are the *structure functions*. They satisfy the modified Jacobi identity,

**Comment 2.1**. The infinitesimal theory of Lie groups is based on the associativity of the operation, *a*·(*b*·*c*) = (*a*·*b*)·*c*. In a similar way, the quasiassociativity identity,

yield the infinitesimal theory of smooth loops [66].

### 2.1. Examples

*Example* 2.1. **Loop 〈 Q**(ℝ/ℤ),

***〉**. The analytical loop

*Q*(ℝ/ℤ) is a loop with the binary operation given by [74, 75]

where *f*(*x*) = (1 − cos 2π*x*)/4.

*Example* 2.2. **Loop Q**ℂ. Let ℂ be the complex plane. It becomes the loop Qℂ, if we define the following nonassociative operation:

where the bar denotes complex conjugation. The associator *l*_{(ζ, η)} is found to be

The basis of left quasi-invariant vectors and 1-forms is given by

The computation of the commutator results in:

and for the structure functions we find that ${C}_{12}^{1}=-\eta $, ${C}_{12}^{2}=\stackrel{\u0304}{\eta}$.

One can show that the metric based on the left fundamental basis forms, *ds*^{2} = θ^{1}θ^{2}, is a left invariant metric. As one can see, this is also the Fubini-Studi metric:

**Comment 2.2**. The loop Qℂ is locally isomorphic to the Bol loop *QS*^{2} [16, 30]. The latter is related to the two-sphere *S*^{2}, which admits a natural quasigroup structure induced by the operation (8). The isomorphism is established by using the stereographic projection from the south pole of the sphere: ζ = *e*^{iφ} tan(θ/2).

*Example* 2.3. **Loop QSU(2)**. Consider a set of the unitary matrices

*QSU*(2) ⊂

*SU*(2). An arbitrary element

*U*

_{η}∈

*QSU*(2) is given by

We introduce the binary (nonassociative) operation on *QSU*(2) as follows:

where

The matrix loop *QSU*(2) forms the nonassociative representation of the loop *Q*ℂ ≃ *QS*^{2}. Indeed, writing the result of multiplication of two elements, *U*_{η}, *U*_{ζ} ∈ *QSU*(2), as

we find (see example 2.2)

*Example* 2.4. **Loop QH^{2}**. Let

*D*⊂ ℂ be the open unit disk,

*D*= {ζ ∈ ℂ : |ζ| < 1}. Then, the set of complex numbers with the operation [58, 59, 64, 65]

forms the two-sided loop *QH*^{2}. The associator, *l*_{(ζ, η)}, on *QH*^{2} is defined by

The basis of left quasi-invariant vectors and 1-forms is given by

The computation of the commutator yields

Using the left fundamental basis forms, we obtain the hyperboloic (left invariant) metric

The loop *QH*^{2} is isomorphic to the geodesic loop of two-dimensional Lobachevsky space. The isomorphism is established by using the relation ζ = *e*^{iφ} tanh(θ/2), where (θ, φ) are inner coordinates on the upper part of two-sheeted hyperboloid, *H*^{2}.

*Example* 2.5. **Loop QSU(1, 1)**. Consider a set of the unitary matrices

*QSU*(1, 1) ⊂

*SU*(1, 1). An arbitrary element,

*U*

_{η}∈

*QSU*(1, 1), is given by

where *D* ⊂ ℂ is the open unit disk, *D* = {ζ ∈ ℂ : |ζ| < 1}. We introduce the nonassociative operation as follows:

where

The loop *QSU*(1, 1) forms the nonassociative representation of the loop *QH*^{2}. Indeed, writing Equation (26) as

we find (see Example 2.4)

*Example* 2.6. **Loop QH**_{ℝ} (*De Sitter loop*). Let us consider the algebra of quaternions, H_{${C}$}, over a complex field ℂ(1, i),

The multiplication of quaternions is defined by the following rules for (*i, j, k*):

For a given quaternion, *q* = α + β*i* + γ*j* + δ*k*, the quaternionic conjugation is defined as

Further, we restrict ourselves to consideration of the set of quaternions H_{ℝ} ⊂ H_{${C}$}:

with the norm defined as ∥ζ∥^{2} = ζζ ^{+} = (ζ^{0})^{2} − (ζ^{1})^{2} − (ζ^{3})^{2} − (ζ^{4})^{2}.

We introduce the binary operation

where / denotes right division and *K* is constant. The set of quaternions H_{ℝ} with the binary operation (77) forms the loop QH_{ℝ}. The associator is given by

**Comment 2.3**. For *K* > 0 the loop QH_{ℝ} is isomorphic to the geodesic loop of spacetime of constant positive curvature (*de Sitter space*) [16].

*Example* 2.7. **Loop** $\mathrm{\text{Q}}{\text{H}}_{\text{\mathbb{R}}}^{{F}}$. Consider the set of quaternions:

where *a*(ζ^{0}) is a smooth function. We introduce a binary operation ${L}_{x}y=x*y,x,y\in {{\text{H}}}_{\text{\mathbb{R}}}^{{F}}$ as

where *x* = ζ^{0} + *a*(ζ^{0})**ζ**, *y* = η^{0} + *a*(η^{0})**η**, and / denotes the right division. It is easy to show that *e* = 0 is the neutral element, so that *L*_{e}*x* = *x*. The associator is defined by

where we set *z* = ξ^{0} + *a*(ξ^{0})**ξ**. The set of quaternions, ${\text{H}}_{\text{\mathbb{R}}}^{{F}}$ with the binary operation *, forms a loop, which will be denoted as $\mathrm{\text{Q}}{\text{H}}_{\text{\mathbb{R}}}^{{F}}$.

*Example* 2.8. **Moufang loop**. A loop *Q* is a Moufang loop if, for all *a, b, c* ∈ *Q*, the following identities hold:

In particular, the sphere *S*^{7}, identified with the unit octonions, may be considered a smooth Moufang loop with the binary operation defined by the octonionic algebra. Other examples of the analytical Moufang loops associated with octonions are the spaces *S*^{3} × *R*^{4} and *S*^{7} × *R*^{7}.

## 3. Nonassociative Geometry in Brief

In this section we outline the algebraic foundations of *nonassociative geometry* (For details see [29, 30, 76–78]).

**Definition 3.1**. Let 〈𝔐, ·, *e*〉 be a partial *groupoid* with a binary operation (*x, y*) ↦ *x*·*y* and a neutral element *e, x*·*e* = *e*·*x* = *x*. Consider 𝔐 to be a smooth manifold (at least *C*^{1}-smooth) and the operation of multiplication (at least *C*^{1}-smooth) to be defined in some neighborhood *U*_{e}. Then, 〈𝔐, ·, *e*〉 is called a *partial loop on* 𝔐.

*Remark* 3.1. The operation of multiplication is locally invertible. This implies that in some neighborhood of the neutral element *e* one has

where *x*·*y* = *L*_{x}*y* = *R*_{y}*x*,

The vector fields *A*_{j} and *B*_{j} defined on *U*_{e} as

are called the *left basic fundamental fields* and *right basic fundamental fields*, respectively.

The solution of the equation

can be written as *f*(*t*) = Exp *tX* and defines the *exponential map*,

Employing the exponential map, one can define a new operation, *tx* = Exp(*t*Exp^{−1}*x*), called the *left canonical unary operation*.

**Definition 3.2**. A smooth loop 〈𝔐, ·, *e*〉 equipped with its canonical left unary operations is called the *left canonical preodule* 〈𝔐, ·, (*t*)_{t ∈ ℝ}, *e*〉.

Introducing one more operation,

we obtain the *canonical left prediodule* of a loop, denoted as 〈𝔐, ·, +, (*t*)_{t ∈ ℝ}, *e*〉.

**Definition 3.3**. A canonical left preodule (prediodule) is called the *left odule (diodule)* if the *monoassociativity* property, *tx*·*ux* = (*t* + *u*)*x* is satisfied.

**Definition 3.4**. Let 𝔐 be a smooth manifold and *L*:(*x, y, z*) ∈ 𝔐 ↦ *L*(*x, y, z*) ∈ 𝔐 be a smooth partial ternary operation, such that *x _{ȧ}y* =

*L*(

*x, a, z*) defines in some neighborhood of the point

*a*the loop with the neutral

*a*. Then, the pair 〈𝔐,

*L*〉 is called a

*loopuscular structure (manifold)*.

**Definition 3.5**. A smooth manifold 𝔐 with a smooth partial ternary operation *L* and smooth binary operations ω_{t}:(*a, b*) ∈ 𝔐 × 𝔐 ↦ ω_{t}(*a, b*) = *t*_{a}*b* ∈ 𝔐, (*t* ∈ ℝ), such that *x _{ȧ}y* =

*L*(

*x, a, y*) and

*t*

_{a}

*z*= ω

_{t}(

*a, z*) determine in some neighborhood of an arbitrary point

*a*the odule with the neutral element

*a*, is called a

*left odular structure (manifold)*〈𝔐,

*L*, (ω

_{t})

_{t ∈ ℝ}〉.

Let 〈𝔐, *L*, (ω_{t})_{t ∈ ℝ}〉 and 〈𝔐, *N*, (ω_{t})_{t ∈ ℝ}〉 be odular structures. Then, we define a *diodular structure (manifold)* as 〈𝔐, *L, N*, (ω_{t})_{t ∈ ℝ}〉. If, in addition, ${x}_{\stackrel{+}{a}}y=N(x,a,y)$ and *t*_{a}*x* = ω_{t}(*a, x*) define a vector space, the diodular structure is called a *linear diodular structure*.

A diodular structure is said to be *geodiodular* if the following identities hold:

**Definition 3.6**. Let 𝔐 be a *C*^{k}-smooth (*k* ≥ 3) affinely connected manifold, with the following operations given on 𝔐 (locally):

where Exp_{x} is the exponential map at the point *x*, and ${\tau}_{x}^{a}$ is the parallel transport along the geodesic connecting the points *a* and *x*. This construction equips a manifold 𝔐 with the linear geodiodular structure, called a *natural linear geodiodular structure of an affinely connected manifold* (𝔐, ∇). In Figure 1, the main operation (48) – (50) characterizing nonassociative geometry are illustrated.

**Figure 1**. Main operations describing nonassociative geometry: **(A)** Multiplication of a “vector” (geodesic) by a scalar (parameter *t*); **(B)** Parallel translation of the geodesic (*ay*) along the geodesic (*ax*); **(C)** Addition of “vectors" (geodesics on the manifold 𝔐); **(D)** Parallel translation of the geodesic (*ax*) along the geodesic triangle (*abc*).

*Remark* 3.2. Any *C*^{k}-smooth (*k* ≥ 3) manifold with affine connection can be considered as a geodiodular structure. A geodiodular space ${{\text{M}}}_{a}^{+}=\langle \U0001d510{,}_{\stackrel{+}{a}},a,{({t}_{a})}_{t\in \text{\mathbb{R}}}\rangle $ plays the role of a tangent space.

**Definition 3.7**. Let 〈𝔐, *L*〉 be a loopuscular structure of a smooth manifold 𝔐. Then the formula

*Y* being a vector field in the neighborhood of a point *a*, defines the *tangent affine connection*. In coordinates, the components of the affine connection are

**Definition 3.8**. Let 〈𝔐, *L*〉 be a loopuscular structure. Then,

is called the *elementary holonomy*.

The elementary holonomy describes the parallel transport along a geodesic triangle path (see Figure 1). Thus, it can be considered as some integral curvature. In the smooth case, differentiating ${({h}_{(x,y)}^{a})}^{i}$ at the point *a* ∈ 𝔐, we obtain (up to a numerical factor) the curvature tensor as

*Remark* 3.3. The curvature of a nonassociative space implies non-trivial elementary holonomy, *h*^{a}(*x, y*) ≠ id. If *h*^{a}(*x, y*) = id for any *a, x, y* ∈ 𝔐, we obtain a space with absolute parallelism. If in addition, the relationship ${x}_{\stackrel{+}{a}}y={x}_{\stackrel{\xb7}{a}}y$ holds, that means the absence of torsion, and the space 𝔐 becomes a classical affine space (see chapter 3 in [30]).

Employing (54), one can show that the elementary holonomy satisfies the so-called *odular Bianchi identities*:

These identities can be considered as a non-local form of the usual Bianchi identities. In the linear approximation, from (56) we obtain the Bianchi identities in their conventional form.

### 3.1. Examples

*Example* 3.1. **Two-dimensional sphere ${S}_{R}^{2}$**. The two-sphere ${S}_{R}^{2}$ with a radius *R* can be equipped by a natural loop structure, that may be described as follows [58, 59, 64, 65]. Let ℂ be a complex plane and ζ, η ∈ ℂ. Let us consider the loop formed by the complex numbers with the non-associative multiplication (left translation) given by

This loop is isomorphic to the local loop related to ${S}_{R}^{2}$. The isomorphism between points of the sphere, ${S}_{R}^{2}$, and of the complex plane ℂ can be established by using the stereographic projection from the south pole of the sphere, writing ζ = *R* tan(θ/2)*e*^{iφ}.

The computation of the associator yields

One can show that for the sphere the elementary holonomy is determined by the associator, ${h}_{(\zeta ,\eta )}={l}_{(\zeta ,{L}_{\zeta}^{-1}\eta )}$ [30, 76]. We obtain

The left invariant metric on ${S}_{R}^{2}$ can be defined as

The geodesic distance between two points on ${S}_{R}^{2}$ is related to the left invariant metric as follows:

where

For *d*_{12}/*R* ≪ 1 we obtain *d*_{12} ≈ ℓ_{12}. Using spherical coordinates, one can show that (61) can be recast as the cosine rule in spherical trigonometry,

where θ_{12} = *d*_{12}/*R*, ${\theta}_{1}=2{tan}^{-1}\left|{\xi}_{1}\right|$, ${\theta}_{2}=2{tan}^{-1}\left|{\xi}_{2}\right|$ and φ_{12} = φ_{1} − φ_{2}.

Let ξ = ζ + *dζ*. Then, (60) leads to

the well-known expression for the element of length of ${S}_{R}^{2}$. Computation of the curvature tensor yields

and using Equation (59) we obtain

which is consistent with Equation (55).

*Example* 3.2. **Lobachevsky space**. Using the isomorphism established in Example 2.4 between the loop *QH*^{2} and Lobachevsky space, one can consider the upper part the two-sheeted unit hyperboloid, *H*^{2}, as a nonassociative space with the left translations defined by

In terms of the inner coordinates (θ, φ) on *H*^{2}, an arbitrary complex number ζ ∈ *D* can be written as ζ = *e*^{iφ} tanh(θ/2).

The computation of the elementary holonomy yields

We introduce the left invariant metric on *H*^{2} as

The Poincaré disk model is associated to the hyperbolic metric *d*(ζ, ξ), assigning to each pair of points ζ, ξ ∈ *D* the distance [79]

Thus, the hyperbolic distance between two points on *H*^{2} is related to the left invariant metric on *D* [see Equation (67)] as follows: cos *d* = 1 + ℓ^{2}/2, where we set ℓ = *g*(ζ, ξ). For *d* ≪ 1 we obtain *d* ≈ ℓ.

*Example* 3.3. **De Sitter spacetime**. De Sitter space, being a spacetime of constant positive curvature with topology *S*^{3} × ℝ, can be treated as a hyperboloid embedded in five-dimensional Minkowski space [80, 81],

with metric given by

The induced metric of the de Sitter space can be written as

The de Sitter loop, QH_{ℝ}, introduced in section 2, is isomorphic to the geodesic loop of de Sitter space. It has a natural geodiodular structure induced by the quaternionic algebra. The left invariant metric on QH_{ℝ} can be defined as follows [16]:

Letting ξ = ζ + *dζ*, this leads to the de Sitter metric (71):

For symmetric spaces the elementary holonomy is determined by the associator, ${h}_{(\zeta ,\eta )}\xi ={l}_{(\zeta ,{L}_{\zeta}^{-1}\eta )}\xi $ [30, 76]. Taking this into account, we obtain

Applying (55), we obtain the curvature tensor of de Sitter space:

**Comment 3.1**. De Sitter space can be obtained as a solution of the diodular Einstein's equations with a “Λ-term” [15].

*Example* 3.4. **Friedmann-Robertson-Walker spacetime**. The metric of the general Friedmann-Robertson-Walker (FRW) model may be written as [80, 81]

where Σ(*r, K*) = sin *r, r*, or sinh*r*, respectively, if *K* = 1, 0, or −1. It can be transformed to the form

which will be used below.

We introduce a binary operation in the neighborhood of the point $\tau \in {H}_{\text{\mathbb{R}}}^{{F}}$ as

where / denotes right division (see Example 2.7), *x* = ζ^{0} + *a*(ζ^{0})**ζ** and *y* = η^{0} + *a*(η^{0})**η**, with *e* = τ being the neutral element. The inverse operation is found to be

The associator is given by

where *z* = ξ^{0} + *a*(ξ^{0})**ξ**. The set of quaternions, ${\text{H}}_{\text{\mathbb{R}}}^{{F}}$, with the binary operation * forms a loop, that will be denoted as $\mathrm{\text{Q}}{\text{H}}_{\text{\mathbb{R}}}^{{F}}$. It is easy to show that the quaternionic algebra induces a natural geoodular structure on $\mathrm{\text{Q}}{\text{H}}_{\text{\mathbb{R}}}^{{F}}$.

Finally, we define the left invariant metric on $\mathrm{\text{Q}}{\text{H}}_{\text{\mathbb{R}}}^{{F}}$ as

where ζ = ζ^{0} + *a*(ζ^{0})**ξ** and ξ = ξ^{0} + *a*(ξ^{0})**ξ**. If we let **ξ** =**ζ** + *d***ζ** and ξ^{0} = ξ^{0} + *dτ*, then (80) yields the FRW metric:

## 4. Nonassociative Discrete Geometry

Since in case of discrete spaces we lack a smooth structure, we must use only algebraic structures to study them. Thus, it is quite natural to employ nonassociative geometry as the adequate algebraic approach (for details see [15–18]).

Consider a nonassociative finite space ${M}$ = 〈𝔐, *L, N*, (ω_{t})_{t ∈ ℝ}〉 equipped with (partial) operations: $L(x,y,z)={x}_{\stackrel{\xb7}{y}}z$, $N(x,y,z)={x}_{\stackrel{+}{y}}z$ and ω_{t}(*x, y*) = *t*_{x}*y*(*t* ∈ ℝ). This means that

• $\langle \U0001d510{,}_{\stackrel{\xb7}{a}},{({t}_{a})}_{t\in \text{\mathbb{R}}}\rangle $ is an odule with a neutral *a* ∈ *M*, for any *a* ∈ 𝔐,

• $\langle \U0001d510{,}_{\stackrel{+}{a}},{({t}_{a})}_{t\in \text{\mathbb{R}}}\rangle $ is a n-dimensional vector space (with zero element *a* ∈ 𝔐),

• The following geoodular identities are valid:

A curve in discrete space can be introduced as an ordered set of points, γ = (*a*_{1}, *a*_{2}, ⋯*a*_{s}) and, further, we can define the vector *y* (of ${{M}}_{{a}_{s}}^{+}$) parallel to a vector *x* (of ${{M}}_{{a}_{1}}^{+}$) along this curve γ as

where ${{M}}_{a}^{+}=\langle \U0001d510{,}_{\stackrel{+}{a}},a,{({t}_{a})}_{t\in \text{\mathbb{R}}}\rangle $. The curvature of a discrete space manifests as a non-trivial elementary holonomy, *h*^{a}(*x, y*) ≠ id,

Given a nonassociative space (in other words, a finite affinely connected space) one can enrich it by adding a metric structure. Specifically, we can define additionally the left invariant metric, ${g}_{b}({L}_{b}^{a}x,{L}_{b}^{a}y)={g}_{a}(x,y)$[17].

In what follows we consider some examples of nonassociative discrete spaces, including de Sitter space and FRW spacetime.

### 4.1. Lobachevsky Space

Here we consider a simple discrete model based on the loop *QH*^{2}. Let ℌ_{n} ⊂ *D* be a triangulated discrete topological space, with *n* being the number of vertices. Then ℌ_{n} becomes a nonassociative space if we define the partial nonassociative operation preserving the geometry of finite space, ℌ_{n},

The “metric" and elementary holonomy are defined by

Figure 2 depicts the discrete Lobachevsky space ℌ_{n} for *n* = 5, 112.

**Figure 2. (Left)** The discrete space ℌ_{n} in the Poincaré disk with 5,112 vertices. **(Right)** Three-dimensional representation of the same space.

### 4.2. The Nonassociative Discrete Geometry of *S*^{2}

Let *n* points be distributed randomly on the surface of the unit sphere, *S*^{2}. Consider the geodesic triangulation of the sphere, with the 2*D*-simplex being a geodesic triangle. We assign to each surface vertex, *p*, the polar coordinates (θ_{p}, φ_{p}). Without loss of generality, assume that the central vertex [i.e., the neutral element of the loop *QS*(2)] is located at the north pole. Using stereographic projection, we obtain

The non-associative (partial) operation (86) now takes the form

Writing ζ_{pq} as

one can rewrite (90) in spherical coordinates as follows:

where *l*(ζ_{p}, ζ_{q}) stands for the associator

The computation of the elementary holonomy yields:

We define the left-invariant metric on *S*^{2} as

**Comment 4.1**. The finite loop contains hidden information concerning the geometry of the sphere. The spherical symmetry is determined by the relation between the elementary holonomy and the associator, ${h}_{({\zeta}_{p},{\zeta}_{q})}={l}_{({\zeta}_{p},{L}_{{\zeta}_{p}}^{-1}{\zeta}_{q})}$. The smooth sphere can be regarded as a result of passage to the limit of *n* → ∞.

### 4.3. Discrete de Sitter Space

Consider a finite set $M={\mathbb{Z}}_{n}^{4}=\left\{p=({p}^{\mu})|{p}^{\mu}\in {\mathbb{Z}}_{n},\mu =0,\dots ,3,n\in \mathbb{N}\right\}$ where ℤ_{n} = {*p* = −*n*, …, *n*} are integers. A partial loop $Q{\text{H}}_{{\mathbb{Z}}_{n}^{4}}$is defined as a set of quaternions

with the indefinite norm $\left|\right|{\zeta}_{p}|{|}^{2}={\ell}^{2}{p}^{\mu}{p}_{\mu}$ and the binary operation ${\text{H}}_{{\mathbb{Z}}_{n}^{4}}\times {\text{H}}_{{\mathbb{Z}}_{n}^{4}}\mapsto {\text{H}}_{{\mathbb{Z}}_{n}^{4}}$ defined by

We define a partial geodiodular finite space ${M}$_{n} at the neutral zero as follows:

• $Q{\text{H}}_{{\mathbb{Z}}_{n}^{4}}$ is the odule with the multiplication operation induced by the quaternionic algebra over ℤ,

• ${{M}}^{+}={\text{H}}_{{\mathbb{Z}}_{n}^{4}}$ is the vector space induced by the quaternionic algebra over ℤ.

The left-invariant metric and elementary holonomy are given by

Taking *q* = *p* + **δ**, |**δ**| ≪ *n*, we obtain ζ_{q} = ζ_{p} + Δζ_{q}, where $\Delta {\zeta}_{q}=\ell ({\delta}^{0}+\mathrm{\text{i}}({\delta}^{1}i+{\delta}^{2}j+{\delta}^{3}k))$. Then, the metric (98) can be recast as

and

as ℓ → 0, *n* → ∞. A similar consideration of the elementary holonomy yields

where $\Delta ({\zeta}_{p},{\zeta}_{q},{\zeta}_{m})={\zeta}_{m}{\zeta}_{q}^{+}{\zeta}_{p}-{\zeta}_{q}{\zeta}_{p}^{+}{\zeta}_{m}$. Approaching the limit ℓ → 0 and *n* → ∞, one recovers the curvature tensor of de Sitter space:

Thus, smooth de Sitter space could be regarded as the passage to the limit of ℓ → 0, *n* → ∞ and ℓ*n* = const. We see that, to some extent, the information about the geometry of de Sitter spacetime is hidden in the finite loop $Q{\text{H}}_{{\mathbb{Z}}_{n}^{4}}$.

### 4.4. Discrete Friedmann-Robertson-Walker Spacetime

Consider a finite set $M={\mathbb{Z}}_{n}^{3}=\left\{p=({p}^{i})|{p}^{i}\in {\mathbb{Z}}_{n},i=1,2,3,n\in \mathbb{N}\right\}$, where ℤ_{n} = {*p* = −*n*, …, *n*} are integers. We denote by $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ a partial loop defined as a set of quaternions,

with the norm $\parallel {\zeta}_{\mathbf{\text{p}}}{\parallel}^{2}={\ell}^{2}{p}^{i}{p}_{i}=-{\ell}^{2}\left({\left({p}^{1}\right)}^{2}+{\left({p}^{2}\right)}^{2}+{\left({p}^{3}\right)}^{2}\right)$ and the binary operation ${\text{H}}_{{\mathbb{Z}}_{n}^{3}}\times {\text{H}}_{{\mathbb{Z}}_{n}^{3}}\mapsto {\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ given by

A partial geodiodular finite space ${M}$_{n} at the neutral element *e* = 0 is defined as

• $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ is the odule with multiplication induced by the quaternionic algebra over ℤ,

• ${{M}}_{n}^{+}={\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ is the vector space induced by the quaternionic algebra over ℤ.

Let ${\text{H}}_{{\mathbb{Z}}_{n}^{3}\times {\text{\mathbb{T}}}_{n}}^{{F}}$ be a set of quaternions,

where 𝕋_{n} = {τ_{n}: *n* ∈ ℕ} and τ_{n} is a discrete parameter. We define the partial loop $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}\times {\text{\mathbb{T}}}_{n}}^{{F}}\simeq \mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}\times {\text{\mathbb{T}}}_{n}$ as the loop with the operation * given by

For the associator we have

We introduce the nonassociative discrete FRW spacetime as ${{M}}^{{F}}={\bigcup}_{n\in \mathbb{N}}{{M}}_{n}^{{F}}$, where a partial geodiodular finite space ${{M}}_{n}^{{F}}$ consists of

• $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}\times {\text{\mathbb{T}}}_{n}}^{{F}}$, the odule with multiplication induced by the quaternionic algebra,

• ${{M}}_{n}^{+}={\text{H}}_{{\mathbb{Z}}_{n}^{3}}\times {\text{\mathbb{T}}}_{n}$, the vector space with structure induced by the quaternionic algebra over ℤ,

• the partial binary operation for the given ${x}_{p}^{n}\in {{M}}_{n}^{{F}}$ and ${x}_{q}^{{n}^{\prime}}\in {{M}}_{{n}^{\prime}}^{{F}}$ is defined as

Thus, τ_{n} can be interpreted as a “discrete time” in our model.

The computation of the elementary holonomy yields

One can introduce the discrete left invariant metric on ${M}$^{${F}$} in the same way as in Equation (98),

Finally, we obtain the metric for FRW spacetime in the continuous limit,

## 5. Emergence of Spacetime

We are now ready to propose a discrete model of spacetime based on nonassociative geometry. In our approach we consider a triangulated topological space, where the vertices are “points" of discrete spacetime, and form the elements of the loop with a partial operation. The edges are identified as the “geodesics" connecting vertices. In contrast to the Regge model (for details and a review, see [82–85]), in nonassociative geometry the curvature does not reside at the vertices, but is defined by the elementary holonomy associated with 2-faces of simplices.

Our approach is based on the following assumptions:

• Vertices of the simplices are considered as partially ordered “atoms" of discrete spacetime.

• The distance between two neighboring atoms can not be less than the fundamental length, ℓ_{p}.

• Interaction between atoms of spacetime, being nonlocal and nonassociative, defines the spacetime geometry.

• Spacetime geometry is encoded in the nonassociative structure of spacetime.

• The foliation of spacetime consists of space-like slices joined by time-like edges.

• Time is quantized and the evolution of spacetime geometry is governed by a random/stochastic process.

• The spacetime dimension may be dynamical.

Let 𝔐_{n} be a triangulated topological space, with *n* being the number of vertices. We identify the elements *x*_{p} ∈ 𝔐_{n} as the vertices of simplices and introduce a partial nonassociative operation, *x*_{p}**x*_{q}. Since the introduced operation is partial, for some elements *x*_{p}, *x*_{q} ∈ 𝔐_{n} the result of multiplication may not belong to 𝔐_{n}: ${x}_{p}*{x}_{q}\mapsto {x}_{p}^{\prime}\notin {\U0001d510}_{n}$. Adding the point ${x}_{p}^{\prime}$ to the set 𝔐_{n} as a new vertex, we obtain 𝔐_{n+1}. This leads to the partial ordering

and the nonassociative discrete spacetime, 𝔐 = ⋃_{n∈ℕ} 𝔐_{n}, becomes a partially ordered finite set or *causal set* ([20, 86], see section 6 for details). The evolution of spacetime geometry is interpreted in terms of the sequence: $\cdots \to {\U0001d510}_{n}\to {\U0001d510}_{n+1}\to \cdots \to {\U0001d510}_{{n}^{\prime}}\dots \phantom{\rule{0.3em}{0ex}}$.

**Comment 5.1**. Albeit the fundamental length, ℓ_{p}, is not necessarily equal to the Planck length, we will take it to be so from this point forward. The known problem of “many-fingered time” [4–8] is overcome in our model as follows. In contrast to CDT, in our approach edges are allowed to be longer than the Planck length. This implies that the local “time” between two slices depends on the ‘coordinates’ (vertices), and each realization of the random process yields different results. Thus, history becomes foliation dependent (i.e., the many-fingered time problem). However, actual physical quantities are defined by probabilities, which are foliation independent. Since in our approach the evolution of the spacetime geometry is governed by a stochastic dynamical process, we avoid the foliation dependence problem.

To illustrate our approach, below we consider a cosmological toy model of (2+1)-D discrete spacetime and a discrete FRW cosmological model.

### 5.1. Toy Model of (2+1)-D Discrete Spacetime

In this section we consider a simple model based on the non-associative compact space 𝔖(τ_{n}) ≃ *Q𝔖*(τ_{n}), obtained by triangulation of the partial loop, *Q𝔖*(τ_{n}), defined as follows. An arbitrary element, ${\zeta}_{p}^{n}\in Q\U0001d516({\tau}_{n})$, takes the form ${\zeta}_{p}^{n}={R}_{n}{\zeta}_{p}$, where *R*_{n} = *R*(τ_{n}) is the “radius” of 𝔖(τ_{n}), τ_{n} is the discrete intrinsic time, a quantum of time is Δτ_{n} = τ_{n+1} − τ_{n} = 1 in Planck units, and ζ_{p} ∈ ℂ is the *p*-th complex number.

The partial nonassociative operation preserving the geometry of spacetime is given by

We assume that evolution of the discrete spacetime geometry is a stochastic process. The radius *R*_{n} is defined by a random walk with step ℓ_{p} and reflecting barrier *R*_{0} at *n* = 0 [87–89]. After *n* steps and for *n* ≫ 1, the expectation value of the scale factor will grow as

The computation yields

In Figure 3 we present the results of our numerical simulations. As one can see (left panel) the analytical formula (115) is in good agreement with our numerical results. We find that, for large *n*, the average area of each 2D-simplex converges to the area of the equilateral triangle, $A=\sqrt{3}/4\approx 0.433$. (see Figure 3, right panel.)

**Figure 3. (Left)** The average radius, 〈*R*〉, as a function of iteration number, *n*, for a 10^{3} × 10^{4} trial (runs × iterations; error bars are omitted for clarity). The red curve corresponds to Equation (116). **(Right)** The average simplex area as a function of iteration number, *n*, for a 20 × 15. It slowly approaches the area of the equilateral triangle ≈0.43.

Assume that 𝔖(τ_{n}) has *m* points. The map *Q𝔖*(τ_{n}) ↦ *Q𝔖*(τ_{n+1}) defined by Equation (113) yields *m*^{2} points belonging to *Q𝔖*(τ_{n+1}). To obtain 𝔖(τ_{n+1}) one should triangulate the set *Q𝔖*(τ_{n+1}). During the triangulation only the points with distance ℓ ≥ ℓ_{p} between closest neighbors will survive. After triangulation, we are left with the two space-like slices 𝔖(τ_{n}) and 𝔖(τ_{n+1}). Next, we triangulate the set 𝔖(τ_{n})⋃𝔖(τ_{n+1}), joining the vertices of the spatial slices by time-like edges. Repeating this procedure, we obtain a partially ordered discrete spacetime, 𝔖 = ⋃_{n}𝔖(τ_{n}): 𝔖(τ_{0})≺ ⋯ ≺𝔖(τ_{n})≺𝔖(τ_{n+1})≺ … .

In Figures 4–7 we present the results of numerical simulations showing a particular realization of the random process. As an initial state at the moment τ_{0} = 0, we take a discrete space 𝔖(0) defined by the tetrahedron formed by 2-dimensional equilateral simplices of the same size, with distance ℓ_{p} between neighboring vertices and ${R}_{0}=\sqrt{3/8}{\ell}_{p}$. As the number of iterations increases, the triangulation of 𝔖(τ_{n}) becomes finer and its elements approach equilateral 2*D*-simplices (see Figure 4). As one can see, with growing *n* the simplicial complex is approximated by *S*^{2}. The results of our numerical simulations illustrating the evolution of spacetime geometry are presented in Figure 5. In Figure 6 the triangulation of the spacetime, 𝔖 = ⋃_{n} 𝔖(τ_{n}) for *n* = 5 is presented. The spatial foliation of the spacetime geometry is depicted in Figure 7. The correspondence between simplicial complexes in the complex plane, ℂ, and 𝔖(τ_{n}) is shown on the right side of the figure. The surfaces τ_{n} = const and τ_{n+1} = const represent adjacent two-dimensional spaces. The timelike edges are omitted for clarity.

**Figure 7**. Spacetime evolution after *n* = 15 steps. Timelike links between vertices of neighboring slices are omitted for clarity.

### 5.2. Discrete Spacetime: FRW Cosmological Model

Here we consider a cosmological discrete spacetime model based on the discrete nonassociative FRW geometry described in section 4.4. We start with the partial loop $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ based on the set of quaternions

Let ${\zeta}_{p}^{n}=R\left({\tau}_{n}\right){\zeta}_{p}$ be an element of the loop $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$, where τ_{n} is a discrete intrinsic time, *R*(τ_{n}) = *a*(τ_{n})ℓ_{p} with *a*(τ_{n}) being the scale factor. We define the partial binary operation in $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$ as follows:

where “/” denotes the right division of the quaternions.

Let 𝔐(τ_{n}) be a triangulated topological space obtained by triangulation of $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n}^{3}}$. Assume that 𝔐(τ_{n}) has *m* points. Then the map 𝔐(τ_{n})↦𝔐(τ_{n+1}) defined by Equation (118) yields *m*^{2} points belonging to $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n+1}^{3}}$. However, only the points with distance ℓ ≥ ℓ_{p} between closest neighbors survive during the triangulation of $\mathrm{\text{Q}}{\text{H}}_{{\mathbb{Z}}_{n+1}^{3}}$. After triangulation, we obtain two space-like slices, 𝔐(τ_{n}) and 𝔐(τ_{n+1}). The next step is to triangulate the set 𝔐(τ_{n})⋃𝔐(τ_{n+1}), joining the vertices of spatial slices by time-like edges. Repeating this procedure, we obtain a partially ordered discrete spacetime: 𝔐(τ_{0}) ≺ ⋯ ≺ 𝔐(τ_{n}) ≺ 𝔐(τ_{n+1}) ≺ … . Since some information is lost during triangulation in the form of “missing” vertices (those that would have been closer than ℓ_{p}), the evolution is irreversible. Thus, in our approach the arrow of time, being related to the partial ordering, appears naturally [18].

We assume that evolution of the discrete spacetime geometry is a stochastic process. The scale factor, *R*_{n}, is defined by a random walk with a step ℓ_{p} and reflecting barrier *R*_{0} at *n* = 0. After *n* steps, the average scale factor will grow as $\langle {R}_{n}\rangle ={R}_{0}+\sqrt{2n/\pi}{\ell}_{p}$ (see Figure 8, right panel).

**Figure 8. (Left)** Initial 3D simplex 𝔐(τ_{0}). **(Right)** The average radius as a function of iteration number, *n*, for a 200 × 25 trial (runs × iterations; error bars are omitted for clarity).

After ~10^{8} steps in Planck units of time, which corresponds to pre-inflationary time, ${t}_{i}~1{0}^{-35}\mathrm{\text{sec}}$ (time of the beginning of inflation), the size of the universe in our model would be ${R}_{pr}~1{0}^{-29}\mathrm{\text{cm}}$. The size of the pre-inflationary universe in the standard FRW-model can be estimated using the relationship between the redshift and the cosmological scale factor, *z*_{i} = *a*(*t*)/*a*(*t*_{i})−1, where *t* denotes the actual age of the universe. To estimate the redshift, we use the ΛCDM model. For the dark energy dominated era the redshift is given by Ryden [90] and Frieman et al. [91]

where ${\Omega}_{\Lambda}=\Lambda /(3{H}_{0}^{2})$. This expression is quite accurate for $t>{t}_{m}=10\mathrm{\text{Myr}}(1{0}^{7}\mathrm{\text{y}})=3.15\xb71{0}^{14}\mathrm{\text{sec}}$. For a radiation-dominated era, we have $a(t)\propto \sqrt{t}$, and for a matter-dominated universe, *a*(*t*) ∝ *t*^{2/3}. During inflation the universe is blown up to enormous size, and the changes in the cosmological scale factor are of order $a({t}_{i})/a({t}_{f})~{e}^{60}\approx 1{0}^{26}$, where *t*_{f} denotes the time of the end of inflation. Combining all results, we obtain

Inserting ${t}_{f}\approx 1{0}^{-33}\mathrm{\text{s}}$, ${t}_{r}\approx 1{0}^{14}\mathrm{\text{s}}$, *t* ≈ 4.35·10^{17}s (age of the universe), Ω_{Λ} = 0.7 and ${H}_{0}=2.33\xb71{0}^{-18}{\mathrm{\text{s}}}^{-1}$ into this expression, we obtain ${z}_{i}\simeq 1{0}^{52}$. The pre-inflationary size of the universe can be estimated as *R*_{i} = *R*(*t*)/*z*_{i}.

If we choose *R*(*t*) as the observable radius of the universe, *R*(*t*) ~ 10^{28}cm, then the radius of the pre-inflationary universe is ${R}_{i}~1{0}^{-24}\mathrm{\text{cm}}$. Comparing this result with the prediction of our model, we find ${R}_{i}/{R}_{pr}=1{0}^{5}$. Thus, starting with only one initial mini-universe, we cannot reach an appropriate size of the universe at the beginning of the inflationary epoch. To overcome this difficulty, one can assume that spacetime would consist of many small fluctuating simplicial regions.

These fluctuations may form mini-clusters, as illustrated in Figure 9, created in different regions and may grow as “mini-universes” (Figure 10). We consider all these mini-universes as a multiverse [92, 93]. After ~10^{7}−10^{8} steps from the very beginning, the mini-universes grow enough and may merge in different regions of space to form pre-inflationary universes (bubbles) with appropriate average size ~ 10^{28} cm. This multiverse state can be metastable and, as a result of the first-order phase transition, the nucleation of bubbles may occur and lead to inflation. Indeed, if the duration of the phase transition is short enough, bubble nucleation of the new phase may have an exponential character [94–96]. Thus, in our model the post-inflationary Universe is the result of first-order inflation, occurring during a strong first-order phase transition [97–99] (Figure 11).

**Figure 10**. “Emergent spacetime”: formation of mini-universe (bubble) from 200 randomly distributed initial 3D-simplices.

**Comment 5.2**. An alternative scenario of arising spacetime structure may be formation of 3D-clusters from 2D-simplicial complexes as illustrated in Figure 9. This implies that the 2D-simplices are more fundamental than 3D-simplices, and 3D volume effectively is created from random foam-like structures. In this scenario the spacetime dimension is dynamical and may vary depending on fusion of simplicial complexes.

## 6. Discussion and Conclusions

The preceding examples show how the framework of nonassociative geometry may describe the discrete structure of spacetime. We deal with points (“atoms”) of discrete nonassociative spacetime. In our approach, at the Planckian scale the standard concept of spacetime is replaced by the nonassociative discrete structure. At large spacetime scales this “looks like” a smooth manifold, and the information about the geometry of the continuous spacetime is hidden in the structure of the discrete nonassociative spacetime.

In our approach it is assumed that the “universe” is created from a few initial points. This leads to a new understanding of the “big-bang”: the initial singularity does not exist. Spacetime evolution is described by a stochastic process preserving spacetime geometry. As a result, we obtain a discrete spacetime with many of the attributes of continuum spacetime: spacelike slices, “many-fingered time,” and causal structure. Note that, while the history of spacetime is foliation dependent, the observables are foliation independent since they are defined by probabilities.

If we admit that the fluctuations can occur in different regions, then foam-like structures can be created and form a spatial foam (Figure 9). This implies that the dimension of spacetime is dynamical and 3D-clusters are created from 2D-simplicial complexes like a volume can be created from foam-like structures. This scenario is very speculative, but deserves further study.

Existence of the fundamental length as a minimal admissible length scale in the Universe implies that the emergence of space in the expanding cosmological models is a continuous process. This point of view is supported by the models with eternal inflation. In these models the inflationary phase of the Universe's expansion lasts forever in at least some regions of the Universe [93, 97].

We would like to compare our approach with relevant models of discrete spacetime proposed in CS [19–21] and CDT [14, 22–28]. Both theories attempt, as their names suggest, to preserve the causality present in continuous spacetime while also constructing a discrete underlying structure. We briefly recount the main features of these approaches to compare and contrast them with our model.

CS proposes that the microstructure of spacetime is a partially ordered set (causal set or causet). A *causal set* is a partially ordered set of elements ${C}$ (sometimes also abbreviated a *poset*) together with the relation ≺, called “precedes,” such that 1) if *x* ≺ *y* and *y* ≺ *z*, then *x* ≺ *z*, ∀*x, y, z* ∈ ${C}$; 2) if *x* ≺ *y* and *y* ≺ *x*, then *x* = *y*, ∀*x, y, z* ∈ ${C}$; 3) for any pair of elements *x* and *z* in ${C}$, the set of other elements {*y*|*x* ≺ *y* ≺ *z*} between *x* and *z* is finite. Thus, the partial order encodes information about the causal structure of spacetime.

These axioms establish the transitivity, non-circularity, and finiteness of causal sets. Thus, the elements of ${C}$ lie as nodes along branches of links between them, with each element representing an abstraction of a fundamental unit of spacetime. The fundamental hypothesis of CS theory is that spacetime is described by a causet made of elements of Planckian size in the quantum regime. The order of the set gives rise to the causality present in macroscopic spacetime, whereas the number of elements determines the spacetime's volume. Together, these two features combine to produce the macroscopic geometry of spacetime.

Like CS, CDT is also concerned with preserving causality at a fundamental scale but is more explicitly geometrical. Formally, CDT seeks a discrete version of Feynman path integral [22, 28],

where the integration is performed over all Lorentzian geometries with topology Σ × [0, 1], with Σ being a compact, connected three-dimensional manifold. The fundamental building blocks of CDT are equilateral 4-simplices of side length *a*, that is, four-dimensional generalizations of equilateral tetrahedra. These simplices are arranged such that their spacelike edges form spacelike surfaces separated by their timelike edges. More explicitly, the partition function, *Z*_{a} is given by

where *S*[*T*] is a suitable gravitational action, and the sum is over the chosen set of abstract triangulations of ${M}$.

So far, CDT has relied on Monte Carlo techniques to build fluctuating spacetimes out of simplices. It should also be stressed that the minimum length cutoff at *a* is *not* identified with the Planck length, but is used as part of a limiting procedure to obtain a continuum spacetime. We refer the interested reader to the cited bibliography and references therein for promising results and developments in both CS and CDT, which would be out of scope for our present discussion.

The description of discrete spacetime by nonassociative geometry is both similar and distinct from those by CS and CDT. Like these two approaches, nonassociative geometry incorporates causality as a fundamental axiom of the model. A given simplicial spacetime is in a precedence relation to the one that follows it, 𝔐(τ_{n}) ≺ 𝔐(τ_{n+1}). A temporal direction is established and can be represented by timelike edges between nodes belonging to times τ_{n} and τ_{n+1}. Since vertices in 𝔐(τ_{n}) determine those in 𝔐(τ_{n+1}) but, due to the minimum length restriction, not vice versa, a definite arrow of time is established. However, unlike the stochastic sequential growth used in CS, where different branches of spacetime grow one element at a time in each discrete time interval, nonassociative geometry elements are created simultaneously in a non-local manner. And unlike CDT, timelike links in nonassociative geometry can have variable length. This yields the partial ordering of spacetimes,

and the nonassociative discrete spacetime, 𝔐 = ⋃_{n∈ℕ} 𝔐(τ_{n}), becomes a finite poset as described in section 5. The evolution of spacetime geometry is interpreted in terms of a sequence: $\cdots \to \U0001d510({\tau}_{n})\to \U0001d510({\tau}_{n+1})\to \cdots \to \U0001d510({\tau}_{{n}^{\prime}})\dots \phantom{\rule{0.3em}{0ex}}$. Thus, in our approach the arrow of time appears naturally. The other difference with CDT is that, in contrast to CDT, in our model the curvature, being defined by the elementary holonomy associated with the faces of 2D-simplices, is nonlocal.

Simplicial cells are fundamental components of nonassociative discrete spacetimes, as in CDT. However, our numerical implementation of nonassociative geometry explicitly adopts the Planck length, ℓ_{p}, as a minimum cutoff length, unlike CDT. Also unlike CDT, nonassociative geometry does not restrict simplices to be equilateral, although the realization we present does converge to this scenario in the limit of high iteration numbers (i.e., large values of time). Another difference with CDT is that spacetime evolution in nonassociative geometry, while having a certain stochastic character, does not rely on Monte Carlo techniques to generate fluctuations in spacetimes, since the “random” component in the evolution resides completely in the random-walk behavior of *R*.

In conclusion, we have shown how one can describe continuum and discrete spacetimes in the framework of the nonassociative geometry. The first feature of our model of discrete spacetime is the obvious absence of an initial singularity. The other important feature is a natural emergence of causal structure and an arrow of time. The irreversibility of evolution is caused by a random process and information loss due to existence of the minimal fundamental Planck length. Certainly, the loopuscular structures used in section 5 are given in an *ad hoc* manner. The open question is then how to describe more general models, using the sum over all possible geometries as in the CDT approach. This is the subject of the future study.

## Author Contributions

AN and HM designed the research and wrote the text. HM wrote the code and prepared figures. All authors reviewed the manuscript.

## Funding

This research was supported by CONACyT Network Project No. 294625 Agujeros Negros y Ondas Gravitatorias.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors acknowledge the support from the CONACyT.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2019.00032/full#supplementary-material

## References

1. Ashtekar A, Stachel J. *Conceptual Problems of Quantum Gravity. Einstein Studies 2*. Boston, MA: Birkhüser (1991).

2. Gross D, Henneaux M, Sevrin A. The quantum structure of space and time. In: *Proceedings of the 23rd Solvay Conference on Physics*. Brussels (2007).

3. Murugan J, Weltman A, Ellis GFR editors. *Foundations of Space and Time: Reflections on Quantum Gravity*. Cambridge: Cambridge University Press (2012).

4. Isham CJ. Canonical quantum gravity and the question of time. In:. Ehlers J, Friedrich H, editors. *Canonical Gravity: From Classical to Quantum: Proceedings of the 117th WE Heraeus Seminar Held at Bad Honnef, Germany, 13–17 September 1993*. Berlin: Springer Berlin Heidelberg (1994). p. 150–69.

5. Kuchar K. *Time and Interpretations of Quantum Gravity.* London: Imperial College of Science and Technology (1992).

6. Rovelli C. Time in quantum gravity: an hypothesis. *Phys Rev D.* (1991) **43**:442–56. doi: 10.1103/PhysRevD.43.442

7. Kiefer C. *Quantum Gravity.* 2nd ed. International Series of Monographs on Physics. New York, NY: Oxford University Press (2007).

9. Anderson E. *The Problem of Time: Quantum Mechanics Versus General Relativity*. 1st ed. Cham: Springer International Publishing AG (2017).

10. Butterfield J, Isham C. Spacetime and the philosophical challenge of quantum gravity. In: Callender C, Huggett N, editors. *Physics Meets Philosophy at the Planck Scale.* Cambridge: Cambridge University Press (2004). p. 33–89.

11. Oriti D (editor). *Approaches to Quantum Gravity: Toward a New Understanding of Space, Time and Matter.* New York, NY: Cambridge University Press (2009).

12. Majid S, Polkinghorne J, Penrose R, Taylor A, Connes A, Heller M. *On Space and Time*. New York, NY: Cambridge University Press (2012).

13. Ambjørn J. Discrete quantum gravity. In: Oriti D, editor *Approaches to Quantum Gravity: Toward a New Understanding of Space, Time and Matter.* New York, NY: Cambridge University Press (2009).

14. Loll R. Discrete approaches to quantum gravity in four dimensions. *Living Rev Relat.* (1998) **1**:13. doi: 10.12942/lrr-1998-13

15. Nesterov AI, Sabinin LV. Non-associative geometry and discrete structure of spacetime. *Comment Math Univ Carolin.* (2000) **41**:347–58.

16. Nesterov AI, Sabinin LV. Nonassociative geometry: towards discrete structure of spacetime. *Phys Rev D.* (2000) **62**:081501. doi: 10.1103/PhysRevD.62.081501

17. Sabinin L. Nonassociative geometry and discrete space-time. *Int J Theor Phys.* (2001) **40**:351–8. doi: 10.1023/A:1003712200867

18. Nesterov AI, Sabinin LV. Nonassociative geometry: Friedmann-Robertson-Walker spacetime. *IJGMMP.* (2006) **3**:1481–91. doi: 10.1142/S0219887806001806

19. Dowker F. Causal sets and discrete spacetime. *AIP Confer Proc.* (2006) **861**:79–88. doi: 10.1063/1.2399569

20. Bombelli L, Lee J, Meyer D, Sorkin RD. Space-time as a causal set. *Phys Rev Lett.* (1987) **59**:521–4. doi: 10.1103/PhysRevLett.59.521

21. Dowker F, Zalel S. Evolution of universes in causal set cosmology. *Comptes Rendus Phys.* (2017) **18**:246–53. doi: 10.1016/j.crhy.2017.03.002

22. Ambjørn J, Carfora M, Marzuoli A. *The Geometry of Dynamical Triangulations.* Berlin: Springer Berlin Heidelberg (1997).

23. Ambjørn J, Loll R. Non-perturbative lorentzian quantum gravity, causality and topology change. *Nuclear Phys. B.* (1998) **536**:407–34. doi: 10.1016/S0550-3213(98)00692-0

24. Loll R. The emergence of spacetime or quantum gravity on your desktop. *Class Quant Grav.* (2008) **25**:114006. doi: 10.1088/0264-9381/25/11/114006

25. Ambjørn J, Jurkiewicz J, Loll R. Reconstructing the universe. *Phys Rev D.* (2005) **72**:064014. doi: 10.1103/PhysRevD.72.064014

26. Ambjørn J, Jordan S, Jurkiewicz J, Loll R. Quantum spacetime, from a practitioner's point of view. *AIP Confer Proc.* (2013) **1514**:60–6. doi: 10.1063/1.4791726

27. Glaser L, Loll R. CDT and cosmology. *Comptes Rendus Phys.* (2017) **18**:265–74. doi: 10.1016/j.crhy.2017.04.002

28. Ambjørn J, Görlich A, Jurkiewicz J, Loll R. Nonperturbative quantum gravity. *Phys Rep.* (2012) **519**:127–210. doi: 10.1016/j.physrep.2012.03.007

29. Sabinin L. Odules as a new approach to a geometry with a connection. *Soviet Math. Dokl.* (1977) **18**:515.

31. Gogberashvili M. Octonionic geometry. *Adv Appl Cliff Algeb.* (2005) **15**:55–66. doi: 10.1007/s00006-005-0003-2

32. Gogberashvili M, Sakhelashvili O. Geometrical applications of split octonions. *Adv Math Phys.* (2015) **2015**:14. doi: 10.1155/2015/196708

33. Günaydin M, Gürsey F. Quark structure and octonions. *J Math Phys.* (1973) **14**:1651–67. doi: 10.1063/1.1666240

34. Günaydin M, Nicolai H. Seven-dimensional octonionic yang-mills instanton and its extension to an heterotic string soliton. *Phys Lett B.* (1995) **351**:169–72. doi: 10.1016/0370-2693(95)00375-U

35. Strobl T. Algebroid yang-mills theories. *Phys Rev Lett.* (2004) **93**:211601.1–4. doi: 10.1103/PhysRevLett.93.211601

36. S Ramgoolam. Towards Gauge theory for a class of commutative and non-associative fuzzy spaces. *JHEP.* (2004) **0403**:034. doi: 10.1088/1126-6708/2004/03/034

37. de Medeiros P, Ramgoolam S. Non-associative gauge theory and higher spin interactions. *JHEP.* (2005) **03**:072. doi: 10.1088/1126-6708/2005/03/072

38. Mylonas D, Schupp P, Szabo RJ. Non-geometric fluxes, quasi-hopf twist deformations, and nonassociative quantum mechanics. *J Math Phys.* (2014) **55**:122301. doi: 10.1063/1.4902378

39. Dzhunushaliev V. Nonassociativity, supersymmetry, and hidden variables. *J Math Phys.* (2008) **49**:042108. doi: 10.1063/1.2907868

40. Dzhunushaliev V. Hidden nonassociative structure in supersymmetric quantum mechanics. *Ann Phys.* (2010) **522**:382–8. doi: 10.1002/andp.201000014

41. Dzhunushaliev V. Observables and unobservables in a non-associative quantum theory. *J Gener Lie Theor Appl.* (2008) **2**:269–72. doi: 10.4172/1736-4337.1000148

42. Dzhunushaliev V. Nonassociative generalization of supersymmetry. *Adv Appl Cliff Algeb.* (2017) **27**:303–10. doi: 10.1007/s00006-015-0580-7

43. Jackiw R. Three-cocycle in mathematics and physics. *Phys Rev Lett.* (1985) **54**:159–62. doi: 10.1103/PhysRevLett.54.159

44. Boulware DG, Deser S, Zumino B. Absence of 3-cocycles in the Dirac Monopole Problem. *Phys Lett B.* (1985) **153**:307–10. doi: 10.1016/0370-2693(85)90554-4

45. Grossman B. A 3-cocycle in quantum mechanics. *Phys Lett B.* (1985) **152**:93–7. doi: 10.1016/0370-2693(85)91146-3

46. Grossman B. Three-cocycle in quantum mechanics. II. *Phys Rev D.* (1985) **33**:2922–9. doi: 10.1103/PhysRevD.33.2922

47. Wu YS, Zee A. Cocycles and magnetic monopole. *Phys Lett B.* (1985) **152**:98–102. doi: 10.1016/0370-2693(85)91147-5

48. Carey AL. The origin of three-cocycles in quantum field theory. *Phys Lett B.* (1987) **194**:267–70. doi: 10.1016/0370-2693(87)90540-5

49. Nesterov AI, Aceves de la Cruz F. Magnetic monopoles with generalized quantization condition. *Phys Lett A.* (2002) **302**:253–60. doi: 10.1016/S0375-9601(02)01172-6

50. Nesterov AI. Three-cocycles, nonassociative gauge transformations and Dirac's monopole. *Phys Lett A.* (2004) **328**:110–5. doi: 10.1016/j.physleta.2004.06.024

51. Bruno P. Non-quantized Dirac monopoles and strings in the Berry phase of anisotropic spin systems. *Phys Rev Lett.* (2004) **93**:247202.1–4. doi: 10.1103/PhysRevLett.93.247202

52. Nesterov AI. Nonassociativity, Dirac monopoles and Aharonov-Bohm effect. *IJGMMP.* (2007) **4**:717–26. doi: 10.1142/S0219887807002259

53. Jo SG. Commutators in an anomalous non-abelian gauge theory. *Phys Lett B.* (1985) **163**:353–9. doi: 10.1016/0370-2693(85)90295-3

54. Niemi AJ, Semenoff GW. Quantum holonomy and the chiral gauge anomaly. *Phys Rev Lett.* (1985) **55**:927–30. doi: 10.1103/PhysRevLett.55.927

55. Semenoff GW. Nonassociative electric fields in chiral gauge theory: an explicit construction. *Phys Rev Lett.* (1988) **60**:680–3. doi: 10.1103/PhysRevLett.60.1590.2

56. Treiman S, Jackiw R, Zumino B, Witten E. *Currenl Algebra and Anomalies.* Singapore: World Scientific (1985).

58. Nesterov AI, Stepanenko VA. On methods of nonassociative algebra in geometry and physics. *L. V. Kirensky Inst Phys.* (1986) 48.

59. Nesterov AI. *Methods of Nonassociative Algebra in Physics*. Dr. Sci. Dissertation, Krasnoyarsk State University, Krasnoyarsk (1989).

60. Ungar AA. The abstract complex Lorentz transformation group with real metric. I. Special relativity formalism to deal with the holomorphic automorphism group of the unit ball in any complex Hilbert space. *J Math Phys.* (1991) **35**:1408–26. doi: 10.1063/1.530597

61. Ungar AA. Thomas precession and its associated grouplike structure. *Am J Phys.* (1991) **59**:824–34. doi: 10.1119/1.16730

62. Nesterov AI. Quasigroups, asymptotic symmetries and conservation laws in general relativity. *Phys Rev D.* (1997) **56**:R7498–502. doi: 10.1103/PhysRevD.56.R7498

63. Nesterov AI. On angular momentum of gravitational radiation. *Phys Lett A.* (1998) **250**:55–61. doi: 10.1016/0375-9601(90)90228-G

64. Nesterov AI. Principal *Q*-bundles. In: Costa R, Cuzzo Jr H, Grishkov A, Peresi LA, editors, *Non Associative Algebra and Its Applications.* New York, NY: Marcel Dekker (2000). p. 247–57.

65. Nesterov AI. Principal loop bundles: toward nonassociative gauge theories. *Int J Theor Phys.* (2001) **40**:339–50. doi: 10.1023/A:1003760016797

66. Nesterov AI. Smooth loops and fiber bundles: theory of principal Q-bundles. *IJGMMP.* (2009) **6**:77–97. doi: 10.1142/S0219887809003412

68. Sabinin L. On the equivalence of categories of loops and homogeneous spaces. *Soviet Math Dokl.* (1972) **13**:970.

69. Sabinin LV. Smooth quasigroups and loops: forty-five years of incredible growth. *Comment Math Univ Carolin.* (2000) **41**:377–400.

70. Kikkawa M. On local loops in affine manifolds. *Sci Hiroshima Univ Ser A-I Math.* (1964) **28**:199–207.

73. Batalin IA. Quasigroup constuction and first class constraints. *J Math Phys.* (1981) **22**:1837–56. doi: 10.1063/1.525155

74. Hudson SN. Lie loops with invariant uniformities. *Trans Amer Math Soc.* (1965) **115**:417–32. doi: 10.2307/1994278

75. Hudson SN. Lie loops with invariant uniformities II. *Trans Amer Math Soc.* (1965) **118**:526–3. doi: 10.2307/1993976

77. Sabinin LV. Differential geometry and quasigroups. *Proc Inst Math Siberian Branch Ac Sci USSR.* (1989) **14**:208.

79. Bianconi G, Rahmede C. Emergent hyperbolic network geometry. *Sci Rep.* (2017) **7**:41974. doi: 10.1038/srep41974

80. Hawking S, Ellis G. *The Large Scale Structure of Space-Time.* Cambridge: Cambridge University Press (1975).

81. Kramer D, Stephani H, MacCallum M, Herlt E. *Exact Solutions of Einstein's Equations.* Berlin: VEB Deutscher Verlag (1980).

82. Regge T. General relativity without coordinates. *Il Nuovo Cimento (1955-1965).* (1961) **19**:558–71. doi: 10.1007/BF02733251

83. Sorkin R. Time-evolution problem in Regge calculus. *Phys Rev D.* (1975) **12**:385–96. doi: 10.1103/PhysRevD.12.385

84. Williams RM, Tuckey PA. Regge calculus: a brief review and bibliography. *Class Quant Grav.* (1992) **9**:1409–22. doi: 10.1088/0264-9381/9/5/021

85. Williams R. Recent progress in Regge calculus. *Nucl Phys B Proc Suppl.* (1997) **57**:73–81. doi: 10.1016/S0920-5632(97)00355-1

87. Chandrasekhar S. Stochastic problems in physics and astronomy. *Rev Mod Phys.* (1943) **15**:1–89. doi: 10.1103/RevModPhys.15.1

88. Lehner G. One dimensional random walk with a partially reflecting barrier. *Ann Math Stat.* (1963) **34**:405–12.

89. Blom G, Holst L, Sandell D. *Problems and Snapshots From the World of Probability.* New York, NY: Springer-Verlag (1994).

91. Frieman JA, Turner MS, Huterer D. Dark energy and the accelerating universe. *Annu Rev Astron Astrophys.* (2008) **46**:385–432. doi: 10.1146/annurev.astro.46.060407.145243

92. Wilczek F. Multiversality. *Class Quant Grav.* (2013) **30**:193001. doi: 10.1088/0264-9381/30/19/193001

93. Linde A. A brief history of the multiverse. *Rep Progr Phys.* (2017) **80**:022001. doi: 10.1088/1361-6633/aa50e4

94. Turner MS, Weinberg EJ, Widrow LM. Bubble nucleation in first-order inflation and other cosmological phase transitions. *Phys Rev D.* (1992) **46**:2384–403. doi: 10.1103/PhysRevD.46.2384

95. Mégevand A, Ramírez S. Bubble nucleation and growth in very strong cosmological phase transitions. *Nucl Phys B.* (2017) **919**:74–109. doi: 10.1016/j.nuclphysb.2017.03.009

96. Mégevand A, Ramírez S. Bubble nucleation and growth in slow cosmological phase transitions. *Nucl Phys B.* (2018) **928**:38–71. doi: 10.1016/j.nuclphysb.2018.01.012

97. Guth AH. Inflation and eternal inflation. *Phys. Rep.* (2000) **333-4**:555–74. doi: 10.1016/S0370-1573(00)00037-5

98. Guth AH. Inflationary universe: a possible solution to the horizon and flatness problems. *Phys Rev D.* (1981) **23**:347–56. doi: 10.1103/PhysRevD.23.347

Keywords: nonassociative geometry, discrete spacetime, quasigroups, smooth loops, quantum gravity

Citation: Nesterov AI and Mata H (2019) How Nonassociative Geometry Describes a Discrete Spacetime. *Front. Phys*. 7:32. doi: 10.3389/fphy.2019.00032

Received: 07 December 2018; Accepted: 22 February 2019;

Published: 18 March 2019.

Edited by:

Oleg N. Kirillov, Northumbria University, United KingdomReviewed by:

Alberto Molgado, Universidad Autónoma de San Luis Potosí, MexicoDouglas Alexander Singleton, California State University, Fresno, United States

Copyright © 2019 Nesterov and Mata. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Alexander I. Nesterov, nesterov@cencar.udg.mx