Spinning Toroidal Brane Cosmology; A Classical and Quantum Survey

We construct a cosmological model based on a free particle model which is constrained on an embedded toroidal brane, with a general rotation around a specific axis in the bulk space. Some related issues such as the rotation axis of the brane, the presence of gravitomagnetic background and its relation to the general angular velocity of the brane, and its quantum mechanics and related issues such as minimal length and minimal momentum of the quantum model in the $\mathbb{T}^3$ brane are studied. Also, some cosmological features such as the constraint which is imposed upon the toroidal universe by the isotropy and homogeneity conditions, the corresponding Hubble law, and accelerating expansion for the spinning toroidal model without considering a cosmological constant are also studied.


Introduction
Global properties of the space is not restricted by general relativity, and a universe with a specific local geometry may have various global topologies [1]. Cosmological models such as Friedmann-Robertson-Walker (FRW) models allow the universe to have the same average local properties everywhere, while there would be a freedom to describe a quite different space for the large scale. This is due to the fact that by having a negative or vanishing local curvature, one can not necessarily conclude that our universe is infinite, since the universe still can have a finite volume due to the global topological multi-connectivity feature of its space [2,3].
On the other hand, flat or hyperbolic FRW solutions can describe multi-connected universes. Hence, within the inflationary scenarios, quantum fluctuations can produce compact spaces of constant curvature, both from flat spaces [4], and curved ones [5,6]. This inspires us to think about a multi-connected finite universe, which the simplest one is the torus [7,8].
But, this is not the whole story. The topology of the universe is not the only thing which is missed in standard cosmological models. Particularly, the spinning universe around the specific axis is not physically desirable [9], and it would challenge inflation theory that disagrees the isotropy and homogeneity features of the universe.
On the other hand, some observational evidences such as the temperature anisotropy of the cosmic microwave background (CMB) radiation [10], the polarization of quasars [11][12][13][14], the acceleration of the cosmic expansion [15][16][17][18][19][20], and angular distribution of the fine-structure constant [21][22][23], may indicate preferred directions in the sky. An important observational evidence obtained by studying the Spiral galaxies with a well-defined handedness in a particular region of the sky shows that the universe may spin on a specific axis [24,25].
In this article, we model a spinning toroidal universe and study its features, while the non-rotating toroidal universe is studied before [26]. This would be important since a spinning universe has not been seriously studied via a cosmological theory.
The tool which is used here is the geometrical embedding approach, where the available configuration space for the test particle is embedded in an extended one, and the corresponding physical and geometrical results can be explained by the formalism of constrained systems.
It is worth mentioning that the most approaches dealing with embedding procedures are concentrated on the tools of Riemannian and pseudo-Riemannian geometries and their corresponding induced metrics. Although this work is based on a similar procedure, instead of inducing space(time) metric, phase space metric (symplectic tensor) is induced [27]. Now, by adding the spin to the embedded space in R 4 , we continue the work [26], where ADD extra dimensions problem [29,30], and making a particle model without referring to general relativity [31,32], is concerned. In that article, from the induced metric on a test particle's phase space, a gauge symmetry structure (symplectic tensor with gauge degrees of freedom) is made, in a way where the extra dimensions of the extended space changed to physical ones (or at least gauged ones).
Here, we follow the different line, i.e. we seek the induced phase space metric on the spinning T 3 in R 4 to make a classical and quantum brane cosmology.
This article is organized in the following way. In the first two sections, we make our spinning toroidal universe and extract its symplectic structure. The third section is dedicated to the features which satisfy the isotropy and homogeneity conditions of the spinning toroidal universe. In section four, we deal with quantum cosmological investigations, and quantum mechanics of the model and some of corresponding features such as minimal length are studied. Last section is devoted to the cosmic dynamics, where the Hubble law and the expansion acceleration of the model is calculated. Some phenomenological and numerical investigations are also suggested to verify the model.

