On Some Mathematical Connections between the Cyclic Universe, Inflationary Universe, p-Adic Inflation, p-Adic Cosmology and Various Sectors of Number Theory

Abstract

This paper is a review, a thesis, of some interesting results that have been obtained in various research concerning the “brane collisions in string and M-theory” (Cyclic Universe), p-adic inflation and p-adic cosmology. In Section 2, we have described some equations concerning cosmic evolution in a Cyclic Universe. In Section 3, we have described some equations concerning the cosmological perturbations in a Big Crunch/Big Bang space-time, the M-theory model of a Big Crunch/Big Bang transition and some equations concerning the solution of a braneworld Big Crunch/Big Bang Cosmology. In Section 4, we have described some equations concerning the generating ekpyrotic curvature perturbations before the Big Bang, some equations concerning the effective five-dimensional theory of the strongly coupled heterotic string as a gauged version of N=1 five-dimensional supergravity with four-dimensional boundaries, and some equations concerning the colliding branes and the origin of the Hot Big Bang. In Section 5, we have described some equations regarding the “null energy condition” violation concerning the inflationary models and some equations concerning the evolution to a smooth universe in an ekpyrotic contracting phase with w>1 . In Section 6, we have described some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-adic string theory. In Section 7, we have described various equations concerning the p-adic minisuperspace model, zeta strings, zeta nonlocal scalar fields and p-adic and adelic quantum cosmology. In Section 8, we have shown various and interesting mathematical connections between some equations concerning the p-adic inflation, the p-adic quantum cosmology, the zeta strings and the brane collisions in string and M-theory. Furthermore, in each section, we have shown the mathematical connections with various sectors of Number Theory, principally the Ramanujan’s modular equations, the Aurea Ratio and the Fibonacci’s numbers.

Share and Cite:

Nardelli, M. (2024) On Some Mathematical Connections between the Cyclic Universe, Inflationary Universe, p-Adic Inflation, p-Adic Cosmology and Various Sectors of Number Theory. Journal of Modern Physics, 15, 1869-1958. doi: 10.4236/jmp.2024.1511080.

1. Introduction

The aim of this paper is to describe the possible mathematical connections between various equations concerning some sectors of Number Theory, cosmic evolution in a Cyclic Universe, cosmological perturbations in a Big Bang/Big Crunch space-time concerning the M-theory model and that of the Brane World Cosmology, some equations concerning the perturbations of the Ekpyrotic curvature before the Big Bang, some equations concerning 5-dimensional Supergravity with 4-dimensional boundaries and some equations concerning the collision of branes and the origin of the Hot Big Bang. We also describe the mathematical connections between some equations concerning the “Zero Energy Condition” violation inherent in inflationary models, the approximate inflationary solutions rolling away from the unstable maximum of p-Adic String Theory and various equations concerning the p-Adic Minisuperspace model, zeta strings and p-Adic and Adelic Quantum Cosmology. As far as Number Theory is concerned, we will describe the possible connection between some equations concerning cosmic evolution in a Cyclic Universe, the Golden Ratio and some formulas connected to it and the Ramanujan Functions connected to the “modes” that correspond to the physical vibrations of a superstring and Bosonic Strings. Some possible connections between some equations concerning cosmological perturbations in a Big Crunch/Big Bang space-time and M-theory model of a Big Crunch/Big Bang transition and mathematical connection with the Aurea ratio and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring. Some possible connections between some equations concerning the solution of a braneworld Big Crunch/Big Bang Cosmology, the Golden Ratio and some formulas related to it and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring and Bosonic Strings. Some possible connections between some equations concerning the generating ekpyrotic curvature perturbations before the Big Bang and the Golden Ratio and some formulas related to it. Some possible connections between some equations concerning the effective five-dimensional theory of the strongly coupled heterotic string as a gauged version of N=1 five-dimensional supergravity with four-dimensional boundaries, the Golden Ratio and some formulas related to it and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring and Bosonic Strings. Some possible connections between some equations concerning the colliding Branes and the Origin of the Hot Big Bang, Aurea ratio and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring. Some possible connections between some equations concerning the evolution to a smooth universe in an ekpyrotic contracting phase with w>1 , the Golden Ratio and some formulas related to it and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring and Bosonic Strings. Some possible connections between some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-Adic String Theory, Aurea ratio and Ramanujan Functions related to the “modes” that correspond to the physical vibrations of a superstring. Some possible connections between an equation of the p-Adic Minisuperspace model and the Ramanujan Function related to the physical vibrations of the bosonic strings. We also show the mathematical connection between some values of the zeta strings and zeta nonlocal scalar fields, the Golden Ratio and some formulas related to it. Finally, we will show the possible mathematical connections between some equations related to the different sectors described above.

2. On Some Equations Concerning Cosmic Evolution in a Cyclic Universe

The action for a scalar field coupled to gravity and a set of fluids ρ i in a homogeneous, flat Universe, with line element d s 2 = a 2 ( τ )( N 2 d τ 2 +d x 2 ) is:

S= d 3 xdτ [ N 1 ( 3 a 2 + 1 2 a 2 ϕ 2 )N( ( aβ ) 4 Σ i ρ i + a 4 V( ϕ ) ) ] . (1)

We use τ to represent conformal time and primes to represent derivatives with respect to τ . N is the lapse function. The background solution for the scalar field is denoted ϕ( τ ) , and V( ϕ ) is the scalar potential.

The equations of motion for gravity, the matter and scalar field ϕ are straightforwardly derived by varying (1) with respect to a,N and ϕ , after which N may be set equal to unity. Expressed in terms of proper time t , the Einstein equations are:

H 2 = 8πG 3 ( 1 2 ϕ ˙ 2 +V+ β 4 ρ R + β 4 ρ M ) , (2)

a ¨ a = 8πG 3 ( ϕ ˙ 2 V+ β 4 ρ R + 1 2 β 4 ρ M ) , (3)

where a dot is a proper time derivative.

With regard the trajectory in the ( a 0 , a 1 ) -plane, the Friedmann constraint reads:

a 0 12 a 1 2 = 4 3 [ ( aβ ) 4 ρ+ 1 16 ( a 0 2 a 1 2 ) 2 V( ϕ 0 ) ] , (4)

Now, we solve the equations of motion immediately before and after the bounce.

Before the bounce there is a little radiation present since it has been exponentially diluted in the preceding quintessence-dominated accelerating phase. Furthermore, the potential V( ϕ ) becomes negligible as ϕ runs off to minus infinity. The Friedmann constraint reads ( a a ) 2 = 1 6 ϕ 2 , and the scalar field equation, ( a 2 ϕ )=0 , where primes denote conformal time derivatives. The general solution is:

ϕ= 3 2 ln( A H 5 ( in )τ ), a=A e ϕ/ 6 =A A H 5 ( in )τ ,

a 0 =A( λ+ λ 1 A H 5 ( in )τ ), a 1 =A( λ λ 1 A H 5 ( in )τ ) , (5)

where λ e ϕ / 6 . We choose τ=0 to be the time when a vanishes, so that τ<0 before collision. A is an integration constant which could be set to unity by rescaling space-time coordinates but it is convenient not to do so. The Hubble constants as defined in terms of the brane scale factors are a 0 / a 0 2 and a 1 / a 1 2 which at τ=0 take the values + λ 3 H 5 (in) and λ 3 H 5 (in) respectively. Re-expressing the scalar field as a function of proper time t= adτ , we obtain:

ϕ= 2 3 ln( 3 2 H 5 ( in )t ) . (6)

The integration constant H 5 ( in )<0 has a natural physical interpretation as a measure of the contraction rate of the extra-dimension. We remember that when the brane separation is small, one can use the usual formula for Kaluza-Klein theory,

d s 5 2 = e 2 3 ϕ d s 4 2 + e 2 2 3 ϕ d y 2 , (7)

where d s 4 2 is the four-dimensional line element, y is the fifth spatial coordinate which runs from zero to L , and L is a parameter with the dimension of length. Thence, we have that:

H 5 d L 5 Ld t 5 d( e 2 3 ϕ ) d t 5 = 2 3 ϕ ˙ e 3 2 ϕ , (8)

where L 5 L e 2 3 ϕ is the proper length of the extra dimension, L is a parameter with dimensions of length, and t 5 is the proper time in the five-dimensional metric,

d t 5 a e 1 6 ϕ dτ= e 1 6 ϕ dt , (9)

with t being FRW proper time. Notice that a shift ϕ can always be compensated for by a rescaling of L . As the extra dimension shrinks to zero, H 5 tends to a constant, H 5 ( in ) .

Immediately after the bounce, scalar kinetic energy dominates and H 5 remains nearly constant. The kinetic energy of the scalar field scales as a 6 and radiation scales as a 4 , so the former dominates at small a . It is convenient to re-scale a so that it is unity at scalar kinetic energy-radiation equality, t r , and denote the corresponding Hubble constant H r . The Friedmann constraint in Equation (4) then reads:

( a ) 2 = 1 2 H r 2 ( 1+ a 2 ) , (10)

and the solution is:

ϕ= 3 2 ln( 2 5 3 τ H 5 2 3 ( out ) H r 1 3 ( H r τ+ 2 3 2 ) ),a= 1 2 H r 2 τ 2 + 2 H r τ . (11)

The brane scale factors are:

a 0 a( λ 1 e ϕ/ 6 +λ e ϕ/ 6 )=A[ λ( 1+ H r τ 2 3 2 )+ λ 1 2 1 6 H r 1 3 H 5 2 3 ( out )τ ] ,

a 1 a( λ 1 e ϕ/ 6 +λ e ϕ/ 6 )=A[ λ( 1+ H r τ 2 3 2 ) 2 1 6 λ 1 H r 1 3 H 5 2 3 ( out )τ ] . (12)

Here, the constant A= 2 1 6 ( H r / H 5 ( out ) ) 1 3 has been defined so that we match a 0 and a 1 to the incoming solution given in (5). As for the incoming solution, we can compute the Hubble constants on the two branes after collision. They are ± λ 3 H 5 ( out )+ 2 5 3 λ 1 H r 2 3 H 5 1 3 on the positive and negative tension branes respectively. For H r < 2 5 2 λ 3 H 5 , the case of relatively little radiation production, immediately after collision a 0 is expanding but a 1 is contracting. Whereas for H r > 2 5 2 λ 3 H 5 , both brane scale factors expand after collision. If no scalar potential V( ϕ ) were present, the scalar field would continue to obey the solution (11), converging to:

ϕ C = 2 3 ln( 2 5 2 H 5 ( out ) H r ) . (13)

This value is actually larger than ϕ for H r < H 5 λ 3 2 5 2 , the case of weak production of radiation. However, the presence of the potential V( ϕ ) alters the expression (13) for the final resting value of the scalar field. As ϕ crosses the potential well travelling in the positive direction, H 5 is reduced to a renormalized value H ^ 5 ( out )< H 5 ( out ) , so that the final resting value of the scalar field can be smaller than ϕ . If this is the case, then a 1 never crosses zero, instead reversing to expansion shortly after radiation dominance. If radiation dominance occurs well after ϕ has crossed the potential well, Equation (13) provides a reasonable estimate for the final resting value, if we use the corrected value H ^ 5 ( out ) . The dependence of (13) is simply understood: while the Universe is kinetic energy dominated, a grows at t 1 3 and ϕ increases logarithmically with time. However, when the Universe becomes radiation dominated and a t 1 2 , Hubble damping increases and ϕ converges to the finite limit above.

With regard Equations (8)-(13), we note the following connections with number theory:

2 5/3 = 32 3 =3.174802104 Φ 16/7 + Φ 26/7 =3.171

2 3/2 = 8 =2.828427125 Φ 13/7 + Φ 14/7 =2.826

2 5/3 =0.314980262 Φ 20/7 + Φ 40/7 =0.3168

2 1/6 = 2 6 =1.122462048 Φ 1/7 + Φ 43/7 =1.1231

2 5/2 = 32 =5.656854249 Φ 25/7 + Φ 37/7 =5.6553

Note that, 32=8×4=24+8 , where 8 and 24 are the “modes” that correspond to the physical vibrations of a superstring and the physical vibrations of the bosonic strings.

Here, we have used the following expression: Φ n/7 , with Φ= 5 +1 2 =1.618033987 that is the Aurea ratio, n is a natural number and 7 are the compactified dimensions of the M-Theory.

Using the following potential:

V( ϕ )= V 0 ( 1 e cϕ )F( ϕ ) , (14)

we consider the motion of ϕ back and forth across the potential well. V may be accurately approximated by V 0 e cϕ . For this pure exponential potential, there is a simple scaling solution:

a( t )= | t | p ,V= V 0 e cϕ = p( 13p ) t 2 ,p= 2 c 2 , (15)

which is an expanding or contracting Universe solution according to whether t is positive or negative. A the end of the expanding phase of the cyclic scenario, there is a period of accelerated expansion which makes the Universe empty, homogeneous and flat, followed by ϕ rolling down the potential V( ϕ ) into the well. After ϕ has rolled sufficiently and the scale factor has begun to contract, the Universe accurately follows the above scaling solution down the well until ϕ encounters the potential minimum. Let us consider the behaviour of ϕ under small shifts in the contracting phase. In the background scalar field equation and the Friedmann equation, we set ϕ= ϕ B +δϕ and H= H B +δH , where ϕ B and H B are the background quantities given from (15). To linear order in δϕ , one obtains:

δ ϕ ¨ + 1+3p t δ ϕ ˙ 13p t 2 δϕ=0 , (16)

with two linearly independent solutions, δϕ t 1 and t 13p , where p1 . In the contracting phase, the former solution grows as t tends to zero. However, this solution is simply an infinitesimal shift in the time to the Big Crunch: δϕ ϕ ˙ . We next the incoming and outgoing collision velocity, which we have parameterized as H 5 ( in ) and H 5 ( out ) . Within the scaling solution (15), we can calculate the value of incoming velocity by treating the prefactor of the potential F( ϕ ) in Equation (14) as a Heaviside function which is unity for ϕ> ϕ min and zero for ϕ< ϕ min , where ϕ min  is the value of ϕ at the minimum of the potential. We compute the velocity of the field as it approaches ϕ min and use energy conservation at the jump in V to infer the velocity after ϕ min is crossed. In the scaling solution, the total energy as ϕ approaches ϕ min from the right is 1 2 ϕ ˙ 2 +V= 3 p 2 t 2 , and this must equal the total energy 1 2 ϕ ˙ 2 evaluated for ϕ just to the left of ϕ min . Hence, we find that ϕ ˙ = 6 p/t = 6p V min / ( 13p ) at the minimum and, according to Equation (8),

H 5 ( in ) 8 c | V min | 1 2 e 3 2 ϕ min 16 c 2 . (17)

Note that from Equation (17), we obtain:

[ H 5 ( in ) ] 2 8 c 2 | V min | ( e 3 2 ϕ min ) 2 16 c 2

where the number 8 is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

At the bounce, this solution is matched to an expanding solution with:

H 5 ( out )=( 1+χ ) H 5 ( in )>0 , (18)

where χ is a small parameter which arises because of the inelasticity of the collision. We shall simply assume a small positive χ is given, and follow the evolution forwards in time. Since χ is small, the outgoing solution is very nearly the time reverse of the incoming solution as ϕ starts back across the potential well after the bounce: the scaling solution is given in (15), but with t positive. We can treat χ as a perturbation and use the solution in Equation (16) discussed above, δϕ t 1 and t 13p . One can straightforwardly compute the perturbation in δ H 5 in this growing mode by matching at ϕ min as before. One finds δ H 5 = 12χ H 5 B / c 2 where H 5 B is the background value, at the minimum. Beyond this point, δ H 5 grows as t 6 /c e 3 2 ϕ for large c , whereas in the background scaling solution H 5 decays with ϕ as e ( 3 2 c/2 )ϕ . The departure occurs when the scalar field has attained the value,

ϕ Dep = ϕ min + 2 c ln c 2 12χ ,| V | ( 12χ c 2 ) 2 | V min | . (19)

As ϕ passes beyond ϕ Dep the kinetic energy overwhelms the negative potential and the field passes onto the plateau V 0 with H 5 nearly constant and equal to:

H ^ 5 ( out )χ ( c 2 12χ ) 6 c H 5 ( in ) , (20)

until the radiation, matter and vacuum energy become significant and H 5 is then damped away to zero. Note that we can rewrite Equation (20) as follow:

H ^ 5 ( out )χ ( c 2 12χ ) 6 c ×( 8 c V min 1 2 e 3 2 ϕ min 16 c 2 ) (21)

