Parametric equations of the expansion of a closed universe with torsion
aa r X i v : . [ g r- q c ] M a r Parametric equations of the expansion of a closed universe with torsion
Prastik Mohanraj
Engineering and Science University Magnet School,500 Boston Post Road, West Haven, Connecticut 06516, USA ∗ Using the Friedmann equations for a closed universe with spin and torsion, parametric equationsfor the scale factor and time are derived for the early universe.
I. INTRODUCTION
Analogous to the Doppler effect perceived with sound waves, where approaching sound waves increase in pitch andreceding sound waves decrease in pitch, the Doppler effect can be perceived with electromagnetic radiation, or light,where approaching waves increase in frequency and receding waves decrease in frequency. An increase in frequencycorresponds to ”bluer” light, making approaching waves experience blueshift. Similarly, receding waves decrease infrequency, producing ”redder” light, and thus experiencing redshift.When galaxies are observed in the sky, over time the light coming from them becomes redder in color, thusthese galaxies are experiencing redshift [1]. In fact, all galaxies in the universe except the galaxies in the neighbor-hood of our own Milky Way appear to experience redshift. As a result, it can be concluded that every galaxy isreceding from our perspective. The recession of galaxies is described by Hubble’s law, which dictates that the speedof recession of a point in space, relative to any observer, is directly proportional to the distance between the twopoints. This law, together with the assumption that the universe is homogeneous and isotropic at largest scales (thecosmological principle), leads to a conclusion that the space itself is stretching by a constant factor over time, thusincreasing distance between two points and causing proportional recession between these points. It can be picturedas how two points on a balloon seem to be moving away from each other when the balloon expands, although thepoints themselves are not experiencing motion, but the surface area on the balloon itself is increasing.Since all these points are receding from one another, in that the space between these points is rapidly increasing, theuniverse is expanding. If the universe is expanding, it must have been much smaller, hotter, and denser in the distantpast. Thus, it is likely that the universe started from an infinitesimally small size, and later undergone rapid ex-pansion to create the universe’s size we observe today. This idea is the foundation of the standard Big Bang cosmology.There are currently three models of a homogeneous and isotropic universe: positively curved universe, nega-tively curved universe, and universe with zero curvature [2]. The 3D plane of the universe becomes curved aroundthe 4D hyperspace of the universe in one of these three models, much like the 2D surface of a balloon wraps aroundthe 3D interior space of the balloon. A negatively curved universe is open, saddle shaped, with a lower energy densitythan the critical energy density. A universe with zero curvature is a flat, Euclidean space with its energy densityequal to the critical energy density. A positively curved universe is closed and spherical, with higher energy densitythan the critical energy density. In this paper, the positively curved, hyper-spherical universe model will be analyzed.For a closed universe, the scale factor a ( t ) is is a function of the cosmic time t , acting as the radius of the4D hyper-sphere whose 3D hyper-surface is the closed universe itself. Hubble’s law states that the rate of expansionof the universe obeys ˙ a = H · a , where H is the Hubble parameter and a dot denotes the derivative with respect to t .The critical energy density is equal to ǫ c = (3 H )(8 πG ) . Units will be used such that the speed of light c = 1.This paper establishes the novel finding of a closed universe with spin and torsion in the Einstein-Cartantheory of gravity, which extends general relativity to include spin of elementary particles. At extremely highdensities, torsion manifests itself as a gravitational repulsive force, which avoids the Big-Bang singularity and replacesthe Big Bang with a bounce that never reaches a singularity. We found the solution to the equation which describesthe expansion of the universe right after the bounce, in the torsion-dominated era. ∗ Electronic address: [email protected]
In Section II, we will see the Friedmann equations and their forms in the different stages of the universe’sevolution. In Section III, we will examine the matter-dominated universe stage of evolution, and present the cycloidas a solution to this evolutionary stage. In Section IV, we will edit the cycloid to form the trochoid, and presentthe novel solution of using a trochoid for the torsion-dominated era of the universe. In Section V, we will applythe torsion-dominated trochoid model to the radiation-dominated model, where the torsion factor has cancelled out,and finally, in Section VI, we will explore the last stage of evolution, the dark-energy-dominated universe, with theFriedmann equations. Section VII provides a conclusion and summary of the analyses presented here.We will now introduce the Friedmann equations, a background of cycloids, which will be crucial to the fol-lowing analyses, and solve for parametric equations in the different eras of the universe, ultimately reaching theparametric equations for the beginning of the universe during the torsion-dominated era.
II. FRIEDMANN EQUATIONS
The Friedmann equations are an extension of the general theory of relativity, applied to the universe on a cosmologicalscale. These equations can be analyzed to develop equations of scale factor as a function of time. The easiestapproach is to develop parametric equations of scale factor and time, in terms of another variable. For reasons laterdescribed, these parametric equations will be derived as a result of similarities between the Friedmann equations fora positively curved universe and the derivation of parametric equations for cycloidal/trochoidal motion.The first Friedmann equation can be written as [2, 3]( ˙ a ) + k = 13 κǫa . (1)Here, ǫ is the energy density of matter in the universe, dominated by relativistic light energy in the early universeand non-relativistic energy of massive particles in the later universe. Also, κ = 8 πG . In this analysis, the universewill be defined as having positive curvature, thus k = +1. The second Friedmann equation can be written as [4]( ˙ ǫa ) + p ( ˙ a ) = f ( a, ǫ ) , (2)where p is the pressure of matter. The term f ( a, ǫ ) is effectively equal to 0, or is adiabatic, for all stages of theuniverse except for the earliest stages where quantum production of particles increases entropy.For relativistic matter, p = ǫ , so Equation (2) can be modified to yield a ˙ ǫ + 4 ˙ aǫ = 0, thus ǫ ∝ a , thus( ˙ a ) + 1 = 2 Aa (3)for the radiation-dominated era of the universe, where the positive constant A ∝ κ , or A ∝ πG . Similarly, fornon-relativistic matter, p = 0, so Equation (2) can be modified to yield a ˙ ǫ + 3 ˙ aǫ = 0, thus ǫ ∝ a , thus( ˙ a ) + 1 = 2 Ca (4)for matter-dominated era of the universe, where the positive constant C ∝ κ , or C ∝ πG . The usage of coefficients2 A and 2 C as opposed to A and C is simply for convenience.These equations, both with the term ( ˙ a ) + 1, resemble the parametric equations of x and y coordinates of acycloid/trochoid graphed on a Cartesian plane. As a result, scale factor a and time t will be solved for in terms ofparametric equations with a common variable θ , similar to parametric equations of cycloids/trochoids with variable θ .We will now see the comparison between cycloidal/trochoidal behavior and these forms of the Friedmann equations. III. MATTER-DOMINATED UNIVERSE: CYCLOID
If we define a point on the circumference of a circle, and roll that circle along the plane in linear motion, the pathtraced by the individual positions of that point create the graph of a cycloid.A cycloid can be imagined as the graph of the continuous positions of a point (0 ,
0) on a circle of radius r ,centered at (0 , r ), as the circle rotates along the x axis in the positive x direction. If the circle rotates along the x axis such that the center of the circle moves a distance of l units from (0 , r ) to ( l, r ), then the circle rotates clockwiseby an angle of lr = θ radians to the vertical. As a result, the x coordinate of the point (0 ,
0) increases by the distancetraveled by the center, or l = rθ units, and decreases by the clockwise horizontal distance traveled from the circle’svertical diameter, or r sin( θ ). So, the x coordinate can be represented as x = r ( θ − sin( θ )) . (5)Additionally, the y coordinate of the point (0 ,
0) increases by the difference between the radius and clockwise verticaldistance from the circle’s horizontal diameter, or r − r cos( θ ). So, the y coordinate can be represented as y = r (1 − cos( θ )) . (6)Trochoids come as a variation of the cycloid theme, and help us reach the parametric equations for the beginning of theuniverse. Here, the impact of the cycloid and trochoid is the most important in determining the parametric equations.We can replace the y coordinate with a , the scale factor, the x coordinate with t , time, and the constant r with C from Equation (4). This alteration to Equations (5) and (6) from the cycloid model, then we obtain thepossible parametric equations for a and t , in terms of variable θ , as t = C ( θ − sin( θ )) , (7) a = C (1 − cos( θ )) . (8)Using these equations, we can produce the following equations for derivatives with respect to θdadθ = C sin( θ )) ,dtdθ = C (1 − cos( θ )) . Breaking up the expression ( ˙ a ) + 1 into the equivalence( ˙ a ) + 1 = (cid:0) dadθ (cid:1) (cid:0) dtdθ (cid:1) + 1and subsequently plugging in the above derivatives gives( ˙ a ) + 1 = ( C sin( θ )) ( C (1 − cos( θ ))) + 1 = sin ( θ )(1 − cos( θ )) + 1 = (sin ( θ )) + (1 − cos( θ )) (1 − cos( θ )) , ( ˙ a ) + 1 = (sin ( θ )) + (1 − θ ) + cos ( θ ))(1 − cos( θ )) = sin ( θ ) + cos ( θ ) + 1 − θ )(1 − cos( θ )) . Using the Pythagorean Identity sin ( θ ) + cos ( θ ) = 1, we obtain( ˙ a ) + 1 = 1 + 1 − θ )(1 − cos( θ )) = 2(1 − cos( θ ))(1 − cos( θ )) = 21 − cos( θ ) = 2 CC (1 − cos( θ )) . Noting that the denominator of the last expression in the above equation equals Equation (8), we derive [2, 3]( ˙ a ) + 1 = 2 Ca .
We have just derived Equation (4), which is the Friedmann equation for the non-relativistic, matter-dominatedera of the universe’s expansion, using the parametric equations of a cycloid. It turns out the cycloid is a perfectsolution to the model given by the Friedmann equation for the non-relativistic, matter-dominated era. A note-worthyidentity that can be realized from the equations for a and t is that the equations for dtdθ and a are identical, thus a = dtdθ .A problem that arises with the cycloid solution is that there exists periodic points where a = 0, at the mini-mum cusps in a cycloid graph. Thus, all of the matter and energy in the universe is being compressed into aninfinitesimal point, indicative of infinite density. The Einstein-Cartan theory of gravity solves this problem. Thistheory extends general relativity by removing its symmetry constraint on the affine connection, which changes theconservation law for the orbital angular momentum into the conservation law for the orbital plus spin angularmomentum [5, 6]. The Einstein-Cartan theory thus admits the exchange between the orbital and spin components ofangular momentum (spinorbit interaction) in curved space-time, which is a result of the Dirac equation in relativisticquantum mechanics. At extremely high densities, torsion is a gravitationally repulsive interaction eliminatingsingularities [7–11]. We will see this phenomenon applied in the next section, where a novel solution is presented. IV. TORSION-DOMINATED UNIVERSE: TROCHOID
Alterations will be performed to Equation (3), since the torsion-dominated era at the beginning of the universe wouldhave had no matter present, only pure energy and radiation. The torsion factor is inversely proportional to a , givingthe Friedmann equation for this era [4, 9–12] as( ˙ a ) + 1 = 2 Aa − ( A − B ) a . (9)Notice how having A ' B results in the torsion factor making a significant impact on the function only for incrediblysmall values of a , which occurred in the tiny moments following the beginning of time. Here, instead of a singularitybeing formed as in the Big Bang Theory, a bounce is created at the minimum value of the universe’s size, sincethe function never reaches the value of 0. This creates an alternative Big Bounce. The torsion factor approaches 0rapidly after the Big Bounce, so the term Aa makes the only impact following this event. As a result, after the BigBounce occurs, the torsion factor can be eliminated from the model without significant alterations in the results.Now, we can alter the cycloid function to reflect this change, in that we can make the cycloid function nothave a equal to 0 at any time t . Instead of Equations (7) and (8) having the same coefficient for all expanded terms,we can make each term have a different coefficient. The new parametric equations are then t = Aθ − B sin( θ ) , (10) a = A − B cos( θ ) . (11)Note that as long as positive constants A > B , the function for a will always take a value greater than 0. Giventhe condition of A ' B from the above alteration of the Friedmann equation, a will reach very close to 0, but neveractually meet it. As a result, a singularity never forms with this model, creating the Big Bounce model discussedearlier. Note how the equivalence a = dtdθ is preserved in these new alterations.Here, we have formed the basic model of a trochoid, a closely related function to cycloids. This form of tro-choids assumes only positive values, which is required for a real universe that never reaches singularities. A slightvariation on the above trochoidal equation for a by writing the equation in a form yields the equation a = A + B cos(2 θ ) . (12)It was hypothesized that Equation (12) is a parametric model of the Friedmann equation from Equation (9), in thatequation (12) is a solution to the torsion-dominated era of the universe’s expansion. This hypothesis yields a = ( A + B cos(2 θ )) ,dadθ = − B sin( θ ) cos( θ ) p A + B cos(2 θ ) . Squaring the above equation, and using the half-angle identities sin( θ ) = q − cos(2 θ )2 and cos( θ ) = q θ )2 yields (cid:0) dadθ (cid:1) = B (1 − cos (2 θ )) A + B cos(2 θ ) . Plugging these above results into ( ˙ a ) + 1 = (cid:0) dadθ (cid:1) (cid:0) dtdθ (cid:1) + 1 = 2 Aa − ( A − B ) a from Equation (9) yields (cid:0) B (1 − cos (2 θ )) A + B cos(2 θ ) (cid:1)(cid:0) dtdθ (cid:1) + 1 = 2 AA + B cos(2 θ ) − ( A − B )( A + B cos(2 θ )) = A + 2 AB cos(2 θ ) + B ( A + B cos(2 θ )) , (cid:0) B (1 − cos (2 θ )) A + B cos(2 θ ) (cid:1)(cid:0) dtdθ (cid:1) = ( A + 2 AB cos(2 θ ) + B ) − (( A + B cos(2 θ )) )( A + B cos(2 θ )) , (cid:0) B (1 − cos (2 θ )) A + B cos(2 θ ) (cid:1)(cid:0) dtdθ (cid:1) = ( A + 2 AB cos(2 θ ) + B ) − ( A + 2 AB cos(2 θ ) + B cos (2 θ ))( A + B cos(2 θ )) = B (1 − cos (2 θ ))( A + B cos(2 θ )) , (cid:0) dtdθ (cid:1) = 1 A + B cos(2 θ ) , (cid:0) dtdθ (cid:1) = A + B cos(2 θ ) , (13) t = Z (cid:0)p A + B cos(2 θ ) (cid:1) dθ. (14)The solution to this equation involves the usage of an elliptic integral of the second kind, denoted as E ( θ, ψ ) = Z θ q − ψ sin ( φ ) dφ. (15)The solution to Equation (14) is t = √ A + B · E ( θ, ψ ) , (16)where ψ = BA + B .The hypothesis for a ( θ ) and the resulting equation for t ( θ ) will now be checked in the original Equation (9)to determine accuracy. We wish to show that the usage of Equation (12) and Equation (16), when substituted into( ˙ a ) + 1, yield Aa − ( A − B ) a , representing the original Friedmann equation from Equation (9).Plugging in the above equations for (cid:0) dtdθ (cid:1) and (cid:0) dadθ (cid:1) into the expanded form of ( ˙ a ) + 1 yields( ˙ a ) + 1 = (cid:0) dadθ (cid:1) (cid:0) dtdθ (cid:1) + 1 = (cid:0) B (1 − cos (2 θ )) A + B cos(2 θ ) (cid:1)(cid:0) A + B cos(2 θ ) (cid:1) + 1 = B (1 − cos (2 θ )) (cid:0) A + B cos(2 θ ) (cid:1) + 1 = B (1 − cos (2 θ )) + (cid:0) A + B cos(2 θ ) (cid:1) (cid:0) A + B cos(2 θ ) (cid:1) , ( ˙ a ) + 1 = A + 2 AB cos(2 θ ) + B (cid:0) A + B cos(2 θ ) (cid:1) = A + 2 AB cos(2 θ ) + A − A + B (cid:0) A + B cos(2 θ ) (cid:1) , ( ˙ a ) + 1 = 2 A ( A + B cos(2 θ )) (cid:0) A + B cos(2 θ ) (cid:1) − ( A − B ) (cid:0) A + B cos(2 θ ) (cid:1) = 2 AA + B cos(2 θ ) − ( A − B ) (cid:0) A + B cos(2 θ ) (cid:1) . Substituting Equation (12) into the above equation yields( ˙ a ) + 1 = 2 Aa − ( A − B ) a . We have derived Equation (9), the Friedmann equation representing the expansion of the universe when it wasinfinitesimally small and infinitesimally young, from parametric equations for a and t given by Equations (12) and(16). So, we haven proven that in the beginning conditions of the positively curved universe, satisfied by the equation( ˙ a ) + 1 = Aa − ( A − B ) a , scale factor a and time t can be represented by parametric equations in terms of θ , as a = A + B cos(2 θ ) and t = √ A + B · E ( θ, ψ ), where ψ = BA + B . These parametric equations constitute a novelsolution to the Friedmann equation for the torsion-dominated early universe. V. RADIATION-DOMINATED UNIVERSE
Within a small amount of time into the initial expansion of the universe, the scale factor a surges massively upward,causing the torsion factor to rapidly approach 0. The torsion-dominated Friedmann equation from Equation (9) givesway to the original relativistic, radiation-dominated universe Friedmann equation from Equation (3), when the torsionfactor ( A − B ) a equals 0. This occurs when the numerator equals 0, when positive constants A and B are equal toeach other. Then, Equation (12) turns into a = A (1 + cos(2 θ )) . With the substitution of the cosine-function double-angle-identity, cos(2 θ ) = 2 cos ( θ ) −
1, we obtain a = 2 A cos ( θ ) , (17) a = √ A cos( θ ) , (18) dadθ = −√ A sin( θ ) , (19) (cid:0) dadθ (cid:1) = 2 A sin ( θ ) . (20)Substituting the result A = B into Equation (14) yields t = Z (cid:0)p A (1 + cos(2 θ )) (cid:1) dθ. Again, substituting the cosine-function double-angle-identity, cos(2 θ ) = 2 cos ( θ ) −
1, yields t = Z (cid:0)p A cos ( θ ) (cid:1) dθ = Z (cid:0) √ A cos( θ ) (cid:1) dθ = √ A sin( θ ) , (21) dtdθ = √ A cos( θ ) , (22) (cid:0) dtdθ (cid:1) = 2 A cos ( θ ) . (23)Plugging in the above equations for (cid:0) dtdθ (cid:1) and (cid:0) dadθ (cid:1) into the expanded form of ( ˙ a ) + 1 yields( ˙ a ) + 1 = (cid:0) dadθ (cid:1) (cid:0) dtdθ (cid:1) + 1 = 2 A sin ( θ )2 A cos ( θ ) + 1 = tan ( θ ) + 1 . Using the Pythagorean Identity tan ( θ ) + 1 = sec ( θ ) then yields( ˙ a ) + 1 = sec ( θ ) = 1cos ( θ ) = 2 A A cos ( θ ) . Substituting Equation (17) then yields ( ˙ a ) + 1 = 2 Aa . We have derived Equation (3), the Friedmann equation representing the expansion of the relativistic, radiation-dominated universe, during the time period following significant impact of the torsion factor. Note how the equivalence a = dtdθ is still preserved in this new set of solutions for the radiation-dominated era. VI. DARK ENERGY UNIVERSE
One further analysis can be put forth for completion of all possible analyses of the Friedmann equations, that beingthe current, dark-energy-dominated universe. This exists in the matter-dominated time period, so Equation (4) willbe altered to add in the dark energy factor. This factor is directly proportional to a , giving the equation( ˙ a ) + 1 = 2 Ca + Da (24)Further into the evolution of the universe, as in our current time period and the future, the scale factor is an extremelylarge number. This means that Ca approaches 0, and makes little contribution. The function converges to( ˙ a ) = Da , ˙ a = √ Da,a = e t √ D . The equation converges to the exponential function of increasing scale factor that is reflected in our currentcosmological observations, and the universe increases in size indefinitely.An important analysis resulting from this model is that the value of the constant D must be greater thansome threshold value, which depends on C , in order for the universe to start acceleration. VII. SUMMARY
We have developed novel parametric solutions for scale factor and time for the torsion-dominated era of the universe,near the Big Bounce, given positive closed curvature. The Friedmann equation used for this specific analysis is( ˙ a ) + 1 = 2 Aa − ( A − B ) a , with parametric solutions as a = A + B cos(2 θ ) ,t = √ A + B ∗ E ( θ, ψ ) , where ψ = BA + B .This solution can be used with the condition A = B to solve for the parametric equations of the relativistic,radiation-dominated era of the universe, at which point the torsion factor added to the Friedmann equation zeroesout. Moreover, the inspiration for such parametric solutions comes from the basis of cycloidal/trochoidal behaviorof a positively curved universe, given the parametric equations of the matter-dominated, non-relativistic era of theuniverse being equal to those of a cycloidal function.Important conclusions that can be made are that there may have been a ”bounce” period near the beginningof the universe, at which point the universe had never formed a singularity, but rather ”bounced” in a state extremelyclose to infinite density. This model suggests that the universe never had a Big Bang, or a definite beginning oftime for that matter, and the beginning of time cannot be traced back to the most recent absolute minimum sizeof the universe, due to the ”bounce” period described here. Rather, the Big Bounce acted as an inflection pointbetween two successive stages of evolution. Also, the models demonstrate that the coefficients of the Friedmannequation forms become irrelevant at extremely large values of scale factor and expansion rate in a universe governedby dark energy, making the fate of such closed universes to inevitably expand forever. An important result isthat the closed universe model follows cycloidal/trochoidal behavior, as is represented with the scenarios whereparametric equations involving correlated trigonometric functions were found as solutions to the Friedmann equations.An important finding that comes up in these analyses is that if the Big Bounce had indeed occurred, thenthe previous evolutionary stage of the universe, prior to the most recent Big Bounce, must not have been able toenter the dark energy phase of evolution and increase in size exponentially and indefinitely, thus reentering bounceperiods. Since we can rule out the possibility that dark energy can never overtake ordinary matter’s gravitationalpull because of scientific evidence against it, the model suggests that the influences of dark energy must increasefor each successive developmental stage of the universe. This means that the effects of dark energy increase oversuccessive evolutionary courses of the universe, each separated by a bounce period. In this present stage, the universemust have reached, for the first time and possibly only time, an evolutionary stage preceded by a bounce period butovertaken by dark energy and expanding indefinitely forever.These analyses remain critical for further research in cosmological evolution and the closed universe model.Overall, the elegance of the Friedmann equations in this setting are highly demonstrated by these analyses. Acknowledgments
I thank Professor Nikodem Poplawski of the University of New Haven for supervising and guiding me with my researchproject, and for his continued support in this work. I also thank Professor Joseph Kolibal for helping me initiate myresearch and providing me with helpful advice throughout my project. [1] J. Rich,
Fundamentals of Cosmology (Springer, 2001).[2] E. A. Lord,
Tensors, Relativity and Cosmology (McGraw-Hill, 1976).[3] L. D. Landau and E. M. Lifshitz,
The Classical Theory of Fields (Pergamon, 1975).[4] N. Pop lawski, Astrophys. J. 832, 96 (2016).[5] T. W. B. Kibble, J. Math. Phys. , 212 (1961); D. W. Sciama, Rev. Mod. Phys. , 463 (1964).[6] F. W. Hehl and J. D. McCrea, Found. Phys. , 267 (1986).[7] W. Kopczy´nski, Phys. Lett. A , 219 (1972); A. Trautman, Nature (Phys. Sci.) , 7 (1973).[8] F. W. Hehl, P. von der Heyde, and G. D. Kerlick, Phys. Rev. D , 1066 (1974); F. W. Hehl, P. von der Heyde, G. D.Kerlick, and J. M. Nester, Rev. Mod. Phys. , 393 (1976).[9] B. Kuchowicz, Gen. Relativ. Gravit. , 511 (1978).[10] M. Gasperini, Phys. Rev. Lett. , 2873 (1986).[11] N. J. Pop lawski, Phys. Lett. B , 181 (2010); N. Pop lawski, Phys. Rev. D , 107502 (2012).[12] S. Desai and N. J. Pop lawski, Phys. Lett. B , 183 (2016).[13] R. K. Pathria, Nature , 298 (1972); V. P. Frolov, M. A. Markov, and V. F. Mukhanov, Phys. Lett. B ; D. A. Eassonand R. H. Brandenberger, J. High Energ. Phys. , 024 (2001); N. J. Pop lawski, Phys. Lett. B687