Spinning toroidal universe
From brane cosmological point of view, this model is built in the form of a T 3 brane which is embedded in the Euclidean space R 4 . Using the theory of constrained systems, we know how to surf the toroidal brane as the boundary of a 3-dimensional universe 1 [26]. Thus, we only need to add the general spin feature, ω, to the T 3 universe, and consequently to our test particle, living on such a brane.
Now, we define the velocity vector of the test particle, living on the spinning T 3 in R 4 as, The direction of the v ω should be chosen in such a way, not to enter a fundamental or coordinate singularity to the kinetic energy (or the Hamiltonian) function, and at the same time, add the fundamental constant ω (or a new Casimir of the corresponding algebra) to the model. For instance, adding the spin velocity term in the form of v ω =ψr sin ϕω, which results to the term v ψ = r sin ϕψ +v ω , enters a coordinate singularity in Hamiltonian, which is not desirable. We ignore such an addition to avoid singularity problems. Whereas, one can add the angular velocity (spin) ω = ωψ to the test particle's spin, (ω +ψ), without entering any singularity to the Hamiltonian. Hence, the test particle's Lagrangian with unit mass living on such a universe is explained as, L = 1 2 (ṙ 2 + r 2φ2 + r 2 (ω +ψ) 2 sin 2 ϕ + r 2χ2 sin 2 ϕ sin 2 ψ), (1.4) and corresponding velocities are, Having Legendre Transformation, the energy function of the test particle on the embedded brane which is spinning aroundψ, can be obtained as, Here, the term in parenthesis includes a removable regular coordinate singularity, although the added term would be remained secure.
By redefining corresponding momenta, the Hamiltonian (1.6) can be redefined as the Hamiltonian of a free massive particle which is minimally coupled to a gauge potential as, where, G =ψωr sin ϕ, and V (r, ϕ) = − 1 2 r 2 ω 2 sin 2 ϕ, are 1-form and scalar gravitomagnetic potential respectively. The negative nature of this potential, which is similar to the confinement potentials for quarks in QCD, confines the test particle on T 3 .
The direct tensor product of this vector potential and a suitable differential operator, d I , one can obtain the background magnetic 2-form of the model as, where, Therefore, B IJ = ω(r ⊗ψ sin ϕ +φ ⊗ψ cos ϕ), (1.9) will not have any coordinate singularity. Putting the gravitomagnetic pseudovector components in an array, indicates that the presence of ω produces a scalar and a magnetic vector potential, which can be interpreted as the magnetic background resulting the constant noncommutativity, which is consistent to the noncommutativity, obtained from the states of string theory [28]. This magnetic background changes the symplectic structure of the phase space of this model to the the symplectic structure of a noncommutative one, where a part of this noncommutativity is controlled via ω, and the other part, as it has been shown in [26], is due to the Casimirs of the model which are two radii of T 3 , embedded in the Euclidean space R 4 . Also, as we see that the components of the magnetic 2-form depends on a cosmological coordinate, ϕ, with the spatial directions as,r ⊗ψ andφ ⊗ψ.
If we neglect the direction along ψ by imposing the constraints, B can be written in a suitable form as, B = ω sin ϕr + ω cos ϕφ, (1.11) with the norm of, B = ω, which is the product of fundamental spin and unit mass, as the gravitational charge current. Thus, one can suggest that due to the spin of the universe, from the R 4 observer's point of view, for every massive particle, m, there is a gravitomagnetic monopole with the vector charge g = m ω [41]. From T 3 observer's point of view 2 , this gravitomagnetic monopole can be measured as g vis .
2 Symplectic structure of spinning toroidal universe By studying the symplectic structure, the dynamics of the particle's motion, the geometry of the phase space, symmetry generators' Lie algebra of the test particle's motion, and consequently symmetry generators' algebra of phase and configuration spaces are accessible.
In the R 4 observer's cosmological coordinates, the brane of our toroidal universe is described with the following constraint, where, r 2 = x 2 + y 2 + z 2 + s 2 , and ς + , and ς − are 3-torus' radii.
The constraint chain structure of the model includes one more constraint. In other words, the effect of the spin in the phase space appears due to the change in the Hamiltonian as, Here, Poisson bracket matrix of constraints, which is the fundamental part of the symplectic structure of the model, does not include ω.
Therefore, in comparison to the non-spinning toroidal universe [26], the spin parameter of the universe appears linearly in the symplectic structure. If we can be assured of two basic cosmological principles, i.e. the isotropy and homogeneity features of the universe, one can introduce the cosmic R 4 observer as the embedded T 3 observer either. This can be done by calculating Dirac brackets and using coordinate reduction methods.
One should notice that if two radii of our toroidal universe is considered to be nearly equal, then our model can satisfy the homogeneity condition by approximation. Also, if we can find the non-unique spin direction, including some free parameters, we can think that our model can be estimated as an isotropic universe.
As we know, by considering two constraints of the phase and configuration spaces, Dirac brackets of basic components of phase space as z µ := (r, ϕ, ψ, χ) ⊕ (p r , p ϕ , p ψ , p χ ) can be obtained as, Now, by imposing the constraints and calculating Dirac brackets, one can redefine the following couple of variables, where can be ignored to make a 3-dimensional coordinate, i.e. r, ϕ and χ, for our 3dimensional test particle. This procedure imports the constraints into the dynamics of the system consistently, i.e. the reduced physical phase space is not separated from the R 4 ⊕ R 4 initial phase space. The initial phase space is important, since it determines the limited domain available for the remained parameters of the brane, such as r, and also it determines the metric (and the symplectic structure) of the smaller phase space. For instance, in initial phase space R 4 ⊕ R 4 , Poisson brackets matrix of basic parameters is obtained as the following symplectic matrix, which indicates the canonical structure of phase and configuration spaces. But, after reducing the phase space to T 3 ⊕ M, and inducing the symplectic structure (phase space metric), we obtain a new J µν (z) which depends on the phase space variables. Due to special form of the constraints, Φ 1 and Φ 2 , and their Poisson brackets matrix, we will see that the geometry of the phase space is noncommutative and non-canonical, although the geometry of the configuration space remains commutative.
Using (2.4), we derive non-vanishing Dirac brackets as, Using the canonical Hamiltonian and the obtained symplectic structure, one can calculate every dynamical variable such as (r,ṙ), and then move on the constrained surface (T 3 hypersurface in configuration space and R 3 hypersurface in phase space).