Also, this equation is related with the number 8, i.e. with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π 2 4 w ϕ w ( it x ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

and with the number 12 (12 = 24/2) that is related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π t 2 w dx e π 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

The time spent to the left of the potential well ( ϕ< ϕ min ) is essentially identical in the incoming and outgoing stages for χ1 , namely:

| t min | c 3 2| V min | . (22)

For the outgoing solution, when ϕ has left the scaling solution but before radiation domination, the definition Equation (8) may be integrated to give the time since the Big Bang at each value of ϕ ,

t( ϕ )= dϕ ϕ ˙ = 2 3 dϕ e 3 2 ϕ H ^ 5 ( ϕ ) 2 3 e 3 2 ϕ H ^ 5 ( out ) . (23)

Also, this equation can be rewritten as follow:

t( ϕ )= dϕ ϕ ˙ = 2 3 d ϕ e 3 2 ϕ H ^ 5 (ϕ) 2 3 e 3 2 ϕ / χ ( c 2 12χ ) 6 c ×( 8 c V min | 1 2 e 3 2 ϕ min  16 c 2 ) . (24)

The time in Equation (23) is a microphysical scale. The corresponding formula for the time before the Big Crunch is very different. In the scaling solution (15), one has for large c :

t( ϕ )= 2 | V min | e c ( ϕ ϕ min )/2 c = 6 e c ( ϕ ϕ min )/2 c 2 | t min | . (25)

The large exponential factor makes the time to the Big Crunch far longer than the time from the Big Bang, for each value of ϕ . This effect is due to the increase in H 5 after the bounce, which, in turn, is due to the positive value of χ . As the scalar field passes beyond the potential well, it runs onto the positive plateau V 0 . The value of H 5 ( out ) is nearly cancelled in the passage across the potential well, and is reduced to H ^ 5 given in Equation (20). Once radiation domination begins, the field quickly converges to the large t (Hubble-damped) limit of Equation (11), namely:

ϕ C = 2 3 ln( 2 5 2 H ^ 5 ( out )/ H r ) , (26)

where H r is the Hubble radius at kinetic-radiation equality. Also, Equation (26) can be rewritten as follow:

ϕ C = 2 3 ln[ 2 5 2 χ ( c 2 12χ ) 6 c ( 8 c V min 1 2 e 3 2 ϕ min 16 c 2 )/ H r ] . (27)

The dependence is obvious: the asymptotic value of ϕ depends on the ratio of H ^ 5 ( out ) to H r . Increasing H ^ 5 ( out ) pushes ϕ further, likewise lowering H r delays radiation domination allowing the logarithmic growth of ϕ in the kinetic energy dominated phase to continue for longer.

The solution of the scalar field equation is, after expanding Equation (11) for large τ , converting to proper time t= a( τ )dτ and matching,

ϕ ˙ 3 H r a 3 ( t ) a 3 0 t dt a 3 V ,ϕ , (28)

where as above, we define a( t ) to be unity at kinetic-radiation equal density. We have that ϕ may reach its maximal value ϕ max and turn around during the radiation, matter or quintessence dominated epoch. For example, ϕ max is reached in the radiation era, if, from Equation (28),

t max t m 10 4 ( t r t m ) 1 5 ( V V ,ϕ ( ϕ C ) ) 2 5 <1 , (29)

where t m is the time of matter domination.

For turn around in the matter era, we require:

3× 10 4 ( t r t m ) 1 6 ( V V ,ϕ ( ϕ C ) ) 1 3 30 . (30)

Finally, if the field runs to very large ϕ C , so that V ,ϕ / V( ϕ C ) c e c ϕ C is exponentially small, then ϕ only turns around in the quintessence-dominated era.

For our scenario to be viable, we require there to be a substantial epoch of vacuum energy domination (inflation) before the next Big Crunch. The number of e-foldings N e of inflation is given by usual slow-roll formula,

N e = dϕ V V ,ϕ e c ϕ c c 2 , (31)

for our model potential. For example, if we demand that the number of baryons per Hubble radius be diluted to below unity before the next contraction, which is certainly over-kill in guaranteeing that the cyclic solution is an attractor, we set e 3 N e 10 80 , or N e 60 . This is easily fulfilled if ϕ C is of order unity Planck units. Hence, Equation (31) can be rewritten as follow:

N e = dϕ V V ,ϕ e c ϕ c c 2 60 . (32)

With regard Equations (29)-(30) and (32), we have the following mathematical connections with the Aurea ratio:

Φ 49/7 + Φ 41/7 =29.03444185+0.059693843=29.0941

Φ 8/7 + Φ 14/7 =0.576974982+0.381966011=0.95894

Φ 10.33/7 =0.4914670835;arcsin( 0.4914670835 ) 180 π =29.437054

Φ 2/7 =0.8715438560;arccos( 0.8715438560 ) 180 π =29.361456

Φ 49/7 + Φ 50/7 =29.03444185+31.10060654=60.135048

Φ 10.33/7 =0.4914670835;arccos( 0.4914670835 ) 180 π =60.562946

Φ 2/7 =0.8715438560;arcsin( 0.8715438560 ) 180 π =60.638544

From the formulae given above, we can also calculate the maximal value ϕ C in the cyclic solution: for large c and for t r χ 1 t min , it is:

ϕ C ϕ min 2 3 ln( χ t r t min ) , (33)

where we used H r 1 t r , the beginning of the radiation-dominated epoch. From Equation (33), we obtain:

t r t min 1 χ ( c 2 N e | V min | V 0 ) 3 2 c 2 . (34)

This equation provides a lower bound on t r . The extreme case is to take | V min |1 . Then, using V 0 10 120 , c10 , N e 60 , we find t r 10 25 seconds. In this case the maximum temperature of the Universe is ≈1010 GeV. This is not very different to what one finds in simple inflationary models.

We have shown that a cyclic universe solution exists provided we are allowed to pass through the Einstein-frame singularity according to the matching conditions, Equations (17) and (18).

Specifically, we assumed that H 5 ( out )=( 1+χ ) H 5 ( in ) where χ is a non-negative constant, corresponding to branes whose relative speed after collision is greater than or equal to the relative speed before collision. Our argument showed that, for each χ0 , there is a unique value of H 5 ( out ) that is perfectly cyclic. Now we show that an increase in velocity is perfectly compatible with energy and momentum conservation in a collision between a positive and negative tension brane, provided a greater density of radiation is generated on the negative tension brane.

We shall assume that all other extra dimensions and moduli are fixed, and the bulk space-time between the branes settles down to a static state after the collision. We shall take the densities of radiation on the branes after collision as being given. By imposing Israel matching in both initial and final states, as well as conservation of total energy and momentum, we shall be able to completely fix the state of the outgoing branes and in particular the expansion rate of the extra dimension H 5 ( out ) , in terms of H 5 ( in ) . The initial state of empty branes with tensions T and T , and with corresponding velocities v + <0 and v >0 obeys:

T 1 v + 2 =T 1 v 2 ; E tot = T 1 v + 2 T 1 v 2 ; P tot = T v + 1 v + 2 T v 1 v 2 (35)

The first equation follows from Israel matching on the two branes as the approach, and equating the kinks in the brane scale factors. The second and third equations are the definitions of the total energy and momentum. The three equations in (35) imply that the incoming, empty state has v + = v , E tot =0 and that the total momentum is:

P tot = TL H 5 ( in ) 1 1 4 ( L H 5 ( in ) ) 2 <0 , (36)

where we identify v + v with the contraction speed of the fifth dimension, | L H 5 ( in ) | . For Equation (17), we can rewrite Equation (36) also as follow:

P tot =TL( 8 c | V min | 1 2 e 3 2 ϕ min  16 c 2 ) 1 1 1 4 [ L( 8 c | V min | 1 2 e 3 2 ϕ min 16 c 2 ) ] 2 <0 . (37)

The corresponding equations for the outgoing state are easily obtained, by replacing T with T+ ρ + T + for the positive tension brane, and T with T+ ρ T for the negative tension brane, assuming the densities of radiation produced at the collision on each brane, ρ + and ρ respectively, are given from a microphysical calculation, and are both positive.

Writing v ± ( out )=tanh( θ ± ) , where θ ± are the associated rapidities, one obtains two solutions:

sinh θ + = 1 2 T ( | P tot |+ | P tot | 1 ( T + 2 T 2 ) ); sinh θ =± 1 2 T + ( | P tot | | P tot | 1 ( T + 2 T 2 ) ) . (38)

where T + T+ ρ + , T T ρ with ρ + and ρ the densities of radiation on the positive and negative tension branes respectively, after collision. Both ρ + and ρ are assumed to be positive. In the first solution, with signs (−+), the velocities of the positive and negative tension branes are the same after the collision as they were before it. In the second, with signs (−−), the positive tension brane continues in the negative y direction but the negative tension brane is also moving in the negative y direction. The corresponding values for v ± (out) and V are:

v + ( out )= | P tot |+ | P tot | 1 ( T + 2 T 2 ) P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2 , v ( out )=± | P tot | | P tot | 1 ( T + 2 T 2 ) P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2

V= P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2 | P tot | ( T + 2 + T 2 )/ ( T + 2 T 2 ) + | P tot | 1 ( T + 2 T 2 ) , orV= P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2 | P tot |+ | P tot | 1 ( T + 2 + T 2 ) , (39)

where the first solution for V holds for the (−+) case, and the second for the (−) case. We are interested in the relative speed of the branes in the outgoing state, since that gives the expansion rate of the extra dimension, v + ( out )+ v ( out )=L H 5 ( out ) , compared to their relative speed 2 v + =L H 5 ( in ) in the incoming state. We find in the (−+) solution,

| H 5 ( out ) H 5 ( in ) |= v + ( out ) v ( out ) 2 v + = P tot 2 +4 T 2 P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2 , (40)

and in the (−−) solution,

| H 5 ( out ) H 5 ( in ) |= T + 2 T 2 P tot 2 P tot 2 +4 T 2 P tot 2 +2( T + 2 + T 2 )+ P tot 2 ( T + 2 T 2 ) 2 . (41)

with P tot given by (36) in both cases. We note that we can rewrite the relation mentioned above, i.e. v + ( out )+ v ( out )=L H 5 ( out ) also as follow:

v + ( out )+ v ( out )=L H 5 ( out )=L[ ( 1+χ )( 8 c | V min | 1 2 e 3 2 ϕ min 16 c 2 ) ] . (42)

At this point, we need to consider how the densities of radiation ρ + and ρ depend on the relative speed of approach of the branes. At very low speeds, | L H 5 ( in ) |1 , one expects the outer brane collision to be nearly adiabatic and an exponentially small amount of radiation to be produced. The (−+) solution has the speeds of both branes nearly equal before and after collision: we assume that it is this solution, rather than the (−) solution which is realised in this low velocity limit. As | L H 5 ( in ) | is increased, we expect ρ + and ρ to grow. Now, if we consider ρ + and ρ to be both P tot T , then the second term in the denominator dominates. If more radiation is produced on the negative tension brane, ρ > ρ + , then:

| H 5 ( out ) H 5 ( in ) |( 1+χ )[ 1+ ( ρ ρ + ) 2T ] (43)

and so χ is small and positive. This is the condition necessary to obtain cyclic behaviour. Conceivably, the brane tension can change from T to T =Tt at collision. Then, we obtain:

( 1+χ )[ 1+ ( ρ ρ + +2t ) 2T ] (44)

For the (−+) solution, we can straightforwardly determine an upper limit for | H 5 ( out )/ H 5 ( in ) |( 1+χ ) . Consider, for example, the case there the brane tension in unchanged at collision, t=0 . The expression in (40) gives | H 5 ( out )/ H 5 ( in ) | as a function of T + , T and P tot . It is greatest, at fixed T and P tot , when T + =T , its smallest value. For P tot 2 < T 2 , it is maximized for T 2 = T 2 P tot 2 , and equal to 1+ P tot 2 4 T 2 when equality holds. For P tot T 2 , it is maximized when T =0 , its smallest value, and P tot 2 =2 T 2 , when it is equal to 4 3 =1.154700538 . This is more than enough for us to obtain the small values of χ needed to make the cyclic scenario work. A reduction in brane tension at collisions t>0 further increases the maximal value of the ratio. To obtain cyclic behaviour, we need χ to be constant from bounce to bounce. That is, compared to the tension before collision, the fractional change in tension and the fractional production of radiation must be constant.

We note that for 4 3 =1.154700538 , we have the following mathematical connections with the Aurea ratio:

4 3 =1.154700538= Φ 1/7 + Φ 22/7 =0.933565132+0.220384833=1.1539499

4 3 =1.154700538= Φ 7/7 + Φ 9/7 =0,6180339887+0.5386437257=1.156677713

(For the various analyzed equations, see Ref. [1]).

3. On Some Equations Concerning Cosmological Perturbations in a Big Crunch/Big Bang Space-Time and M-Theory Model of a Big Crunch/Big Bang Transition

We consider a positive or negative tension brane with cosmological symmetry but which moves through the five-dimensional bulk. The motion through the warped bulk induces expansion or contraction of the scale factor on the brane. The scale factor on the brane obeys a “modified Friedmann” equation,

H ± 2 =± 1 3 M 5 3 L ρ ± + ρ ± 36 M 5 6 K b ± 2 + e b ± 4 , (45)

where ρ ± is the density (not including the tension) of matter or radiation confined to the brane, b ± is the brane scale factor, and H ± is the induced Hubble constant on the positive (negative) tension brane. Choosing conformal time on each brane, and neglecting the ρ 2 terms, Equations (45) become:

b + 2 =+ 1 3 M 5 3 L ρ + b + 4 K b + 2 +e, b 12 = 1 3 M 5 3 L ρ b 4 K b 2 +e . (46)

where prime denotes conformal time derivative. The corresponding acceleration equations for b + and b , from which e disappears, are derived by differentiating Equations (46) and using d( ρ b 4 )= b 3 ( ρ3P )db with P being the pressure of matter or radiation on the brines. We now show that these two equations can be derived from a single action provided we equate the conformal times on each brane. Consider the action,

S= dt N d 3 x[ 3 M 5 3 L( N 2 b + 2 K b + 2 ) ρ + b + 4 +3L( N 2 b 12 K b + 2 ) ρ b 4 ] , (47)

where N is a lapse function introduced to make the action time reparameterization invariant. Varying with respect to b ± and then setting N=1 gives the correct acceleration equations for b + and b following from (46). These equations are equivalent to (46) up to two integration constants.

We rewrite the action (47) in terms of a four-dimensional effective scale factor a and a scalar field ϕ , defined by:

b + =acosh( ϕ 6 ), b =asinh( ϕ 6 )

Clearly, a and ϕ transform as a scale factor and as a scalar field under rescalings of the spatial coordinates x . To interpret ϕ more physically, note that for static branes the bulk space-time is perfect Anti-de Sitter space with line element d Y 2 + e 2Y/L ( d t 2 +d x 2 ) . The separation between the branes is given by:

d=Lln( a + a )=Lln[ coth ( ϕ 6 ) ]

so d tends from zero to infinity as ϕ tends from minus infinity to zero. In terms of a and ϕ , the action (47) becomes:

S= dt d 3 x [ 3 M 5 3 L( a ˙ 2 K a 2 )+ 1 2 a 2 ϕ ˙ 2 ]+ S m , (48)

which is recognized as the action for Einstein gravity with line element a 2 ( t )( d t 2 + γ ij d x i d x j ) , γ ij being the canonical metric on H 3 , S 3 or E 3 with curvature K , and a minimally coupled scalar field ϕ . The matter action S m is conventional, except that the scale factor appearing is not the Einstein-frame scale factor but instead b + =acosh( ϕ/ 6 ) and b =asinh( ϕ/ 6 ) on the positive and negative tension branes respectively.

Now, we wish to make use of two very powerful principles. The first is the assertion that even in the absence of symmetry, the low energy modes of the five-dimensional theory should be describable with a four-dimensional effective action. The second is that since the original theory was coordinate invariant, the four-dimensional effective action must be coordinate invariant too. Since the five-dimensional theory is local and causal, it is reasonable to expect these properties in the four-dimensional theory. If furthermore the relation between the four-dimensional induced metrics on the branes and the four-dimensional fields is local, then covariance plus agreement with the above results forces the relation to be:

g μν + = ( cosh( ϕ/ 6 ) ) 2 g μν 4d , g μν = ( sinh( ϕ/ 6 ) ) 2 g μν 4d . (49)

When we couple matter to the brane metrics, these expressions should enter the action for matter confined to the positive and negative tension branes respectively. Likewise, we can from (48) and covariance immediately infer the effective action for the four-dimensional theory:

S= d 4 x g ( M 4 2 2 R 1 2 ( μ ϕ ) 2 )+ S m [ g ]+ S m + [ g + ] , (50)

where we have defined the effective four-dimensional Planck mass M 4 2 = ( 8π G 4 ) 1 = M 5 3 L . The two brane geometries are determined according to Formulae (49), and the background solution relevant post-collision is assumed to consist of two flat, parallel branes with radiation densities ρ ± . The corresponding four-dimensional effective theory has radiation density ρ r , and a massless scalar field with kinetic energy density ρ ϕ . The four-dimensional Friedmann equation in conformal time then reads:

a 2 = 1 3 ( ρ r a 4 + ρ ϕ a 4 )4 A 4 ( r 4 + A 4 a 2 ) , (51)

where we have defined the constants A 4 and r 4 , and used the fact that the massless scalar kinetic energy ρ ϕ a 6 . The solution to (51) and the massless scalar field equation ( a 2 ϕ )=0 is:

a 2 =4 A 4 τ( 1+ r 4 τ ) , ϕ= 3 2 ln( A 4 τ 1+ r 4 τ ) , (52)

From these solutions, we reconstruct the scale factors on the branes according to (49), obtaining:

b ± =1± A 4 τ+ r 4 τ , (53)

so we see that with the choice of normalization for the scale factor a made in (51), the brane scale factors are unity at collision. We may now directly compare the predictions (53) with the exact five-dimensional solution, equating the terms linear in τ to obtain:

A 4 =( 1/L )( 1+ L 2 ( r + r ) 12 )tanh( y 0 /2 ) , r 4 = L( r + + r ) 12tanh( y 0 /2 ) , (54)

where y 0 is the rapidity associated with the relative velocity of the branes at collision V=tanh( y 0 ) and r ± is the value of the radiation density ρ ± on each brane at collision. Thence, Equation (53) can be rewritten also:

b ± =1±( 1/L )( 1+ L 2 ( r + r ) 12 )tanh( y 0 /2 )τ+ L( r + + r ) 12tanh( y 0 /2 ) τ . (55)

Furthermore, we define the fractional density mismatch on the two branes as:

f= r + r r + + r , (56)

so that we have:

r + r = 12f r 4 L tanh( y 0 /2 ) . (57)

Now, we describe the perturbations of the brane-world system in terms of the four-dimensional effective theory. We shall now describe the scalar perturbations, in longitudinal gauge with a spatially flat background where the scale factor and the scalar field are given by (52). The perturbed line element is:

d s 2 = a 2 ( τ )( 1+2Φ )d τ 2 +( 12Ψ )d x 2 . (58)

Since there are no anisotropic stresses in the linearized theory, we have Φ=Ψ .

A complete set of perturbation equations consists of the radiation fluid equations, the scalar field equation of motion and the Einstein momentum constraint:

δ r = 4 3 ( k 2 v r 3 Φ ) v r = 1 4 δ r +Φ ( δϕ ) +2 ( δϕ ) = k 2 ( δϕ )+4 ϕ Φ

Φ +Φ= 2 3 a 2 ρ r v r + 1 2 ϕ ( δϕ ) , (59)

where primes denote τ derivatives, δ r is the fractional perturbation in the radiation density, v r is the scalar potential for its velocity, i.e. v r = v r , δϕ is the perturbation in the scalar field, and from (52), we have the background quantities:

a a = 1+2 r 4 τ 2τ( 1+ r 4 τ ) and 2 3 ϕ = 1 τ( 1+ r 4 τ )

We are interested in solving these equations in the long wavelength limit, | kτ |1 . Solving all the above equations for δlna , one finds:

δ i 1+ w i δ 1+w ,i=1,,N , (60)

for adiabatic perturbations. The components of the background energy density in the four-dimensional effective theory are scalar kinetic energy, with w ϕ =1 , and radiation, with w r = 1 3 . It follows that for adiabatic perturbations, at long wavelengths, we must have:

δ ϕ 3 2 δ r . (61)

In longitudinal gauge, the fractional energy density perturbation and the velocity potential perturbation in the scalar field (considered as a fluid with w=1 ) are given by:

δ ϕ =2( ( δϕ ) ϕ Φ ), v ϕ = δϕ ϕ . (62)

From Equations (59) above (and using ϕ a 2 ), it follows that:

( δ ϕ 3 2 δ r ) =2 k 2 ( v r δϕ ϕ ) . (63)

Maintaining the adiabaticity condition (61) up to order ( kτ ) 2 then requires that the fractional velocity perturbations for the scalar field and the radiation should be equal: v r δϕ/ ϕ . Expressing the radiation velocity in terms of δϕ , the momentum constraint then yields:

δϕ ( 1+ 2 3 ρ r ρ ϕ ) 1 ( 2( Φ +Φ ) ϕ ) , (64)

where ρ ϕ = 1 2 ϕ 2 a 2 .

The above equations may be used to determine the leading terms in an expansion in | kτ | of all the quantities of interest about the singularity. We shall choose to parameterize the expressions in terms of the parameters describing the comoving energy density perturbation, ε m = 2 3 2 k 2 Φ , which has the following series expansion about τ=0 :

ε m = ε 0 D( τ )+ ε 2 E( τ ) , (65)

where ε 0 and ε 2 are arbitrary constants, and

D( τ )=12 r 4 τ 1 2 k 2 τ 2 ln| kτ |+,E( τ )= τ 2 + (66)

For adiabatic perturbations, we obtain:

δ ϕ = ε 0 ( 9 4 k 2 τ 2 3 8 ln| kτ |+ 1 4 3 4 r 4 2 k 2 )+ ε 2 3 4 k 2 +O( τ,τln| kτ | ) , (67)

v ϕ = ε 0 ( 3 4 k 2 τ ( 1 r 4 τ ) )+O( τ,τln| kτ | ) ,

δ r = 2 3 δ ϕ +O( τ 2 , τ 2 ln| kτ | ) , v r = v ϕ +O( τ,τln| kτ | ) ,

Φ= ε 0 ( 3 8 k 2 τ 2 + 3 28 ln| kτ |+ 3×5 8 r 4 2 k 2 ) ε 2 3 8 k 2 +O( τ,τln| kτ | ) , (68)

δϕ 6 = ε 0 ( 3 8 k 2 τ 2 ( 12 r 4 τ )+ 1 28 ln| kτ |+ 1 8 + 13 8 r 4 2 k 2 ) ε 2 1 8 k 2 +O( τ,τln| kτ | ) , (69)

ζ 4,M = 1 2 k 2 ε 2 + ε 0 ( 1 8 k 2 ( k 2 +16 r 4 2 )+ 1 4 ln| kτ | )+O( τ,τln| kτ | ) , (70)

where ζ 4,M is the curvature perturbation on comoving slices introduced by Mukhanov.

With regard the (69), we note that is possible the following mathematical connection with the Aurea ratio:

1 6 =0.40824829= Φ 13/7 = ( 5 +1 2 ) 13/7 =0.40914770.409

Furthermore, in Equations (68)-(69), there is the number 8, that is related to the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π t 2 w dx e π 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] ,

and that 2, 3, 5, 8 and 13 are Fibonacci’s numbers.

Now, we consider the propagation of metric perturbations through a collision of tensionless branes where the background space-time is precisely C / Z 2 × R 3 . The form that we take for the five-dimensional cosmological background metric is:

d s 2 = n 2 ( t,y )( d t 2 + t 2 d y 2 )+ b 2 ( t,y ) δ ij d x i d x j , (71)

and we write the most general scalar metric perturbation about this as:

d s 2 = n 2 ( t,y )( ( 1+2Φ )d t 2 2Wdtdy+ t 2 ( 12Γ )d y 2 2 i αd x i dt + 2 t 2 i βdyd x i )+ b 2 ( t,y )( ( 12Ψ ) δ ij 2 i j χ )d x i d x j . (72)

For perturbations on C × R 3 it is straightforward to find a gauge in which the metric takes the form:

d s 2 =( 1+ 4 3 k 2 χ )( d t 2 + t 2 d y 2 )+( ( 1 2 3 k 2 χ ) δ ij +2 k i k j χ )d x i d x j , (73)

and χ satisfies a massless scalar equation of motion on M C × R 3 . To be precise, the gauge is:

α=β=0,Γ=ΦΨ k 2 χ,Φ= 2 3 k 2 χ,Ψ= 1 3 k 2 χ,W=0 . (74)

Notice that the non-zero variables can all be related to χ according to:

( Γ,Φ,Ψ )=( 2 3 ,+ 2 3 ,+ 1 3 ) k 2 χ . (75)

We shall, henceforth, refer to these as the “Milne ratio conditions”. Furthermore, imposing the Z 2 symmetry, we obtain Neumann boundary conditions on χ ,

χ ( y ± )=0 , (76)

where y ± =± y 0 /2 are the location of the two Z 2 fixed points. In the model space-time, the lowest energy mode for χ is y-independent and has the asymptotic form:

χ( t,y )=Q+Pln| kt | , (77)

with Q and P being arbitrary constants. We have the following relations:

Q out = Q in +2( γln2 ) P in , P out = P in . (78)

These relations are sufficient to determine the metric fluctuations after the bounce. We are only interested in the long-wavelength part of the spectrum, and, for the cases of interest, P is suppressed by k 2 compared to Q . As a result, we obtain the approximate matching rule:

Q out = Q in , P out = P in . (79)

The key conditions (74) through (76) are satisfied precisely for all time in a compactified Milne mod Z 2 background.

Now, we wish to use the four-dimensional effective (moduli) theory to infer the boundary data for the five-dimensional bulk perturbations. In any four-dimensional gauge, the four-dimensional metric perturbation h μν and scalar field perturbation δϕ determined the induced metric perturbations on the branes via Formulae (49):

h μν ± = h μν +2 ( ln Ω ± ) ,ϕ δϕ g μν , (80)

where Ω + =cosh( ϕ/ 6 ) and Ω =sinh( ϕ/ 6 ) and the metric perturbations are fractional, i.e. δ g μν = a 2 h μν , δ g μν ± = b ± 2 h μν ± . This formula is particularly easy to use in five-dimensional longitudinal gauge. This gauge may always be chosen, and it is completely gauge fixed. In this gauge, the five-dimensional metric takes the form:

d s 2 = n 2 ( t,y )( ( 1+2 Φ L )d t 2 2 W L dtdy+ t 2 ( 12 Γ L )d y 2 ) + b 2 ( t,y )( ( 12 Ψ L ) δ ij )d x i d x j . (81)

In the absence of anisotropic stresses, the brane trajectories are unperturbed in this gauge. An immediate consequence is that the four-dimensional longitudinal gauge scalar perturbation variables Φ ± and Ψ ± describing perturbations of the induced geometry on each brane:

d s ± 2 = b ± 2 ( τ ± )( ( 1+2 Φ ± )d τ ± 2 +( 12 Ψ ± )d x 2 ) , (82)

are precisely the boundary values of the five-dimensional longitudinal gauge perturbations Φ ± Φ L ( y ± ) and Ψ ± Ψ L ( y ± ) . Using (80) and (82), we find for the induced perturbations:

Φ + = Φ 4 + 1 6 tanh( ϕ 6 )δϕ , Ψ + = Φ 4 1 6 tanh( ϕ 6 )δϕ ,

Φ = Φ 4 + 1 6 coth( ϕ 6 )δϕ , Ψ = Φ 4 1 6 coth( ϕ 6 )δϕ . (83)

The brane conformal times may be expressed in terms of t by integrating,

τ ± = 0 t dt q( t, y ± ) , (84)

where qb/n . So, for example, the boundary value of the bulk metric perturbation Φ L on the positive tension brane is given explicitly by:

Φ L ( t, y + )= Φ 4 ( q ( t, y + ) 1 dt ) + 1 6 tanh( ϕ ( q ( t, y + ) 1 dt )/ 6 )δϕ( q ( t, y + ) 1 dt ), (85)

where y + is the location of the positive tension brane.

Also, in these Equations (83)-(85), we have the following connection with the Aurea ratio:

1 6 =0.40824829= Φ 13/7 = ( 5 +1 2 ) 13/7 =0.40914770.409 .

Now, using (67) and the following equation:

Ψ( t )= ζ 4,M ( t )+ ε 0 tanh( y 0 /2 ) 32 k 2 L 2 cosh 2 ( y 0 /2 ) [ 18( y 0 sinh( y 0 ) ) L 2 ( r + r )( 3 y 0 +sinh( y 0 ) ) ]+O( ρ ± 2 L 2 ,t,tln| kt | ) (86)

to find Q and P before and after the bounce for all components of the metric perturbations and matching according to the rule given in Equation (79) results in ζ 4,M inheriting two separate scale-invariant long wavelength contributions in the post-singularity state. The first occurs as a direct consequence of the sign change in (79), and is independent of the amount of radiation generated at the singularity. The second is proportional to the difference in the densities of the radiation on the two branes. At leading order in velocities, we have:

Δ ζ 4,M = 3 64 ε 0 k 2 L 2 ( V in 4 + V out 4 ) ( r + r ) ε 0 V out 2 32 k 2 +O( r ± V 3 , V 5 L 2 , ρ ± 2 L 2 ) , (87)

where V in and V out are the relative velocities of the branes before and after collision. Note that since P ε 0 , matching P is in fact equivalent to matching ε 0 across the collision. In terms of four-dimensional parameters including r 4 given in (57) defining the abundance of the radiation and the fractional density mismatch f defined in (56), we find again at leading order in velocities:

Δ ζ 4,M = 3 64 ε 0 k 2 L 2 ( V in 4 + V out 4 ) 3 ε 0 28 k 2 f r 4 V out 3 L . (88)

This is the final result, relevant to tracking perturbations across the singularity in the ekpyrotic and cyclic models.

The result for the long wavelength curvature perturbation amplitude in the four-dimensional effective theory, propagated into the hot Big Bang after the brane collision is:

ζ M = 9 ε 0 16 k 2 L 2 tanh( θ/2 ) cosh 2 ( θ/2 ) ( θsinhθ ) 3 ε 0 V coll 4 64 k 2 L 2 (89)

where θ is the rapidity corresponding to the relative speed V coll of the branes at collision, and the second formula assumes V coll is small. L is the bulk curvature scale, and ε 0 16 k 2 has a scale invariant power spectrum.

With regard Equations (88) and (89), we have the following mathematical connections with the Aurea ratio and the Fibonacci’s numbers: 64=48+16=24+24+8+8 ; 64= 8 2 =34+21+8+1 that are Fibonacci’s numbers; 24=3×8 ; 64=( 3×8 )+( 3×8 )+8+8 , with 3 and 8 that are Fibonacci’s numbers; 8 Φ 53/7 306.342224 .

Furthermore, in Equation (88), 2, 3 and 8 are Fibonacci’s numbers and 8 is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

The d+1 -dimensional space-time we consider is a direct product of d1 -dimensional Euclidean space, R d1 , and a two-dimensional time-dependent space-time known as compactified Milne space-time, or M C . The line element for M C × R d1 is thus:

d s 2 =d t 2 + t 2 d θ 2 +d x 2 ,0θ θ 0 ,<t< , (90)

where x are Euclidean coordinates on R d1 , θ parameterizes the compact dimension and t is the time. The compact dimension may either be a circle, in which case we identify θ with θ+ θ 0 , or a Z 2 orbifold in which case we identify θ with θ+2 θ 0 and further identify θ with 2 θ 0 θ . The fixed points θ=0 and θ= θ 0 are then interpreted as tensionless Z 2 -branes approaching at rapidity θ 0 , colliding at t=0 to re-emerge with the same relative rapidity. The orbifold reduction is the case of prime interest in the ekpyrotic/cyclic models, originally motivated by the construction of heterotic M-theory from eleven-dimensional supergravity. In these models, the boundary branes possess nonzero tension. However, the tension is a subdominant effect near t=0 and the brane collision is locally well-modelled by M C × R d1 .

Now consider a string loop of radius R in M-theory frame. Its mass M is 2πR times the effective string tension μ 2 L , where L is the size of the extra dimension. The effective Einstein-frame gravitational coupling is given by κ d 2 = κ d+1 2 /L . The gravitational potential produced by such a loop in d spacetime dimensions is:

Φ= κ d 2 M ( d2 ) A d2 R d3 (91)

where A D is the area of the unit D -sphere, A D = 2π D+1 2 / Γ( ( D+1 )/2 ) . Specializing to the case of interest, namely 2-branes in eleven-dimensional M-theory, the tension μ 2 is related to the eleven-dimensional gravitational coupling by a quantization condition relating to the four-form flux, reading:

μ 2 3 = 2 π 2 / ( n κ 11 2 ) (92)

with n an integer. Equations (91) and (92) then imply that the typical gravitational potential around a string loop is:

Φ= 105 64π μ 2 2 R 6 n ( μ 2 2 R 6 n ) 1 θ 0 2 /n (93)

up to numerical factors.

With regard (93), we note that are possible the following mathematical connections with the Aurea ratio and the Fibonacci’s numbers:

105=89+16=55+34+8+8;105=21×5 64= 8 2 =34+21+8+1=48+16=24+24+8+8=( 3×8 )+( 3×8 )+8+8

