Arxiv:2008.02136V1 [Gr-Qc] 4 Aug 2020 ∗ O H Pntno Ne Ed Oterpeetto of Spin) with Representation ﬂuid the (Ideal ﬂuid to Spin Leads [11]

Home , Black hole cosmology

arXiv:2008.02136v1 [gr-qc] 4 Aug 2020 ∗ o h pntno nE ed oterpeetto of spin) with representation fluid the (ideal fluid to spin leads [11]. a as EC matter in fermionic law tensor the conservation spin the spin the of the for expansion in multipole quadratic The are tensor. that terms matter additional con- of tensor acquires Levi-Civita energy–momentum symmetric can the which the EC in with nection, Consequently, GR as the matter. rewritten to fermionic proportional be be of to tensor [6– out tensor turns spin tensor torsion torsion the for The to action by respect 8]. the with determined varying matter is from and gravity torsion obtained con- equations Instead, not field is [10]. the [9], zero connection to an- affine strained the the is of which part tensor, torsion tisymmetric the con- that This requires interaction. sistency Dirac spin-orbit parti- total the the with allows Dirac consistent the that a be equation rela- for must of spacetime general law curved momentum in in conservation angular cle spin) as The and scalar, (orbital [6–8]. Ricci (GR) and the tivity Sciama to field by gravitational proportional the expanded for is theory, density spacetime Lagrangian this equips the In Einstein– Kibble, which classical gravity torsion. the of with theory by (EC) provided Cartan is singularities ing avoid however, hole must, black universe singularity. A new gravitational a a [4]. to universe bridge baby a a creating to paradox: goes Genesis, information information Hole hole an black Black the or solve we Cosmology naturally which Hole could hypothesis, Black in- This as the to horizon. of refer the because at hole black redshift the finite outside could observed dynamics be subsequent not Con- and of [3]. formation par- side bridge its the Einstein–Rosen sequently, other to an connected the through hole, on universe black ent the exist of would horizon universe event the a the universe another Such in in might bounce existing [2–5]. hole a universe black from parent closed a universe of The baby interior a as [1]. originated closed have be may universe [email protected] Introduction. h ipetadms aua ehns o prevent- for mechanism natural most and simplest The rvttoa olpeo udwt oso noaunivers a into torsion with fluid a of collapse Gravitational ahccelre nsz hntepeiu yl,utli re it until originated cycle, therefore bounces. have previous T might several the universe have then Our matter. may size indefinitely. horizon in of event larger amounts the cycle enormous d each of produce side production other Particle and the pa inflation shear. Quantum of over bounce. dominate period Einste nonsingular to gravitation the a torsion that and with show helps metric it We Tolman replaces source. the and a use larity as We fluid spin hole. relativistic black a into sion ecnie rvttoa olpeo peial symmet spherically a of collapse gravitational consider We eetosrain ugs htour that suggest observations Recent eateto ahmtc n hsc,Uiest fNwHa New of University Physics, and Mathematics of Department 0 otnPs od etHvn T056 USA 06516, CT Haven, West Road, Post Boston 300 ioe Pop lawski Nikodem fasi udaegvnby given are pressure fluid [18]. and spin density a al. of et energy Hehl effective by the discovered a Consequently, was gives tensor which energy–momentum value, the nonzero for in averaging terms even spin macroscopic occur the because can of singularity spins a oriented of randomly Kopczyński avoidance and was The [13], which Trautman aligned, [14]). [12], are Hehl fermions by of discovered the spins a when by 15–17] in [6, described singularity universe cosmological Friedmann–Lemaˆıtre–Robertson–Walker isotropic a metric (FLRW) and of homogeneous formation the vent α oriae hti ohsnhoosadcmvn [16]. comoving of and system synchronous a pressure both in of is described absence that ini- be the coordinates can is In collapse that hole. a black such fluid gradients, a spin into rest a at of tially sphere symmetric spherically theory field quantum in correc- [27]. diagrams) radiative Feynman of (loop divergence tions com- ultraviolet may and is it the Furthermore, [25] eliminate tensor [26]. universe acceleration cosmic spin the present the the in the explain asymmetry also if nonsin- may matter-antimatter Torsion A occurs [24]. cos- antisymmetric also that the pletely [23]. of [5] bounce radiation observations inflation gular background Planck exponential microwave the of mic with period consistent is finite a [22]. universe erate zero closed is new, energy total a whose as black 21] expands 5, a nonsingu- then [4, a in and reaches singularity bounce matter collapsing a lar the of Instead, the formation that hole. violate gravity the repulsive (1) prevent as in act may and coupling condition spin-torsion correc- energy strong negative the densities, the from density, high tions nuclear extremely EC than At and greater GR. neglected much be to can reduces torsion effectively of effects the densities, where n pressure, and oso a eeaegaiainlrplinadpre- and repulsion gravitational generate can Torsion nti ril,w osdrgaiainlclas fa of collapse gravitational consider we article, this In gen- can bounce a after production particle Quantum = κ rmabakhole. black a from ce h omlgclsz n expands and size cosmological the aches ( ǫ lrplino oso rvnsasingu- a prevents torsion of repulsion al ~ ∗ i peeo udwt pnadtor- and spin with fluid a of sphere ric rn xaso a eeaeafinite a generate can expansion uring c and ) tcepouto uigcontraction during production rticle 2 uhauies soclaoy with oscillatory, is universe a Such / 2[,1–0 with 18–20] [5, 32 nCra edeutoswt a with equations field in–Cartan p n ǫ ersligcoe nvreon universe closed resulting he ˜ f = r h hroyai nrydensity energy thermodynamic the are stenme est ffrin,and fermions, of density number the is ǫ − αn f 2 , ven, nabakhole black a in e p ˜ = κ p 8 = − αn πG/c f 2 , 4 tlower At . (1) 2

0 We use the Tolman metric [28] and the EC field equa- a spin fluid, Tµν = (˜ǫ +p ˜)uµuν pg˜ µν , are: T0 =ǫ ˜, 1 2 3 − tions with a relativistic spin fluid as a source. The non- T1 = T2 = T3 = p˜. The Einstein field equations µ µ − relativistic spin fluid was considered in [29]. We demon- Gν = κTν in this frame of reference are: strate that, after an event horizon forms and if there is G0 = κǫ,˜ G1 = G2 = G3 = κp,˜ G1 =0. (4) no shear, gravitational repulsion of torsion prevents a sin- 0 1 2 3 − 0 gularity and replaces it with a nonsingular bounce. The The covariant conservation of the energy–momentum resulting universe on the other side of the event hori- tensor gives zon is closed and oscillatory with an infinite number of bounces and cycles. Without torsion, a singularity would ′ 2ǫ˜˙ ′ 2˜p be reached and the metric would be described by the in- λ˙ + 2µ ˙ = , ν = , (5) ǫ˜+˜p ǫ˜+˜p terior Schwarzschild solution, which is equivalent to the − − Kantowski–Sachs metric describing an anisotropic uni- where the constants of integration depend on the allowed verse with topology R S2 [30]. Thanks to torsion, the transformations τ τ ′(τ) and R R′(R). × 3 universe in a black hole becomes closed with topology S If the pressure is→ homogeneous (no→ pressure gradients), (3-sphere). then p′ = 0 and p = p(τ). In this case, the second If we include quantum particle production in changing equation in (5) gives ν′ = 0. Therefore, ν = ν(τ) and a ′ gravitational fields [31], then two effects appear. Dur- transformation τ τ (τ) can bring ν to zero and g00 = ing contraction, particle production with torsion act to- eν to 1. The system→ of coordinates becomes synchronous gether to reverse gravitational attraction generated by [16]. Defining r(τ, R)= eµ/2 turns (2) into shear. During expansion, this production can generate a finite period of inflation and produce enormous ds2 = c2dτ 2 eλ(τ,R)dR2 r2(τ, R)(dθ2 +sin2θ dφ2). (6) amounts of matter. Consequently, each cycle is larger − − and longer then the previous cycle [5, 32]. The number The Einstein equations (3) reduce to of bounces and cycles is finite because the universe even- −λ e ′′ ′ ′ ′ 1 tually reaches a size at which the cosmological constant κǫ˜ = (2rr + r 2 rr λ )+ (rr˙λ˙ +r ˙2 + 1), becomes dominant and expands indefinitely. − r2 − r2 1 − ′ Gravitational collapse of a homogeneous sphere. For κp˜ = ( e λr 2 +2rr¨ +r ˙2 + 1), a spherically symmetric gravitational field in spacetime − r2 − −λ filled with an ideal fluid, the geometry is given by the e ′′ ′ ′ r˙λ˙ 1 2¨r 2κp˜ = (2r r λ )+ + λ¨ + λ˙ 2 + , Tolman metric [16, 28]: − − r − r 2 r ′ ′ 2 ν(τ,R) 2 2 λ(τ,R) 2 µ(τ,R) 2 2 2 2r ˙ λr˙ =0. (7) ds = e c dτ e dR e (dθ +sin θ dφ ), − − − (2) Integrating the last equation in (7) gives where ν, λ and µ are functions of a time coordinate τ and a radial coordinate R. We can still apply coordi- r′2 ′ ′ eλ = , (8) nate transformations τ τ (τ) and R R (R) without 1+ f(R) changing the form of the→ metric (2).→ The components of the Einstein tensor corresponding to (2) that do not where f is a function of R satisfying a condition 1+f > 0 vanish identically are [16, 28]: [16]. Substituting (8) into the second equation in (7) 2 2 ′2 ′ ′ −ν 2 gives 2rr¨ +r ˙ f = κpr˜ , which is integrated to 0 −λ ′′ 3µ µ λ e ˙ µ˙ − − G0 = e µ + + λµ˙ + − 4 − 2 2 2 2 F (R) κ 2 − r˙ = f(R)+ pr˜ dr, (9) +e µ, r − r Z −λ ′2 2 1 e µ ′ ′ −ν µ˙ ν˙ 3µ ˙ where F is a positive function of R. Substituting (8) into G1 = + µ ν + e µ¨ + − 2 2 − 2 4 the third equation in (7) does not give a new relation. −µ +e , Substituting (8) into the first equation in (7) and using −ν (9) gives 2 3 e ˙ ˙ ¨ ˙ 2 2 G2 = G3 = (λν˙ +µ ˙ ν˙ λµ˙ 2λ λ 2¨µ µ˙ ) − 4 − − − − − ′ −λ F (R) e ′′ ′ ′′ ′ ′ ′ ′ ′ ′ ′ κ(˜ǫ +p ˜)= . (10) (2ν + ν 2 +2µ + µ 2 µ λ ν λ + µ ν ), r2r′ − 4 − − −λ Combining (9) and (10) gives 1 e ′ ′ ˙ ′ ′ G0 = (2µ ˙ +µµ ˙ λµ µν˙ ), (3) 2 − − R κ ′ r˙2 = f(R)+ ǫr˜ 2r dR. (11) where a dot denotes differentiation with respect to cτ and r Z a prime denotes differentiation with respect to R. 0 In the comoving frame of reference, the spatial com- Every particle in a collapsing fluid sphere is represented ponents of the four-velocity uµ vanish. Accordingly, the by a radial coordinate R that ranges from 0 (at the cen- nonzero components of the energy–momentum tensor for ter of the sphere) to R0 (at the surface of the sphere). 3

If the mass of the sphere is M, then the Schwarzschild which has a form of the first law of thermodynamics 2 radius rg = 2GM/c of the black hole that forms from for the energy density and pressure (1) [5]. If we as- the sphere is equal to [16] sume that the spin fluid is composed by an ultrarela- 4 tivistic matter in kinetic equilibrium, then ǫ = h⋆T , R0 ′ 3 r = κ ǫr˜ 2r dR. (12) p = ǫ/3, and nf = hnfT , where T is the temperature g 2 4 ~ 3 Z0 of the fluid, h⋆ = (π /30)(gb + (7/8)gf)kB/( c) , and 2 3 3 hnf = (ζ(3)/π )(3/4)gfk /(~c) [5, 20]. For standard- Spinless dustlike sphere. Before considering gravita- B ′ model particles, g = 29 and g = 90. Since p = 0, the tional collapse of a sphere composed of a spin fluid, it is b f temperature does not depend on R: T = T (τ). Substi- instructive to consider spinless dust, for which the pres- tuting these relations into (20) gives sure vanishes and thusp ˜ = 0. Substituting (10) into (12) ′ gives r2r T 3 = g(R), (21) r = F (R ) F (0) = F (R ), (13) g 0 − 0 where g is a function of R. Putting this equation into which determines the value of R0. If f < 0, then (9) has (11) gives a solution R κ ′ F F r˙2 = f(R)+ (h T 4 αh2 T 6) r2r dR. (22) r = (1 + cos η), τ τ (R)= (η + sin η), r ⋆ − nf Z −2f − 0 2( f)3/2 0 − (14) Equations (21) and (22) give the function r(τ, R), which where η is a parameter and τ0(R) is a function of R [16, with (8) gives λ(τ, R). The integration of (22) also con- 28]. Choosing tains the initial value τ0(R). The metric (6) depends thus on three arbitrary functions: f(R), g(R), and τ (R). f(R)= sin2 R, F (R)=2a sin3 R, τ (R) = const. 0 0 0 We seek a solution of (21) and (22) as − (15) gives [16] f(R)= sin2 R, r(τ, R)= a(τ) sin R, (23) − r = a sin R(1 + cos η) τ τ = a (η + sin η), (16) 0 − 0 0 where a(τ) is a nonnegative function of τ. This choice is where a0 is a constant. Initially, at τ = τ0 and η = 0, analogous to a dust sphere: the first equation in (15), the the sphere is at rest:r ˙ = 0. Clearly, a singularity r = 0 first equation in (16), and (18). Accordingly, (21) gives is reached for all particles in a finite time. The value 3 3 2 of a0 and R0 can be determined from (13), the second a T sin R cos R = g(R), (24) equation in (15), and the first equation in (16) for η = 0: in which separation of the variables τ and R leads to 1/2 3 1/2 rg r0 sin R0 = , a0 = , (17) g(R) = const sin2 R cos R, a3T 3 = const. (25) r0 4rg · where r0 = r(0, R0) is the initial radius of the sphere. An Consequently, we find event horizon for the entire sphere forms when r(τ, R )= 0 ˙ r , that is, at cos(η/2) = sin R . T H g 0 aT = a0T0, + =0, (26) Substituting (15) and (16) into (8) gives eλ(τ,R) = T c a2(1 + cos η)2. If we define 0 where a0 = a(0), T0 = T (0), and H = ca/a˙ is the Hubble parameter. Substituting (23) into (22) gives a(τ)= a0(1 + cos η), (18)

2 κ 4 2 6 2 then the square of an infinitesimal interval in the interior a˙ +1= (h⋆T αhnfT )a . (27) of a collapsing dust (6) turns into [16] 3 − ds2 = c2dτ 2 a2(τ)dR2 a2(τ) sin2 R(dθ2 + sin2θ dφ2). Using (26) in (27) yields − − (19) 4 4 2 6 6 2 κ h⋆T0 a0 αhnfT0 a0 The initial value of a is equal to 2a0. This metric has a a˙ = 1+ 2 4 . (28) form of the closed FLRW metric and describes a part of − 3 a − a a closed universe with 0 R R0. Substituting (23) into (8) gives eλ(τ,R) = a2. Conse- Spin-fluid sphere. ≤ ≤ We now proceed to the main part of quently, the square of an infinitesimal interval in the in- the article and consider gravitational collapse of a sphere terior of a collapsing spin fluid (6) is also given by (19). composed of a spin fluid to demonstrate the formation µ/2 Substituting the initial values a(0) = 2a0 (17) and of a nonsingular universe. Substituting r = e and (8) a˙ (0) = 0 into (27), in which the second term on the right- into the first equation in (5) gives 2 3 4 hand side is negligible, gives Mc = (4π/3)r0h⋆T0 . This d ′ d ′ relation indicates the equivalence of mass and energy of a (˜ǫr2r )+˜p (r2r )=0, (20) dτ dτ fluid sphere with radius r0 and determines T0. An event 4 horizon for the entire sphere forms when r(τ, R0) = rg, nomenologically given by 1/2 which is equivalent to a = (rg r0) . Equation (28) has 1 d(√ gn ) βH4 two turning points,a ˙ = 0, if [20] − f = , (31) c√ g dt c4 3 ~4 4 − r0 3πG hnf 2 > l , (29) 6 4 2 3 Planck where g = a sin R sin θ is the determinant of the met- rg 8h⋆ ∼ ric tensor in− (19) and β is a nondimensional production which is satisfied for astrophysical systems that form rate [5]. With particle production, the second equation black holes. in (26) turns into Avoidance of singularity. Equation (28) can be solved analytically in terms of an elliptic integral of the second T˙ H βH3 = 3 3 1 . (32) kind [20], giving the function a(τ) and then r(τ, R) = T c 3c hnfT − a(τ) sin R. The value of a never reaches zero because as a decreases, the right-hand side of (21) becomes negative, Particle production changes the power law nf(a): contradicting the left-hand side. The change of the sign −(3+δ) 1/2 n a , (33) occurs when a< (rgr0) , that is, after the event horizon f ∼ forms. Consequently, all particles with R> 0 fall within where δ varies with τ. Putting this relation into (31) the event horizon but never reach r = 0 (the only particle gives at the center is the particle that is initially at the center, with R = 0). A singularity is therefore avoided. Nonzero δ aδa˙ 3. (34) values of a in (19) give finite values of T and therefore ∼− finite values of ǫ, p, and nf. During contraction,a ˙ < 0 and thus δ > 0. The term The resulting universe on the other side of the event 2 −6−2δ 2 −6 nf a grows faster than σ a and a singularity horizon has a closed geometry (constant positive curva- is avoided.∼ Particle production and∼ torsion act together ture). The quantity a(τ) is the scale factor of this uni- to reverse the effects of shear, generating a nonsingular verse. The universe is oscillatory: the value of a oscillates bounce. The dynamics of the nonsingular, relativistic between the two turning points. The value of R0 does not universe in a black hole is described by equations (27) change. A turning point at whicha> ¨ 0 is a bounce, and 3 1/2 and (32), with the initial conditions a(0) = (r0/rg) a turning point at whicha< ¨ 0 is a crunch. The universe anda ˙(0) = 0, that give the functions a(τ) and T (τ). The has therefore an infinite number of bounces and crunches, shear would enter the right-hand side of (27) as an addi- and each cycle is alike. tional positive term that is proportional to a−4. When 4 The Raychaudhuri equation for a congruence of the universe becomes nonrelativistic, the term h⋆T in geodesics without four-acceleration and rotation is (27) changes into a positive term that is proportional to 2 2 µ ν − dθ/ds = θ /3 2σ Rµν u u , where θ is the expan- a 1. The cosmological constant enters (27) as a positive − 2 − − sion scalar, σ is the shear scalar, and Rµν is the Ricci term that is proportional to a2. tensor [8]. For a spin fluid, the last term in this equation Particle production increases the maximum size of the is equal to κ(˜ǫ+3˜p)/2. Consequently, the necessary and − scale factor that is reached at a crunch. Consequently, sufficient condition for avoiding a singularity in a black the new cycle is larger and lasts longer then the previous 2 hole is κ(˜ǫ +3˜p)/2 > 2σ . For a relativistic spin fluid, cycle. According to (17), R0 is given by p = ǫ/3,− this condition is equivalent to 3 rg 2 2 sin R0 = , (35) 2καnf > 2σ + κǫ. (30) a(0) Without torsion, the left-hand side of (31) would be ab- where a(0) is the initial scale factor that is equal to the sent and this inequality could not be satisfied, resulting maximum scale factor in the first cycle. Since the maxi- in a singularity. Torsion therefore provides a necessary mum scale factor in the next cycle is larger, the value of condition for preventing a singularity. In the absence of sin R0 decreases. As cycles proceed, R0 approaches π. shear, this condition is also sufficient. Inflation and end of oscillations. During contraction, The presence of shear opposes the effects of torsion. H is negative and the temperature T increases. Dur- The shear scalar σ2 grows with decreasing a like a−6, ing expansion, if β is too big, then the right-hand side of 2 ∼ which is the same power law as that for nf . Therefore, if (32) could become positive. In this case, the temperature the initial shear term dominates over the initial torsion would grow with increasing a, which would lead to eter- term in (31), then it will dominate at later times during nal inflation [5]. Consequently, there is an upper limit contraction and a singularity will form. To avoid a sin- to the production rate: the maximum of the function 2 3 3 3 gularity if the shear is present, nf must grow faster than (βH )/(3c hnfT ) must be lesser than 1. −6 3 3 3 a . Consequently, fermions must be produced in a If (βH )/(3c hnfT ) in (32) increases after a bounce ∼black hole during contraction. to a value that is slightly lesser than 1, then T would Particle production. The production rate of particles become approximately constant. Accordingly, H would in a contracting or expanding universe [31] can be phe- be also nearly constant and the scale factor a would grow 5 exponentially, generating inflation. Since the energy den- a baby universe in a parent black hole existing in another sity would be also nearly constant, the universe would universe. This hypothesis is supported by the presented produce enormous amounts of matter and entropy. Such analysis of gravitational collapse of a spin fluid with tor- an expansion would last until the right-hand side of (32) sion and particle production. A more realistic scenario of drops below 1. Consequently, inflation would last a finite gravitational collapse should involve a fluid sphere that period of time. After this period, the effects of torsion is inhomogeneous and rotating. If the pressure in the weaken and the universe smoothly enters the radiation- sphere is not homogeneous, then the system of coordi- dominated expansion, followed by the matter-dominated nates cannot be comoving and synchronous [16, 34]. Con- expansion. sequently, ν and the temperature would depend on R and If the universe during expansion does not reach a criti- the equations of the collapse and the subsequent dynam- cal size at which the cosmological constant is significant, ics of the universe would be more complicated. If the then it recollapses to another bounce and starts a new sphere were rotating, then further complications would oscillation cycle [33]. The new cycle is larger and longer appear [35] and the angular momentum of the forming then the previous cycle [5, 32]. After a finite series of cy- Kerr black hole would be another parameter in addition cles, the universe reaches the critical size which prevents to the mass [36]. Nevertheless, the general character of the next contraction and and enters the cosmological- the effects of torsion and particle production in avoid- constant-dominated expansion, during which it expands ing a singularity and generating a bounce in a black hole indefinitely. The value of R0 asymptotically tends to π, would still be valid. which is the maximum value of R in a closed isotropic universe given by (19). The last bounce is the big bang. This work was funded by the University Research Final remarks. Our universe might have been born as Scholar program at the University of New Haven.

[1] W. Handley, arXiv:1908.09139; E. Di Valentino, A. Mel- Mathematical Society, 1927); E. Schrödinger, Space-time chiorri, and J. Silk, Nature Astron. 4, 196 (2020). Structure (Cambridge University Press, 1954); J. A. [2] I. D. Novikov, J. Exp. Theor. Phys. Lett. 3, 142 (1966); Schouten, Ricci-Calculus (Springer-Verlag, 1954). R. K. Pathria, Nature 240, 298 (1972); V. P. Frolov, M. [10] F. W. Hehl and J. D. McCrea, Found. Phys. 16, 267 A. Markov, and V. F. Mukhanov, Phys. Lett. B 216, (1986); N. Poplawski, arXiv:1304.0047. 272 (1989); Phys. Rev. D 41, 383 (1990); L. Smolin, [11] K. Nomura, T. Shirafuji, and K. Hayashi, Prog. Theor. Class. Quantum Grav. 9, 173 (1992); S. Hawking, Black Phys. 86, 1239 (1991). Holes and Baby Universes and other Essays (Bantam [12] F. W. Hehl, Abh. Braunschw. Wiss. Ges. 18, 98 (1966). Dell, 1993); W. M. Stuckey, Am. J. Phys. 62, 788 (1994); [13] A. Trautman, Bull. Acad. Polon. Sci., Serie Sci. Math. D. A. Easson and R. H. Brandenberger, J. High Energ. Astr. Phys. 20, 185 (1972); Symp. Math. 12, 139 (1973); Phys. 06, 024 (2001); J. Smoller and B. Temple, Proc. Nature Phys. Sci. 242, 7 (1973). Natl. Acad. Sci. USA 100, 11216 (2003). [14] W. Kopczyński, Phys. Lett. A 39, 219 (1972); W. [3] N. J. Poplawski, Phys. Lett. B 687, 110 (2010); N. Kopczyński, Phys. Lett. A 43, 63 (1973). Poplawski, arXiv:1910.10819; arXiv:1912.02173. [15] A. Friedmann, Z. Phys. A 10, 377 (1922); G. Lemaˆıtre, [4] N. J. Poplawski, Phys. Lett. B 694, 181 (2010); Phys. Ann. Soc. Sci. Bruxelles A 53, 51 (1933); H. P. Robert- Lett. B 701, 672 (2011). son, Astrophys. J. 82, 284 (1935); A. G. Walker, Proc. [5] N. Poplawski, Astrophys. J. 832, 96 (2016); Int. J. Mod. London Math. Soc. 42, 90 (1937). Phys. D 27, 1847020 (2018). [16] L. D. Landau and E. M. Lifshitz, The Classical Theory [6] E. A. Lord, Tensors, Relativity and Cosmology (McGraw- of Fields (Pergamon, 1975). Hill, 1976). [17] V. A. Fock, The Theory of Space, Time and Gravitation [7] D. W. Sciama, Proc. Camb. Phil. Soc. 54, 72 (1958); (Macmillan, 1964); P. A. M. Dirac, General Theory of T. W. B. Kibble, J. Math. Phys. 2, 212 (1961); D. W. Relativity (Wiley, 1975). Sciama, in Recent Developments in General Relativity, p. [18] F. W. Hehl, P. von der Heyde, and G. D. Kerlick, Phys. 415 (Pergamon, 1962); Rev. Mod. Phys. 36, 463 (1964); Rev. D 10, 1066 (1974). Rev. Mod. Phys. 36, 1103 (1964); F. W. Hehl and B. K. [19] I. S. Nurgaliev and W. N. Ponomariev, Phys. Lett. B Datta, J. Math. Phys. 12, 1334 (1971); F. W. Hehl, Gen. 130, 378 (1983). Relativ. Gravit. 4, 333 (1973); Gen. Relativ. Gravit. 5, [20] G. Unger and N. Poplawski, Astrophys. J. 870, 78 (2019). 491 (1974); F. W. Hehl, P. von der Heyde, G. D. Kerlick, [21] B. Kuchowicz, Gen. Relativ. Gravit. 9, 511 (1978); M. and J. M. Nester, Rev. Mod. Phys. 48, 393 (1976); V. de Gasperini, Phys. Rev. Lett. 56, 2873 (1986); Y. N. Sabbata and M. Gasperini, Introduction to Gravitation Obukhov and V. A. Korotky, Class. Quantum Grav. 4, (World Scientific, 1985); V. de Sabbata and C. Sivaram, 1633 (1987); N. J. Poplawski, Gen. Relativ. Gravit. 44, Spin and Torsion in Gravitation (World Scientific, 1994); 1007 (2012). N. J. Poplawski, Phys. Lett. B 690, 73 (2010); Phys. [22] N. J. Poplawski, Class. Quantum Grav. 31, 065005 Lett. B 727, 575 (2013). (2014); N. Poplawski, Mod. Phys. Lett. A 33, 1850236 [8] N. J. Poplawski, Classical Physics: Spacetime and Fields, (2018). arXiv:0911.0334. [23] S. Desai and N. J. Poplawski, Phys. Lett. B 755, 183 [9] L. P. Eisenhart, Non-Riemannian Geometry (American (2016). 6

[24] N. Poplawski, Phys. Rev. D 85, 107502 (2012); J. (1971); Phys. Rev. D 3, 2546 (1971); Y. B. Zeldovich and Magueijo, T. G. Zlosnik, and T. W. B. Kibble, Phys. Rev. A. A. Starobinskii, J. Exp. Theor. Phys. Lett. 26, 252 D 87, 063504 (2013); J. L. Cubero and N. J. Poplawski, (1977); V. A. Beilin, G. M. Vereshkov, Y. S. Grishkan, Class. Quantum Grav. 37, 025011 (2020). N. M. Ivanov, V. A. Nesterenko, and A. N. Poltavtsev, [25] N. J. Poplawski, Phys. Rev. D 83, 084033 (2011). J. Exp. Theor. Phys. 51, 1045 (1980). [26] N. Poplawski, Gen. Relativ. Gravit. 46, 1625 (2014). [32] J. D. Barrow and M. P. D¸abrowski, Mon. Not. Roy. As- [27] N. Poplawski, Found. Phys., in press (2020). tron. Soc. 275, 850 (1995); J. D. Barrow and C. Ganguly, [28] R. A. Tolman, Proc. Natl. Acad. Sci. USA 20, 169 (1934); Int. J. Mod. Phys. D 26, 1743016 (2017). J. R. Oppenheimer and H. Snyder, Phys. Rev. 56, 455 [33] H. Bondi, Cosmology (Cambridge University Press, (1939). 1960); J. D. North, The Measure of the Universe (Claren- [29] M. Hashemi, S. Jalalzadeh, and A. H. Ziaie, Eur. Phys. don Press, 1965). J. C 75, 53 (2015). [34] E. M. Lifshitz and I. M. Khalatnikov, J. Exp. Theor. [30] R. Kantowski and R. K. Sachs, J. Math. Phys. 7, 443 Phys. 12, 108 (1961). (1966); R. W. Brehme, Am. J. Phys. 45, 423 (1977); N. [35] A. G. Doroshkevich, Y. B. Zel’dovich, and I. D. Novikov, Poplawski, arXiv:2007.11556. J. Exp. Theor. Phys. 22, 122 (1966). [31] L. Parker, Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. [36] R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963). 183, 1057 (1969); Y. B. Zeldovich, J. Exp. Theor. Phys. Lett. 12, 307 (1970); L. Parker, Phys. Rev. D 3, 346