Isotropic spinning torus universe
In previous part, we found the spin direction of the 3-dimensional universe. Since this direction is located in the coordinates of a 4-dimensional universe, a 3-dimensional observer living on the toroidal universe cannot perceive this direction completely. Thus, in order to find the portion of the direction of the spin, which is seen by the 3-dimensional observer, we follow the following root.
As it is mentioned, after obtaining the dynamical equations and algebraic calculations, one can move on the constrained surface of Φ 1 and Φ 2 , and omit the couple variable (ψ, p ψ ). But, one should notice that such omitted variables cannot be considered as completely nonphysical ones, since the invisible spin of the universe is directed toward ψ. Hence, this direction and its corresponding variable should not be ignored completely.
On the other hand, if we want to ignore a canonical couple for the sake of the second order constraints Φ 1 and Φ 2 , one should have them symplected in the formalism of constrained systems. Hence, to sacrifice a canonical couple to reduce the configuration space from R 4 to T 3 , we redefine the hypersurface T 3 , embedded in R 4 as, Hence, the visible direction of such a universe tangents to this hypersurface and the direction in the configuration space (the coordinate part of the phase space), which is invisible for the T 3 observer, is perpendicular to this hypersurface asn invis = ∇Φ 1 |∇Φ 1 | , As we see,χ andφ are not included in (3.2), which indicates that those directions and their corresponding angles remain in the reduced T 3 .
On the other hand, from the R 4 observer's point of view, the spin direction of T 3 is directed towardψ, and we expected that the directionψ to be invisible for the T 3 observer. But, we see that the T 3 observer can also detect the spin of the universe and somehow its direction.
From R 4 observer's point of view, the invisible direction for T 3 can be obtained by solving the following relation, Therefore, the portion of the spin (according to the spherical coordinate R 4 ) which is seen by the real observer is obtained as, As we see, in the above direction, there is still a small portion forψ, which means from the R 4 observer's point of view, the T 3 observer can see the spin vectorΩ = ωψ. Nevertheless, the presence of free and arbitrary parameters u ϕ , u χ in (3.4), indicates the fact that the embedded T 3 universe is isotropic according to its corresponding observer.