and 3, 5, 8, 21, 34, 55 and 89 are Fibonacci’s numbers. Furthermore, we have that:

105 64 =1.6406251.640; Φ 5/7 + Φ 21/7 = ( 5 +1 2 ) 5/7 + ( 5 +1 2 ) 21/7 =1.646251.646

Thence, the gravitational potential on the scale of the loops is of order θ 0 2 and therefore is consistently small for small collision rapidity. Since the mean separation of the loops when they are produced is of order their size R , this potential Φ is the typical gravitational potential throughout space. Multiplying the tt component of the background metric (90) by 1+2Φ and redefining t , we conclude that the outgoing metric has an expansion rapidity of order θ 0 ( 1+C θ 0 2 ) with C a constant of order unity. We conclude that for small θ 0 the gravitational back-reaction due to string loop productions is small.

On Some Equations Concerning the Solution of a Braneworld Big Crunch/Big Bang Cosmology

We shall employ a coordinate system in which the five-dimensional line element for the background takes the form:

d s 2 = n 2 ( t,y )( d t 2 + t 2 d y 2 )+ b 2 ( t,y )d x 2 , (94)

where y parameterizes the fifth dimension and x i ,i=1,2,3 , the three non-compact dimensions. Cosmological isotropy excludes dtd x i or dyd x i terms, and homogeneity ensures n and b are independent of x . The t,y part of the background metric may then be taken to be conformally flat.

We find it simplest to work in coordinates in which the brane locations are fixed but the bulk evolves. The bulk metric is therefore given by (94), with the brane locations fixed at y=± y 0 for all time t . The five-dimensional solution then has to satisfy both the Einstein equations and the Israel matching conditions on the branes. The bulk Einstein equations read G μ v =Λ δ μ ν , where the bulk cosmological constant is Λ=6/ L 2 . Evaluating the linear combinations G 0 0 + G 5 5 and G 0 0 + G 5 5 ( 3/2 ) G i i , we find:

β ,ττ β ,yy + β ,τ 2 β ,y 2 +12 e 2v =0 , v ,ττ v ,yy + 1 3 ( β ,y 2 β ,τ 2 )2 e 2v =0 , (95)

where t/L = e τ ,β3lnb and vln( nt/L ) . The Israel matching conditions on the branes read:

b ,y b = n ,y n = nt L (96)

where all quantities are to be evaluated at the brane locations y=± y 0 .

Now, we express the metric as a series of Dirichlet or Neumann polynomials in y 0 and y , bounded at order n by a constant times y 0 n , such that the series satisfies the Israel matching conditions exactly at every order in y 0 . To implement this, we first change variables from b and n to those obeying Neumann boundary conditions. From (96), b/n is Neumann. Likewise, if we define N( t,y ) by:

nt= 1 N( t,y )y , (97)

then one can easily check that N( t,y ) is also Neumann on the branes. Since N and b/n obey Neumann boundary conditions on the branes, we can expand both in a power series:

, b/n = q 0 ( t )+ n=3 q n ( t ) P n ( y ) , (98)

where P n ( y ) are polynomials:

P n ( y )= y n n n2 y n2 y 0 2 , n=3,4, (99)

satisfying Neumann boundary conditions and each bounded by | P n ( y ) |< 2 y 0 n / ( n2 ) , for the relevant range of y . Note that the time-dependent coefficients in this ansatz may also be expanded as a power series in y 0 . By construction, our ansatz satisfies the Israel matching conditions exactly at each order in the expansion. Substituting the series ansatze (98) into the background Einstein Equations (95), we may determine the solution order by order in the rapidity each order in y 0 , one generically obtains a number of linearly independent algebraic equations, and at most one ordinary differential equation in t .

The first few terms of the solution are:

N 0 = 1 t 1 2 t y 0 2 + 1 24 t( 89 t 2 ) y 0 4 + , N 3 = 1 6 +( 5 72 2 t 2 ) y 0 2 +

and

q 0 =1 3 2 t 2 y 0 2 +( t 2 7 8 t 4 ) y 0 4 + , q 3 =2 t 3 y 0 2 + (100)

With regard Equations (100), we note that are possible the following mathematical connections with the Aurea ratio and the Fibonacci’s numbers:

5 72 =0.06944 Φ 39/7 = ( 5 +1 2 ) 39/7 =0.06849207

7 8 =0.875 Φ 2/7 = ( 5 +1 2 ) 2/7 =0.87154

Furthermore, 72 (=24 × 3) and 8 are connected with the “modes” that correspond to the physical vibrations of the bosonic strings and to the physical vibrations of a superstring by the following Ramanujan functions:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log [ 10+11 2 4 + 10+7 2 4 ,

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

To calculate the affine distance between the branes along a spacelike geodesic we must solve the geodesic equations in the bulk. Let us first consider the situation in Birkhoff-frame coordinates for which the bulk metric is static and the branes are moving. The Birkhoff-frame metric takes the form:

d s 2 =d Y 2 N 2 ( Y )d T 2 + A 2 ( Y )d x 2 , (101)

where for Schwarzschild-AdS with a horizon at Y=0 ,

A 2 ( Y )= cosh( 2Y/L ) cosh( 2 Y 0 /L ) , N 2 ( Y )= cosh( 2 Y 0 /L ) cosh( 2Y/L ) [ sinh( 2Y/L ) sinh( 2 Y 0 /L ) ] 2 . (102)

At T=0 , the Y -coordinate of the branes is represented by the parameter Y 0 ; their subsequent trajectories Y ± ( T ) can then be determined by integrating the Israel matching conditions, which read tanh( 2 Y ± /L )=± 1 V ± 2 , where V ± = ( d Y ± / dT )/ N( Y ± ) are the proper speeds of the positive- and negative-tension branes respectively. From this, it further follows that Y 0 is related to the rapidity y 0 of the collision by tanh y 0 =sech( 2 Y 0 /L ) .

For the purpose of measuring the distance between the branes, a natural choice is to use spacelike geodesics that are orthogonal to the four translational Killing vectors of the static bulk, corresponding to shifts in x and T . Taking the x and T coordinates to be fixed along the geodesic then, we find that Y ,λ is constant for an affine parameter λ along the geodesic. To make the connection to our original brane-static coordinate system, recall that the metric function b 2 ( t,y )= A 2 ( Y ) , and thus:

Y ,λ 2 = ( b b ,t t ,λ +b b ,y y ,λ ) 2 b 4 θ 2 = n 2 ( t ,λ 2 + t 2 y ,λ 2 ) , (103)

where we have introduced the constant θ=tanh y 0 =V/c . Adopting y now as the affine parameter, we have:

0=( b ,t 2 b 2 + n 2 ( b 4 θ 2 ) ) t ,y 2 +2 b ,t b ,y b 2 t ,y +( b ,y 2 b 2 n 2 t 2 ( b 4 θ 2 ) ) , (104)

where t is to be regarded now as a function of y . We can solve this equation order by order in y 0 using the series ansatz:

t( y )= n=0 c n y n , (105)

where the constants c n are themselves series in y 0 . Using the series solution for the background geometry given in Equations (100), and imposing the boundary condition that t( y 0 )= t 0 , we obtain:

c 0 = t 0 + t 0 y 0 2 2 2 t 0 2 y 0 3 + ( t 0 +36 t 0 3 ) y 0 4 24 t 0 2 ( 1+5 t 0 2 ) y 0 5 +( t 0 720 + 17 t 0 3 4 +4 t 0 5 ) y 0 6 t 0 2 ( 13+250 t 0 2 +795 t 0 4 ) y 0 7 60 +O( y 0 8 ) (106)

c 1 =2 t 0 2 y 0 2 +( 5 t 0 2 3 +5 t 0 4 ) y 0 4 8 t 0 3 y 0 5 +( 91 t 0 2 180 + 23 t 0 4 6 + 53 t 0 6 4 ) y 0 6 +O( y 0 7 ) (107)

c 2 = t 0 2 t 0 ( 1+6 t 0 2 ) y 0 2 4 + t 0 2 y 0 3 ( t 0 48 2 t 0 3 +4 t 0 5 ) y 0 4 + ( t 0 2 +23 t 0 4 ) y 0 5 2 +O( y 0 6 ) (108)

c 3 = 5 t 0 2 y 0 2 3 t 0 2 ( 25+201 t 0 2 ) y 0 4 18 +O( y 0 5 ) (109)

c 4 = 5 t 0 24 +( 5 t 0 48 + 7 t 0 3 4 ) y 0 2 5 t 0 2 y 0 3 12 +O( y 0 4 ) (110)

c 5 = 61 t 0 2 y 0 2 60 +O( y 0 3 ) (111)

c 6 = 61 t 0 720 +O( y 0 2 ) (112)

c 7 =0+O( y 0 ). (113)

Substituting t 0 = x 0 / y 0 and y=ω y 0 , we find x( ω )= x 0 / y 0 +O( y 0 ) , i.e. to lowest order in y 0 , the geodesics are trajectories of constant time lying solely along the ω direction. Hence, in this limit, the affine and metric separation of the branes, defined with the following equation:

d m =L 1 1 x 4 1ω x 4 dω +O( y 0 2 )=Lln( 1+ x 4 1 x 4 )+O( y 0 2 ) , (114)

must necessarily agree. To check this, the affine distance between the branes is given by:

d a I = y 0 y 0 n t 2 t 2 dy =2 t 0 y 0 + ( t 0 +5 t 0 3 ) y 0 3 3 4 t 0 2 y 0 4 + ( t 0 10 t 0 3 +159 t 0 5 ) y 0 5 60 2( t 0 2 +30 t 0 4 ) y 0 6 3 + ( t 0 +31115 t 0 3 5523 t 0 5 +12795 t 0 7 ) y 0 7 2520 +O( y 0 8 ), (115)

which to lowest order in y 0 reduces to:

d a L =2 x 0 + 5 x 0 3 3 + 53 x 0 5 20 + 853 x 0 7 168 +O( x 0 8 )+O( y 0 2 ), (116)

in agreement with the series expansion of (114).

We obtain also the following equation:

δ d m L = 1 1 n x 4 ξ 40 ( x 4 ) b 2 e 1 2 x 4 2 dω = 2 x 4 ξ 40 ( x 4 ) ( 1 x 4 2 ) 2 = 1 ( x 4 2 1 ) [ 8 3 k ˜ x 4 ( A 0 J 0 ( k ˜ x 4 )+ B 0 Y 0 ( k ˜ x 4 ) )4( A 0 J 1 ( k ˜ x 4 )+ B 0 Y 1 ( k ˜ x 4 ) ) ]. (117)

To evaluate the perturbation δ d a in the affine distance between the branes, consider:

δ g μν x ˙ μ x ˙ v dλ = 1 2 dλ g ρσ x ˙ ρ x ˙ σ ( δ g μν x ˙ μ x ˙ ν + g μν,κ δ x κ x ˙ μ x ˙ ν +2 g μν x ˙ μ δ x ˙ ν ) =[ x ˙ v δ x v g ρσ x ˙ ρ x ˙ σ ]+ 1 2 δ g μv x ˙ μ x ˙ ν g ρσ x ˙ ρ x ˙ σ dλ (118)

where dots indicate differentiation with respect to the affine parameter λ , and in going to the second line we have integrated by parts and made use of the background

geodesic equation x ¨ σ = 1 2 g μν,σ x ˙ μ x ˙ ν and the constraint g μν x ˙ μ x ˙ ν =1 . If the

endpoints of the geodesics on the branes are unperturbed, this expression is further simplified by the vanishing of the surface term. Converting to coordinates where t 0 = x 0 / y 0 and y=ω y 0 , to lowest order in y 0 the unperturbed geodesics lie purely in the ω direction, and so the perturbed affine distance is identical to the following perturbed metric distance:

δ d m L = y 0 y 0 nt Γ L dy = 1 1 n x 4 ξ 4 b 2 e 1 2 x 4 2 dω . (119)

Explicitly, we find:

δ d a L = 2( B+A t 0 2 ) y 0 t 0 [ B( 4+3 t 0 2 ) 12 t 0 + A( t 0 +9 t 0 3 ) 3 ] y 0 3 +( 4B+4A t 0 2 ) y 0 4 [ B( 2+2169 t 0 2 +135 t 0 4 )+2A t 0 2 ( 1+1110 t 0 2 +375 t 0 4 ) 120 t 0 ] y 0 5 +[ 4A t 0 2 ( 1+42 t 0 2 )B( 4+57 t 0 2 ) 6 ] y 0 6 [ B( 4+88885 t 0 2 +952866 t 0 4 +28875 t 0 6 )+4A t 0 2 ( 1152481 t 0 2 +293517 t 0 4 +36015 t 0 6 ) 10080 t 0 ] y 0 7 +O( y 0 8 ), (120)

which, substituting t 0 = x 0 / y 0 and dropping terms of O( y 0 2 ) , reduces to:

δ d a L = 2 B ˜ x 0 2A x 0 B ˜ 4 x 0 3A x 0 3 9 8 B ˜ x 0 3 25 4 A x 0 5 275 96 B ˜ x 0 5 343 24 A x 0 7 +O( x 0 8 ), (121)

where B ˜ =B y 0 2 . Also, this expression is in accordance with the series expansion of (119). However, the perturbed affine and metric distance do not agree at O( y 0 2 ) .

With regard Equations (116) and (121), we have the following mathematical connections with Aurea ratio and with Fibonacci’s numbers:

5 3 =1.666 Φ 5/7 + Φ 20/7 = ( 5 +1 2 ) 5/7 + ( 5 +1 2 ) 20/7 =1.663054757

53 20 =2.65 Φ 12/7 + Φ 14/7 = ( 5 +1 2 ) 12/7 + ( 5 +1 2 ) 14/7 =2.663697

853 168 =5.077381 Φ 21/7 + Φ 16/7 + Φ 10/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 16/7 + ( 5 +1 2 ) 10/7 =5.0718

9 8 =1.125 Φ 7/7 + Φ 10/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 10/7 =1.12089

25 4 =6.25 Φ 28/7 Φ 7/7 = ( 5 +1 2 ) 28/7 ( 5 +1 2 ) 7/7 =6.236067977

275 96 =2.864583 Φ 14/7 + Φ 20/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 20/7 =2.870901

343 24 =14.2916 Φ 35/7 + Φ 17/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 17/7 =14.3078

Furthermore, we have that:

53=21+342 ; 853=2+8+233+610 ; where 2, 8, 21, 34, 233 and 610 are Fibonacci’s numbers;

9= 3 2 ; 25=34 3 2 ; 275=233+34+8 ; 343=233+34+8+13+55 ; where 3, 8, 13, 34, 55 and 233 are Fibonacci’s numbers.

(For the various analyzed equations, see Ref. [2]-[5]).

4. On Some Equations Concerning the Generating Ekpyrotic Curvature Perturbations before the Big Bang

With regard the ekpyrotic perturbations including gravity, we consider the action for N decoupled fields interacting only through gravity:

d 4 x g [ 1 2 R 1 2 i=1 N ( ϕ i ) 2 i=1 N V i ( ϕ i ) ] , (122)

where we have chosen units in which 8πG M Pl 2 =1 . In a flat Friedmann-Robertson-Walker background with line element d s 2 =d t 2 + a 2 ( t )d x 2 , the scalar field and Friedmann equations are given by:

ϕ ¨ i +3H ϕ ˙ i + V i, ϕ i =0 (123)

and

H 2 = 1 3 [ 1 2 i ϕ ˙ i 2 + i V i ( ϕ i ) ] , (124)

where H= a ˙ /a and V i, ϕ i = V i / ϕ i with no summation implied. Another useful relation is:

H ˙ = 1 2 i ϕ ˙ i 2 . (125)

If all the fields have negative exponential potentials V i ( ϕ i )= V i e c i ϕ i then as is well-known, the Einstein-scalar equations admit the scaling solution:

a= ( t ) p , ϕ i = 2 c i ln( A i t ), V i = 2 A i 2 c i 2 ,p= i 2 c i 2 . (126)

Thus, if c i 1 for all i , we have a very slowly contracting universe with p1 .

We focus on the entropy perturbation since this is a local, gauge-invariant quantity, and on the case of only two scalar fields. The entropy perturbation equation:

δ ¨ s+[ k 2 + ϕ ˙ 2 2 V , 1 ϕ 1 2 ϕ ˙ 1 ϕ ˙ 2 V , ϕ ϕ ϕ 2 + ϕ ˙ 1 2 V , ϕ 2 ϕ 2 ϕ ˙ 1 2 + ϕ ˙ 2 2 +3 ( ϕ ˙ 2 V , ϕ 1 ϕ ˙ 1 V , ϕ 2 ϕ ˙ 1 2 + ϕ ˙ 2 2 ) 2 ]δs=0 (127)

in flat spacetime is replaced by:

δ ¨ s+3H δ ˙ s+[ k 2 a 2 + ϕ ˙ 2 2 V , ϕ 1 ϕ 1 2 ϕ ˙ 1 ϕ ˙ 2 V , ϕ 1 ϕ 2 + ϕ ˙ 1 2 V , ϕ 2 ϕ 2 ϕ ˙ 1 2 + ϕ ˙ 2 2 +3 ( ϕ ˙ 2 V , ϕ 1 ϕ ˙ 1 V , ϕ 2 ϕ ˙ 1 2 + ϕ ˙ 2 2 ) 2 ]δs = 4 k 2 θ ˙ a 2 ϕ ˙ 1 2 + ϕ ˙ 2 2 Φ. (128)

For simplicity we will focus attention on straight line trajectories in scalar field space. Since θ ˙ =0 , the entropy perturbation is not sourced by the Newtonian potential Φ and we can solve the equations rather simply. We shall assume that the background solution obeys scaling symmetry so that ϕ ˙ 2 =γ ϕ ˙ 1 . Denoting τ derivatives with primes, and introducing the re-scaled entropy field:

δS=a( τ )δs , (129)

Equation (128) becomes:

δ S +( k 2 a a + a 2 V ,ϕϕ )δS=0 . (130)

The crucial term governing the spectrum of the perturbations is then:

τ 2 ( a a V ,ϕϕ a 2 ) . (131)

When this quantity is approximately 2, we will again get nearly scale-invariant perturbations. It is customary to define the quantity:

ε 3 2 ( 1+w ) ϕ ˙ 1 2 + ϕ ˙ 2 2 2 H 2 = ( 1+ γ 2 ) ϕ ˙ 2 2 H 2 . (132)

In the background scaling solution,

ε= c 2 2( 1+ γ 2 ) . (133)

We proceed by evaluating the quantity in (131) in an expansion in inverse powers of ε and its derivatives with respect to N , where N=ln( a/ a end ) , where a end is the value of a at the end of the ekpyrotic phase. Note that N decreases as the fields roll downhill and the contracting ekpyrotic phase proceeds. We obtain the first term in (131) by differentiating (125), obtaining:

a a =2 H 2 a 2 ( 1 1 2 ε ) . (134)

The second term in (131) is found by differentiating (132) twice with respect to time and using the background equations and the definition of N . We obtain:

a 2 V ,ϕϕ = a 2 H 2 ( 2 ε 2 6ε 5 2 ε ,N )+O( ε 0 ) . (135)

Finally, need to express ( a /a )=aH in terms of the conformal time τ . From (134), we obtain:

= 2 ( 1ε ) (136)

which integrates to:

1 = 0 τ dτ ( ε1 ) . (137)

Now, inserting 1= d( τ )/ dτ under the integral and using integration by parts we can re-write this as:

1 =ετ( 1 1 ε ( ετ ) 1 0 τ ε τdτ ) . (138)

Using the same procedure once more, the integral in this expression can be written as:

( ετ ) 1 0 τ ε τdτ = ε τ ε ( ετ ) 1 0 τ d dτ ( ε τ )τdτ . (139)

Now using the fact that ε = ε ,N , and that to leading order in 1/ε , can be replaced by its value in the scaling solution (with constant ε ), τ= ε 1 , we can re-write the second term on the righthand side as:

( ετ ) 1 0 τ d dτ ( ε τ )τdτ = ( ετ ) 1 0 τ dτ 1 ε ( ε ,N ε ) ,N , (140)

which shows that this term is of order 1/ ε 2 and can thus be neglected. Altogether, we obtain:

1 = 0 τ dτ ( ε1 )ετ( 1 1 ε ε ,N ε 2 ) . (141)

Using (134) and (135) with (141), we can calculate the crucial term entering the entropy perturbation equation,

τ 2 ( a a V ,ϕϕ a 2 )=2( 1 3 2ε + 3 4 ε ,N ε 2 ) . (142)

The deviation from scale-invariance in the spectral index of the entropy perturbation is then given by:

n s 1= 2 ε ε ,N ε 2 . (143)

The first term on the right-hand side is the gravitational contribution, which, being positive, tends to make the spectrum blue. The second term is the non-gravitational contribution, which tends to make the spectrum red.

Now defining to be the curvature perturbation on comoving spatial slices, for N scalar fields with general Kahler metric g ij ( ϕ ) on scalar field space, the linearized Einstein-scalar field equations lead to:

R ˙ = H H ˙ ( g ij D 2 ϕ i D t 2 s j k 2 a 2 Ψ ) , (144)

where the N1 entropy perturbations:

s i =δ ϕ i ϕ ˙ i g jk ( ϕ ) ϕ ˙ j δ ϕ k g lm ( ϕ ) ϕ ˙ l ϕ m (145)

are just the components of δ ϕ i orthogonal to the background trajectory, and the operator D 2 / D t 2 is just the geodesic operator on scalar field space. Things simplify because the scalar field space is flat, so the metric is g ij = δ ij , and D/ Dt reduces to an ordinary time derivative. Considering only two scalar fields, we have:

s 1 = ϕ ˙ 2 δs/ ϕ ˙ 1 2 + ϕ ˙ 2 2 , s 2 =+ ϕ ˙ 1 δs/ ϕ ˙ 1 2 + ϕ ˙ 2 2 . (146)

For a straight line trajectory in field space, the right-hand side of (145) vanishes even if the entropy perturbation is nonzero.

We assume that the scalar field bounce occurs after the ekpyrotic potentials are turned off, so that the universe is kinetic-dominated from the 4 d point-of-view. The scalar field trajectory is ϕ ˙ 2 = γ ˜ ϕ ˙ 1 , for t< t b , and ϕ ˙ 2 = γ ˜ ϕ ˙ 1 , for t> t b , with ϕ ˙ 1 constant and negative in the vicinity of the bounce. The bounce leads to a delta function on the right-hand side of (144),

D 2 ϕ 2 D t 2 =δ( t t b )2 ϕ ˙ 2 ( t b + ) , (147)

where t b is the time of the bounce of the negative-tension brane. As can be readily seen from (144), if the entropy perturbations already have acquired a scale-invariant spectrum by the time t b , then the bounce leads to their instantaneous conversion into curvature perturbations with precisely the same long wavelength spectrum. We can estimate the amplitude of the resulting curvature perturbation by integrating Equation (145) using (147). Since we have assumed the universe is kinetic-dominated at this time, H=1/ ( 3t ) . Since the entropy perturbation:

δs ( ϕ ˙ 1 δ ϕ 2 ϕ ˙ 2 δ ϕ 1 )/ ϕ ˙ 1 2 + ϕ ˙ 2 2 (148)

is canonically normalized, its spectrum is given by:

δ ϕ 2 = k 2 dk 4 π 2 1 k 3 t 2 (149)

up to non-scale invariant corrections. This expression only holds as long as the ekpyrotic behaviour is still underway: the ekpyrotic phase ends at a time t end approximately given by | V min |=2/ ( c 2 t end 2 ) . After t end , the entropy perturbation obeys δ ¨ s+ t 1 δ ˙ =0 , which has the solution δs=A+Bln( t ) .

Matching this solution to the growing mode solution t 1 in the ekpyrotic phase, one finds that by t b the entropy grows by an additional factor of 1+ln( t end / t b ) . Employing the Friedmann equation to relate ϕ ˙ 2 = γ ˜ ϕ ˙ 1 to H , putting everything together and restoring the Planck mass, we find for the variance of the spatial curvature perturbation in the scale-invariant case,

2 = c 2 | V min  | 3 π 2 M Pl 2 γ ˜ 2 ( 1+ γ ˜ 2 ) 2 ( 1+ln( t end / t b ) ) 2 dk k dk k Δ 2 ( k ) (150)

for the perfectly scale-invariant case. Notice that the results depend only logarithmically on t b : the main dependence is on the minimum value of the effective potential and the parameter c . Observations on the current Hubble horizon indicate Δ 2 ( k )2.2× 10 9 . Ignoring the logarithm in (150), this requires c | V min | 1 2 10 3 M Pl , or approximately the GUT scale. This is of course entirely consistent with the heterotic M-theory.

With regard Equations (150), we have the following mathematical connections with the Aurea ratio:

1 3 π 2 =0.03377 1 25 Φ 16/7 = 1 25 ( 5 +1 2 ) 16/7 =0.03329 . If we take M Pl =0.4340 , then we have that 1 3 π 2 M Pl 2 =0.179307948 Φ 25/7 = ( 5 +1 2 ) 25/7 =0.1793145665 . Furthermore, we obtain:

arcsin( 0.1793145665 ) 180 π +arccos( 0.1793145665 ) 180 π =10.3289+79.670290=55+34+1,

with 34 and 55 Fibonacci’s numbers.

If the entropic perturbations are suddenly converted to curvature perturbations, the curvature perturbations inherit the spectral tilt given in (143). We now begin by re-expressing Equation (143) in terms of N , the number of e-folds before the end of the ekpyrotic phase (where dN=( ε1 )N and ε1 ):

n s 1= 2 ε dlnε dN . (151)

This expression is identical to the case of the Newtonian potential perturbations, except that the first term has the opposite sign. In this expression, ε( N ) measures the equation of state during the ekpyrotic phase, which must decrease from a value much greater than unity to a value of order unity in the last N e-folds. If we estimate ε N α , then the spectral tilt is:

n s 1 2 N α α N (152)

Here, we see that the sign of the tilt is sensitive to α . For nearly exponential potentials ( α1 ) , the spectral tilt is n s 1+1/N 1.02 , slightly blue, because the first term dominates. However, there are well-motivated examples in which the equation of state does not decrease linearly with N . We have introduced α to parameterize these cases. If α>0.14 , the spectral tilt is red. For example, n s =0.97 for α2 . These examples represent the range that can be achieved for the entropically induced curvature perturbations in the simplest models, roughly 0.97< n s <1.02 .

For comparison, if we use the same estimating procedure for the Newtonian potential fluctuations in the cyclic model (assuming they converted to curvature fluctuations before the bounce through 5 d effects), we obtain 0.95< n s <0.97 . This range agrees with the estimate obtained by an independent analysis based on studying inflaton potentials directly.

With regard the values of n s , i.e. 0.95, 0.97 and 1.02, we have the following mathematical connections with the Aurea ratio:

0.950.9552046220= Φ 0.67/7 = ( 5 +1 2 ) 0.67/7 ;

0.970.9773457024= Φ 0.33/7 = ( 5 +1 2 ) 0.33/7 ;

1.021.0231794109= Φ 0.33/7 = ( 5 +1 2 ) 0.33/7 ;

A second way of analyzing the spectral tilt is to assume a form for the scalar field potential. Consider the case where the two fields have steep potentials that can be modelled as V( ϕ 1 )= V 0 e cdϕ and ϕ ˙ 2 =γ ϕ ˙ 1 . Then, Equation (143) becomes:

n s 1= 4( 1+ γ 2 ) c 2 M Pl 2 4 c ,ϕ c 2 , (153)

where we have used the fact that c( ϕ ) has the dimensions of inverse mass and restored the factors of Planck mass. The presence of M Pl clearly indicates that the first term on the right is a gravitational term. It is also the piece that makes a blue contribution to the spectral tilt. The second term is the non-gravitational term and agrees precisely with the following flat space-time result:

n s 1=4 c ,ϕ c 2 , (154)

although the agreement is not at all obvious at intermediate steps of the calculation. For a pure exponential potential, which has c ,ϕ =0 , the non-gravitational contribution is zero, and the spectrum is slightly blue. For plausible values of c=20 and γ=1/2 , say, the gravitational piece is about one percent and the spectral tilt is n s 1.01 , also consistent with our earlier estimate.

We note that c=20 , is related with the Aurea ratio by the following mathematical formula:

[ Φ 35/7 + Φ 14/7 ] 3 2 =[ ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 14/7 ] 3 2 11.090+2.618=13.708 3 2 =20.562.

In the cyclic model, the steepness of the potential must decrease as the field rolls downhill in order that the ekpyrotic phase comes to an end, which corresponds to c ,ϕ >0 . If c( ϕ ) changes from some initial value c ¯ 1 to some value of order unity at the end of the ekpyrotic phase after ϕ changes by an amount Δϕ , then c ,ϕ c ¯ / Δϕ . When c is large, the non-gravitational term in Equation (153) typically dominates and the spectral tilt is a few per cent towards the red.

For example, suppose c ϕ β and c( ϕ )dϕ 125 ; then, the spectral tilt is:

n s 1=0.03 β 1+β , (155)

which corresponds to 0.97< n s <1 for positive 0<β< , in agreement with our earlier estimate. With regard the value 125, we have the following mathematical connection with the Fibonacci’s numbers:

125=2+5+8+21+34+55

We note that negative potentials of this type with very large values of c have been argued to arise naturally in string theory. Our expression for the spectral tilt of the entropically induced curvature spectrum can also be expressed in terms of the customary “fast-roll” parameters:

ε ¯ ( V V ,ϕ ) 2 = 1 c 2 , η ¯ ( V V ,ϕ ) ,ϕ . (156)

Note that ε ¯ =1/ ( 2( 1+ γ 2 )ε ) . Then, the spectral tilt is:

n s 1= 4( 1+ γ 2 ) M Pl 2 ε ¯ 4 η ¯ . (157)

This result can be compared with the spectral index of the time-delay (Newtonian potential) perturbation, where the corresponding formula is:

n s 1= 4 M Pl 2 ε ¯ 4 η ¯ . (158)

Here, the first term is again gravitational, but it has the opposite sign of the gravitational contribution to the entropically induced fluctuation spectrum. So, the tilt is typically a few per cent redder. Finally, for inflation, the spectral tilt is:

n s 1=6ε+2η (159)

where the result is expressed in terms of the slow-roll parameters ε( 1/2 ) ( M Pl V ,ϕ /V ) 2 and η M Pl 2 V ,ϕϕ /V . Here, we have revealed the factors of M Pl to illustrate that both inflationary contributions are gravitational in origin. This gives the same range for n s as the Newtonian potential perturbations in the cyclic model.

4.1. On Some Equations Concerning the Effective Five-Dimensional Theory of the Strongly Coupled Heterotic String as a Gauged Version of N = 1 Five-Dimensional Supergravity with Four-Dimensional Boundaries

We will now briefly review the effective description of strongly coupled heterotic string theory as 11-dimensional supergravity with boundaries given by Horava and Witten. The bosonic part of the action is of the form:

S= S SG + S YM (160)

where S SG is the familiar 11-dimensional supergravity:

S SG = 1 2 κ 2 M 11 g [ R+ 1 24 G IJKL G IKL + 2 1728 ε I 1 I 11 C I 1 I 2 I 3 G I 4 I 7 G I 8 I 11 ] (161)

and S YM are the two E 8 Yang-Mills theories on the orbifold planes explicitly given by:

S YM = 1 8π κ 2 ( κ 4π ) 2/3 M 10 ( 1 ) g { tr ( F ( 1 ) ) 2 1 2 tr R 2 } 1 8π κ 2 ( κ 4π ) 2/3 M 10 ( 2 ) g { tr ( F ( 2 ) ) 2 1 2 tr R 2 } (162)

Here, F IJ ( i ) are the two E 8 gauge field strengths and C IJK is the 3-form with field strength G IJKL =24 [I C JKL] . In order for the above theory to be supersymmetric as well as anomaly free, the Bianchi identity for G should receive a correction such that:

( dG ) 11 JKL ¯ = 1 2 2 π ( κ 4π ) 2/3 { J ( 1 ) δ( x 11 )+ J ( 2 ) δ( x 11 πρ ) } JKKL ¯ (163)

where the sources are given by:

J ( i ) =( tr F ( i ) F ( i ) 1 2 trRR ) . (164)

With regard Equations (161) and (163), we have the following mathematical connections with the Aurea ratio and the Ramanujan modular equations:

2 =1.414213562 Φ 5/7 = ( 5 +1 2 ) 5/7 =1.4101875 ;

1728=4324=24243= 24 2 3 ;

432=306.342Hz1.4101875817 ; 432=2 3 3 8 ;

1 2 2 π =0.112539539 Φ 31.67/7 = ( 5 +1 2 ) 31.67/7 =0.1133912969

We note that 2, 3 and 8 are Fibonacci’s numbers, while 24 is related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

While the standard embedding of the spin connection into the gauge connection:

tr F ( 1 ) F ( 1 ) =trRR (165)

leads to vanishing source terms in the weakly coupled heterotic string Bianchi identity, in the present case, one is left with non-zero sources ±trRR on the two hyperplanes. As a result, the antisymmetric tensor field G and, hence, the second term in the gravitino supersymmetry variation:

δ Ψ I = D I η+ 2 288 ( Γ IKLM 8 g IJ Γ KLM ) G JKLM η+ (166)

do not vanish.

With regard the mathematical connections with the Aurea ratio and Fibonacci’s numbers, we note that:

432=306.342 Φ 5/7 =306.342 ( 5 +1 2 ) 5/7 =306.3421.4101875817 ;

432288=144 and 144 is a Fibonacci’s number. Furthermore, 288=2412 and 24 are related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

Now, let us start with the zeroth order metric:

d s 11 2 = η μν d x μ d x ν + R 0 2 ( d x 11 ) 2 + V 0 1/3 Ω AB d x A d x B , (167)

where Ω AB is a Calabi-Yau metric with Kahler form ω a b ¯ =i Ω a b ¯ . (Here, a and b ¯ are holomorphic and anti-holomorphic indices). To keep track of the scaling properties of the solution, we have introduced moduli V 0 and R 0 for the Calabi-Yau volume and the orbifold radius, respectively. To order κ 2/3 , the metric can be written in the form:

d s 11 2 =( 1+ b ^ ) η μν d x μ d x ν + R 0 2 ( 1+ γ ^ ) ( d x 11 ) 2 + V 0 1/3 ( Ω AB + h AB )d x A d x B (168)

where the functions b ^ , γ ^ and h AB depend on x 11 and the Calabi-Yau coordinates. Furthermore, G ABCD and G ABC11 receive a contribution of order κ 2/3 from the Bianchi identity source terms. The general explicit form of the corrections are:

b ^ = 2 3 R 0 V 0 2/3 α( | x 11 | πρ/2 ) , γ ^ = 2 2 3 R 0 V 0 2/3 α( | x 11 | πρ/2 ) (169)

h AB = 2 3 R 0 V 0 2/3 α( | x 11 | πρ/2 ) Ω AB , G ABCD = 1 6 α ε ABCD EF ω EF ε( x 11 ) (170)

with

α= 1 8 2 πv ( κ 4π ) 2/3 X ωtr R ( Ω ) R ( Ω ) , v= X Ω . (171)

Here, ε( x 11 ) is the step function which is +1 (−1) for x 11 positive (negative).

With regard Equations (169), (170) and (171), we have the following mathematical connections with the Aurea ratio:

2 3 =0.47140452 Φ 11/7 = ( 5 +1 2 ) 11/7 =0.469451 ;

2 2 3 =0.942809041 Φ 7/7 + Φ 35/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 35/7 =0.618034+0.090170=0.708204;

0.708204 4 3 =0.944272 ;

1 6 =0.1666 Φ 26/7 = ( 5 +1 2 ) 26/7 =0.1674018269 ;

1 8 2 π =0.0281348 1 25 Φ 18.33/7 = 1 25 ( 5 +1 2 ) 18.33/7 =0.02835642 .

In the five-dimensional space M 5 of the reduced theory, the orbifold fixed planes constitute four-dimensional hypersurfaces which we denote by M 4 ( i ) ,i=1,2 . There will be an E 6 gauge field A μ ( 1 ) accompanied by gauginos and gauge matter fields on the orbifold plane M 4 ( 1 ) . We will set these gauge matter fields to zero in the following. The field content of the orbifold plane M 4 ( 2 ) consists of an E 8 gauge field A μ ( 2 ) and the corresponding gauginos. In addition, there is another important boundary effect which results from the non-zero internal gauge field and gravity curvatures. More precisely, note that:

X Ω tr F AB ( 1 ) F ( 1 )AB = X Ω tr R AB R AB =16 2 πv ( 4π κ ) 2/3 α, F AB ( 2 ) =0 . (172)

In view of the boundary action (162), it follows that we will retain cosmological type terms with opposite signs on the two boundaries. Note that the size of those terms is set by the same constant α , given by Equation (171), which determines the magnitude of the non-zero mode. The boundary cosmological terms are another important ingredient in reproducing the 11-dimensional background as a solution of the five-dimensional theory.

We can perform the Kaluza-Klein reduction on the metric:

d s 11 2 = V 2/3 g αβ d x α d x β + V 1/3 Ω AB d x A d x B . (173)

The complete configuration for the antisymmetric tensor field that we use in the reduction is given by:

C αβγ , G αβγδ =24 [α C βγδ] , C αAB = 1 6 A α ω AB ,

G αβAB = αβ ω AB , αβ = α A β β A α ,

C ABC = 1 6 ξ ω ABC , G αABC = α ξ ω ABC (174)

and the non-zero mode is:

G ABCD = α 6 ε ABCD   EF ω EF ε( x 11 ) , (175)

where α was defined in Equation (171).

We can now compute the five-dimensional effective action of Horava-Witten theory. Using the field configuration (172)-(175), we find from the action (160)-(161) that:

S 5 = S grav + S hyper + S bound (176)

where

S grav = 1 2 κ 5 2 M 5 g [ R+ 3 2 αβ αβ + 1 2 ε αβγδδ A α βγ s ] (177)

S hyper = 1 2 κ 5 2 M 5 g [ 1 2 V 2 α V α V+2 V 1 α ξ α ξ ¯ + 1 24 V 2 G αβγγ G αβγδ + 2 24 ε αβγγδ G αβγδ ( i( ξ ε ξ ¯ ξ ¯ ε ξ )+2α A ε )+ 1 3 V 2 α 2 ] (178)

S bound = 1 2 κ 5 2 { 2 2 M 4 ( 1 ) g V 1 α +2 2 M 4 ( 2 ) g V 1 α } 1 16π α GUT i=1 2 M 4 ( i ) g Vtr F μν ( i ) 2 . (179)

In this expression, we have now dropped higher-derivative terms. The 4-form field strength G αβγδ is subject to the Bianchi identity:

( dG ) 11μvρσ = κ 5 2 4 2 π α GUT { J ( 1 ) δ( x 11 )+ J ( 2 ) δ( x 11 πρ ) } μνρσ (180)

which follows directly from the 11-dimensional Bianchi identity (163). The currents J ( i ) have been defined in Equation (164). The five-dimensional Newton constant κ 5 and the Yang-Mills coupling α GUT are expressed in terms of 11-dimensional quantities as:

κ 5 2 = κ 2 v , α GUT = κ 2 2v ( 4π κ ) 2/3 . (181)

Since we have compactified on a Calabi-Yau space, we expect the bulk part of the above action to have eight preserved supercharges and, therefore, to correspond to minimal N=1 supergravity in five dimensions. Accordingly, let us compare the result (177)-(179) to the known N=1 supergravity matter theories in five dimensions. In these theories, the scalar fields in the universal hypermultiplet parameterize a quaternionic manifold with coset structure M Q =SU( 2,1 )/SU( 2 )×U( 1 ) . Hence, to compare our action to these we should dualize the three-form C αβγ to a scalar field σ by setting (in the bulk):

G αβγδ = 1 2 V 2 ε αβγδε ( ε σi( ξ ε ξ ¯ ξ ¯ ε ξ )2α A ε ) . (182)

Then, the hypermultiplet part of the action (178) can be written as:

S hyper = v 2 κ 2 M 5 g [ h uv α q u α q v + 1 3 V 2 α 2 ] (183)

where q u =( V,σ,ξ, ξ ¯ ) . The covariant derivative α is defined as α q u = α q u +α A α k u with k u =( 0,2,0,0 ) . The sigma model metric h uv = u v K Q can be computed from the Kahler potential:

K Q =ln( S+ S ¯ 2C C ¯ ),S=V+ξ ξ ¯ +iσ,C=ξ . (184)

Consequently, the hypermultiplet scalars q u parameterize a Kahler manifold with metric h uv . It can be demonstrated that k u is a Killing vector on this manifold.

To analyze the supersymmetry properties of the solutions shortly to be discussed, we need the supersymmetry variations of the fermions associated with the theory (176). They can be obtained either by a reduction of the 11-dimensional gravitino variation (166) or by generalizing the known five-dimensional transformations by matching onto gauged four-dimensional N=2 theories. It is sufficient to keep the bosonic terms only. Both approaches lead to:

δ ψ α i = D α ε i + 2 i 8 ( γ α βγ 4 δ α β γ γ ) J ˜ βγ ε i 1 2 V 1/2 ( α ξ ( τ 1 i τ 2 ) j i α ξ ¯ ( τ 1 +i τ 2 ) j i ) ε j 2 i 96 V ε α βγδδ G βγγδ ( τ 3 ) j i ε j 2 12 α V 1 ε( x 11 ) γ α ( τ 3 ) j i ε j

δ ζ i = 2 48 V ε αβγδδ G αβγγ γ ε ε i i 2 V 1/2 γ α ( α ξ ( τ 1 i τ 2 ) j ) 2 + α ξ ¯ ( τ 1 +i τ 2 ) j i ) ε j + i 2 V 1 γ β β V ε i i 2 α V 1 ε( x 11 ) ( τ 3 ) j i ε j (185)