Toward the spin direction and the fundamental constant ω
As we said before, for R 4 observer, the spin of the embedded T 3 is, (3.5) and the observer in toroidal universe sees, It should be noticed thatr andψ are R 4 unit vectors. Thus, no one can claim that the spin direction and its value are completely observed by T 3 observer. Hence, to obtain the exact value, observed by the T 3 observer, one should go on the constraints surface, i.e. ψ should be solved from Φ 1 = 0. Then, by replacing the obtained ψ in the (3.6), in the homogeneous limit of the universe, ǫ → 0, and for far or near distances (IR or UV), we would have, For the homogeneous limit of the universe, (ǫ → 0), and far distances, ς = 1, (r−2) → 0, the spin of the universe is obtained as, Since the observer on T 3 , does not have any liberty to choose theφ direction, one can have the average volume over the space as, and vanish its dependence to the ϕ, where, the star sign, ⋆, indicates the spatial average in far distances, ((r − 2) → 0), and α is a numerical coefficient obtained via the spatial average volume integration, It is interesting to see the effect of inhomogeneity of the universe, ǫ, on the fundamental constant of the universe's spin (or its Casimir), ω in (3.10). In other words, ω can be regarded as the bare spin in primordial field theories, without having the spacetime interaction or renormalization, and Ω (T 3 ) vis ⋆ is the observable and renormalized spin.
Similar to QED where mass and charge of the electron is corrected by the coefficient α QED = 1 137 , in our model the intrinsic spin of the universe is also corrected via ǫ. This correspondence (or similarity) may intuit us to have a deeper understanding of the results of this model.

Homogeneous torus universe
In algebraic geometry presentation of the universe with topology of T 3 , from the R 4 observer's point of view, the radial component, r, is finite.  Table 1: a correspondence (not an exact duality in a perfect sense) between a typical field theory and present particle model Also, for the observer in R 4 , the UV region 0 ≤ r < ς + − ς − , especially the location r = 0 is not available by the T 3 observer, and without imposing a fine tuning condition the overall universe will not be homogeneous. In other words, our visible universe with minimal length ǫ, can be considered as a homogeneous T 3 universe, with infinitesimally unequal radii, and the parameter ǫ controls the inhomogeneity. Hence, we introduce a special scheme to introduce the 3 dimensional brane as a T 3 shaped universe which is deviation from R 4 with two almost equal radii, ς + = ς and ς − = (1 − ǫ)ς. As we mentioned, the geometrical deformation and inhomogeneity of such universe is controlled by ǫ.