where τ i are the Pauli spin matrices. Thence, we see that the relevant five-dimensional effective theory for the reduction of Horava-Witten theory is a gauged N=1 supergravity theory with bulk and boundary potentials.

The theory (176) has all of the prerequisites necessary for such a three-brane solution to exist. Generally, in order to have a ( D2 ) -brane in a D -dimensional theory, one needs to have a ( D1 ) -form field, or, equivalently, a cosmological constant. This cosmological term is provided by the bulk potential term in the action (176). From the viewpoint of the bulk theory, we could have multi three-brane solutions with an arbitrary number of parallel branes located at various places in the x 11 direction. As is well known, however, elementary brane solutions have singularities at the location of the branes, needing to be supported by source terms. The natural candidates for those source terms are the boundary actions. Given the anomaly-cancellation requirements, this restricts the possible solutions to those representing a pair of parallel three-branes corresponding to the orbifold planes. It is clear that in order to find a three-brane solution, we should start with the Ansatz:

d s 5 2 =a ( y ) 2 d x μ d x ν η uv +b ( y ) 2 d y 2 (186)

V=V( y )

where a and b are functions of y= x 11 and all other field vanish. The general solution for this Ansatz, satisfying the equations of motion derived from the action (176), is given by:

a= a 0 H 1/2 ,b= b 0 H 2 ,V= b 0 H 3 andH= 2 3 α| y |+ c 0 (187)

where a 0 , b 0 and c 0 are constants. We note that the boundary source terms have fixed the form of the harmonic function H in the above solution. Without specific information about the sources, the function H would generically be glued together from an arbitrary number of linear pieces with slopes ± 2 3 α . The edges of each piece would then indicate the location of the source terms. The necessity of matching the boundary sources at y=0 and πρ , however, has forced us to consider only two such linear pieces, namely y[ 0,πρ ] and y[ πρ,0 ] . These pieces are glued together at y=0 and πρ . Therefore, we have:

a= a 0 H 1/2 ,b= b 0 H 2 ,V= b 0 H 3 andH= 2 3 α| y |+ c 0 (188)

which shows that the solution represents two parallel three-branes located at the orbifold planes. We stress that this solution solves the five-dimensional theory (176) exactly, whereas the original deformed Calabi-Yau solution was only an approximation to order κ 2/3 . It is straightforward to show that the linearized version of (187), that is, the expansion to first order in α=O( κ 2/3 ) , coincides with Witten’s solution (168)-(169)-(170) upon appropriate matching of the integration constants. Hence, we have found an exact generalization of the linearized 11-dimensional solution. We still have to check that our solution preserves half of the supersymmetries. When g αβ and V are the only non-zero fields, the supersymmetry transformations (185) simplify to:

δ ψ α i = D α ε i 2 12 αε( y ) V 1 γ α ( τ 3 ) j i ε j δ ζ i = i 2 V 1 γ β β V ε i i 2 αε( y ) V 1 ( τ 3 ) j i ε j . (189)

With regard Equations (179), (180), (182), (185) and (187)-(189), we have the following mathematical connections with Aurea ratio and with the Ramanujan’s modular equations:

1 2 =0.707106781 Φ 5/7 = ( 5 +1 2 ) 5/7 =0.709125 ;

2 24 =0.058925565 1 25 Φ 7.67/7 = 1 25 ( 5 +1 2 ) 7.67/7 =0.0590213661 ;

2 2 =2.828427125 Φ 15/7 = ( 5 +1 2 ) 15/7 =2.8043399 ;

1 4 2 =0.17677 Φ 25.33/7 = ( 5 +1 2 ) 25.33/7 =0.17525 ;

1 4 2 π =0.056269769 1 25 Φ 8.33/7 = 1 25 ( 5 +1 2 ) 8.33/7 =0.05640 ;

2 8 =0.176776695 Φ 25/7 = ( 5 +1 2 ) 25/7 =0.179314 ;

2 12 =0.11785113 Φ 31/7 = ( 5 +1 2 ) 31/7 =0.118708 ;

2 48 =0.029462782 1 25 Φ 18/7 = 1 25 ( 5 +1 2 ) 18/7 =0.0290137 ;

2 96 =0.014731391 1 25 Φ 28/7 = 1 25 ( 5 +1 2 ) 28/7 =0.0145898 ;

2 3 =0.47140452 Φ 11/5 = ( 5 +1 2 ) 11/5 =0.469451 .

Furthermore, the number 8, 12, 24, 48 and 96 are connected with the “modes” that correspond to the physical vibrations of a superstring and to physical vibrations of the bosonic strings by the following Ramanujan functions:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