Quantum cosmology investigations
To study the infinitesimal structure of space by a test particle which is localized in a small region, quantum mechanics is a suitable tool. Although the test particle's coordinates which are used here is from the R 4 observer's view, transforming to the particle's local coordinates would be an easy task.
Therefore, the origin can be chosen by the local observer in T 3 , with the uncertainty of the order of ςǫ. This observer can put his test particle in the origin and investigate the fine structure of the space, by studying its corresponding quantum mechanics. Also, he can find the proper estimate for the parameters ǫ, ς and ω of the model [43].
Obtaining the symplectic algebra of the model, J µν (z), the Heisenberg algebra of the main quantum observables in the UV limit can be obtained as, where, J µν (0) is the constant matrix with symplectic structure, and J µν (r) is the first order approximation of the symplectic structure for ǫ = 0, which will be done only for the dimensionless radial componentr.
Using Dirac brackets, the observer may make the quantum mechanics for the test particle to clearly describe the structure of the phase space of the model in quantum scale. Hence, by a good approximation, the Poisson structure J µν can be quantized.
By scaling the radius component of the model, r, and its corresponding momenta, p r , with the specific length ς of the universe, one can have the quantized version of the model.
[r,p r ] = iℏ 1 −r 2 (1 + ǫ) sin 2 ϕ , Comparing commutators above with the common structure of quantum mechanics variables in pure 3-dimensional space of R 3 , one can see how the embedding procedure changes the algebra of basic variables of the model. Also, we can see that the inhomogeneity of T 3 is entered to the above algebra.
One can also compare these brackets with the non-rotating torus universe, previously modelled in [26]. We can see that except three last commutators, including ω, all commutators are the same as before. Moreover, one can compare these results with similar type of Dirac brackets for particle on a torus knot in [44].
Moreover, although the test particle is confined on a 3-dimensional space, due to the embedding procedure, the above mentioned quantum commutators includeψ andp ψ which can be regarded as hidden variables (like Einstein hidden variables) that effect the quantum state of the particle.
On the other hand, as we see ω, which indicates the spin of T 3 , is appeared in three commutators. Two of them include ω multiplied by the great radius of the universe, ς, which can be ignored by a renormalizing procedure which vanish infinities. In other words, we claim that those terms including ω does not have any observable effect, but since they also includep ψ , one can say that via those commutators, the momentum is injected form the hidden dimension of the universe to T 3 and spins the brane.
Nevertheless, the ω is appeared multiplied by ǫ in the commutator of variables (p r , p ϕ ), and hence it can employ an observable quantum mechanical effect. Also, all above mentioned commutators can be interpreted in such a way where due to the embedding T 3 and its spin in R 4 , the Planck constant is corrected. This claim can be investigated in high energy (tiny distances) regime. By the way, these forms of commutators in this model, inspire us to investigate IR/UV mixing [33].

Minimal length in T 3 universe
Like several works [34][35][36][37][38][39], we see that in quantum commutators of our model (4.3), which are deviated from the common Heisenberg algebra, a minimal length is appeared. Hence, one could expand this symplectic algebra according to the dimension parameter of minimal length ǫ, by calculating them in quantum regime of small lengths,r = r ς , and investigate such a claim 3 .
We think that the amount of inhomogeneity and deformation in this universe, ǫ , can be estimated with the help of the cosmic data. Also, the fundamental length of this model can be regarded as the Hubble length, as R H ≈ 2ςN , where N is the winding number on the periodic T 3 [40].
The above-mentioned basic commutators help us to study the dynamics of the quantum test particle in a small localized region around the observer. This means that the test particle here acts as an elementary particle, like an electron, to test the effect of the topology and spin of the universe on the fine structure of the space.
On the other hand, as we will see further, this test particle can be regarded as a far universe (from the R 4 observer that is coincided to the origin of T 3 observer), where studying its dynamics may help us to understand the effect of the topology and spin of the universe on cosmic parameters.
The advantage of choosing the spherical coordinates to calculate basic operators in quantum Hilbert space is that one can easily obtain those commutators in the infinitesimal quantum distance approximations, by expanding the results around r = 0, where we can compare our results to the minimal lengths, and minimal momentum quantum mechanics 4 [37][38][39]. Although models dealing with such phenomena, consider different issues such as correction and renormalization of Planck scale in high energy physics and their effects in quantum mechanics in curved spacetime, we consider the topology of the universe and the extrinsic geometry imposed via embedment, and also the spin of the universe to have the basic commutators changed.
Thus, the quantum mechanics obtained from our model which is expressed by the basic commutators algebra (4.3), includes two fundamental constants, ς and ω, and can be compared to Amelino-Commelia models of doubly special relativity, where it is based on the change in the particles' energy-momentum dispersion in small scales and high energies [45][46][47]. The symmetry properties and the algebra of symmetry generators of our model can be compared with [48][49][50][51]. Moreover, our results can be compared to the models including extra fundamental parameters rather than ℏ and c, such as those models studying the Lorentz symmetry below Planck scale [52][53][54].

Omitting the canonical couple (ψ,p ψ )
After quantizing a physical states in the formalism of constrained systems, solving the Schrödinger equation solely can not completely determine all quantum states as H | phys >= E | phys >. This happens due to the presence of gauge symmetry which is imported by first class constraints. Hence, one should have them in a more limited way of constrained equations asF C | phys >= 0, whereF C is the operator form of first class constraints [26].
Hence, from the equivalence sets of states, first class constraints extract a complete set in the operator form. The operator form of constraints effects the states obtained from Hamiltonian, similar to the way of choosing a specific gauge and breaking the gauge symmetry.
In second class systems, we reduce the basic variables (and also quantum observable states), and with the help of those reduced variables and their algebra we find a proper representation for basic observables. Then we obtain the reduced Hamiltonian, H red , defined according to the reduced observables, and rewrite the Schrödinger equation as, H red | phys >= E | phys >, and then we can have the states.
The above mentioned method can be used only if by calculating Dirac brackets one can omit a canonical couple, such as two second class constraints as For this model, one can omit the canonical couple (ψ,p ψ ) by imposing the constraints, although their effect on the model would not vanish, e.g. the interaction ofp ψ with other basic observables indicate that this operator conveys the intrinsic spin of the universe to the observable part of T 3 . Therefore, being embedded in R 4 does not only effects the T 3 Hamiltonian, and it also effects the universe's observable representations. In other words, although this embedding is explained by second class constraints, it is also a gauge symmetry, and (ψ,p ψ ) fixes this symmetry [56][57][58]. 5 As we see in (4.3), ω is not included in the commutator of (r,p r ), and it is only appeared in the commutators of a phase space variable and one of angular momenta. Therefore one can conclude that although we got rid of operators (ψ,p ψ ) by using constrained relations, their effect to the algebra is not negligible, since they have (speciallyp ψ has) a nonvanishing algebra with the basic observables of T 3 . Hence, ω can be regarded as a new Casimir which is entered to the algebra of observables of the model.

Cosmic Dynamics
Now, we can study the dynamics of a cosmic object, such as a galaxy, to detect the local and global structure of the space, by studying its Hamilton function and symplectic structure, and reading off the Hamilton equations of motion for far distance object as, r = {r, H} * ,ṗ r = {p r , H} * . (5.1) On a universe, where the T 3 observer does not sense any deviation or angular rotation, the test particle's energy is negligible in comparison to the parameters of the universe. As a mater of fact, the total energy of the universe is constructed via its geometrical structure, ς , and movement, ω. Hence, we consider that the energy of a sample universe, such as our test particle, is almost zero, From the first Hamilton equations (5.1), one can obtain the relation between distance and velocity, (the Hubble law), and from the second one the acceleration of the universe is obtained.ṙ = v cosmos r = p r (r, ϕ),r = a cosmos r (r, p r , ϕ).