24= 4[ antilog 0 cosπtx w coshπx e π 2 4 w ϕ w ( it x ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ]

We note also that 8 is a Fibonacci’s number.

The Killing spinor equations δ ψ α i =0,δ ζ i =0 are satisfied for the solution (187) if we require that the spinor ε i is given by:

ε i = H 1/4 ε 0 i , γ 11 ε 0 i = ( τ 3 ) j i ε 0 j (190)

where ε 0 i is a constant symplectic Majorana spinor. This shows that we have indeed found a BPS solution preserving four of the eight bulk supercharges.

4.2. On Some Equations Concerning the Colliding Branes and the Origin of the Hot Big Bang

We have derived the five-dimensional effective action of heterotic M-theory in the precedent subsection in 3. Now, we shall use a simplified action describing gravity g γδ , the universal “breathing” modulus of the Calabi-Yau three-fold ϕ , a four-form gauge field A γδδζ with field strength =dA and a single bulk M5-brane. It is given by:

S= M 5 3 2 5 d 5 x g ( R 1 2 ( ϕ ) 2 3 2 e 2ϕ 2 5! ) 3 i=1 3 α i M 5 3 4 ( i ) d 4 ξ ( i ) ( h ( i ) e ϕ ε μvκ 4! A γδδζ μ X ( i ) γ V X ( i ) δ κ X ( i ) ε λ X ( i ) ζ ), (191)

where γ,δ,ε,ζ=0,,4;μ,ν,=0,,3 . The space-time is a five-dimensional manifold M 5 with coordinates x γ . The four-dimensional manifolds M 4 ( i ) ,i=1,2,3 are the visible, hidden, and bulk branes respectively, and have internal coordinates ξ ( i ) μ and tension α i M 5 3 . Note that α i has dimension of mass. If we denote α 1 α , α 2 αβ , and α 3 β , then the visible brane has tension α M 5 3 , the hidden brane ( αβ ) M 5 3 , and the bulk brane β M 5 3 . It is straightforward to show that the tension of the bulk brane, β M 5 3 , must always be positive. Furthermore, one can easily deduce that the tension on the visible brane, α M 5 3 , can be either positive or negative. We will take α>0 , so that the tension on the visible brane is negative. Furthermore, we will choose β such that αβ>0 , that is, the tension of the hidden brane is positive. The tensor h μν ( i ) is the induced metric on M 4 ( i ) . The functions X ( i ) γ ( ξ ( i ) μ ) are the coordinates in M 5 of a point on M 4 ( i ) with coordinates ξ ( i ) μ . In other words, X ( i ) γ ( ξ ( i ) μ ) describe the embedding of the branes into M 5 . The BPS solution of Lukas, Ovrut and Waldram is then given by:

d s 2 =D( y )( N 2 d τ 2 + A 2 d x 2 )+ B 2 D 4 ( y )d y 2 ; e ϕ =B D 3 ( y )

0123Y =α A 3 N B 1 D 2 ( y )  for  y<Y

1123v =( αβ ) A 3 N B 1 D 2 ( y )   fory>Y (192)

where

D( y )=αy+C    fory<Y

D( y )=( αβ )y+C+βY   fory>Y (193)

and A,B,C,N and Y are constants. Note that A,B,C,N are dimensionless and Y has the dimension of length. The visible and hidden boundary branes are located at y=0 and y=R , respectively, and the bulk brane is located at y=Y , 0YR . We assume that C>0 , so that the curvature singularity at D=0 does not fall between the boundary branes. Note that y=0 lies in the region of smaller volume while y=R lies in the region of larger volume. Note that inserting the solution of the four-form equation of motion into Equation (191) yields precisely the bulk action with charge α in the interval 0yY and charge α+β in the interval YyR . The formulation of the action Equation (191) using the four-form A is particularly useful when the theory contains bulk branes, as is the case in ekpyrotic theory.

The following equation:

| δ k |=4 a c H c ( 1 2 H ˙ c a c H c 2 )| Δτ( k ) | , (194)

expresses the density perturbation in terms of the time delay at the time of collision, Δτ( k ) . If we consider the exponential potential V=v e mαY , then Equation (194) yields:

| δ k | 4 m 2 α 2v mC+2 | Δτ( k ) | . (195)

Now, we compute the spectrum of quantum fluctuations of the brane δ Y k and use the result to compute the time delay, Δτ( k ) .

For the calculation of quantum fluctuations, it is sufficient to work at the lowest order in β/α . Without loss of generality, we can therefore set A=N=1 . In that case, the bulk brane Lagrangian is given by:

L β =3β M 5 3 B[ 1 2 D ( Y ) 2 η μν μ Y ν YV( Y ) ] . (196)

Note that this agrees with L β given in the following equation:

L β = 3β M 5 3 A 3 B N [ 1 2 D 2 ( Y ) Y ˙ 2 N 2 V( Y ) ] , (197)

when we set A=N=1 and spatial gradients of Y to zero. Let us first consider the spatially homogeneous motion of the brane which will be described by Y 0 ( τ ) . It is governed by the following equation of motion:

1 2 D ( Y 0 ) 2 Y ˙ 0 2 +V( Y 0 )=E , (198)

where E is a constant. Equation (198) is simply the statement that the energy E of the bulk brane is conserved to this order in β/α . Since we have chosen the visible brane to lie at y=0 and the hidden universe to lie at y=R , we focus on the branch Y ˙ <0 in which case the bulk brane moves towards the visible brane. The solution to Equation (198) is then given by:

( τ )= 0 Y 0 D( Y )d Y 2( EV( Y ) ) (199)

with τ0 , and with the collision occurring at τ=0 . Let us now consider fluctuations around the background solution Y 0 ( τ ) . Namely, if Y= Y 0 ( τ )+δY( τ, x ) , with δY( τ, x ) Y 0 ( τ ) , we can expand the action to quadratic order in δY :

L fuc 1 2 D 0 2 [ δ Y ˙ 2 + ( ( δY ) ) 2 ] +[ α 2 D 0 2 ( V 0 E )α D 0 1 d V 0 d Y 0 + 1 2 d 2 V 0 d Y 0 2 ] ( δY ) 2 , (200)

where we have used Equation (198), and where we have introduced D 0 D( Y 0 ) and V 0 V( Y 0 ) for simplicity. The key relation is the fluctuation equation as derived from the action (200):

x 2 d 2 f k d x 2 [ a ¨ pert a pert τ 2 x 2 ] f k =0;x| k |( τ ) , (201)

where f k D 0 δ Y k and where a pert is defined by:

a ¨ pert a pert D 0 3 ( α d V 0 d Y 0 D 0 d 2 V 0 d Y 0 2 ) . (202)

The fluctuation Equation (201), can be compared with the corresponding equation for the perturbations of a scalar field with no potential and minimally coupled to an FRW background with scale factor a( τ ) :

δ ϕ ¨ k +2 a ˙ a δ ϕ ˙ k + k 2 δ ϕ k =0 . (203)

Defining f k =aδ ϕ k ^ , Equation (203) becomes:

x 2 d 2 f k d x 2 [ a ¨ a τ 2 x 2 ] f k =0 . (204)

Let us now discuss the Hubble horizon for the perturbations. Recall that in usual 4 d cosmology (see Equation (204)), we have:

x=k( τ )=( k a )a( τ )= k phys a( τ ) k phys H 1 , (205)

where H 1 a 2 / a ˙ is the Hubble radius as derived from the scale-factor a . By definition, a mode is said to be outside the Hubble horizon when its wavelength is larger than the Hubble radius. From Equation (205), we see that this occurs when x<1 . Therefore, a mode with amplitude f k crosses outside the horizon when xO( 1 ) . Similarly, in our scenario, we can write:

x=k( τ )= k phys D 0 1/2 ( τ ) k phys H pert 1 , (206)

where k phys =k/ D 0 1/2 . The role of the Hubble radius is replaced by:

H pert 1 D 0 1/2 ( τ )= D 0 1/2 0 Y 0 D( Y )dY 2( EV( Y ) ) , (207)

which is to be thought of as an effective Hubble radius for the perturbations. So, the length scale at which amplitudes freeze depends on a B (rather than a pert ), but the amplitude itself, as derived from Equation (201), depends on a pert . The feature of two different scale factors is novel aspect of ekpyrotic scenario. With regard the comparison to inflationary cosmology, we have that in inflation, the wavelengths are stretched superluminally while the horizon is nearly constant. In the ekpyrotic scenario, the wavelengths are nearly constant while the horizon shrinks. We can obtain a spectrum which is scale-invariant. Writing the equation for the perturbations in the form of Equation (201) is useful since one can read off from it the spectral slope of the power spectrum. It is determined by the value of ( a ¨ pert / a pert ) τ 2 . In particular, one obtains a scale-invariant spectrum if ( a ¨ pert / a pert ) τ 2 =2 when the modes observed on the CMB cross outside the horizon.

Combining Equations (199) and (202), we find:

a ¨ pert a pert τ 2 = D 0 3 ( α d V 0 d Y 0 D 0 d 2 V 0 d Y 0 2 ) [ 0 Y D( Y )d Y 2( EV( Y ) ) ] 2 . (208)

The spectrum will be scale-invariant if the right hand side of Equation (208) equals 2 when the modes of interest cross outside the horizon. Thence, we have:

a ¨ pert a pert τ 2 = D 0 3 ( α d V 0 d Y 0 D 0 d 2 V 0 d Y 0 2 ) [ 0 Y D( Y )d Y 2( EV( Y ) ) ] 2 =2 . (209)

With regard Equation (209), we have the following mathematical connections with the Aurea ratio:

Φ 14/7 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 14/7 =2.618034+0.381966=3; 2 3 3=2;

Φ 21/7 + Φ 7/7 + Φ 28/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 28/7 =4.236068+1.618034+0.145898=6; 1 3 6=2.

We note also that 2 is prime number and Fibonacci’s number.

Thence we can rewrite Equation (209) also as follow:

a ¨ pert a pert τ 2 = D 0 3 ( α d V 0 d Y 0 D 0 d 2 V 0 d Y 0 2 ) [ 0 Y D( Y )d Y 2( EV( Y ) ) ] 2 = 2 3 [ Φ 14/7 + Φ 14/7 ]= 2 3 [ ( 5+1 2 ) 14/7 + ( 5 +1 2 ) 14/7 ] = 2 3 ( 2.618034+0.381966 )=2;

We can add a potential V( Y ) of the form that might result from the exchange of wrapped M2branes. We would like to think of V as the potential derived from the superpotential W for the modulus Y in the 4 d low energy theory. Typically, superpotentials for such moduli are of exponential form, for example,

W e cY , (210)

where c is a positive parameter with dimension of mass. The corresponding potential is constructed from W and the Kahler potential K according to the usual prescription:

V= e K/ M pl 2 [ K ij D i W D j W ¯ 3 M pl 2 W W ¯ ] . (211)

where D i =/ ϕ i + K i / M pl 2 is the Kahler covariant derivative, K i = K/ ϕ i , K ij = 2 K/ ϕ i ϕ j and a sum over each superfield ϕ i is implicit. Equations (210) and (211) imply that V decays exponentially with Y . Here, it will suffice to perform the calculation using a simple exponential potential, namely:

V( Y )=v e mαY , (212)

where v and m are positive, dimensionless constants. Note that, in the case where the potential is generated by the exchange of wrapped M2-branes, the parameter m is of the form m= c T 3 ν/α , where c is a constant, T 3 is the tension of the M2-brane, and v is the volume of the curve on which it is wrapped. The perturbation modes of interest are those which are within the current Hubble horizon. As the wavelengths corresponding to those modes passed outside the effective Hubble horizon on the moving bulk-brane, the amplitudes became fixed. Scale invariance will require mD>1 during this period. We know that, if the potential V is negligible compared to E , the spectrum of fluctuations is not scale-invariant. Hence, we consider the limit where | E || V 0 | . This condition, as seen from the equation of motion for Y 0 , Equation (198), is satisfied if Y ˙ 0 =0 initially, or, equivalently, if the bulk brane begins nearly at rest. For the brane to be nearly at rest, one must have | E || V 0 | initially. As the brane traverses the fifth dimension, | V | increases exponentially, whereas E is constant. Hence, the condition | E || V 0 | is automatically satisfied. The bulk brane beginning nearly at rest is precisely what we expect for a nearly BPS initial state. Applying the condition | E |<| V 0 | , Equation (199) reduces to:

τ 2 1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) , (213)

where we have neglected the endpoint contribution at Y=0 . On the other hand, Equation (202) gives:

a ¨ pert a pert = m 2 α 2 v e mα Y 0 D 0 2 ( 1+ 1 m D 0 ) . (214)

Combining the above two expressions, we obtain:

D 0 2 m 2 α 2 v e mα Y 0 =( 1+ 1 m D 0 ), a ¨ pert a pert ; τ 2 =2( 1+ 1 m D 0 ), a ¨ pert a pert ( 1 2 m D 0 )

a ¨ pert a pert τ 2 =2( 1+ 1 m D 0 )( 1 2 m D 0 ). (215)

The right hand side of Equation (215) is approximately equal to 2 in the limit of large m D 0 . Thence, we have:

a ¨ pert a pert τ 2 =2( 1+ 1 m D 0 )( 1 2 m D 0 )=2 . (216)

Also here, with regard the numerical result of Equation (216), we have the following mathematical connections with Aurea ratio:

Φ 14/7 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 14/7 =2.618034+0.381966=3; 2 3 3=2;

Φ 21/7 + Φ 7/7 + Φ 28/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 28/7 =4.236068+1.618034+0.145898=6; 1 3 6=2.

Hence, the exponential potential of Equation (212) results in a nearly scale-invariant spectrum of perturbations provided that | E |<| V 0 | and m D 0 >1 are satisfied when modes pass outside the effective Hubble horizon.

We next compute the perturbation amplitude, by using Equation (201) to calculate | Δ Y k | . The conditions | E || V 0 | and m D 0 1 must be satisfied when wavelengths pass outside the horizon. These conditions can be relaxed once the mode is well outside the horizon. In the limit that m D 0 >1 when the relevant modes cross outside the horizon, Equation (201) reduces to:

x 2 d 2 f k d x 2 [ 2 x 2 ] f k =0 , (217)

with solution:

f k = x 1/2 ( C 1 ( k ) J 3/2 ( x )+ C 2 ( k ) J 3/2 ( x ) ) , (218)

where J ±3/2 are Bessel functions. The coefficients C 1 ( k ) and C 2 ( k ) are fixed by requiring that modes well-within the horizon (i.e. x1 ) be Minkowskian vacuum fluctuations, that is:

f k = 1 6kβ M 5 3 B e ikτ forx1 . (219)

Using this initial condition, we find the following amplitude for modes outside the horizon (with x1 ):

Δ f k k 3/2 f k ( 2π ) 3/2 = i ( τ ) ( 2π ) 3/2 6β M 5 3 B . (220)

Substituting Equation (213) and using f k = D 0 δ Y k , we find:

Δ Y k = mα 2 ( 2π ) 3/2 3β M 5 3 B v e mα Y 0 D 0 2 . (221)

Finally, we define the time-delay Δτ( k ) by:

| Δτ( k ) |=| Δ Y k Y ˙ 0 |= m 2 α 16 π 3/2 3β M 5 3 B ( 2 m D 0 ) , (222)

where we have used the equation of motion for Y 0 , Equation (198). Note that the time-dependence of Δτ( k ) is mild, a necessary condition for the validity of the time-delay formalism. The factor of m D 0 mD( Y 0 ( τ ) )  is to be evaluated at time τ when a given mode crosses outside the horizon during the motion of the bulk brane. Let D k denote the value of D 0 at horizon crossing for mode k . Since horizon crossing occurs when x=1 , or, equivalently, when ( τ )= k 1 , Equation (213) gives:

D k 2 m log( m 2 α 2k v e mC 2 ) . (223)

Substituting Equations (222) and (223) into Equation (195), we find:

| δ k |= α 2 m 4 2v 4 π 3/2 3β M 5 3 B ( mC+2 ) ( 2 m D k ) . (224)

This expression for | δ k | increases gradually with increasing k , corresponding to a spectrum tilted slightly towards the blue. The blue tilt is due to the fact that, in this example, D is decreasing as the brane moves. That is, the spectral index,

n s 1+ dlog | δ k | 2 dlogk 1+ 4 m D k , (225)

exceeds unity. The current CMB data constrains the spectral index to lie in the range about 0.8< n s <1.2 . Therefore, for our results to be consistent with experiments, we must have:

m D k >20 , (226)

a constraint that is easily satisfied.

With regard the value s of n s , i.e. 0.8 and 1.2, we have the following mathematical connections with Aurea ratio:

Φ 7/7 = ( 5 +1 2 ) 7/7 =1.6180339871.618034;1.618034 3 4 =1.2135255 ;

Φ 0 + Φ 42/7 + Φ 56/7 = ( 5 +1 2 ) 0 + ( 5 +1 2 ) 42/7 + ( 5 +1 2 ) 56/7 =1+0.055728+0.021286=1.077014;1.077014 3 4 =0.8077605.

For the value of Equation (226), we have the following mathematical connections:

Φ 35/7 + Φ 21/7 + Φ 14/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 14/7 =11.090170+4.236068+0.381966=15.708204 4 3 =20.944272.

Thence, we obtain:

m D k >20 = Φ 35/7 + Φ 21/7 + Φ 14/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 14/7 =11.090170+4.236068+0.381966=15.708204 4 3 =20.944272.

We note that 20.94427221 ; 21=13+8 , where 8 and 13 are Fibonacci’s numbers. Furthermore, the number 8 is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

Furthermore, we consider the power-law potential:

V( Y )=vD ( Y ) q =v ( αY+C ) q , (227)

where v>0 and q<0 are constants. In this case, Equation (208) gives:

a ¨ pert a pert τ 2 2 ( 1 2 q ) ( 1 4 q ) 2 2 (228)

for | q |>1 . Hence, a power-law potential can also lead to a nearly scale-invariant spectrum provided that its exponent is sufficiently large. We can straightforwardly extend our analysis to an arbitrary potential V( Y ) . Let us suppose that V( Y ) satisfies:

| D( Y ) dV dY |α| V( Y ) |,| D( Y ) d 2 V d Y 2 |α| dV dY | . (229)

Then, Equation (208) reduces to:

a ¨ pert a pert τ 2 2( V V V 2 ) . (230)

Hence, the conditions for scale invariance are Equations (229) as well as:

V V V 2 1 . (231)

Also, for the numerical value of Equations (228) and (231), we have the following mathematical connections with the Aurea ratio:

Φ 14/7 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 14/7 =2.618034+0.381966=3; 2 3 3=2;

Φ 21/7 + Φ 7/7 + Φ 28/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 28/7 =4.236068+1.618034+0.145898=6; 1 3 6=2;

Φ 14/7 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 14/7 =2.618034+0.381966=3; 1 3 3=1.

If g ¯ μν is the unperturbed, homogeneous metric (see Equation (192) with A and N functions of time), the perturbed 5d metric can be written as:

g μν = g ¯ μν + A 2 ( t )D( y,t ) h μν ( x ,t ), (232)

where μ,v=0,,3 . We can treat the tensor perturbations h μν as functions of x and t only. We are interested in the tensor perturbations which satisfy the conditions: h 0μ =0 , h j i =0 , and i h ij =0 . The perturbed 5 d Einstein action to quadratic order is:

S fluct T M 5 3 2 d 5 x g R= M 5 2 8 d 4 x a 2 ( h ˙ v μ h ˙ μ v i h v μ i h μ v ) (233)

where the second expression is obtained by integrating over y . The tensor action is analogous to the scalar action given in Equation (200). From the action, we can derive the tensor analogue of the scalar fluctuation equation of motion, Equation (201):

x 2 d 2 f k T d x 2 [ a ¨ a τ 2 x 2 ] f k T =0 , (234)

where

h μ v d 3 k ( 2π ) 3 ε v μ h k ( τ ) (235)

and

f k T a h k (236)

The critical difference between this tensor equation and the scalar fluctuation equation, Equation (201), is that the effective scale factor a pert in Equation (201) has been replaced by a .

We introduced a potential to insure that a pert led to a nearly scale-invariant spectrum, ( a ¨ pert / a pert ) τ 2 2 . However, a( τ ) in the tensor equation is approximately constant (recall that a= ( B I 3 ( 0 ) M 5 ) 1/2 +O( β/α ) ) . Consequently, the root mean square tensor fluctuation amplitude:

| Δ h k | k 3/2 h k ( 2π ) 3/2 k ( 2π ) 3/2 . (237)

is not scale-invariant.

With regard Equation (237), we have obtained the following mathematical connections with the Aurea ratio:

( 2π ) 3/2 =15.74960995 ;

Φ 32/7 + Φ 40/7 = ( 5 +1 2 ) 32/7 + ( 5 +1 2 ) 40/7 =15.7500424315.75 ;

Φ 19/7 = ( 5 +1 2 ) 19/7 =0.2708618458; arcsin( 0.2708618458 ) 180 π =15.715558;

Φ 28/7 + Φ 14/7 + Φ 0 = ( 5 +1 2 ) 28/7 + ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 0 =10.472136 ;

10.472136 3 2 =15.708204 .

The cyclic story can be described in terms of an ordinary four-dimensional field theory, which can be obtained by taking the long wavelength limit of the brane picture. The distance between branes becomes a moduli (scalar) field ϕ . The interbrane interaction is replaced by a scalar field potential, V( ϕ ) . The different stages in the cyclic model in the brane picture are in one-to-one correspondence to the motion of the scalar field along the potential. Then, the action S describing gravity, the scalar field ϕ , and the matter-radiation fluid is:

S= d 4 x g [ 1 16πG 1 2 ( ϕ ) 2 V( ϕ )+ β 4 ( ϕ ) ρ R ] , (238)

where g is the determinant of the Friedmann-Robertson-Walker metric g μν , G is Newton’s constant and R is the Ricci scalar.

The β factor has the property that β as a0 such that aβ constant. The revised solution to the equation of motion is ρ R 1/ ( aβ ) 4 which approaches a constant as a0 . The energy, once thinned out during the dark energy dominated phase, remains thinned out at the bounce. The β -factor simply reflects the fact that the extra-dimension collapses but our three-dimensions do not. As a result, entropy produced during one cycle is not concentrated at the crunch and does not contribute significantly to the entropy density at the beginning of the next cycle. Hence, cycles can continue for an arbitrarily long time and there is no practical way of distinguishing one cycle from the next.

If the cyclic model can be described in terms of ordinary field theory, then it may seem surprising that it is possible to generate a nearly scale invariant spectrum density perturbations. There are actually three distinct ways of producing a nearly scale-invariant spectrum, and that inflation represents only one of them. The three ways can be characterized by:

w ( 1 2 ϕ ˙ 2 V ) ( 1 2 ϕ ˙ 2 +V ) ,

the effective equation of state of the scalar field. Case I is where w1 and the universe is expanding, the example of inflation. Case II is a contracting universe with w0 . Case III is a contracting universe with w1 , that is the situation that applies in the cyclic model.

What is required to obtain w1 ? From the expression for w , it is apparent that this is only possible if the potential is negative. In particular, for a negative exponentially steep potential Vexp( cϕ ) , the solutions to the equation of motion have a scaling solution in which ϕ ˙ 2 / 2V is constant and approximately −1. Consequently, w is much greater than unity and nearly constant. The generation of fluctuations for w1 can be understood heuristically by examining the perturbed Klein-Gordon equation:

δ ϕ k =( k 2 + a a + V ,ϕϕ )δ ϕ k (239)

where ϕ( x,t ) has been expanded in Fourier components δ ϕ k ( t ) with wavenumber k and prime is derivative with respect to conformal time η . The a /a term is due to gravitational expansion, and the last term is due to the self-interaction of the scalar field. This equation applies equally to inflation and to cyclic models. The cyclic model corresponds to the limit where the gravity term is negligible and, instead, the perturbation equation is driven by the potential term. For the negative exponential potential, for example, the scaling solution corresponds to V ϕ,ϕ 2/ η 2 .

We have defined that w1 . We take the following values: w=4.97 and w=4.23 . We obtain the following mathematical connections with the Aurea ratio:

Φ 7/7 + Φ 7/7 + Φ 21/7 + Φ 63/7 =1.618034+0.618034+0.236068+0.013156=2.485292;

2.4852922=4.970584( forΦ= 5 +1 2 ) ;

Φ 21/7 = Φ 3 = ( 5 +1 2 ) 21/7 =4.236067977; Φ 21/7 = Φ 3 = ( 5 +1 2 ) 3 =0.236067977;

arcsin( 0.236067977 ) 180 π =13.654585;arcsin( 0.23 ) 180 π =13.2913 .

We note that 3, 21 and 13 are Fibonacci’s numbers.

For inflation, the most stringent constraints are on the flat part of the potential, the range of the inflaton field where the density perturbations are generated. The constraints are commonly expressed as bounds on two “slow-roll” parameters:

ε= ( V V ) 2 1 (240)

and

η= V V 1 (241)

For the cyclic model, the analogous constraints are on the steep portion of the potential where perturbations are generated. The constraints can be expressed in terms of two “fast-roll” parameters:

ε ¯ = ( V V ) 2 1 (242)

and

η ¯ =1 V V ( V ) 2 1 . (243)

The first constraint forces the slope to be steep and the second fixes the curvature, where each applies to the range of ϕ where the fluctuations are generated that are within the horizon today. The result is that the constraints in the two models are remarkably similar.

We note that if we take for η ¯ =0.090170 , we have the following mathematical connections with the Aurea ratio:

Φ 21/7 + Φ 49/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 49/7 =0.236068+0.034442=0.270510;

0.270510 1 3 =0.090170 ;

Φ 35/7 = ( 5 +1 2 ) 35/7 =0.090169944 ; arcsin( 0.090169944 ) 180 π =5.1733845

that is a Fibonacci’s number.

(For the various analyzed equations, see Ref. [6]-[8]).

5. On Some Equations Concerning the “Null Energy Condition” (NEC) Violation regarding the Inflationary Models

The metric of the higher dimensional theory is -flat (RF) or -flat up to a conformal factor (CRF):

d s 2 = e 2Ω ( d t 2 + a ¯ 2 ( t )d x 2 )+ g mn d y m d y n , (244)

where the x are the non-compact spatial dimensions; y{ y m } are the extra dimensions; a ¯ ( t ) is the usual FRW scale factor; and

g mn ( t,y )= e 2 Ω ¯ g ¯ mn (245)

where g ¯ mn has Ricci (scalar) curvature =0 , as evaluated in the compact dimensions. We call the metric R -flat (RF) if Ω ¯ = const. and conformally R -flat (CRF) if Ω( t,y )= Ω ¯ ( t,y ) .

Now we develop some basic relations that make it possible to detect easily when a higher dimensional theory is forced to violate the NEC.

To describe a spatially-flat FRW spacetime after dimensional reduction, the metric g mn ( t,y ) and warp function Ω( t,y ) must be functions of time t and extra-dimensional coordinates y m only. We parameterize the rate of change of g mn using quantities ξ and σ mn defined by:

1 2 d g mn dt = 1 k ξ g mn + σ mn (246)

where g mn σ mn =0 and where ξ and σ are functions of time and the extra dimensions.

The space-space components of the energy-momentum tensor are block diagonal with a 3×3 block describing the energy-momentum in the three non-compact dimensions and k×k block for the k compact directions. The 00 component is the higher dimensional energy density ρ .

Associated with the two blocks of space-space components of T IJ are two trace averages:

p 3 1 3 γ 3 μv T μv and p k 1 k γ k mn T mn , (247)

where γ 3,k are respectively the 3×3 and k×k blocks of the higher dimensional space-time metric. Violating the NEC means that T MN n M n N <0 for at least one null vector n M and at least one spacetime point. We find simple methods for identifying a subset of cases where the NEC must be violated. For this purpose, the following two lemmas are very useful:

Lemma 1: If ρ+ p 3 or ρ+ p k is less than zero for any space-time point, then the NEC is violated.

The second lemma utilizes the concept of A-averaged quantities:

Q A = ( Q e AΩ g d k y )/ ( e AΩ g d k y ) ; (248)

that is, quantities averaged over the extra dimensions with weight factor e AΩ where, for simplicity, we restrict ourselves to constant A . Using the fact that the weight function in the A -average is positive definite, a straightforward consequence is:

Lemma 2: If ρ+ p 3 A <0 or ρ+ p k A <0 for any A and any { t,x } , then the NEC must be violated.