(5.3)
Obtaining p r from (5.2), we encounter several situations. In near distances, v cosmos r ∼ r 2 , and for distances near to the boundary of periodic universe T 3 , v cosmos r ∼ r − 3 2 . Due to some observational evidences 6 , extreme points, i.e. r ≈ 0 and r ≈ 2ς are put away. Hence, one can calculate around a median point, such as r ≈ ς. Thus, for the velocity of the test particle with respect to the first order of ς, we have, Also, for the acceleration to the zeroth order of ς, we obtain, a cosmos = 2ςω 2 sin 2 φ (2 cos 2φ − 5) 3 ((4 − 64ǫ) cos 2φ + 64ǫ + 2 cos 4φ − 33). (5.5) The angle ϕ is a dynamical variable, but since we consider the universe in a particular time, we can think of a constant angle, having an average volume integrationF = The constant term in (5.6), which includes the T 3 universe parameters ς, ω and ǫ, is appeared due to the translation velocity of T 3 universe, from the R 4 observer. The second term, can be regarded as the Hubble law. As we see, the Hubble constant in this model is independent of the periodic radius of this universe, 2ς, and to this order of approximation is only obtained via the spin of the universe (i.e. no spin, no Hubble law).
Here, the Hubble constant is renormalized up to the first order of ǫ. This means that if we desire not to have a blue shift for the distances greater than the half of the radius of the universe, ǫ < 0.45. Also, the positive acceleration limit, indicates a better limit of ǫ < 0.37.
It is worth mentioning that the Hubble law gets a nonlinear form in the extreme point r = 2ς and near boundaries as v cosmos r ∼ (r − 2ς) −3 2 , which is another testimony for our periodic T 3 universe, embedded in R 4 .
In r ∼ ς, (5.6) and (5.7) can be written in the form of its two characteristics of the model, i.e. ς and ω as 7 , It is evident that in distances smaller than ς, the correction of the Hubble law gets smaller than the common form, and the acceleration decreases infinitesimally from a constant value.
Nevertheless, to study the expansion of the space for a distance δr ≪ ς, (e.g. r = 0.01ς) we have, δa ≈ − δr ς (2.83ǫ + 5.08)a ch . (5.11) Hence, the far universe decelerates, and this to much less than a percent of the constant acceleration we obtained before (5.7).
Therefore, the universe near the origin experiences the greater acceleration than a far one, and form the observer's point of view the final state of the universe seems to have accelerating expansion.
Also, using the relation (5.10), for r = 0, we can obtain the classical big bang acceleration as, a bigbang = a ch (6.25 − 0.28ǫ). (5.12) But, what about the acceleration near the boundary, (r = 2ς)? As one can see, the universe experiences a big negative acceleration as, Here, v ch and a ch indicate the characteristic velocity and acceleration respectively.
which happens due to the periodicity of T 3 universe, and can be investigated in further works related to cyclic models [60][61][62].
As we see, (5.13) contains an infinite and singular term, which will be thrown away, and a negative one, which returns the far universe to the permitted area of r < 2ς, where the second half period of a periodic universe appears.
This return and having a negative acceleration in boundaries, can be regarded as a big crunch (in r ≈ 2ς) that happens in a periodic universe after a big bang (in r = 0).
Thus, with the help of embedding approach and constrained formalism, without using field theory or general relativity and its corrections such as f (R) gravity, we obtained a big crunch phase for our model, and all that is done only by embedding a T 3 in R 4 , without considering a potential or a cosmological constant [64].

Some possible phenomenological and numerical investigations
Our cosmological model here is based upon the particle nature (not the field) and its interactions (in quantum scale of a spinless, chargeless and without any other quantum feature of quantum particles, containing only mass, unit inertia, and in large scales as a far universe with the negligible energy in comparison to the energy of the universe) in a background space (not spacetime).
Since we did not use any general relativistic and field theory formalism, no cosmological constant is entered to our model. Anyway, our test particle (which is an observing far universe) experiences an acceleration for the universe. As we showed, we obtained a linear expansion of a positive acceleration over the median point of a periodic universe. Fitting this acceleration with observational data, one can tune a ch and find its proper limit.
Moreover, since we are in a positive acceleration era, the positivity condition imposes a limit on the inhomogeneity parameter of the universe. Also, we have another limit which is imposed on ǫ by the Hubble law of our model. Fitting this parameter with observational data imposes a limit on v ch , according to the half radius of the periodic universe. One should notice that none of those fittings will provide an absolute value for ς.
But, if we assume that there exist a fine tuning between two radii of the T 3 , the quantum structure of spacetime introduces a fundamental minimal length as ς + − ς − = ǫς. Hence, one can estimate the order of magnitude of ς by studying the algebra of quantum commutators in (4.3), which will be important in at the order of Planck length.
Thus, by extracting the minimal length of noncommutative algebra of quantum commutators in our model, and comparing them with the quantum gravity specific length (Planck length), one can find a proper limit of ς. In addition, the above mentioned algebra which includes ǫω also provides a minimal momentum, that may help us to impose limits for ς, ǫ and ω.
Also, studying the quantum harmonic oscillator of the model may help us to relate the model to the CMB fluctuations.