To illustrate the utility of A -averaging, we introduce the CRF metric into the higher-dimensional Einstein equations, and then try to express terms dependent on a ¯ in terms of the 4 d effective scale factor using the relation a( t ) e ϕ/2 a ¯ ( t ) , where:

e ϕ k e 2Ω g d k y (249)

and is the 4+k -dimensional Planck length. The 4 d effective scale factor, a( t ) , obeys the usual 4 d Friedmann equations:

( a ˙ a ) 2 = 1 3 ρ 4d (250)

( a ˙ a ) 2 +2 a ¨ a = p 4d (251)

Note that the 4 d effective energy density ρ 4d and pressure p 4d are generally different from ρ and p 3 in the higher dimensional theory if the warp factor is non-trivial. Then, using the Einstein equations, we obtain:

e ϕ e 2Ω ( ρ+ p 3 ) A =( ρ 4d + p 4d ) k+2 2k ξ A 2 k+2 2k ( ξ ξ A ) 2 A σ 2 A (252)

e ϕ e 2Ω ( ρ+ p k ) A = 1 2 ( ρ 4d +3 p 4d )+2( A 4 1 ) k+2 2k ( ξ ξ A ) 2 A k+2 2k ξ A 2 σ 2 A +[ 5+ 10 k +k+A( 3+ 6 k ) ] e 2Ω ( Ω ) 2 A + k+2 2k 1 a 3 d dt ( a 3 ξ A ). (253)

There is a range where:

4A 105k+ k 2 3k6 A * , (254)

which is the case for 13k3 (for CRF). Some theorems below rely on choosing A=2 ; for this value to be within the range given in Equation (254), it is necessary that 8k3 . Since this includes the relevant string and M-theory models, we will implicitly assume this range of k for CRF models. We note that 3, 8 and 13 are Fibonacci’s numbers.

The two relations in Equation (252) can be rewritten:

e ϕ e 2Ω ( ρ+ p 3 ) A = ρ 4d ( 1+w ) k+2 2k ξ A 2 +non-positivetermsforallA (255)

e ϕ e 2Ω ( ρ+ p k ) A = 1 2 ρ 4d ( 1+3w )+ k+2 2k 1 a 3 d dt ( a 3 ξ A )+non-positivetermsforsomeA, (256)

where the values of A that make the last term non-positive are those that are in the range in Equation (254). Recall that w represents the ratio of the total 4 d effective pressure p 4d to the total 4 d effective energy density ρ 4d .

On the left hand side of Equations (255) and (256), both ϕ and A depend on the warp factor, Ω, but the combination is invariant under shifts ΩΩ+C , where C is a constant. Furthermore, the combination tends to have a weak dependence on Ω. For example, if ρ+ p k is homogeneous in { y m } , the left hand side reduces to K( ρ+ p k ) , where the dimensionless coefficient K is not very sensitive to Ω or A ; in particular,

K= k I( A+2 )/ I( A )I( 2 )

where

I( A ¯ ) e A ¯ Ω g d k y (257)

In this notation, the k -dimensional volume of the compact space is V k =I( 0 ) ; then, K is equal to k / V k , a coefficient which is strictly less than unity. Similarly, if ρ+ p k is smooth and Ω has a sharp maximum on some subspace of dimension m and volume v m , then the left hand side of Equation (256) is O( 1 )( m / v m ) ( ρ+ p k ) max , where ( ρ+ p k ) max is the value of ρ+ p k evaluated on the subspace where Ω is maximal.

If the NEC is violated, it must be violated in the compact dimensions; it must be violated strongly ( w k significantly below the minimally requisite value for NEC violation); and the violation in the compact dimensions must vary with time in a manner that precisely tracks the equation-of-state in the 4 d effective theory. The magnitude of the NEC violation is proportional to ρ 4d according to Equation (256), which is roughly 10100 times greater during the inflationary epoch than during the present dark energy dominated epoch. Hence, the source of NEC violation for inflation must be different and 10100 stronger.

The fact that NEC violation is required to have inflation in theories with extra dimensions is unexpected since this was not a requirement in the original inflationary models based on four dimensions only. Curiously, a criticism raised at times about models with bounces from a contracting phase to an expanding phase, such as the ekpyrotic and cyclic alternatives to inflationary cosmology, is that the bounce requires a violation of the NEC (or quantum gravity corrections to GR as the FRW scale factor a( t )0 that serve the same function).

If is true that the violation of the condition NEC (condition of null energy) is required for the inflationary universe model and for the cyclic universe model, then it is possible that for the cyclic model the acceleration and initial exponential expansion of the inflationary phase, is equivalent to the collision between the two Brane-worlds and to the consequent acceleration of the expansion of space immediately after the Big Bang. This could be the explanation of the various cosmological and mathematical connections between the two models. Then, the inflation and the Big Bang would be only phases of the cyclic universe. Every cycle has its phase of Big Bang and its phase of inflation.

On Some Equations Concerning the Evolution to a Smooth Universe in an Ekpyrotic Contracting Phase with w > 1

With regard the evolution to a smooth universe in an ekpyrotic contracting phase with w>1 , we find that the ekpyrotic smoothing mechanism is robust in the sense that the ratio to the proper volume of the smooth region to the mixmaster-like region grows exponentially fast along time slices of constant mean curvature.

In this system, the spacetime is described in terms of a coordinate system ( t, x i ) and a tetrad ( e 0 , e α ) where both the spatial coordinate index i and the spatial tetrad index α go from 1 to 3. Choose e 0 to be hypersurface orthogonal with the relation between tetrad and coordinates of the form e 0 = N 1 t and e α = e α   i i where N is the lapse and the shift is chosen to be zero. Choose the spatial frame { e α } to be Fermi propagated along the integral curves of e 0 . The commutators of the tetrad components are decomposed as follows:

[ e 0 , e α ]= u ˙ α e 0 ( H δ α β + σ α β ) e β

[ e α , e β ]=( 2 a [α δ β] γ + ε αβδ n δ γ ) e γ (258)

where n αβ is symmetric, and σ αβ is symmetric and trace free. The scale invariant tetrad variables are defined by 0 e 0 /H and α e α /H while scale invariant versions of the other gravitational variables are given by:

{ E α i , Σ αβ , A α , N αβ } { e α i , σ αβ , a α , n αβ }/H (259)

Note that the relation between the scale invariant tetrad variables and the coordinate derivatives is:

0 = N 1 t

α = E α i i (260)

where N=NH is the scale invariant lapse. The matter model is a scalar field ϕ with potential V of the form:

V( ϕ )= V 0 e cϕ , (261)

where V 0 and c are positive constants. The scale invariant matter variables are given by:

W= 0 ϕ

S α = α ϕ (262)

V ¯ =V/ H 2 . (263)

The time coordinate t is chosen so that

e t =3H . (264)

Note that this means that the surfaces of constant time are constant mean curvature surfaces. Note also that the singularity is approached as t . Due to Equation (264) the scale invariant lapse satisfies an elliptic equation:

α α N+2 A α α N+N( 3+ Σ αβ Σ αβ + W 2 V ¯ )=3 . (265)

We note that 3 is a Fibonacci’s number. Furthermore, we have the following mathematical connection with the Aurea ratio:

Φ 14/7 + Φ 0 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 0 + ( 5 +1 2 ) 14/7 =2.618034+1+0.381966=4

4 3 4 =3 .

The gravitational quantities E α i , A α , N αβ and Σ αβ satisfy the following hyperbolic evolution equations:

t E α i = E α i N( E α i + Σ α β E β i ) (266)

t A α = A α + 1 2 Σ α β β N α N+N( 1 2 β Σ α β A α Σ α β A β ) (267)

t N αβ = N αβ ε γδ(α Σ δ β) γ N+N( N αβ +2 N (α γ Σ β)γ ε γδ(α γ Σ δ β) ) (268)

t Σ αβ = Σ αβ + <α β> N+ A <α β> N+ ε γδ(α N β) δ γ N +N [ 3 Σ αβ <α A β> 2 N <α γ N β>γ + N γ γ N αβ + ε γδ(α ( γ N β) δ 2 A γ N β) δ )+ S <α S β> ] (269)

Here, parentheses around a pair of indices denote the symmetric part, while angle brackets denote the symmetric trace-free part. The equations of motion for the matter variables are as follows:

t ϕ=NW (270)

t S α = S α +W α N+N[ α W( S α + Σ α β S β ) ] (271)

t W=W+ S α α N+N( α S α 3W2 A α S α V ¯ ϕ ) (272)

In addition, the variables are subject to the vanishing of the following constraint quantities:

( C com ) λi = ε αβλ [ α E β i A α E β i ] N λγ E γ i (273)

( e J ) γ = α N αγ + ε αβγ α A β 2 A α N αγ (274)

( C C ) α = β Σ α β 3 Σ α β A β ε αβγ N βδ Σ δ γ W S α (275)

C G =1+ 2 3 α A α A α A α 1 6 N αβ N αβ + 1 12 ( N γ γ ) 2 1 6 Σ αβ Σ αβ 1 6 W 2 1 6 S α S α 1 3 V ¯ (276)

( C S ) α = S α α ϕ (277)

With regard the value 12 of Equation (276), we have the mathematical connection with the following Ramanujan’s modular equation:

The number 12 = 24/2, is related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

For simplicity, we choose the initial conformal metric to be flat and ( x,y,z ) to be the usual cartesian coordinates for that metric, and we choose the spatial triad to lie along those spatial directions. Thus, the scale free spatial triad becomes:

E α i = H 1 ψ 2 δ α i . (278)

It then follows from Equation (258) that:

A α =2 ψ 1 E α i i ψ

N αβ =0 . (279)

The shear is essentially the trace-free part of the extrinsic curvature, and the constraint equations simplify for a particular rescaling of the trace-free part of the extrinsic curvature with the conformal factor. We therefore introduce the quantity Z αβ defined by:

Σ αβ = ψ 6 Z αβ . (280)

Similar considerations apply to the matter variables, leading us to define the quantity Q given by:

W= ψ 6 Q . (281)

Here, we will specify Q,ϕ and a part of Z ik and solve the constraint equations for the conformal factor ψ and the rest of Z ik . From Equation (275) and our ansatz for the scale invariant variables we obtain:

i Z ik =Q k ϕ . (282)

In the vacuum case, this equation simply becomes the conditions that Z ik is divergence-free, which is in turn simply an algebraic condition on the Fourier coefficients of Z ik . Note that since Σ αβ must be trace-free, so must Z ik . A simple, but still fairly general divergence-free and trace-free Z ik is the following:

Z ik =( b 2 κ 0 κ a 1 cosx+ b 1 a 2 cosx 0 a 2 cosx b 1 b 2 a 1 cosx ) , (283)

where κ, a 1 , a 2 , b 1 and b 2 are constants. We still keep this divergence-free part of Z ik but now add to it a piece that has a non-zero divergence. We simply specify the Fourier coefficients of ϕ and Q via:

Q( x,t=0 )= f 1 H cos( m 1 x+ d 1 ) , ϕ( x,t=0 )= f 2 cos( m 2 x+ d 2 ) (284)

where f 1 , m 1 , d 1 , f 2 , m 2 and d 2 are constants. This turns Equation (282) into an algebraic equation for the Fourier coefficients of this non-zero divergence piece of Z ik which we then solve. Now imposing Equation (276) our ansatz yields:

i i ψ=( 3 4 H 2 1 4 V ) ψ 5 1 8 ( i ϕ i ϕ )ψ 1 8 ( Q 2 + Z ik Z ik ) H 2 ψ 7 , (285)

which is solved for the conformal factor ψ using the numerical methods.

With regard Equation (285), the number 8 is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

The constraint Equations (273) and (274) are automatically satisfied by this ansatz. We then satisfy Equation (277) by using the given value of ϕ to compute the initial value of S α .

Now we show results from a single example that demonstrates the generic behaviour: evolution from a highly inhomogeneous, anisotropic universe with significant curvature at the initial time to a universe containing distinct volumes of either smooth, homogeneous w>1 matter dominated regions, or w=1 mixmaster-like regions. Whenever a w1 region forms it grows exponentially fast in proper volume relative to w=1 regions. The particular initial conditions for this example are (283)-(284):

a 1 =0.70, a 2 =0.10,κ=0.01, b 1 =1.80, b 2 =0.15, f 1 =2.00 ,

m 1 =1, d 1 =1.7, f 2 =0.15, m 2 =2, d 2 =1.0 , and V 0 =0.1,c=10 (286)

for the scalar field potential parameters (262).

With regard the values of (286), we take the following: 0.10 1.15 0.70 1.70 1.80 and 10. We have the following mathematical connections with the Aurea ratio:

Φ 35/7 + Φ 49/7 + Φ 63/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 49/7 + ( 5 +1 2 ) 63/7 =0.090170+0.034442+0.013156=0.137767;0.137767 3 4 =0.103326

Φ 28/7 + Φ 42/7 + Φ 56/7 + Φ 84/7 = ( 5 +1 2 ) 28/7 + ( 5 +1 2 ) 42/7 + ( 5 +1 2 ) 56/7 + ( 5 +1 2 ) 84/7 =0.145898+0.055728+0.021286+0.003106=0.226018 2 3 =0.150679

Φ 7/7 + Φ 21/7 + Φ 35/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 35/7 =0.618034+0.236068+0.090170=0.944272 3 4 =0.708204

Φ 7/7 + Φ 7/7 + Φ 42/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 42/7 =1.618034+0.618034+0.055728=2.291796 3 4 =1.718847

Φ 7/7 + Φ 21/7 + Φ 49/7 + Φ 84/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 49/7 + ( 5 +1 2 ) 84/7 =0.618034+0.236068+0.034442+0.003106 =0.8916492=1.7832980.902=1.80

Φ 35/7 + Φ 14/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 14/7 =11.090170+2.618034=13.708204 3 4 =10.281153.

Now, let χ 1 ( n ) be a complex character to the modulus 5 such that χ 1 ( 2 )=i , and let:

κ= 102 5 2 5 1 =0.284078227 . (287)

The function:

f( s )= 1iκ 2 L( s, χ 1 )+ 1+iκ 2 L( s, χ ¯ 1 ) ,

where L( s,χ )= n=1 χ( n ) n s , is called the Davenport-Heilbronn function and satisfies the Riemann-type equation:

( π 5 ) s/2 Γ( ( s+1 ) 2 )f( s )=g( 1s ) . (288)

We note that 10.281153κ10 . Furthermore:

Φ 14/7 = ( 5 +1 2 ) 14/7 =0.381966 ; 0.381966 3 4 =0.286475 . Thence, we can write also:

κ= 102 5 2 5 1 =0.2840782270.2841 ( 5 +1 2 ) 14/7 3 4 0.2865

It is enlightening to visualize the evolution via the behaviour of the matter ( Ω m ) , shear ( Ω s ) and curvature ( Ω k ) contributions to the normalized energy density, defined as:

Ω m 1 6 W 2 + 1 6 S α S α + 1 3 V ¯

Ω s 1 6 Σ αβ Σ αβ (289)

Ω k 2 3 α A α + A α A α + 1 6 N αβ N αβ 1 12 ( N γ γ γ ) 2 , (290)

where Ω m + Ω s + Ω k =1 by (276).

We note that Equation (290), i.e. 12 = 24/2, is related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

The effective equation of state parameter w takes the following form in Hubble normalized variables:

w= 1 2 W 2 + 1 2 S α S α V ¯ 1 2 W 2 + 1 2 S α S α + V ¯ . (291)

It is evident that at late times the region that has smoothed out and become matter dominated coincides with w1 , whereas the mixmaster-like regime evolves to w=1 . We can calculate the behaviour of the solution in the asymptotic matter dominated region as follows. At late times, all spatial derivatives have become negligible. The constraint (276) then reduces to:

W 2 +2 V ¯ 6 10 , (292)

and slicing condition for N (265) becomes:

3N 3 3 V ¯ , (293)

Furthermore, V ¯ is finite and non-zero. This implies from (262), (265), (266) that asymptotic form:

ϕ( x,t ) ϕ 0 ( x )+ 2t c (294)

and thus W (263) tends to:

W 2 cN . (295)

Combining these relations gives:

Wc , (296)

V ¯ 3 c 2 2 , (297)

N 2 c 2 , (298)

and from (291):

w c 2 /3 1 . (299)

We have that c= 2 =1.414213562 ; Wc and V ¯ =2 . Thence, we have the following mathematical connections with the Aurea ratio:

Φ 7/7 + Φ 21/7 + Φ 35/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 35/7 =0.618034+0.236068+0.090170=0.944272 3 2 =1.416408;

Φ 14/7 + Φ 14/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 14/7 =2.618034+0.381966=3; 2 3 3=2;

Φ 21/7 + Φ 7/7 + Φ 28/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 28/7 =4.236068+1.618034+0.145898=6; 1 3 6=2.

Let S denote the proper spatial volume element associated with the spatial metric h ij of t = const. slices, i.e. S= deth . The fractional change of S with respect to time is:

t lnS= 1 2 h ij t h ij , (300)

which can be written as:

t lnS=3N . (301)

In the asymptotic regime where spatial gradients are negligible, N approaches a constant (293), and thus (301) can be integrated to give:

S m e 6t/ c 2 ,w1 , S v e t ,w=1 (302)

where we have used (293) where w1 , and note that V ¯ 0 when w=1 . Thus, at late times the ratio of the proper volume of matter to mixmaster-like regions of the universe grows as:

= S m dx S v dx e t( 16/ c 2 ) . (303)

Thus, as long as c> 6 (which is equivalent to w>1 ), R as t . We note that 6 =2.449489743 , is related to the following mathematical connection with the Aurea ratio:

Φ 7/7 + Φ 28/7 + Φ 42/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 28/7 + ( 5 +1 2 ) 42/7 =0.618034+0.145898+0.055728=0.8196603=2.458980

(For the various analyzed equations, see Ref. [9] [10]).

6. On Some Equations Concerning the Approximate Inflationary Solutions Rolling Away from the Unstable Maximum of p-Adic String Theory

The action of p-adic string theory is given by:

S= m s 4 g p 2 d 4 x ( 1 2 ϕ p 2 m s 2 ϕ+ 1 p+1 ϕ p+1 ) m s 4 g p 2 d 4 x ( 1 2 ϕ e m p 2 ϕ+ 1 p+1 ϕ p+1 ) (304)

where = t 2 + 2 in the flat space and we have defined:

1 g p 2 1 g s 2 p 2 p1 and m p 2 2 m s 2 lnp . (305)

The dimensionless scalar field ϕ( x ) describes the open string tachyon, m s is the string mass scale and g s is the open string coupling constant. Though the action (304) was originally derived for p a prime number, it appears that it can be continued to any positive integer and even makes sense in the limit p1 . Setting =0 in the action, the resulting potential takes the form:

V= m s 4 g p 2 ( 1 2 ϕ 2 1 p+1 ϕ p+1 ) . (306)

The action (304) is a simplified model of the bosonic string which only qualitatively reproduces some aspects of a more realistic theory. That being said, there are several nontrivial similarities between p-adic string theory and the full string theory.

The field equation that results from (304) is:

e / m p 2 ϕ= ϕ p (307)

We are interested in perturbing around the solution ϕ=1 , which is a critical point of the potential, representing the unstable tachyonic maximum.

One may wonder whether the field theory (304) naively allows for slow roll inflation in the conventional sense. Naively one might expect that for a slowly rolling field the higher powers of in the kinetic term are irrelevant and one may approximate (304) by a local field theory. The action (304) can be rewritten as:

S= d 4 x [ 1 2 χχV( χ )+ ] (308)

where we have defined the field χ as:

χ= χ 0 ϕ

χ 0 = p g s lnp 2( p1 ) m s (309)

and the potential is:

V( χ )= m s 2 lnp χ 2 m s 4 g s 2 p 2 p 2 1 ( χ χ 0 ) p+1 . (310)

In (308), the denotes terms with higher powers of . Thence, Equation (308) can be rewritten also:

S= d 4 x [ 1 2 p g s lnp 2( p1 ) m s ϕ p g s lnp 2( p1 ) m s ϕ m s 2 lnp χ 2 m s 4 g s 2 p 2 p 2 1 ( χ χ 0 ) p+1 + ]. (311)

Working in the context of the action (308) let us consider the slow roll parameters describing the flatness of the potential (310) about the unstable maximum χ= χ 0 . It is straightforward to show that:

M p 2 2 1 V ( χ 0 ) 2 ( V( χ ) χ ) 2 | χ= χ 0 =0 (312)

M p 2 1 V( χ 0 ) 2 V( χ ) χ 2 | χ= χ 0 = 4 g s 2 lnp p 2 1 p 2 ( M p m s ) 2 (313)

With regard the approximate solution for the classical background, we must solve the Friedmann equation:

H 2 = 1 3 M p 2 ρ ϕ (314)

to second order in u . To find the energy density ρ ϕ , we turn to the stress energy tensor for the p-adic scalar field. A convenient expression for T μν is:

T μν = m s 4 2 g p 2 g μν [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ ( α e τ m p 2 ϕ )( α e ( 1τ ) m p 2 ϕ ) ] m s 4 m p 2 g p 2 0 1 dτ ( μ e τ m p 2 ϕ )( ν e ( 1τ ) m p 2 ϕ ). (315)

One may verify that the T μν is symmetric by changing the dummy integration variable τ1τ in the last term. For homogeneous ϕ( t ) , the above expression simplifies, and for T 00 , we find:

ρ ϕ = T 00 = m s 4 2 g p 2 [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ t ( e τ m p 2 ϕ t ( e ( 1τ ) m p 2 ϕ ) ]. (316)

One can evaluate the above expression term by term, keeping up to O( e 2λt ) u 2 . The final result reads:

T 00 = m s 4 2 g p 2 [ 1u( 1+ e μ 1 ) 2[ 1( p+1 )u ] p+1 +u( e μ 1 1 ) ]+O( u 2 ) = m s 4 ( p1 ) 2 g p 2 ( p+1 ) +O( u 2 ). (317)

The O( u ) terms cancel out and matching the coefficients in the Friedmann equation gives us the simple results:

H 0 2 = m s 4 6 M p 2 p1 g p 2 ( p+1 ) (318)

and

H 1 =0 (319)

for zeroth and first order respectively. The O( u 2 ) contribution to T 00 is quite complicated but once we use (319) it simplifies greatly. Matching coefficient at order O( u 2 ) in the Friedmann equation gives:

H 2 = λ m s 4 4 g p 2 m p 2 M p 2 e μ 1 = 1 8 g s 2 p 3 lnp p1 ( m s M p ) 2 λ . (320)

We note that the number 8 in Equation (320) is a Fibonacci’s number and can be connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log [ 10+11 2 4 + 10+7 2 4 .

Because of our sign convention for H r , the fact that H 2 >0 means that the expansion is slowing as ϕ rolls from the unstable maximum, as one would expect in a conventional inflationary model.

We are approximating the background dynamics as de Sitter which amounts to working in the limit u0 , so that:

H 2 H 0 2 = m s 4 6 M p 2 p1 g p 2 ( p+1 ) (321)

ϕ ϕ 0 1 . (322)

We expand the p-adic tachyon field in perturbation theory as:

ϕ( t, x )= ϕ ( 0 ) ( t )+δϕ( t, x )=1+δϕ( t, x ) . (323)

The perturbed Klein-Gordon Equation (307) takes the form:

e / m p 2 δϕ=pδϕ . (324)

One can construct solutions by taking δϕ to be an eigenfunction of the operator. If we choose δϕ to satisfy:

δϕ=+Bδϕ (325)

then this is also a solution to (324) if:

B= m p 2 lnp=2 m s 2 (326)

where in the second equality we have used (306).

For fields which are on-shell (that is, when (325) is solved) the field obeys:

( 1 e / m p 2 )δϕ=( 1 e B/ m p 2 )δϕ=( 1 e B/ m p 2 ) 1 B ( B )δϕ =( 1 e B/ m p 2 ) 1 B δϕ= p1 2 m s 2 δϕ. (327)

Thus, for on-shell fields the kinetic term in the Lagrangian can be written as:

on-shell = m s 4 g p 2 1 2 ϕ( 1 e / m p 2 )ϕ+= m s 4 g p 2 p1 2 m s 2 1 2 ϕϕ+= 1 2 φφ+ (328)

In (328), we have defined the “canonical” field:

φAϕ (329)

where

A m s p 2 g s . (330)

The field φ has a canonical kinetic term in the action, at least while (325) is satisfied. Now, let us return to the task of solving (325), bearing in mind that δφ=Aδϕ is the appropriate canonically normalized field. We write the quantum mechanical solution in term of annihilation/creation operators as:

δφ( t, x )= 1 ( 2π ) 3/2 d 3 k [ a k φ k ( t ) e ikx +h.c. ] (331)

and the mode functions φ k ( t ) are given by:

φ k ( t )= 1 2 π a 3 H 0 e iπ 2 ( v+1/2 ) H v ( 1 ) ( k a H 0 ) (332)

where the order of the Hankel functions is:

v= 9 4 + B H 0 2 = 9 4 + 2 m s 2 H 0 2 (333)

and of course a= e H 0 t . In the second equality in (333) we have used (326) and (306). In writing (331) we have used the usual Bunch-Davies vacuum normalization so that on small scales, ka H 0 , one has:

| φ k | a 1 2k

which reproduces the standard Minkowski space fluctuations. This is the usual procedure in cosmological perturbation theory. On large scales, ka H 0 , the solutions (331) behave as:

| φ k | H 2 k 3 ( k a H 0 ) 3/2 v

which gives a large-scale power spectrum for the fluctuations:

P δφ = ( H 0 2π ) 2 ( k a H 0 ) n s 1

with spectral index:

n s 1=32v

From (333) it is clear that to get an almost scale-invariant spectrum we require m s H 0 . In this limit, we have:

n s 1 4 3 ( m s H 0 ) 2 (334)

which gives a red tilt to the spectrum, in agreement with the latest WMAP data. For n s 0.95 one has m s 0.2 H 0 . Comparing (331) to the corresponding solution in a local field theory we see that the p-adic tachyon field fluctuations evolve as though the mass-squared of the field was 2 m s 2 which may be quite different from the mass scale which one would infer by truncating the infinite series of derivatives: 2 V/ χ 2 ( χ= χ 0 ) .

We note that for n s =0.95 , from Equation (334), we obtain that 2 m s 2 H 0 2 =0.075 . Thence, we have the following mathematical connections with the Aurea ratio:

0.95 Φ 7/7 + Φ 35/7 = ( 5 +1 2 ) 7/7 + ( 5 +1 2 ) 35/7 =0.618034+0.90170=0.708204 4 3 =0.944272

From Equation (333), we have that:

v= 9 4 + 2 m s 2 H 0 2 = 9 4 +0.075 = 2.25+0.075 = 2.325 =1.524795068 ;

1.52480 Φ 0 + Φ 28/7 = ( 5 +1 2 ) 0 + ( 5 +1 2 ) 28/7 =1+0.145898=1.145898 4 3 =1.527864

We note that for Equations (287)-(288)-(289), we have that:

κ 100 = 102 5 2 5 1 =0.00284078227;1.5278640.002840=1.525024

We now want to fix the parameters of the model by comparing to the observed features of the CMB perturbation spectrum. There are three dimensionless parameters, g s ,p and the ratio m s / M p . The important question is whether there is a sensible parameter range which can account for CMB observations, i.e. the spectral tilt and the amplitude of fluctuations. Using (318) in (334), we can relate the tilt to the model parameters via:

| n s 1 |= 8( p+1 ) p 2 ( M p m s ) 2 g s 2 ( m s M p ) 2 = 8( p+1 ) p 2 g s 2 | n s 1 | . (335)

Also, for this equation, we note that the number 8 is a Fibonacci’s number and is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:

8= 1 3 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

Thus, one can have a small tilt while ensuring that the string scale is smaller than the Planck scale, provided that g s 2 /p 1 . Henceforth we will use (335) to determine m s / M p in terms of p, g s , and | n s 1 |0.05 . All the dimensionless parameters in our solution, m s / H 0 ,λ/ m s , ϕ 2 and H 2 / m s , are likewise functions of n s 1,p and g s . From (318) and (335) we see that for p1 ,

m s H 0 6 g p M p m s g s 6 p M p m s 1 2 3| n s 1 | . (336)

It may seem strange to have H exceeding m s since that means the energy density exceeds the fundamental scale, but this is an inevitable property of the p-adic tachyon at its maximum, as shown in Equation (317). This is similar to other attempts to get tachyonic or brane-antibrane inflation from string theory, since the false vacuum energy is just the brane tension which goes like m s 4 / g s .

Next, we determine λ/ m s , where λ is the mass scale appearing in the power series in e λt which provides the ansatz for the background solutions. We consider the following equation for λ in the H 0 m s limit:

e μ 1 = e ( λ 2 +3 H 0 )/ m p 2 =p

The positive root for λ gives:

λ m s n s 1 3 . (337)

In order to fix the amplitude of the density perturbations we consider the curvature perturbation ζ . We assume that:

ζ H φ ˙ δφ

as in conventional inflation models. To evaluate the prefactor H/ φ ˙ we must work beyond zeroth order in the small u expansion. We take ϕ=1u to evaluate the prefactor, even though the perturbation δχ is computed in the limit that ϕ=1 . This should reproduce the full answer up to O( u ) corrections. The prefactor is:

H φ ˙ H 0 Aλu 2 3/2 g s p 1 | n s 1 | 1 u m s 1

With regard Equations (336) and (337), we have the following mathematical connections with the Aurea ratio:

m s H 0 6 g p M p m s g s 6 p M p m s 1 2 3| n s 1 | =0.193649167 ;

0.19365 Φ 28/7 = ( 5 +1 2 ) 28/7 =0.145898 4 3 =0.194531 .

λ m s n s 1 3 =0.129099444 ;

0.12901 Φ 28/7 + Φ 49/7 + Φ 63/7 = ( 5 +1 2 ) 28/7 + ( 5 +1 2 ) 49/7 + ( 5 +1 2 ) 63/7 =0.145898+0.034442+0.013156 =0.193496 2 3 =0.1289970.1290.

We should evaluate u at the time of horizon crossing, t * , defined to be approximately 60 e-foldings before the end of inflation t end , assuming that the energy scale of inflation is high. The inflation ends when u1/ p 1/2 . From Equations (336)-(337), we see that H 0 /λ =2/ | n s 1 | ; therefore, we can write the scale factor a( t ) e H 0 t in the form:

a( t )u ( t ) 2/ | n s 1 | (338)

so that a * = e 60 a end corresponds to:

u * = e 30| n s 1 | u end e 30| n s 1 | 1 p 1/2 (339)

We note that H 0 /λ =2/ | n s 1 | for n s =0.95 is equal to 40. This value can be related with the following mathematical connections with the Aura ratio:

Φ 14/7 + Φ 0 + Φ 35/7 = ( 5 +1 2 ) 14/7 + ( 5 +1 2 ) 0 + ( 5 +1 2 ) 35/7 =2.618034+1+0.090170=3.708204 4 3 =4.944272;

Φ 35/7 + Φ 7/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 7/7 =11.090170+1.618034=12.708204 4 3 =16.944272;

Φ 35/7 + Φ 14/7 = ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 14/7 =11.090170+2.618034=13.708204 4 3 =18.277605;

H 0 /λ =2/ | n s 1 | =40;4.9+16.9+18.2=40 .

The power spectrum of the curvature perturbation is given by:

P ζ = | H φ ˙ | 2 P δφ A ζ 2 ( k a H 0 ) n s 1 (340)

where the amplitude of fluctuations A ζ can now be read off as:

A ζ 2 = 8 3 π 2 g s 2 p e 60| n s 1 | | n s 1 | 3 (341)

Thence, we can rewrite Equation (340) as follows:

P ζ = | H φ ˙ | 2 P δφ 8 3 π 2 g s 2 p e 60| n s 1 | | n s 1 | 3 . (342)

As an example, taking n s 0.95 one can fix the amplitude of the density perturbations A ζ 2 10 10 by choosing:

g s p 0.48× 10 7 . (343)

Setting A ζ 2 = 10 10 and using (341), we obtain an expression for g s in terms of p and | n s 1 | :

g s = 3 π 2 8 p e 30| n s 1 | | n s 1 | 3/2 × 10 5 . (344)

Combining (342) with (335), we also obtain:

m s M p = 3 π 2 p+1 p e 30| n s 1 | | n s 1 |× 10 5 . (345)

The string scale is bounded from above as m s / M p 0.94× 10 6 and that for typical values of p, n s , it is close to m s / M p 0.61× 10 6 . Furthermore, from (342) that g s is unconstrained and that g s , p are not independent parameters.

Now we define the Hubble slow roll parameters ε H , η H by:

ε H 1 2 M p 2 φ ˙ 2 H 2 , ε H η H φ ¨ H φ ˙ . (346)

These are the appropriate parameters to describe the rate of time variation of the inflaton as compared to the Hubble scale. Using the solution ϕ1u (recall that φ=Aϕ , A= m s p/ ( 2 g s ) ) we find that:

ε H 1 2 p+1 p e 60| n1 | | n s 1 | , η H | n s 1 | 2 . (347)

We see that the Hubble slow-roll parameters are small. This means that p-adic tachyon field rolls slowly in the conventional sense. One reaches the same conclusion if one defines the potential slow roll parameters using the correct canonical field, which is φ (329):

M p 2 2 ( 1 V V φ ) 2 | φ=A =0 , M p 2 1 V 2 V φ 2 | φ=A = 1 2 | n 1| . (348)

With regard Equations (347) and (348), we have that:

1 2 | n s 1 |=0.025 ;

Φ 21/7 + Φ 35/7 + Φ 70/7 = ( 5 +1 2 ) 21/7 + ( 5 +1 2 ) 35/7 + ( 5 +1 2 ) 70/7 =0.236068+0.090170+0.008131=0.334369 3 4 =0.250776;

1 25 ×0.250776=0.0250776

On the other hand, consider the potential slow roll parameter which one would naively define using the derivative truncated action (308):

M p 2 2 ( 1 V V χ ) 2 | χ= χ 0 =0 , M p 2 1 V 2 V χ 2 | χ= χ 0 =( p1 lnp ) 1 2 | n s 1 | (349)

where in (349), we have used Equations (313) and (335). We see that (349) can be enormous, though the tachyon field rolls slowly. Taking the largest allowed value of p , p 10 14 , and n s 0.95 , we have M p 2 V 1 | 2 V/ χ 2 | 10 11 . Since large values of p are required if one wants to obtain g s 1 , it follows that it is somewhat natural for p-adic inflation to operate in the regime where the higher derivative corrections play an important role in the dynamics.

(For the various analyzed equations, see Ref. [11]).

7. On Some Equations Concerning p-Adic Minisuperspace Model, Zeta Strings, Zeta Nonlocal Scalar Fields and p-Adic and Adelic Quantum Cosmology

Consider the standard Minkowski signature minisuperspace model of a homogeneous isotropic universe with a cosmological constant λ . The usual parametrization of the metric:

d s 2 = N 2 d t 2 + a 2 d Ω 3 2 (350)

leads to classical solutions which are trigonometrical functions of time. In the p-adic case, we prefer to work with rational functions. We shall use the following ansatz:

d s 2 = N ( t ) 2 q( t ) d t 2 +q( t )d Ω 3 2 (351)

Here, N and a are functions of time and d Ω 3 2 is the metric on the unit 3-sphere. The action for this metric is the same as the corresponding usual case:

S[ q( t ) ]= 1 2 t 1 t 2 dt N( q ˙ 2 4 N 2 λ 2 q1 ) (352)

We assume that the cosmological constant λ is a rational number. The classical equations of motions have the form:

q ¨ =2λ . (353)

The solution of this equation for the boundary conditions:

q( 0 )= q 1 ,q( T )= q 2 , (354)

is the following:

q( t )=λ t 2 +[ q 2 q 1 T λT ]t+ q 1 . (355)

Here, q( t ),p( t ) Q p . The Green function corresponding to the transition from the point q 1 to the point q 2 has the form:

G p ( q 1 , q 2 )= Q p dT K T ( q 1 ,0| q 2 ,T ) (356)

where K T ( q 1 ,0| q 2 ,T ) is the propagator:

K T ( q 1 ,0| q 2 ,T )= χ p ( S ) t dq( t ) . (357)

In the path, integral one integrates over trajectories with the boundary conditions (354). One can perform the Gaussian path integral (357) in the usual way using shifting to the classical solution. One gets:

K T ( q 1 ,0| q 2 ,T )=c( T ) χ p ( S cl ) (358)

where S cl is the action calculated on the trajectories (355).

S cl = S cl ( q 2 , q 1 ,T )= λ 2 T 3 24 +[ λ( q 1 + q 2 )2 ] T 4 + ( q 1 q 2 ) 2 8T . (359)

The factor c( T ) is the same as for a free particle:

c( t )= λ p ( T ) | T | p 1/2 . (360)

Therefore, one has the Green function:

G p ( q 1 , q 2 )= Q p dT λ p ( T ) | T | p 1/2 χ p ( λ 2 T 3 24 +[ λ( q 1 + q 2 )2 ] T 4 + ( q 1 q 2 ) 2 8T ) . (361)

The corresponding wave function has the form:

Ψ( q )= Q p dT λ p ( T ) | T | p 1/2 χ p ( 1 l Pl 2 [ q 2 8T +( λq2 ) T 4 λ 2 T 3 24 ] ) (362)

where we restore the explicit dependence on the Planck length.

We note that the number 24 in Equation (362) can be related to the physical vibrations of the bosonic strings by the following Ramanujan function:

24= 4[ antilog 0 cosπtx w coshπx e π x 2 w dx e π t 2 4 w ϕ w ( it w ) ] 142 t 2 w log[ 10+11 2 4 + 10+7 2 4 ] .

Now let us estimate the integral on T applying the stationary phase approximation. The saddlepoint equation has the form:

S = λ 2 T 2 8 q 2 8 T 2 + λq 4 1 2 =0 (363)

which yields:

T 1,2 2 = ( 1± 1λq ) 2 λ 2 . (364)

As is known, for p1( mod4 ) , there is the square root of −1 in Q p , so we get nontrivial saddle points. For p3( mod4 ) , we have no saddle point at all. To make sense of the saddle points in the case p1( mod4 ) , we should be sure that the square root 1λq also has a sense. For this purpose, we have to assume that | λq | p <1 . The corresponding actions have the form:

S= 1 3λ [ 1± ( 1λq ) 3/2 ] , S= 1 3λ [ 1± ( 1λq ) 3/2 ] , (365)

In order that these expressions are rational, we have to assume that λ is rational as well as that λq= ξ 2 is such a rational that the solution of the equation:

ξ 2 + η 2 =1 (366)

in respect to η is also rational. Let us consider the Euclidean metric:

d s 2 = N ( t ) 2 q( t ) d t 2 +q( t )d Ω 3 2 . (367)

The Euclidean action for this metric is the same as the corresponding action in the usual case:

S[ q( t ) ]= 1 2 t 1 t 2 dt N( q ˙ 2 4 N 2 + λ 2 q1 ) . (368)

We shall prove that it is possible to restore:

Ψ E ( q )= dT exp( 2πi l Pl 2 S E )dq (369)

when l Pl 0 in the corresponding p-adic wave function. Indeed, in the p-adic case for Euclidean metric, we get a basic Green function:

G p ( q 1 , q 2 )= dT λ p ( T ) | T | p 1/2 χ p ( λ 2 T 3 24 +[ λ( q 1 q 2 )2 ] T 4 ( q 1 q 2 ) 2 8T ) . (370)

Now let us estimate the integral on T in:

Q p dT λ p ( T ) | T | 1/2 χ p ( λ 2 T 3 24 +( λq2 ) T 4 q 2 8T ) (371)

applying the stationary phase approximation. Formally, there are the following saddle points:

T= 1 λ [ 1± ( 1λq ) 1/2 ] , T= 1 λ [ 1± ( 1λq ) 1/2 ] , (372)

for

| λq | p <1

with corresponding actions:

S= 1 3λ [ 1± ( 1λq ) 3/2 ] , S= 1 3λ [ 1± ( 1λq ) 3/2 ] . (373)

Note that now these stationary points have sense for all p and for q satisfying (192) according our general formula for the wave function of the universe:

Ψ( G 3 )= am p χ p ( S p )D ( g μν ) p , (374)

we write:

Ψ E ( q )= p=2,3,5, λ p | λ | p ( 1λ q p ( ± 1 3λ ) )| [ ( 1λ q 1 ) 3/2 ±1 ]

= | ( 1λq ) 1/2 ±1 | | λ | exp{ ± 2πi 3λ [ ( 1λ q 1 ) 3/2 ±1 ] }. (375)

7.1. Zeta Strings and Zeta Nonlocal Scalar Fields

The exact tree-level Lagrangian for effective scalar field φ which describes open p-adic string tachyon is:

p = 1 g 2 p 2 p1 [ 1 2 φ p 2 φ+ 1 p+1 φ p+1 ] , (376)

where p is any prime number, = t 2 + 2 is the D-dimensional d’Alambertian and we adopt metric with signature ( ++ ) . Now, we want to show a model which incorporates the p-adic string Lagrangians in a restricted adelic way. Let us take the following Lagrangian:

. (377)

Recall that the Riemann zeta function is defined as:

ζ( s )= n1 1 n s = p 1 1 p s , s=σ+iτ , σ>1 . (378)

Employing usual expansion for the logarithmic function and definition (378) we can rewrite (377) in the form:

L= 1 g 2 [ 1 2 ϕζ( 2 )ϕ+ϕ+ln( 1ϕ ) ] , (379)

where | ϕ |<1 . ζ( 2 ) acts as pseudodifferential operator in the following way:

ζ( 2 )ϕ( x )= 1 ( 2π ) D e ikk ζ( k 2 2 ) ϕ ˜ ( k )dk , k 2 = k 0 2 k 2 >2+ε , (380)

where ϕ ˜ ( k )= e ikx ϕ( x )dx is the Fourier transform of ϕ( x ) .

Dynamics of this field ϕ is encoded in the (pseudo)differential form of the Riemann zeta function. When the d’Alambertian is an argument of the Riemann zeta function we shall call such string a “zeta string”. Consequently, the above ϕ is an open scalar zeta string. The equation of motion for the zeta string ϕ is:

ζ( 2 )ϕ= 1 ( 2π ) D k 0 2 k ˜ 2 >2+ε e ixk ζ( k 2 2 ) ϕ ˜ ( k )dk = ϕ 1ϕ (381)

which has an evident solution ϕ=0 .

For the case of time dependent spatially homogeneous solutions, we have the following equation of motion:

ζ( t 2 2 )ϕ( t )= 1 2π | k 0 |> 2 +ε e i k 0 t ζ( k 0 2 2 ) ϕ ˜ ( k 0 )d k 0 = ϕ( t ) 1ϕ( t ) . (382)

With regard the open and closed scalar zeta strings, the equations of motion are:

ζ( 2 )ϕ= 1 ( 2π ) D e ikk ζ( k 2 2 ) ϕ ˜ ( k )dk = n1 θ n( n1 ) 2 ϕ n , (383)

ζ( 4 )θ= 1 ( 2π ) D e ikk ζ( k 2 4 ) θ ˜ ( k )dk = n1 [ θ n 2 + n( n1 ) 2( n+1 ) θ n( n1 ) 2 1 ( ϕ n+1 1 ) ], (384)

and one can easily see trivial solution ϕ=θ=0 .

The exact tree-level Lagrangian of effective scalar field φ , which describes open p-adic string tachyon, is:

p = m p D g p 2 p 2 p1 [ 1 2 φ p 2 m p 2 φ+ 1 p+1 φ p+1 ] , (385)

where p is any prime number, = t 2 + 2 is the D-dimensional d’Alambertian and we adopt metric with signature ( ++ ) , as above. Now, we want to introduce a model which incorporates all the above string Lagrangians (385) with p replaced by nN . Thence, we take the sum of all Lagrangians n in the form:

L= n=1 + C n n = n=1 + C n m n D g n 2 n 2 n1 [ 1 2 ϕ n 2 m n 2 ϕ+ 1 n+1 ϕ n+1 ] , (386)

whose explicit realization depends on particular choice of coefficients C n , masses m n and coupling constants g n .

Now, we consider the following case:

C n = n1 n 2+h , (387)

where h is a real number. The corresponding Lagrangian reads:

L h = m D g 2 [ 1 2 ϕ n=1 + n 2 m 2 h ϕ+ n=1 + n h n+1 ϕ n+1 ] (388)

and it depends on parameter h . According to the Euler product formula, one can write:

. (389)

Recall that standard definition of the Riemann zeta function is:

ζ( s )= n=1 + 1 n s = p 1 1 p s ,s=σ+iτ,σ>1 , (390)

which has analytic continuation to the entire complex s plane, excluding the point s=1 , where it has a simple pole with residue 1. Employing definition (390), we can rewrite (388) in the form:

L h = m D g 2 [ 1 2 ϕζ( 2 m 2 +h )ϕ+ n=1 + n h n+1 ϕ n+1 ]. (391)

Here, ζ( 2 m 2 +h ) acts as a pseudodifferential operator:

ζ( 2 m 2 +h )ϕ( x )= 1 ( 2π ) D e ixk ζ( k 2 2 m 2 +h ) ϕ ˜ ( k )dk , (392)

where ϕ ˜ ( k )= e ikx ϕ( x )dx is the Fourier transform of ϕ( x ) .

We consider Lagrangian (391) with analytic continuations of the zeta function and the power series n h n+1 ϕ n+1 , i.e.

L h = m D g 2 [ 1 2 ϕζ( 2 m 2 +h )ϕ+AC n=1 + n h n+1 ϕ n+1 ] , (393)

where AC denotes analytic continuation.

Potential of the above zeta scalar field (393) is equal to L h at =0 , i.e.

V h ( ϕ )= m D g 2 ( ϕ 2 2 ζ( h )AC n=1 + n h n+1 ϕ n+1 ) , (394)

where h1 since ζ( 1 )= . The term with ζ -function vanishes at h=2,4,6, . The equation of motion in differential and integral form is:

ζ( 2 m 2 +h )ϕ=AC n=1 + n h φ n , (395)

1 ( 2π ) D R D e ixk ζ( k 2 2 m 2 +h ) ϕ ˜ ( k )dk =AC n=1 + n h ϕ n , (396)

respectively.

Now, we consider five values of h , which seem to be the most interesting, regarding the Lagrangian (393): h=0 , h=±1 , and h=±2 . For h=2 , the corresponding equation of motion now read:

ζ( 2 m 2 2 )ϕ= 1 ( 2π ) D R D e ikk ζ( k 2 2 m 2 2 ) ϕ ˜ ( k )dk = ϕ( ϕ+1 ) ( 1ϕ ) 3 . (397)

This equation has two trivial solutions: ϕ( x )=0 and ϕ( x )=1 . Solution ϕ( x )=1 can be also shown taking ϕ ˜ ( k )=δ( k ) ( 2π ) D and ζ( 2 )=0 in (397).

For h=1 , the corresponding equation of motion is:

ζ( 2 m 2 1 )ϕ= 1 ( 2π ) D R D e ixk ζ( k 2 2 m 2 1 ) ϕ ˜ ( k )dk = ϕ ( 1ϕ ) 2 . (398)

where ζ( 1 )= 1 12 .

The equation of motion (398) has a constant trivial solution only for ϕ( x )=0 .

For h=0 , the equation of motion is:

ζ( 2 m 2 )ϕ= 1 ( 2π ) D R D e ikk ζ( k 2 2 m 2 ) ϕ ˜ ( k )dk = ϕ 1ϕ . (399)

It has two solutions: ϕ=0 and ϕ=3 . The solution ϕ=3 follows from the Taylor expansion of the Riemann zeta function operator:

ζ( 2 m 2 )=ζ( 0 )+ n1 ζ ( n ) ( 0 ) n! ( 2 m 2 ) n , (400)

as well as from ϕ ˜ ( k )= ( 2π ) D 3δ( k ) .

For h=1 , the equation of motion is:

1 ( 2π ) D R D e ikk ζ( k 2 2 m 2 +1 ) ϕ ˜ ( k )dk = 1 2 ln ( 1ϕ ) 2 , (401)

where ζ( 1 )= gives V 1 ( ϕ )= .

In conclusion, for h=2 , we have the following equation of motion:

1 ( 2π ) D R D e ikk ζ( k 2 2 m 2 +2 ) ϕ ˜ ( k )dk = 0 ϕ ln ( 1w ) 2 2w dw . (402)

Since holds equality:

0 1 ln( 1w ) w dw = n=1 1 n 2 =ζ( 2 )

one has trivial solution ϕ=1 in (402).

Now, we want to analyze the following case: C n = n 2 1 n 2 . In this case, from the Lagrangian (386), we obtain:

L= m D g 2 [ 1 2 ϕ{ ζ( 2 m 2 1 )+ζ( 2 m 2 ) }ϕ+ ϕ 2 1ϕ ]. (403)

The corresponding potential is:

V( ϕ )= m D g 317ϕ 24( 1ϕ ) ϕ 2 . (404)

The equation of motion is:

[ ζ( 2 m 2 1 )+ζ( 2 m 2 ) ]ϕ= ϕ[ ( ϕ1 ) 2 +1 ] ( ϕ1 ) 2 . (405)

Its weak field approximation is:

[ ζ( 2 m 2 1 )+ζ( 2 m 2 )2 ]ϕ=0 , (406)

which implies condition on the mass spectrum:

ζ( M 2 2 m 2 1 )+ζ( M 2 2 m 2 )=2 . (407)

From (407) it follows one solution for M 2 >0 at M 2 2.79 m 2 and many tachyon solutions when M 2 <38 m 2 .

With regard the extension by ordinary Lagrangian, we have the Lagrangian, potential, equation of motion and mass spectrum condition that, when C n = n 2 1 n 2 , are:

L= m D g 2 [ ϕ 2 { m 2 ζ( 2 m 2 1 )ζ( 2 m 2 )1 }ϕ+ ϕ 2 2 ln ϕ 2 + ϕ 2 1ϕ ] , (408)

V( ϕ )= m D g 2 ϕ 2 2 [ ζ( 1 )+ζ( 0 )+1ln ϕ 2 1 1ϕ ] , (409)

[ ζ( 2 m 2 1 )+ζ( 2 m 2 ) m 2 +1 ]ϕ=ϕln ϕ 2 +ϕ+ 2ϕ ϕ 2 ( 1ϕ ) 2 , (410)

ζ( M 2 2 m 2 1 )+ζ( M 2 2 m 2 )= M 2 m 2 . (411)

In addition to many tachyon solutions, Equation (411) has two solutions with positive mass: M 2 2.67 m 2 and M 2 4.66 m 2 .

Now, we describe the case of C n =μ( n ) n1 n 2 . Here μ( n ) is the Mobius function, which is defined for all positive integers and has values 1,0,1 depending on factorization of n into prime numbers p . It is defined as follows:

μ( n )={ 0, n= p 2 m ( 1 ) k , n= p 1 p 2 p k , p i p j 1, n=1,( k=0 ) (412)

The corresponding Lagrangian is:

L μ = C 0 0 + m D g 2 [ 1 2 ϕ n=1 + μ( n ) 2 m 2 ϕ+ n=1 + μ( n ) n+1 ϕ n+1 ] (413)

Recall that the inverse Riemann zeta function can be defined by:

1 ζ( s ) = n=1 + μ( n ) n s ,s=σ+it,σ>1 . (414)

Now, (413) can be rewritten as:

L μ = C 0 0 + m D g 2 [ 1 2 ϕ 1 ζ( 2 m 2 ) ϕ+ 0 ( ϕ )dϕ ] , (415)

where M( ϕ )= n=1 + μ( n ) ϕ n =ϕ ϕ 2 ϕ 3 ϕ 5 + ϕ 6 ϕ 7 + ϕ 10 ϕ 11 The corresponding potential, equation of motion and mass spectrum formula, respectively, are:

V μ ( ϕ )= L μ ( =0 )= m D g 2 [ C 0 2 ϕ 2 ( 1ln ϕ 2 ) ϕ 2 0 ϕ ( ϕ )dϕ ] , (416)

1 ζ( 2 m 2 ) ϕ( ϕ ) C 0 m 2 ϕ2 C 0 ϕlnϕ=0 , (417)

1 ζ( M 2 2 m 2 ) C 0 M 2 m 2 +2 C 0 1=0 , | ϕ |1 , (418)

where usual relativistic kinematic relation k 2 = k 0 2 + k 2 = M 2 is used.

Now, we take the pure numbers concerning Equations (407) and (411). They are: 2.79, 2.67 and 4.66. We note that all the numbers are related with Φ= 5 +1 2 , thence with the aurea ratio, by the following expressions:

2.79 Φ 15/7 ;2.67 Φ 13/7 + Φ 21/7 ;4.66 Φ 22/7 + Φ 30/7 . (419)

7.2. p-Adic and Adelic Quantum Cosmology

Adelic quantum cosmology is an application of adelic quantum theory to the universe as a whole. In the path integral approach to standard quantum cosmology starting point is Feynman’s idea that the amplitude to go from one state with intrinsic metric h ij , and matter configuration ϕ on an initial hypersurface Σ , to another state with metric h ij , and matter configuration ϕ on a final hypersurface Σ , is given by a functional integral of χ ( S | g μν ,Φ ] ) over all four-geometries g μν , and matter configurations Φ, which interpolate between the initial and final configurations, i.e.

h ij , ϕ , Σ | h ij ,ϕ, Σ = D ( g μν ) D ( Φ ) χ ( S [ g μν ,Φ ] ) . (420)

The S [ g μν ,Φ ] is the usual Einstein-Hilbert action:

S[ g μv ,Φ ]= 1 16πG ( M d 4 x g ( R2Λ )+2 M d 3 x h K ) 1 2 M d 4 x g [ g μv μ Φ ν Φ+V( Φ ) ] (421)

for the gravitational field and matter fields Φ. In (421), R is scalar curvature of four-manifold M , Λ is cosmological constant, K is trace of the extrinsic curvature K ij at the boundary M of the manifold M . To perform p-adic and adelic generalization, we first make p-adic counterpart of the action (421) using form-invariance under change of real to the p-adic number fields. Then, we generalize (420) and introduce p-adic complex-valued cosmological amplitude:

h ij , ϕ , Σ | h ij ,ϕ, Σ p = D ( g μν ) p D ( Φ ) p χ p ( S p [ g μν ,Φ ] ) . (422)

The space of all 3-metrics and matter field configurations ( h ij ( x ),ϕ( x ) ) on a 3-surface is called superspace (this is the configuration space in quantum cosmology). Superspace is the infinite dimensional one with a finite number of coordinates ( h ij ( x ),ϕ( x ) ) at each point x of the 3-surface. One useful approximation is to truncate the infinite degrees of freedom to a finite number, thereby obtaining some particular minisuperspace model. Usually, one restricts the four-metric to be of the form d s 2 = N 2 ( t )d t 2 + h ij d x i d x j , where N( t ) is the laps function. For such minisuperspaces, functional integrals (420) and (422) are reduced to functional integration over three-metrics, matter configurations and to one usual integral over the laps function. If one takes boundary condition q α ( t 2 )= q 2 α , q α ( t 1 )= q 1 α then integral in (420) and (422), in the gauge N=0 , is a minisuperspace propagator. In this case, it holds:

q 2 α | q 1 α v = dN K v ( q 2 α ,N| q 1 α ,0 ) , (423)

where

K v ( q 2 α ,N| q 1 α ,0 )= D q α χ v ( S v [ q α ] ) (424)

is an ordinary quantum-mechanical propagator between fixed q α in fixed time N . For quadratic classical action S p cl ( q 2 ,N| q 1 ,0 ) , (424) becomes:

K p ( q 2 ,N| q 1 ,0 )= λ p ( 2 S p cl 2 q 2 q 1 ) | 2 S p cl q 2 q 1 | p 1/2 χ p ( S p cl ( q 2 ,N| q 1 ,0 ) ) . (425)

If system has n decoupled degrees of freedom, its p-adic kernel is a product:

K p ( q 2 ,N| q 1 ,0 )= α=1 n λ p ( 2 S p cl 2 q 2 α q 1 α ) | 2 S p cl q 2 α q 1 α | p 1/2 χ p ( S p cl ( q 2 α ,N| q 1 α ,0 ) ) . (426)

p-adic and adelic wave functions of the universe may be found by means of the following equation:

U( t ) ψ αβ ( x )=χ( E α t )ψ ψ αβ ( x ) , (427)

where ψ αβ ( x ) are adelic wave eigenfunctions, E=( E , E 2 ,, E p , ) is the corresponding adelic energy, α=( α , α 2 ,, α p , ) and β=( β , β 2 ,, β p , ) are indicies for energy levels and their degeneration, respectively.

The corresponding adelic eigenstates have the form:

Ψ( q α )= ψ ( q α ) pS ψ p ( q p α ) pS Ω( | q p α | p ) . (428)

A necessary condition to construct an adelic model is existence of the p-adic (vacuum) state Ω( q α | p ) , which satisfies:

| q 1 α | p 1 K p ( q 2 α ,N| q 1 α ,0 )d q 1 α =Ω( | q 2 α | p ) (429)

for all but a finite number of p .

Now, we describe the p-adic and adelic model with cosmological constant in D=3 dimensions. This model has the metric:

d s 2 = σ 2 ( N 2 ( t )d t 2 + a 2 ( t )( d θ 2 + sin 2 θd φ 2 ) ) , (430)

where σ=G . The corresponding v-adic action is:

S v [ a ]= 1 2 0 1 dtN a 2 ( t )( a ˙ 2 N 2 a 2 + 1 a 2 λ ) , (431)

where λ=Λ σ 2 . The Euler-Lagrange equation of motion:

a ¨ N 2 aλ=0

has the solution:

a( t )= 1 2sinh( N λ ) ( ( a 2 a 1 e N λ ) e N λ t +( a 1 e N λ a 2 ) e N λ t ) , (432)

where the boundary conditions are a( 0 )= a 1 , a( 1 )= a 2 . For the classical action, it gives:

S v cl ( a 2 ,N| a 1 ,0 )= 1 2 λ [ N λ +λ( 2 a 1 a 2 sinh( N λ ) a 1 2 + a 2 2 tanh( N λ ) ) ] . (433)

Quantum-mechanical propagator has the form:

K v ( a 2 ,N| a 1 ,0 )= λ v ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | v 1/2 χ v ( S v cl ( a 2 ,N| a 1 ,0 ) ). (434)

Equation (429), in a more explicit form, reads:

Ω( | a 2 | p )= λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 χ p ( N 2 + λ 2tanh( N λ ) a 2 2 ) × | a 1 | p 1 χ p ( λ 2tanh( N λ ) a 1 2 λ sinh( N λ ) a 2 a 1 )d a 1 . (435)

We note that the p-adic Gauss integral over the region of integration | x | p p v is:

| x | p p v χ p ( α x 2 +βx )dx = p v Ω( p v | β | p ), | α | p p 2v 1; (436)

| x | p p v χ p ( α x 2 +βx )dx = λ p ( α ) | 2α | p 1/2 χ p ( β 2 4α )Ω( p v | β 2α | p ), | α | p p 2v >1; (437)

where Ω( u ) is defined as follows:

Ω( u )=1,u1 ; Ω( u )=0,u>1 .

Using (437), for  v=0 , we obtain:

Ω( | a 2 | p )= χ p ( N 2 + λ 2 tanh( N λ ) a 2 2 )Ω( | a 2 | p ) (438)

with condition | λ 2tanh( N λ ) | p >1 , i.e. | N | p <1 . For p2 , left hand side is equal to Ω( | a 2 | p ) if | λ | p 1 holds. Applying also the (436) to (435), we have:

Ω( | a 2 | p )= λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 × χ p ( N 2 + λ a 2 2 2coth( N λ ) )Ω( | λ a 2 sinh( N λ ) | p ). (439)

It becomes an equality if condition | N | p 1 take place.

Thence, we can rewrite Equation (435) as follow:

Ω( | a 2 | p )= λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 × χ p ( N 2 + λ 2tanh( N λ ) a 2 2 ) × | a 1 | p 1 χ p ( λ 2tanh( N λ ) a 1 2 λ sinh( N λ ) a 2 a 1 )d a 1 = λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 × χ p ( N 2 + λ a 2 2 2coth( N λ ) )Ω( | λ a 2 sinh ( N λ ) p ). (440)

The de Sitter minisuperspace model in quantum cosmology is the simplest, nontrivial and exactly soluble model. This model is given by the Einstein-Hilbert action with cosmological term (421) without matter fields, and by Robertson-Walker metric:

d s 2 = σ 2 ( N 2 ( t )d t 2 + a 2 ( t )d Ω 3 2 ) , (441)

where σ 2 = 2G 3π and a( t ) is the scale factor. Instead of (441), we shall use:

d s 2 = σ 2 ( N 2 ( t ) q( t ) d t 2 +q( t )d Ω 3 2 ) . (442)

The corresponding v-adic action for this one-dimensional minisuperspace model is:

S v [ q ]= 1 2 t 1 t 2 dt N( q ˙ 2 4 N 2 λq+1 ) , (443)

where λ= Λ σ 2 3 . The classical equation of motion ( N=1 ) q ¨ =2λ , with the boundary conditions q( 0 )= q 1 and q( T )= q 2 ,( T= t 2 t 1 ) gives:

q( t )=λ t 2 +( q 2 q 1 T λT )t+ q 1 . (444)

After substitution (444) into (443) and integration, one obtains that the classical action is:

S v cl ( q 2 ,T| q 1 ,0 )= λ 2 T 3 24 [ λ( q 1 + q 2 )2 ] T 4 ( q 2 q 1 ) 2 8T . (445)

Since (445) is quadratic in q 2 and q 1 , quantum-mechanical propagator has the form:

K v ( q 2 ,T| q 1 ,0 )= λ v ( 8T ) | 4T | v 1/2 χ v ( S v cl ( q 2 ,T| q 1 ,0 ) ) . (446)

Equation (429) reads:

Ω( | q 2 | p )= λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T ) × | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1 . (447)

(For the various analyzed equations, see Ref. [12]-[16]).

8. Mathematical Connections

Now, we describe some possible mathematical connections. We take Equation (31) of Section 2. We note that can be related with Equations (315), (316) of Section 6, hence we have the following connections:

N e = dϕ V V ,ϕ e c ϕ c c 2 m s 4 2 g p 2 g μν [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ ( α e τ m p 2 ϕ )( α e ( 1τ ) m p 2 ϕ ) ] m s 4 m p 2 g p 2 0 1 dτ ( μ e τ m p 2 ϕ )( ν e ( 1τ ) m p 2 ϕ );

N e = dϕ V V ,ϕ e c ϕ c c 2 m s 4 2 g p 2 [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ t ( e τ m p 2 ϕ ) t ( e ( 1τ ) m p 2 ϕ ) ]. (448)

Thence, there is a mathematical connections between the slow-roll formula regarding the number of e-foldings N e of inflation and the equations of the stress energy tensor for the p-adic scalar field in padic inflation.

Now, we take Equations (48), (50), (117) and (119) of Section 3. We note that can be related with Equation (447) of Section 6, hence we obtain the following connections:

S= dt d 3 x [ 3 M 5 3 L( a ˙ 2 K a 2 )+ 1 2 a 2 ϕ ˙ 2 ]+ S m

λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T )× | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1 ; (449)