Discussion
In this article we modelled a spinning, finite and periodic cosmological universe based on a T 3 , embedded in a R 4 space, where its inhomogeneity and deformation is controlled via the dimensionless parameter ǫ, as the difference of two torus' radii.
The spin vector is chosen in a way not to enter a fundamental or coordinate singularity to the Hamiltonian function, and at the same time adds the fundamental constant angular velocity, ω, to the model. Also, we showed that that the presence of ω produces a scalar and a magnetic vector potential, which can be interpreted as the magnetic background resulting the constant noncommutativity. This magnetic background changes the symplectic structure of the phase space of this model to the the symplectic structure of a noncommutative one, where a part of this noncommutativity is controlled via ω, and the other part, as it has been shown in [26], is due to the Casimirs, i.e. ς 1 and ς 2 , radii of T 3 .
One should notice that we assumed two radii of our toroidal universe to be nearly equal, then our model satisfies the homogeneity condition by approximation. We also showed how the homogeneity of the T 3 universe, with infinitesimally unequal radii, is controlled by the parameter ǫ. Also, we found the non-unique spin direction as a preferable direction in space, including some free parameters, which control the universe's anisotropy. Hence, with the help of that and assuming the isotropic condition of the universe, we obtained the invisible spin direction for T 3 , from the R 4 observer's point of view.
Then, by going to the constrained surface and reducing a couple of phase space coordinates, we imported the constraints into the dynamics of the system consistently.
Afterwards, by using Dirac brackets and quantizing Poisson structure of the model, we clearly describe the structure of the phase space of the model in quantum scale. This is the starting point of making the quantum mechanics as (4.3), which is the noncommutative fundamental quantum commutators of the model, although we did not present any representation for these fundamental operators. Quantum commutators obtained via this procedure show that the embedding procedure changes the algebra of basic variables of the model and the inhomogeneity of T 3 is entered to the obtained algebra.
The form of commutators indicate that ω can employ an observable quantum mechanical effects and in our model the momentum is injected form the hidden dimension of the universe to the T 3 , which spins the brane. This can be interpreted in such a way that due to the embedding T 3 and its spin in R 4 , the Planck constant is corrected. This claim can be investigated in high energy regime and inspires us to investigate IR/UV mixing. We can also see that quantum commutators include a minimal length and minimal momentum, and thus the model can be compared to some models based on doubly special relativity.
Calculating the velocity of the test particle, we obtained the corresponding Hubble law in our model. As we see, the Hubble constant which is normalized up to the first order of ǫ, is independent of the periodic radius of this universe, 2ς, and only depends to the spin of the universe.
Moreover, to study the expansion of the space in different distances we calculated the acceleration and its variation, which shows that the universe near the origin experiences the greater acceleration than a far one, and form the observer's point of view the final state of the universe seems to experience an accelerating expansion. This would be an important result since we obtained the acceleration of the universe without considering a cosmological constant.
Also, we calculated a positive classical acceleration at (r = 0) and a big deceleration near the boundary, (r = 2ς). Hence, one can say that in origin of the model we have a big bang and after a half period in r ≈ 2ς we encounter a big crunch which is the feature of a periodic universe. We hope that this feature of the model can be investigated in further works related to cyclic and bouncing cosmologies.