S= d 4 x g ( M 4 2 2 R 1 2 ( μ ϕ ) 2 )+ S m [ g ]+ S m + [ g + ]

λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T )× | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1 ; (450)

δ d m L = 1 1 n x 4 ξ 40 ( x 4 ) b 2 e 1 2 x 4 2 dω = 2 x 4 ξ 40 ( x 4 ) ( 1 x 4 2 ) 2 = 1 x 4 2 1 [ 8 3 k ˜ x 4 ( A 0 J 0 ( k ˜ x 4 )+ B 0 Y 0 ( k ˜ x 4 ) )4( A 0 J 1 ( k ˜ x 4 )+ B 0 Y 1 ( k ˜ x 4 ) ) ]

λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T )× | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1 ; (451)

δ d m L = y 0 y 0 nt Γ L dy = 1 1 n x 4 ξ 4 b 2 e 1 2 x 4 2 dω

λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T )× | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1 . (452)

Thence, mathematical connections between some equations concerning cosmological perturbations in a Big Crunch/Big Bang space-time and M-theory model of a Big Crunch/Big Bang transition (48)-(50), some equations concerning the solution of a braneworld Big Crunch/Big Bang cosmology (117)-(119) and the equation concerning the de Sitter minisuperspace model in p-adic quantum cosmology (447).

Now, we take Equations (315), (316), (334), (341) and (349) of Section 6. We note that can be related with Equations (151), (153), (155), (213), (223) and (225) of Section 4, hence we obtain the following mathematical connections:

τ 2 1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα γ 0 ( 1 2 m D 0 )

m s 4 2 g p 2 g μν [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m ρ 2 ϕ ) + 1 m p 2 0 1 dτ ( α e τ m ρ 2 ϕ )( α e ( 1τ ) m ρ 2 ϕ ) ]

m s 4 m p 2 g p 2 0 1 dτ ( μ e τ m p 2 ϕ )( ν e ( 1τ ) m p 2 ϕ )

m s 4 2 g p 2 [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ t ( e τ m p 2 ϕ ) t ( e ( 1τ ) m p 2 ϕ ) ]; (453)

n s 1 4 3 ( m s H 0 ) 2 2 ε dlnε dN 0.03 β 1+β

1+ dlog | δ k | 2 dlogk 1+ 4 m ×1/ [ 2 m log( m 2 α 2k v e mC 2 ) ]

1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) ; (454)

A ζ 2 = 8 3 π 2 g s 2 p e 60| n s 1 | | n s 1 | 3 2 ε dlnε dN 4( 1+ γ 2 ) c 2 M Pl 2 4 c ϕ c 2 0.03 β 1+β

1+ dlog | δ k | 2 dlogk 1+ 4 m ×1/ [ 2 m log( m 2 α 2k v e mC 2 ) ]

1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) ; (455)

M p 2 1 V 2 V χ 2 | χ= χ 0 =( p1 lnp ) 1 2 | n s 1 | 2 ε dlnε dN 4( 1+ γ 2 ) c 2 M Pl 2 4 c ,ϕ c 2 0.03 β 1+β

1+ dlog | δ k | 2 dlogk 1+ 4 m ×1/ [ 2 m log( m 2 α 2k v e mC 2 ) ]

1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) . (456)

Thence, there are mathematical connections between some equations concerning the generating ekpyrotic curvature perturbations before the Big Bang, some equations concerning the colliding branes and the origin of the hot Big Bang and some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-adic string theory.

Now, we take Equations (306), (315) and (316) of Section 6 and Equations (381), (396) and (415), of Section 7. We obtain the following mathematical connections:

S= m s 4 g p 2 d 4 x ( 1 2 ϕ p 2 m s 2 ϕ+ 1 p+1 ϕ p+1 ) m s 4 g p 2 d 4 x ( 1 2 ϕ e m p 2 ϕ+ 1 p+1 ϕ p+1 )

m s 4 2 g p 2 g μν [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ ( α e τ m p 2 ϕ )( α e ( 1τ ) m p 2 ϕ ) ] m s 4 m p 2 g p 2 0 1 dτ ( μ e τ m p 2 ϕ )( ν e ( 1τ ) m p 2 ϕ )

m s 4 2 g p 2 [ ϕ e m p 2 ϕ 2 p+1 ϕ p+1 + 1 m p 2 0 1 dτ ( e τ m p 2 ϕ )( e ( 1τ ) m p 2 ϕ ) + 1 m p 2 0 1 dτ t ( e τ m p 2 ϕ ) t ( e ( 1τ ) m p 2 ϕ ) ]

ζ( 2 )ϕ= 1 ( 2π ) D k 0 2 k 2 >2+ε e ikk ζ( k 2 2 ) ϕ ˜ ( k )dk = ϕ 1ϕ

1 ( 2π ) D R D e ikk ζ( k 2 2 m 2 +h ) ϕ ˜ ( k )dk =AC n=1 + n h ϕ n

C 0 0 + m D g 2 [ 1 2 ϕ 1 ζ( 2 m 2 ) ϕ+ 0 ( ϕ )dϕ ] . (457)

Thence, mathematical connections between some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-adic string theory and some equations concerning the zeta strings and the zeta nonlocal scalar fields.

In conclusion, with regard Section 7, we have the following mathematical connections between Equations (440) and (447) and Equation (213) of Section 4:

Ω( | a 2 | p )= λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 × χ p ( N 2 + λ 2tanh( N λ ) a 2 2 )

× | a 1 | p 1 χ p ( λ 2tanh( N λ ) a 1 2 λ sinh( N λ ) a 2 a 1 )d a 1 = λ p ( λ 2sinh( N λ ) ) | λ sinh( N λ ) | p 1/2 × χ p ( N 2 + λ a 2 2 2coth( N λ ) )Ω( | λ a 2 sinh( N λ ) | p )

1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) ;

Ω( | q 2 | p )= λ p ( 8T ) | 4T | p 1/2 χ p ( λ 2 T 3 24 T 2 + λ q 2 T 4 + q 2 2 8T ) × | q 1 | p 1 χ p [ q 1 2 8T +( λT 4 q 2 4T ) q 1 ]d q 1

1 2v [ 0 Y 0 D( Y ) e mα Y /2 d Y ] 2 2 D 0 2 m 2 α 2 v e mα Y 0 ( 1 2 m D 0 ) . (458)

9. Conclusion

In this work, we have described mathematical connections between the slow-roll formula regarding the number of e-foldings N e of inflation and the equations of the stress-energy tensor for the p-adic scalar field in padic inflation. We have described mathematical connections between some equations concerning cosmological perturbations in a Big Crunch/Big Bang space-time and M-theory model of a Big Crunch/Big Bang transition (48)-(50), some equations concerning the solution of a braneworld Big Crunch/Big Bang cosmology (117)-(119) and the equation concerning the de Sitter minisuperspace model in p-adic quantum cosmology (447). We have described mathematical connections between some equations concerning the generating ekpyrotic curvature perturbations before the Big Bang, some equations concerning the colliding branes and the origin of the hot Big Bang and some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-adic string theory. We have described mathematical connections between some equations concerning the approximate inflationary solutions rolling away from the unstable maximum of p-adic string theory and some equations concerning the zeta strings and the zeta nonlocal scalar fields. In conclusion, with regard to Section 8, we have the following mathematical connections between Equations (440) and (447) and Equation (213) of Section 4. Thence, there are mathematical connections between some equations concerning the p-adic quantum cosmology and the fundamental equation concerning the colliding branes and the origin of the hot Big Bang.

Acknowledgements

I would like to thank Paul J. Steinhardt of the Department of Physics of Princeton University for his availability and the very useful references and advice. Furthermore, I would like to thank also Branko Dragovich of the Institute of Physics of Belgrade (Serbia) for his availability and friendship with regard to me and Gianmassimo Tasinato of Heidelberg University for his availability and the useful reference concerning the p-adic Inflation.

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Steinhardt, P.J. and Turok, N. (2002) Cosmic Evolution in a Cyclic Universe. Physical Review D, 65, Article ID: 126003.
https://doi.org/10.1103/physrevd.65.126003
[2] Tolley, A.J., Turok, N. and Steinhardt, P.J. (2004) Cosmological Perturbations in a Big-Crunch/Big-Bang Space-Time. Physical Review D, 69, Article ID: 106005.
https://doi.org/10.1103/physrevd.69.106005
[3] Steinhardt, P.J. and Turok, N. (2004) The Cyclic Model Simplified. arXiv: astro-ph/0404480.
[4] Turok, N., Perry, M. and Steinhardt, P.J. (2004) M Theory Model of a Big Crunch/Big Bang Transition. Physical Review D, 70, Article ID: 106004.
https://doi.org/10.1103/physrevd.70.106004
[5] McFadden, P.L., Turok, N. and Steinhardt, P.J. (2007) Solution of a Braneworld Big Crunch/Big Bang Cosmology. Physical Review D, 76, Article ID: 104038.
https://doi.org/10.1103/physrevd.76.104038
[6] Lehners, J., McFadden, P., Turok, N. and Steinhardt, P.J. (2007) Generating Ekpyrotic Curvature Perturbations before the Big Bang. Physical Review D, 76, Article ID: 103501.
https://doi.org/10.1103/physrevd.76.103501
[7] Lukas, A., Ovrut, B.A., Stelle, K.S. and Waldram, D. (1999) Universe as a Domain Wall. Physical Review D, 59, Article ID: 086001.
https://doi.org/10.1103/physrevd.59.086001
[8] Khoury, J., Ovrut, B.A., Steinhardt, P.J. and Turok, N. (2001) Ekpyrotic Universe: Colliding Branes and the Origin of the Hot Big Bang. Physical Review D, 64, Article ID: 123522.
https://doi.org/10.1103/physrevd.64.123522
[9] Steinhardt, P.J. and Wesley, D. (2008) Dark Energy, Inflation and Extra Dimensions. arXiv: 0811.1614.
[10] Garfinkle, D., Lim, W.C., Pretorius, F. and Steinhardt, P.J. (2008) Evolution to a Smooth Universe in an Ekpyrotic Contracting Phase with w > 1. arXiv: 0808.0542.
[11] Barnaby, N., Biswas, T. and Cline, J.M. (2007) p-Adic Inflation. Journal of High Energy Physics, 4, Article No. 56.
https://doi.org/10.1088/1126-6708/2007/04/056
[12] Aref’eva, I.Ya., Dragovich, B., Frampton, P.H. and Volovich, I.V. (1990) Wave Function of the Universe and p-Adic Gravity.
https://lib-extopc.kek.jp/preprints/PDF/1990/9007/9007373.pdf
[13] Dragovich, B. (2007) Zeta Strings. arXiv: hep-th/0703008.
[14] Dragovich, B. (2008) Zeta Nonlocal Scalar Fields. Theoretical and Mathematical Physics, 157, 1671-1677.
[15] Dragovich, B. and Nesic, L. (2000) p-Adic and Adelic Generalization of Quantum Cosmology. arXiv: gr-qc/0005103.
[16] Lange, C., Nardelli, M. and Bini, G. (2008) Sistema Musicale Aureo Phi(n/7) e connessioni matematiche tra Numeri Primi e “Paesaggio” della Teoria delle Stringhe.
http://eprints.bice.rm.cnr.it/679/1/NardLanBin02.pdf

Copyright © 2025 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.