On Some Mathematical Connections between the Cyclic Universe, Inflationary Universe, p-Adic Inflation, p-Adic Cosmology and Various Sectors of Number Theory ()
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
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
, 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
in a homogeneous, flat Universe, with line element is:
. (1)
We use
to represent conformal time and primes to represent derivatives with respect to
.
is the lapse function. The background solution for the scalar field is denoted
, and
is the scalar potential.
The equations of motion for gravity, the matter and scalar field
are straightforwardly derived by varying (1) with respect to
and
, after which
may be set equal to unity. Expressed in terms of proper time
, the Einstein equations are:
, (2)
, (3)
where a dot is a proper time derivative.
With regard the trajectory in the
-plane, the Friedmann constraint reads:
, (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
becomes negligible as
runs off to minus infinity. The Friedmann constraint reads
, and the scalar field equation,
, where primes denote conformal time derivatives. The general solution is:
,
, (5)
where
. We choose
to be the time when
vanishes, so that
before collision.
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
and
which at
take the values
(in) and
(in) respectively. Re-expressing the scalar field as a function of proper time
, we obtain:
. (6)
The integration constant
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,
, (7)
where
is the four-dimensional line element,
is the fifth spatial coordinate which runs from zero to
, and
is a parameter with the dimension of length. Thence, we have that:
, (8)
where
is the proper length of the extra dimension,
is a parameter with dimensions of length, and
is the proper time in the five-dimensional metric,
, (9)
with
being FRW proper time. Notice that a shift
can always be compensated for by a rescaling of
. As the extra dimension shrinks to zero,
tends to a constant,
.
Immediately after the bounce, scalar kinetic energy dominates and
remains nearly constant. The kinetic energy of the scalar field scales as
and radiation scales as
, so the former dominates at small
. It is convenient to re-scale
so that it is unity at scalar kinetic energy-radiation equality,
, and denote the corresponding Hubble constant
. The Friedmann constraint in Equation (4) then reads:
, (10)
and the solution is:
. (11)
The brane scale factors are:
,
. (12)
Here, the constant
has been defined so that we match
and
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
on the positive and negative tension branes respectively. For
, the case of relatively little radiation production, immediately after collision
is expanding but
is contracting. Whereas for
, both brane scale factors expand after collision. If no scalar potential
were present, the scalar field would continue to obey the solution (11), converging to:
. (13)
This value is actually larger than
for
, the case of weak production of radiation. However, the presence of the potential
alters the expression (13) for the final resting value of the scalar field. As
crosses the potential well travelling in the positive direction,
is reduced to a renormalized value , so that the final resting value of the scalar field can be smaller than
. If this is the case, then
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 . The dependence of (13) is simply understood: while the Universe is kinetic energy dominated,
grows at
and
increases logarithmically with time. However, when the Universe becomes radiation dominated and
, Hubble damping increases and
converges to the finite limit above.
With regard Equations (8)-(13), we note the following connections with number theory:
Note that,
, 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:
, with
that is the Aurea ratio,
is a natural number and 7 are the compactified dimensions of the M-Theory.
Using the following potential:
, (14)
we consider the motion of
back and forth across the potential well.
may be accurately approximated by
. For this pure exponential potential, there is a simple scaling solution:
, (15)
which is an expanding or contracting Universe solution according to whether
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
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
and
, where
and
are the background quantities given from (15). To linear order in
, one obtains:
, (16)
with two linearly independent solutions,
and
, where
. In the contracting phase, the former solution grows as
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
and
. Within the scaling solution (15), we can calculate the value of incoming velocity by treating the prefactor of the potential
in Equation (14) as a Heaviside function which is unity for
and zero for
, where
is the value of
at the minimum of the potential. We compute the velocity of the field as it approaches
and use energy conservation at the jump in
to infer the velocity after
is crossed. In the scaling solution, the total energy as
approaches
from the right is
, and this must equal the total energy
evaluated for
just to the left of
. Hence, we find that at the minimum and, according to Equation (8),
. (17)
Note that from Equation (17), we obtain:
where the number 8 is connected with the “modes” that correspond to the physical vibrations of a superstring by the following Ramanujan function:
At the bounce, this solution is matched to an expanding solution with:
, (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
positive. We can treat
as a perturbation and use the solution in Equation (16) discussed above,
and
. One can straightforwardly compute the perturbation in
in this growing mode by matching at
as before. One finds
where
is the background value, at the minimum. Beyond this point,
grows as
for large
, whereas in the background scaling solution
decays with
as
. The departure occurs when the scalar field has attained the value,
. (19)
As
passes beyond
the kinetic energy overwhelms the negative potential and the field passes onto the plateau
with
nearly constant and equal to:
, (20)
until the radiation, matter and vacuum energy become significant and
is then damped away to zero. Note that we can rewrite Equation (20) as follow:
(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:
and with the number 12 (12 = 24/2) that is related to the physical vibrations of the bosonic strings by the following Ramanujan function:
The time spent to the left of the potential well
is essentially identical in the incoming and outgoing stages for
, namely:
. (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
,
. (23)
Also, this equation can be rewritten as follow:
. (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
:
. (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
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
. The value of
is nearly cancelled in the passage across the potential well, and is reduced to
given in Equation (20). Once radiation domination begins, the field quickly converges to the large
(Hubble-damped) limit of Equation (11), namely:
, (26)
where
is the Hubble radius at kinetic-radiation equality. Also, Equation (26) can be rewritten as follow:
. (27)
The dependence is obvious: the asymptotic value of
depends on the ratio of to
. Increasing pushes
further, likewise lowering
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
and matching,
, (28)
where as above, we define
to be unity at kinetic-radiation equal density. We have that
may reach its maximal value
and turn around during the radiation, matter or quintessence dominated epoch. For example,
is reached in the radiation era, if, from Equation (28),
, (29)
where
is the time of matter domination.
For turn around in the matter era, we require:
. (30)
Finally, if the field runs to very large
, so that
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
of inflation is given by usual slow-roll formula,
, (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
, or
. This is easily fulfilled if
is of order unity Planck units. Hence, Equation (31) can be rewritten as follow:
. (32)
With regard Equations (29)-(30) and (32), we have the following mathematical connections with the Aurea ratio:
From the formulae given above, we can also calculate the maximal value
in the cyclic solution: for large
and for
, it is:
, (33)
where we used
, the beginning of the radiation-dominated epoch. From Equation (33), we obtain:
. (34)
This equation provides a lower bound on
. The extreme case is to take
. Then, using
,
,
, we find
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
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
, there is a unique value of
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
, in terms of
. The initial state of empty branes with tensions
and
, and with corresponding velocities
and
obeys:
(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
and that the total momentum is:
, (36)
where we identify
with the contraction speed of the fifth dimension,
. For Equation (17), we can rewrite Equation (36) also as follow:
. (37)
The corresponding equations for the outgoing state are easily obtained, by replacing
with
for the positive tension brane, and
with
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
, where
are the associated rapidities, one obtains two solutions:
. (38)
where
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
direction but the negative tension brane is also moving in the negative
direction. The corresponding values for
(out) and
are:
(39)
where the first solution for
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,
, compared to their relative speed
in the incoming state. We find in the (−+) solution,
(40)
and in the (−−) solution,
. (41)
with
given by (36) in both cases. We note that we can rewrite the relation mentioned above, i.e.
also as follow:
. (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,
, 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
is increased, we expect
and
to grow. Now, if we consider
and
to be both
, then the second term in the denominator dominates. If more radiation is produced on the negative tension brane,
, then:
(43)
and so
is small and positive. This is the condition necessary to obtain cyclic behaviour. Conceivably, the brane tension can change from
to
at collision. Then, we obtain:
(44)
For the (−+) solution, we can straightforwardly determine an upper limit for
. Consider, for example, the case there the brane tension in unchanged at collision,
. The expression in (40) gives
as a function of
and
. It is greatest, at fixed
and
, when
, its smallest value. For
, it is maximized for
, and equal to
when equality holds. For
, it is maximized when
, its smallest value, and
, when it is equal to
. 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
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
, we have the following mathematical connections with the Aurea ratio:
(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,
, (45)
where
is the density (not including the tension) of matter or radiation confined to the brane,
is the brane scale factor, and
is the induced Hubble constant on the positive (negative) tension brane. Choosing conformal time on each brane, and neglecting the
terms, Equations (45) become:
. (46)
where prime denotes conformal time derivative. The corresponding acceleration equations for
and
, from which
disappears, are derived by differentiating Equations (46) and using
with
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,
, (47)
where
is a lapse function introduced to make the action time reparameterization invariant. Varying with respect to
and then setting
gives the correct acceleration equations for
and
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
and a scalar field
, defined by:
Clearly,
and
transform as a scale factor and as a scalar field under rescalings of the spatial coordinates
. To interpret
more physically, note that for static branes the bulk space-time is perfect Anti-de Sitter space with line element . The separation between the branes is given by:
so
tends from zero to infinity as
tends from minus infinity to zero. In terms of
and
, the action (47) becomes:
, (48)
which is recognized as the action for Einstein gravity with line element
,
being the canonical metric on
or
with curvature
, and a minimally coupled scalar field
. The matter action
is conventional, except that the scale factor appearing is not the Einstein-frame scale factor but instead
and
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:
,
. (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:
, (50)
where we have defined the effective four-dimensional Planck mass
. 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
, and a massless scalar field with kinetic energy density
. The four-dimensional Friedmann equation in conformal time then reads:
, (51)
where we have defined the constants
and
, and used the fact that the massless scalar kinetic energy
. The solution to (51) and the massless scalar field equation
is:
,
, (52)
From these solutions, we reconstruct the scale factors on the branes according to (49), obtaining:
, (53)
so we see that with the choice of normalization for the scale factor
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:
,
, (54)
where
is the rapidity associated with the relative velocity of the branes at collision
and
is the value of the radiation density
on each brane at collision. Thence, Equation (53) can be rewritten also:
. (55)
Furthermore, we define the fractional density mismatch on the two branes as:
, (56)
so that we have:
. (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:
. (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:
, (59)
where primes denote
derivatives,
is the fractional perturbation in the radiation density,
is the scalar potential for its velocity, i.e.
,
is the perturbation in the scalar field, and from (52), we have the background quantities:
We are interested in solving these equations in the long wavelength limit,
. Solving all the above equations for
, one finds:
, (60)
for adiabatic perturbations. The components of the background energy density in the four-dimensional effective theory are scalar kinetic energy, with
, and radiation, with
. It follows that for adiabatic perturbations, at long wavelengths, we must have:
. (61)
In longitudinal gauge, the fractional energy density perturbation and the velocity potential perturbation in the scalar field (considered as a fluid with
) are given by:
. (62)
From Equations (59) above (and using
), it follows that:
. (63)
Maintaining the adiabaticity condition (61) up to order
then requires that the fractional velocity perturbations for the scalar field and the radiation should be equal:
. Expressing the radiation velocity in terms of
, the momentum constraint then yields:
, (64)
where
.
The above equations may be used to determine the leading terms in an expansion in
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,
, which has the following series expansion about
:
, (65)
where
and
are arbitrary constants, and
(66)
For adiabatic perturbations, we obtain:
, (67)
,
,
,
, (68)
, (69)
, (70)
where
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:
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:
,
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
. The form that we take for the five-dimensional cosmological background metric is:
, (71)
and we write the most general scalar metric perturbation about this as:
(72)
For perturbations on
it is straightforward to find a gauge in which the metric takes the form:
, (73)
and
satisfies a massless scalar equation of motion on
. To be precise, the gauge is:
. (74)
Notice that the non-zero variables can all be related to
according to:
. (75)
We shall, henceforth, refer to these as the “Milne ratio conditions”. Furthermore, imposing the
symmetry, we obtain Neumann boundary conditions on
,
, (76)
where
are the location of the two
fixed points. In the model space-time, the lowest energy mode for
is y-independent and has the asymptotic form:
, (77)
with
and
being arbitrary constants. We have the following relations:
,
. (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,
is suppressed by
compared to
. As a result, we obtain the approximate matching rule:
,
. (79)
The key conditions (74) through (76) are satisfied precisely for all time in a compactified Milne
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
and scalar field perturbation
determined the induced metric perturbations on the branes via Formulae (49):
, (80)
where
and
and the metric perturbations are fractional, i.e.
,
. 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:
(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:
, (82)
are precisely the boundary values of the five-dimensional longitudinal gauge perturbations
and
. Using (80) and (82), we find for the induced perturbations:
,
,
,
. (83)
The brane conformal times may be expressed in terms of
by integrating,
, (84)
where
. So, for example, the boundary value of the bulk metric perturbation
on the positive tension brane is given explicitly by:
(85)
where
is the location of the positive tension brane.
Also, in these Equations (83)-(85), we have the following connection with the Aurea ratio:
.
Now, using (67) and the following equation:
(86)
to find
and
before and after the bounce for all components of the metric perturbations and matching according to the rule given in Equation (79) results in
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:
, (87)
where
and
are the relative velocities of the branes before and after collision. Note that since
, matching
is in fact equivalent to matching
across the collision. In terms of four-dimensional parameters including
given in (57) defining the abundance of the radiation and the fractional density mismatch
defined in (56), we find again at leading order in velocities:
. (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:
(89)
where
is the rapidity corresponding to the relative speed
of the branes at collision, and the second formula assumes
is small.
is the bulk curvature scale, and
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:
;
that are Fibonacci’s numbers;
;
, with 3 and 8 that are Fibonacci’s numbers;
.
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:
The
-dimensional space-time we consider is a direct product of
-dimensional Euclidean space,
, and a two-dimensional time-dependent space-time known as compactified Milne space-time, or
. The line element for
is thus:
, (90)
where
are Euclidean coordinates on
,
parameterizes the compact dimension and
is the time. The compact dimension may either be a circle, in which case we identify
with
, or a
orbifold in which case we identify
with
and further identify
with
. The fixed points
and
are then interpreted as tensionless
-branes approaching at rapidity
, colliding at
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
and the brane collision is locally well-modelled by
.
Now consider a string loop of radius
in M-theory frame. Its mass
is
times the effective string tension
, where
is the size of the extra dimension. The effective Einstein-frame gravitational coupling is given by
. The gravitational potential produced by such a loop in
spacetime dimensions is:
(91)
where
is the area of the unit
-sphere,
. Specializing to the case of interest, namely 2-branes in eleven-dimensional M-theory, the tension
is related to the eleven-dimensional gravitational coupling by a quantization condition relating to the four-form flux, reading:
(92)
with
an integer. Equations (91) and (92) then imply that the typical gravitational potential around a string loop is:
(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:
and 3, 5, 8, 21, 34, 55 and 89 are Fibonacci’s numbers. Furthermore, we have that:
Thence, the gravitational potential on the scale of the loops is of order
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
, this potential Φ is the typical gravitational potential throughout space. Multiplying the
component of the background metric (90) by
and redefining
, we conclude that the outgoing metric has an expansion rapidity of order
with
a constant of order unity. We conclude that for small
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:
, (94)
where
parameterizes the fifth dimension and
, the three non-compact dimensions. Cosmological isotropy excludes
or
terms, and homogeneity ensures
and
are independent of
. The
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
for all time
. 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
, where the bulk cosmological constant is
. Evaluating the linear combinations
and
, we find:
,
, (95)
where
and
. The Israel matching conditions on the branes read:
(96)
where all quantities are to be evaluated at the brane locations
.
Now, we express the metric as a series of Dirichlet or Neumann polynomials in
and
, bounded at order
by a constant times
, such that the series satisfies the Israel matching conditions exactly at every order in
. To implement this, we first change variables from
and
to those obeying Neumann boundary conditions. From (96),
is Neumann. Likewise, if we define
by:
, (97)
then one can easily check that
is also Neumann on the branes. Since
and
obey Neumann boundary conditions on the branes, we can expand both in a power series:
,
, (98)
where
are polynomials:
,
(99)
satisfying Neumann boundary conditions and each bounded by
, for the relevant range of
. Note that the time-dependent coefficients in this ansatz may also be expanded as a power series in
. 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
, one generically obtains a number of linearly independent algebraic equations, and at most one ordinary differential equation in
.
The first few terms of the solution are:
,
and
,
(100)
With regard Equations (100), we note that are possible the following mathematical connections with the Aurea ratio and the Fibonacci’s numbers:
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:
,
.
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:
, (101)
where for Schwarzschild-AdS with a horizon at
,
,
. (102)
At
, the
-coordinate of the branes is represented by the parameter
; their subsequent trajectories
can then be determined by integrating the Israel matching conditions, which read
, where
are the proper speeds of the positive- and negative-tension branes respectively. From this, it further follows that
is related to the rapidity
of the collision by
.
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
and
. Taking the
and
coordinates to be fixed along the geodesic then, we find that
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
, and thus:
, (103)
where we have introduced the constant
. Adopting
now as the affine parameter, we have:
, (104)
where
is to be regarded now as a function of
. We can solve this equation order by order in
using the series ansatz:
, (105)
where the constants
are themselves series in
. Using the series solution for the background geometry given in Equations (100), and imposing the boundary condition that
, we obtain:
(106)
(107)
(108)
(109)
(110)
(111)
(112)
(113)
Substituting
and
, we find
, i.e. to lowest order in
, 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:
, (114)
must necessarily agree. To check this, the affine distance between the branes is given by:
(115)
which to lowest order in
reduces to:
(116)
in agreement with the series expansion of (114).
We obtain also the following equation:
(117)
To evaluate the perturbation
in the affine distance between the branes, consider:
(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
and the constraint
. 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
and
, to lowest order in
the unperturbed geodesics lie purely in the
direction, and so the perturbed affine distance is identical to the following perturbed metric distance:
. (119)
Explicitly, we find:
(120)
which, substituting
and dropping terms of
, reduces to:
(121)
where
. Also, this expression is in accordance with the series expansion of (119). However, the perturbed affine and metric distance do not agree at
.
With regard Equations (116) and (121), we have the following mathematical connections with Aurea ratio and with Fibonacci’s numbers:
Furthermore, we have that:
;
; where 2, 8, 21, 34, 233 and 610 are Fibonacci’s numbers;
;
;
;
; 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
decoupled fields interacting only through gravity:
, (122)
where we have chosen units in which
. In a flat Friedmann-Robertson-Walker background with line element
, the scalar field and Friedmann equations are given by:
(123)
and
, (124)
where
and
with no summation implied. Another useful relation is:
. (125)
If all the fields have negative exponential potentials
then as is well-known, the Einstein-scalar equations admit the scaling solution:
. (126)
Thus, if
for all
, we have a very slowly contracting universe with
.
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:
(127)
in flat spacetime is replaced by:
(128)
For simplicity we will focus attention on straight line trajectories in scalar field space. Since
, 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
. Denoting
derivatives with primes, and introducing the re-scaled entropy field:
, (129)
Equation (128) becomes:
. (130)
The crucial term governing the spectrum of the perturbations is then:
. (131)
When this quantity is approximately 2, we will again get nearly scale-invariant perturbations. It is customary to define the quantity:
. (132)
In the background scaling solution,
. (133)
We proceed by evaluating the quantity in (131) in an expansion in inverse powers of
and its derivatives with respect to
, where
, where
is the value of
at the end of the ekpyrotic phase. Note that
decreases as the fields roll downhill and the contracting ekpyrotic phase proceeds. We obtain the first term in (131) by differentiating (125), obtaining:
. (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
. We obtain:
. (135)
Finally, need to express
in terms of the conformal time
. From (134), we obtain:
(136)
which integrates to:
. (137)
Now, inserting
under the integral and using integration by parts we can re-write this as:
. (138)
Using the same procedure once more, the integral in this expression can be written as:
. (139)
Now using the fact that
, and that to leading order in
,
can be replaced by its value in the scaling solution (with constant
),
, we can re-write the second term on the righthand side as:
, (140)
which shows that this term is of order
and can thus be neglected. Altogether, we obtain:
. (141)
Using (134) and (135) with (141), we can calculate the crucial term entering the entropy perturbation equation,
. (142)
The deviation from scale-invariance in the spectral index of the entropy perturbation is then given by:
. (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
scalar fields with general Kahler metric
on scalar field space, the linearized Einstein-scalar field equations lead to:
, (144)
where the
entropy perturbations:
(145)
are just the components of
orthogonal to the background trajectory, and the operator
is just the geodesic operator on scalar field space. Things simplify because the scalar field space is flat, so the metric is
, and
reduces to an ordinary time derivative. Considering only two scalar fields, we have:
. (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 , for
, and , for
, with
constant and negative in the vicinity of the bounce. The bounce leads to a delta function on the right-hand side of (144),
, (147)
where
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
, 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,
. Since the entropy perturbation:
(148)
is canonically normalized, its spectrum is given by:
(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
approximately given by
. After
, the entropy perturbation obeys
, which has the solution
.
Matching this solution to the growing mode solution
in the ekpyrotic phase, one finds that by
the entropy grows by an additional factor of
. Employing the Friedmann equation to relate to
, putting everything together and restoring the Planck mass, we find for the variance of the spatial curvature perturbation in the scale-invariant case,
(150)
for the perfectly scale-invariant case. Notice that the results depend only logarithmically on
: the main dependence is on the minimum value of the effective potential and the parameter
. Observations on the current Hubble horizon indicate
. Ignoring the logarithm in (150), this requires
, 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:
. If we take
, then we have that
. Furthermore, we obtain:
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
, the number of e-folds before the end of the ekpyrotic phase (where
and
):
. (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,
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
e-folds. If we estimate
, then the spectral tilt is:
(152)
Here, we see that the sign of the tilt is sensitive to
. For nearly exponential potentials
, the spectral tilt is
, slightly blue, because the first term dominates. However, there are well-motivated examples in which the equation of state does not decrease linearly with
. We have introduced
to parameterize these cases. If
, the spectral tilt is red. For example,
for
. These examples represent the range that can be achieved for the entropically induced curvature perturbations in the simplest models, roughly
.
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
. This range agrees with the estimate obtained by an independent analysis based on studying inflaton potentials directly.
With regard the values of
, i.e. 0.95, 0.97 and 1.02, we have the following mathematical connections with the Aurea ratio:
;
;
;
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
and
. Then, Equation (143) becomes:
, (153)
where we have used the fact that
has the dimensions of inverse mass and restored the factors of Planck mass. The presence of
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:
, (154)
although the agreement is not at all obvious at intermediate steps of the calculation. For a pure exponential potential, which has
, the non-gravitational contribution is zero, and the spectrum is slightly blue. For plausible values of
and
, say, the gravitational piece is about one percent and the spectral tilt is
, also consistent with our earlier estimate.
We note that
, is related with the Aurea ratio by the following mathematical formula:
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
. If
changes from some initial value
to some value of order unity at the end of the ekpyrotic phase after
changes by an amount
, then . When
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
and
; then, the spectral tilt is:
, (155)
which corresponds to
for positive
, in agreement with our earlier estimate. With regard the value 125, we have the following mathematical connection with the Fibonacci’s numbers:
We note that negative potentials of this type with very large values of
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:
,
. (156)
Note that
. Then, the spectral tilt is:
. (157)
This result can be compared with the spectral index of the time-delay (Newtonian potential) perturbation, where the corresponding formula is:
. (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:
(159)
where the result is expressed in terms of the slow-roll parameters
and
. Here, we have revealed the factors of
to illustrate that both inflationary contributions are gravitational in origin. This gives the same range for
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:
(160)
where
is the familiar 11-dimensional supergravity:
(161)
and
are the two
Yang-Mills theories on the orbifold planes explicitly given by:
(162)
Here,
are the two
gauge field strengths and
is the 3-form with field strength
. In order for the above theory to be supersymmetric as well as anomaly free, the Bianchi identity for
should receive a correction such that:
(163)
where the sources are given by:
. (164)
With regard Equations (161) and (163), we have the following mathematical connections with the Aurea ratio and the Ramanujan modular equations:
;
;
;
;
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:
While the standard embedding of the spin connection into the gauge connection:
(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
on the two hyperplanes. As a result, the antisymmetric tensor field
and, hence, the second term in the gravitino supersymmetry variation:
(166)
do not vanish.
With regard the mathematical connections with the Aurea ratio and Fibonacci’s numbers, we note that:
;
and 144 is a Fibonacci’s number. Furthermore,
and 24 are related to the physical vibrations of the bosonic strings by the following Ramanujan function:
.
Now, let us start with the zeroth order metric:
(167)
where
is a Calabi-Yau metric with Kahler form
. (Here,
and
are holomorphic and anti-holomorphic indices). To keep track of the scaling properties of the solution, we have introduced moduli
and
for the Calabi-Yau volume and the orbifold radius, respectively. To order
, the metric can be written in the form:
(168)
where the functions
and
depend on
and the Calabi-Yau coordinates. Furthermore,
and
receive a contribution of order
from the Bianchi identity source terms. The general explicit form of the corrections are:
,
(169)
,
(170)
with
,
. (171)
Here,
is the step function which is +1 (−1) for
positive (negative).
With regard Equations (169), (170) and (171), we have the following mathematical connections with the Aurea ratio:
;
;
;
.
In the five-dimensional space
of the reduced theory, the orbifold fixed planes constitute four-dimensional hypersurfaces which we denote by
. There will be an
gauge field
accompanied by gauginos and gauge matter fields on the orbifold plane
. We will set these gauge matter fields to zero in the following. The field content of the orbifold plane
consists of an
gauge field
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:
. (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:
. (173)
The complete configuration for the antisymmetric tensor field that we use in the reduction is given by:
,
,
,
,
,
(174)
and the non-zero mode is:
, (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:
(176)
where
(177)
(178)
(179)
In this expression, we have now dropped higher-derivative terms. The 4-form field strength
is subject to the Bianchi identity:
(180)
which follows directly from the 11-dimensional Bianchi identity (163). The currents
have been defined in Equation (164). The five-dimensional Newton constant
and the Yang-Mills coupling
are expressed in terms of 11-dimensional quantities as:
. (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
supergravity in five dimensions. Accordingly, let us compare the result (177)-(179) to the known
supergravity matter theories in five dimensions. In these theories, the scalar fields in the universal hypermultiplet parameterize a quaternionic manifold with coset structure
. Hence, to compare our action to these we should dualize the three-form
to a scalar field
by setting (in the bulk):
. (182)
Then, the hypermultiplet part of the action (178) can be written as:
(183)
where
. The covariant derivative
is defined as
with
. The sigma model metric
can be computed from the Kahler potential:
. (184)
Consequently, the hypermultiplet scalars
parameterize a Kahler manifold with metric
. It can be demonstrated that
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
theories. It is sufficient to keep the bosonic terms only. Both approaches lead to:
(185)
where
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
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
-brane in a
-dimensional theory, one needs to have a
-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
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:
(186)
where
and
are functions of
and all other field vanish. The general solution for this Ansatz, satisfying the equations of motion derived from the action (176), is given by:
(187)
where
and
are constants. We note that the boundary source terms have fixed the form of the harmonic function
in the above solution. Without specific information about the sources, the function
would generically be glued together from an arbitrary number of linear pieces with slopes
. The edges of each piece would then indicate the location of the source terms. The necessity of matching the boundary sources at
and
, however, has forced us to consider only two such linear pieces, namely
and
. These pieces are glued together at
and
. Therefore, we have:
(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
. It is straightforward to show that the linearized version of (187), that is, the expansion to first order in
, 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
and
are the only non-zero fields, the supersymmetry transformations (185) simplify to:
(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:
;
;
;
;
;
;
;
;
;
.
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:
We note also that 8 is a Fibonacci’s number.
The Killing spinor equations
are satisfied for the solution (187) if we require that the spinor
is given by:
(190)
where
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
, the universal “breathing” modulus of the Calabi-Yau three-fold
, a four-form gauge field
with field strength
and a single bulk M5-brane. It is given by:
(191)
where
. The space-time is a five-dimensional manifold
with coordinates
. The four-dimensional manifolds
are the visible, hidden, and bulk branes respectively, and have internal coordinates
and tension
. Note that
has dimension of mass. If we denote
,
, and
, then the visible brane has tension
, the hidden brane
, and the bulk brane
. It is straightforward to show that the tension of the bulk brane,
, must always be positive. Furthermore, one can easily deduce that the tension on the visible brane,
, can be either positive or negative. We will take
, so that the tension on the visible brane is negative. Furthermore, we will choose
such that
, that is, the tension of the hidden brane is positive. The tensor
is the induced metric on
. The functions
are the coordinates in
of a point on
with coordinates
. In other words,
describe the embedding of the branes into
. The BPS solution of Lukas, Ovrut and Waldram is then given by:
(192)
where
(193)
and
and
are constants. Note that
are dimensionless and
has the dimension of length. The visible and hidden boundary branes are located at
and
, respectively, and the bulk brane is located at
,
. We assume that
, so that the curvature singularity at
does not fall between the boundary branes. Note that
lies in the region of smaller volume while
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
and charge
in the interval
. 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:
, (194)
expresses the density perturbation in terms of the time delay at the time of collision,
. If we consider the exponential potential
, then Equation (194) yields:
. (195)
Now, we compute the spectrum of quantum fluctuations of the brane
and use the result to compute the time delay,
.
For the calculation of quantum fluctuations, it is sufficient to work at the lowest order in
. Without loss of generality, we can therefore set
. In that case, the bulk brane Lagrangian is given by:
. (196)
Note that this agrees with
given in the following equation:
, (197)
when we set
and spatial gradients of
to zero. Let us first consider the spatially homogeneous motion of the brane which will be described by
. It is governed by the following equation of motion:
, (198)
where
is a constant. Equation (198) is simply the statement that the energy
of the bulk brane is conserved to this order in
. Since we have chosen the visible brane to lie at
and the hidden universe to lie at
, we focus on the branch
in which case the bulk brane moves towards the visible brane. The solution to Equation (198) is then given by:
(199)
with
, and with the collision occurring at
. Let us now consider fluctuations around the background solution
. Namely, if
, with
, we can expand the action to quadratic order in
:
(200)
where we have used Equation (198), and where we have introduced
and
for simplicity. The key relation is the fluctuation equation as derived from the action (200):
, (201)
where
and where
is defined by:
. (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
:
. (203)
Defining
, Equation (203) becomes:
. (204)
Let us now discuss the Hubble horizon for the perturbations. Recall that in usual 4 d cosmology (see Equation (204)), we have:
, (205)
where
is the Hubble radius as derived from the scale-factor
. 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
. Therefore, a mode with amplitude
crosses outside the horizon when
. Similarly, in our scenario, we can write:
, (206)
where
. The role of the Hubble radius is replaced by:
, (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
(rather than
), but the amplitude itself, as derived from Equation (201), depends on
. 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 . In particular, one obtains a scale-invariant spectrum if when the modes observed on the CMB cross outside the horizon.
Combining Equations (199) and (202), we find:
. (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:
. (209)
With regard Equation (209), we have the following mathematical connections with the Aurea ratio:
We note also that 2 is prime number and Fibonacci’s number.
Thence we can rewrite Equation (209) also as follow:
We can add a potential
of the form that might result from the exchange of wrapped M2branes. We would like to think of
as the potential derived from the superpotential
for the modulus
in the 4 d low energy theory. Typically, superpotentials for such moduli are of exponential form, for example,
, (210)
where
is a positive parameter with dimension of mass. The corresponding potential is constructed from
and the Kahler potential
according to the usual prescription:
. (211)
where
is the Kahler covariant derivative,
,
and a sum over each superfield
is implicit. Equations (210) and (211) imply that
decays exponentially with
. Here, it will suffice to perform the calculation using a simple exponential potential, namely:
, (212)
where
and
are positive, dimensionless constants. Note that, in the case where the potential is generated by the exchange of wrapped M2-branes, the parameter
is of the form
, where
is a constant,
is the tension of the M2-brane, and
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
during this period. We know that, if the potential
is negligible compared to
, the spectrum of fluctuations is not scale-invariant. Hence, we consider the limit where
. This condition, as seen from the equation of motion for
, Equation (198), is satisfied if
initially, or, equivalently, if the bulk brane begins nearly at rest. For the brane to be nearly at rest, one must have
initially. As the brane traverses the fifth dimension,
increases exponentially, whereas
is constant. Hence, the condition
is automatically satisfied. The bulk brane beginning nearly at rest is precisely what we expect for a nearly BPS initial state. Applying the condition
, Equation (199) reduces to:
, (213)
where we have neglected the endpoint contribution at
. On the other hand, Equation (202) gives:
. (214)
Combining the above two expressions, we obtain:
(215)
The right hand side of Equation (215) is approximately equal to 2 in the limit of large
. Thence, we have:
. (216)
Also here, with regard the numerical result of Equation (216), we have the following mathematical connections with Aurea ratio:
Hence, the exponential potential of Equation (212) results in a nearly scale-invariant spectrum of perturbations provided that
and
are satisfied when modes pass outside the effective Hubble horizon.
We next compute the perturbation amplitude, by using Equation (201) to calculate
. The conditions
and
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
when the relevant modes cross outside the horizon, Equation (201) reduces to:
, (217)
with solution:
, (218)
where
are Bessel functions. The coefficients
and
are fixed by requiring that modes well-within the horizon (i.e.
) be Minkowskian vacuum fluctuations, that is:
. (219)
Using this initial condition, we find the following amplitude for modes outside the horizon (with
):
. (220)
Substituting Equation (213) and using
, we find:
. (221)
Finally, we define the time-delay
by:
, (222)
where we have used the equation of motion for
, Equation (198). Note that the time-dependence of
is mild, a necessary condition for the validity of the time-delay formalism. The factor of
is to be evaluated at time
when a given mode crosses outside the horizon during the motion of the bulk brane. Let
denote the value of
at horizon crossing for mode
. Since horizon crossing occurs when
, or, equivalently, when
, Equation (213) gives:
. (223)
Substituting Equations (222) and (223) into Equation (195), we find:
. (224)
This expression for
increases gradually with increasing
, corresponding to a spectrum tilted slightly towards the blue. The blue tilt is due to the fact that, in this example,
is decreasing as the brane moves. That is, the spectral index,
, (225)
exceeds unity. The current CMB data constrains the spectral index to lie in the range about
. Therefore, for our results to be consistent with experiments, we must have:
, (226)
a constraint that is easily satisfied.
With regard the value s of
, i.e. 0.8 and 1.2, we have the following mathematical connections with Aurea ratio:
;
For the value of Equation (226), we have the following mathematical connections:
Thence, we obtain:
We note that
;
, 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:
.
Furthermore, we consider the power-law potential:
, (227)
where
and
are constants. In this case, Equation (208) gives:
(228)
for
. 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
. Let us suppose that
satisfies:
. (229)
Then, Equation (208) reduces to:
. (230)
Hence, the conditions for scale invariance are Equations (229) as well as:
. (231)
Also, for the numerical value of Equations (228) and (231), we have the following mathematical connections with the Aurea ratio:
If
is the unperturbed, homogeneous metric (see Equation (192) with
and
functions of time), the perturbed
metric can be written as:
(232)
where
. We can treat the tensor perturbations
as functions of
and
only. We are interested in the tensor perturbations which satisfy the conditions:
,
, and
. The perturbed 5 d Einstein action to quadratic order is:
(233)
where the second expression is obtained by integrating over
. 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):
, (234)
where
(235)
and
(236)
The critical difference between this tensor equation and the scalar fluctuation equation, Equation (201), is that the effective scale factor
in Equation (201) has been replaced by
.
We introduced a potential to insure that
led to a nearly scale-invariant spectrum, . However,
in the tensor equation is approximately constant (recall that
. Consequently, the root mean square tensor fluctuation amplitude:
. (237)
is not scale-invariant.
With regard Equation (237), we have obtained the following mathematical connections with the Aurea ratio:
;
;
;
.
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,
. 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
describing gravity, the scalar field
, and the matter-radiation fluid is:
, (238)
where
is the determinant of the Friedmann-Robertson-Walker metric
,
is Newton’s constant and
is the Ricci scalar.
The
factor has the property that
as
such that
constant. The revised solution to the equation of motion is
which approaches a constant as
. 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:
the effective equation of state of the scalar field. Case I is where
and the universe is expanding, the example of inflation. Case II is a contracting universe with
. Case III is a contracting universe with
, that is the situation that applies in the cyclic model.
What is required to obtain
? From the expression for
, it is apparent that this is only possible if the potential is negative. In particular, for a negative exponentially steep potential
, the solutions to the equation of motion have a scaling solution in which
is constant and approximately −1. Consequently,
is much greater than unity and nearly constant. The generation of fluctuations for
can be understood heuristically by examining the perturbed Klein-Gordon equation:
(239)
where
has been expanded in Fourier components
with wavenumber
and prime is derivative with respect to conformal time
. The
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
.
We have defined that
. We take the following values:
and
. We obtain the following mathematical connections with the Aurea ratio:
;
.
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:
(240)
and
(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:
(242)
and
. (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
, we have the following mathematical connections with the Aurea ratio:
;
;
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):
, (244)
where the
are the non-compact spatial dimensions;
are the extra dimensions;
is the usual FRW scale factor; and
(245)
where
has Ricci (scalar) curvature
, as evaluated in the compact dimensions. We call the metric
-flat (RF) if
= const. and conformally
-flat (CRF) if
.
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
and warp function
must be functions of time
and extra-dimensional coordinates
only. We parameterize the rate of change of
using quantities
and
defined by:
(246)
where
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
block describing the energy-momentum in the three non-compact dimensions and
block for the
compact directions. The
component is the higher dimensional energy density
.
Associated with the two blocks of space-space components of
are two trace averages:
, (247)
where
are respectively the
and
blocks of the higher dimensional space-time metric. Violating the NEC means that
for at least one null vector
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
or
is less than zero for any space-time point, then the NEC is violated.
The second lemma utilizes the concept of A-averaged quantities:
; (248)
that is, quantities averaged over the extra dimensions with weight factor
where, for simplicity, we restrict ourselves to constant
. Using the fact that the weight function in the
-average is positive definite, a straightforward consequence is:
Lemma 2: If
or
for any
and any
, then the NEC must be violated.
To illustrate the utility of
-averaging, we introduce the CRF metric into the higher-dimensional Einstein equations, and then try to express terms dependent on
in terms of the 4 d effective scale factor using the relation , where:
(249)
and
is the
-dimensional Planck length. The 4 d effective scale factor,
, obeys the usual 4 d Friedmann equations:
(250)
(251)
Note that the 4 d effective energy density
and pressure
are generally different from
and
in the higher dimensional theory if the warp factor is non-trivial. Then, using the Einstein equations, we obtain:
(252)
(253)
There is a range where:
, (254)
which is the case for
(for CRF). Some theorems below rely on choosing
; for this value to be within the range given in Equation (254), it is necessary that
. Since this includes the relevant string and M-theory models, we will implicitly assume this range of
for CRF models. We note that 3, 8 and 13 are Fibonacci’s numbers.
The two relations in Equation (252) can be rewritten:
(255)
(256)
where the values of
that make the last term non-positive are those that are in the range in Equation (254). Recall that
represents the ratio of the total 4 d effective pressure
to the total 4 d effective energy density
.
On the left hand side of Equations (255) and (256), both
and
depend on the warp factor, Ω, but the combination is invariant under shifts
, where
is a constant. Furthermore, the combination tends to have a weak dependence on Ω. For example, if
is homogeneous in
, the left hand side reduces to
, where the dimensionless coefficient
is not very sensitive to Ω or
; in particular,
where
(257)
In this notation, the
-dimensional volume of the compact space is
; then,
is equal to
, a coefficient which is strictly less than unity. Similarly, if
is smooth and Ω has a sharp maximum on some subspace of dimension
and volume
, then the left hand side of Equation (256) is
, where
is the value of
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 (
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
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
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
, 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
and a tetrad
where both the spatial coordinate index
and the spatial tetrad index
go from 1 to 3. Choose
to be hypersurface orthogonal with the relation between tetrad and coordinates of the form
and
where
is the lapse and the shift is chosen to be zero. Choose the spatial frame
to be Fermi propagated along the integral curves of
. The commutators of the tetrad components are decomposed as follows:
(258)
where
is symmetric, and
is symmetric and trace free. The scale invariant tetrad variables are defined by
and
while scale invariant versions of the other gravitational variables are given by:
(259)
Note that the relation between the scale invariant tetrad variables and the coordinate derivatives is:
(260)
where
is the scale invariant lapse. The matter model is a scalar field
with potential
of the form:
, (261)
where
and
are positive constants. The scale invariant matter variables are given by:
(262)
. (263)
The time coordinate
is chosen so that
. (264)
Note that this means that the surfaces of constant time are constant mean curvature surfaces. Note also that the singularity is approached as
. Due to Equation (264) the scale invariant lapse satisfies an elliptic equation:
. (265)
We note that 3 is a Fibonacci’s number. Furthermore, we have the following mathematical connection with the Aurea ratio:
.
The gravitational quantities
and
satisfy the following hyperbolic evolution equations:
(266)
(267)
(268)
(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:
(270)
(271)
(272)
In addition, the variables are subject to the vanishing of the following constraint quantities:
(273)
(274)
(275)
(276)
(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:
.
For simplicity, we choose the initial conformal metric to be flat and
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:
(278)
It then follows from Equation (258) that:
. (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
defined by:
. (280)
Similar considerations apply to the matter variables, leading us to define the quantity
given by:
. (281)
Here, we will specify
and a part of
and solve the constraint equations for the conformal factor
and the rest of
. From Equation (275) and our ansatz for the scale invariant variables we obtain:
. (282)
In the vacuum case, this equation simply becomes the conditions that
is divergence-free, which is in turn simply an algebraic condition on the Fourier coefficients of
. Note that since
must be trace-free, so must
. A simple, but still fairly general divergence-free and trace-free
is the following:
, (283)
where
and
are constants. We still keep this divergence-free part of
but now add to it a piece that has a non-zero divergence. We simply specify the Fourier coefficients of
and
via:
,
(284)
where
and
are constants. This turns Equation (282) into an algebraic equation for the Fourier coefficients of this non-zero divergence piece of
which we then solve. Now imposing Equation (276) our ansatz yields:
, (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:
.
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
.
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
matter dominated regions, or
mixmaster-like regions. Whenever a
region forms it grows exponentially fast in proper volume relative to
regions. The particular initial conditions for this example are (283)-(284):
,
, and
(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:
Now, let
be a complex character to the modulus 5 such that
, and let:
. (287)
The function:
,
where
, is called the Davenport-Heilbronn function and satisfies the Riemann-type equation:
. (288)
We note that
. Furthermore:
;
. Thence, we can write also:
It is enlightening to visualize the evolution via the behaviour of the matter
, shear
and curvature
contributions to the normalized energy density, defined as:
(289)
, (290)
where
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:
.
The effective equation of state parameter
takes the following form in Hubble normalized variables:
. (291)
It is evident that at late times the region that has smoothed out and become matter dominated coincides with
, whereas the mixmaster-like regime evolves to
. 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:
, (292)
and slicing condition for
(265) becomes:
, (293)
Furthermore,
is finite and non-zero. This implies from (262), (265), (266) that asymptotic form:
(294)
and thus
(263) tends to:
. (295)
Combining these relations gives:
, (296)
, (297)
, (298)
and from (291):
. (299)
We have that
;
and
. Thence, we have the following mathematical connections with the Aurea ratio:
Let
denote the proper spatial volume element associated with the spatial metric
of
= const. slices, i.e.
. The fractional change of
with respect to time is:
, (300)
which can be written as:
. (301)
In the asymptotic regime where spatial gradients are negligible,
approaches a constant (293), and thus (301) can be integrated to give:
,
(302)
where we have used (293) where
, and note that
when
. Thus, at late times the ratio
of the proper volume of matter to mixmaster-like regions of the universe grows as:
. (303)
Thus, as long as
(which is equivalent to
),
as
. We note that
, is related to the following mathematical connection with the Aurea ratio:
(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:
(304)
where
in the flat space and we have defined:
and
. (305)
The dimensionless scalar field
describes the open string tachyon,
is the string mass scale and
is the open string coupling constant. Though the action (304) was originally derived for
a prime number, it appears that it can be continued to any positive integer and even makes sense in the limit
. Setting
in the action, the resulting potential takes the form:
. (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:
(307)
We are interested in perturbing around the solution
, 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:
(308)
where we have defined the field
as:
(309)
and the potential is:
. (310)
In (308), the
denotes terms with higher powers of
. Thence, Equation (308) can be rewritten also:
(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
. It is straightforward to show that:
(312)
(313)
With regard the approximate solution for the classical background, we must solve the Friedmann equation:
(314)
to second order in
. To find the energy density
, we turn to the stress energy tensor for the p-adic scalar field. A convenient expression for
is:
(315)
One may verify that the
is symmetric by changing the dummy integration variable
in the last term. For homogeneous
, the above expression simplifies, and for
, we find:
(316)
One can evaluate the above expression term by term, keeping up to
. The final result reads:
(317)
The
terms cancel out and matching the coefficients in the Friedmann equation gives us the simple results:
(318)
and
(319)
for zeroth and first order respectively. The
contribution to
is quite complicated but once we use (319) it simplifies greatly. Matching coefficient at order
in the Friedmann equation gives:
. (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:
.
Because of our sign convention for
, the fact that
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
, so that:
(321)
. (322)
We expand the p-adic tachyon field in perturbation theory as:
. (323)
The perturbed Klein-Gordon Equation (307) takes the form:
. (324)
One can construct solutions by taking
to be an eigenfunction of the operator. If we choose
to satisfy:
(325)
then this is also a solution to (324) if:
(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:
(327)
Thus, for on-shell fields the kinetic term in the Lagrangian can be written as:
(328)
In (328), we have defined the “canonical” field:
(329)
where
. (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
is the appropriate canonically normalized field. We write the quantum mechanical solution in term of annihilation/creation operators as:
(331)
and the mode functions
are given by:
(332)
where the order of the Hankel functions is:
(333)
and of course
. 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,
, one has:
which reproduces the standard Minkowski space fluctuations. This is the usual procedure in cosmological perturbation theory. On large scales,
, the solutions (331) behave as:
which gives a large-scale power spectrum for the fluctuations:
with spectral index:
From (333) it is clear that to get an almost scale-invariant spectrum we require
. In this limit, we have:
(334)
which gives a red tilt to the spectrum, in agreement with the latest WMAP data. For
one has
. 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
which may be quite different from the mass scale which one would infer by truncating the infinite series of derivatives:
.
We note that for
, from Equation (334), we obtain that
. Thence, we have the following mathematical connections with the Aurea ratio:
From Equation (333), we have that:
;
We note that for Equations (287)-(288)-(289), we have that:
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,
and the ratio
. 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:
. (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:
Thus, one can have a small tilt while ensuring that the string scale is smaller than the Planck scale, provided that
. Henceforth we will use (335) to determine
in terms of
, and
. All the dimensionless parameters in our solution,
and
, are likewise functions of
and
. From (318) and (335) we see that for
,
. (336)
It may seem strange to have
exceeding
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
.
Next, we determine
, where
is the mass scale appearing in the power series in
which provides the ansatz for the background solutions. We consider the following equation for
in the
limit:
The positive root for
gives:
. (337)
In order to fix the amplitude of the density perturbations we consider the curvature perturbation
. We assume that:
as in conventional inflation models. To evaluate the prefactor
we must work beyond zeroth order in the small
expansion. We take
to evaluate the prefactor, even though the perturbation
is computed in the limit that
. This should reproduce the full answer up to
corrections. The prefactor is:
With regard Equations (336) and (337), we have the following mathematical connections with the Aurea ratio:
;
.
;
We should evaluate
at the time of horizon crossing,
, defined to be approximately 60 e-foldings before the end of inflation
, assuming that the energy scale of inflation is high. The inflation ends when
. From Equations (336)-(337), we see that
; therefore, we can write the scale factor
in the form:
(338)
so that
corresponds to:
(339)
We note that
for
is equal to 40. This value can be related with the following mathematical connections with the Aura ratio:
.
The power spectrum of the curvature perturbation is given by:
(340)
where the amplitude of fluctuations
can now be read off as:
(341)
Thence, we can rewrite Equation (340) as follows:
. (342)
As an example, taking
one can fix the amplitude of the density perturbations
by choosing:
. (343)
Setting
and using (341), we obtain an expression for
in terms of
and
:
. (344)
Combining (342) with (335), we also obtain:
. (345)
The string scale is bounded from above as
and that for typical values of
, it is close to
. Furthermore, from (342) that
is unconstrained and that
,
are not independent parameters.
Now we define the Hubble slow roll parameters
,
by:
,
. (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
(recall that
,
) we find that:
,
. (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):
,
. (348)
With regard Equations (347) and (348), we have that:
;
On the other hand, consider the potential slow roll parameter which one would naively define using the derivative truncated action (308):
,
(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
,
, and
, we have
. Since large values of
are required if one wants to obtain
, 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:
(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:
(351)
Here,
and
are functions of time and
is the metric on the unit 3-sphere. The action for this metric is the same as the corresponding usual case:
(352)
We assume that the cosmological constant
is a rational number. The classical equations of motions have the form:
. (353)
The solution of this equation for the boundary conditions:
, (354)
is the following:
. (355)
Here,
. The Green function corresponding to the transition from the point
to the point
has the form:
(356)
where
is the propagator:
. (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:
(358)
where
is the action calculated on the trajectories (355).
. (359)
The factor
is the same as for a free particle:
. (360)
Therefore, one has the Green function:
. (361)
The corresponding wave function has the form:
(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:
.
Now let us estimate the integral on
applying the stationary phase approximation. The saddlepoint equation has the form:
(363)
which yields:
. (364)
As is known, for
, there is the square root of −1 in
, so we get nontrivial saddle points. For
, we have no saddle point at all. To make sense of the saddle points in the case
, we should be sure that the square root
also has a sense. For this purpose, we have to assume that
. The corresponding actions have the form:
,
, (365)
In order that these expressions are rational, we have to assume that
is rational as well as that
is such a rational that the solution of the equation:
(366)
in respect to
is also rational. Let us consider the Euclidean metric:
. (367)
The Euclidean action for this metric is the same as the corresponding action in the usual case:
. (368)
We shall prove that it is possible to restore:
(369)
when
in the corresponding p-adic wave function. Indeed, in the p-adic case for Euclidean metric, we get a basic Green function:
. (370)
Now let us estimate the integral on
in:
(371)
applying the stationary phase approximation. Formally, there are the following saddle points:
,
, (372)
for
with corresponding actions:
,
. (373)
Note that now these stationary points have sense for all
and for
satisfying (192) according our general formula for the wave function of the universe:
, (374)
we write:
(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:
, (376)
where
is any prime number,
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:
,
,
. (378)
Employing usual expansion for the logarithmic function and definition (378) we can rewrite (377) in the form:
, (379)
where
.
acts as pseudodifferential operator in the following way:
, (380)
where is the Fourier transform of
.
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:
(381)
which has an evident solution
.
For the case of time dependent spatially homogeneous solutions, we have the following equation of motion:
. (382)
With regard the open and closed scalar zeta strings, the equations of motion are:
, (383)
(384)
and one can easily see trivial solution
.
The exact tree-level Lagrangian of effective scalar field
, which describes open p-adic string tachyon, is:
, (385)
where
is any prime number,
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
replaced by
. Thence, we take the sum of all Lagrangians
in the form:
, (386)
whose explicit realization depends on particular choice of coefficients
, masses
and coupling constants
.
Now, we consider the following case:
, (387)
where
is a real number. The corresponding Lagrangian reads:
(388)
and it depends on parameter
. According to the Euler product formula, one can write:
. (389)
Recall that standard definition of the Riemann zeta function is:
, (390)
which has analytic continuation to the entire complex
plane, excluding the point
, where it has a simple pole with residue 1. Employing definition (390), we can rewrite (388) in the form:
(391)
Here,
acts as a pseudodifferential operator:
, (392)
where is the Fourier transform of
.
We consider Lagrangian (391) with analytic continuations of the zeta function and the power series
, i.e.
, (393)
where
denotes analytic continuation.
Potential of the above zeta scalar field (393) is equal to
at
, i.e.
, (394)
where
since
. The term with
-function vanishes at
. The equation of motion in differential and integral form is:
, (395)
, (396)
respectively.
Now, we consider five values of
, which seem to be the most interesting, regarding the Lagrangian (393):
,
, and
. For
, the corresponding equation of motion now read:
. (397)
This equation has two trivial solutions:
and
. Solution
can be also shown taking
and
in (397).
For
, the corresponding equation of motion is:
. (398)
where
.
The equation of motion (398) has a constant trivial solution only for
.
For
, the equation of motion is:
. (399)
It has two solutions:
and
. The solution
follows from the Taylor expansion of the Riemann zeta function operator:
, (400)
as well as from
.
For
, the equation of motion is:
, (401)
where
gives
.
In conclusion, for
, we have the following equation of motion:
. (402)
Since holds equality:
one has trivial solution
in (402).
Now, we want to analyze the following case:
. In this case, from the Lagrangian (386), we obtain:
(403)
The corresponding potential is:
. (404)
The equation of motion is:
. (405)
Its weak field approximation is:
, (406)
which implies condition on the mass spectrum:
. (407)
From (407) it follows one solution for
at
and many tachyon solutions when
.
With regard the extension by ordinary Lagrangian, we have the Lagrangian, potential, equation of motion and mass spectrum condition that, when
, are:
, (408)
, (409)
, (410)
. (411)
In addition to many tachyon solutions, Equation (411) has two solutions with positive mass:
and
.
Now, we describe the case of
. Here
is the Mobius function, which is defined for all positive integers and has values
depending on factorization of
into prime numbers
. It is defined as follows:
(412)
The corresponding Lagrangian is:
(413)
Recall that the inverse Riemann zeta function can be defined by:
. (414)
Now, (413) can be rewritten as:
, (415)
where
The corresponding potential, equation of motion and mass spectrum formula, respectively, are:
, (416)
, (417)
,
, (418)
where usual relativistic kinematic relation
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
, thence with the aurea ratio, by the following expressions:
. (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
, and matter configuration
on an initial hypersurface
, to another state with metric
, and matter configuration
on a final hypersurface
, is given by a functional integral of
over all four-geometries
, and matter configurations Φ, which interpolate between the initial and final configurations, i.e.
. (420)
The
is the usual Einstein-Hilbert action:
(421)
for the gravitational field and matter fields Φ. In (421), R is scalar curvature of four-manifold
, Λ is cosmological constant,
is trace of the extrinsic curvature
at the boundary
of the manifold
. 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:
. (422)
The space of all 3-metrics and matter field configurations
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
at each point
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
, where
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
,
then integral in (420) and (422), in the gauge
, is a minisuperspace propagator. In this case, it holds:
, (423)
where
(424)
is an ordinary quantum-mechanical propagator between fixed
in fixed time
. For quadratic classical action
, (424) becomes:
. (425)
If system has
decoupled degrees of freedom, its p-adic kernel is a product:
. (426)
p-adic and adelic wave functions of the universe may be found by means of the following equation:
, (427)
where
are adelic wave eigenfunctions,
is the corresponding adelic energy,
and
are indicies for energy levels and their degeneration, respectively.
The corresponding adelic eigenstates have the form:
. (428)
A necessary condition to construct an adelic model is existence of the p-adic (vacuum) state
, which satisfies:
(429)
for all but a finite number of
.
Now, we describe the p-adic and adelic model with cosmological constant in
dimensions. This model has the metric:
, (430)
where
. The corresponding v-adic action is:
, (431)
where
. The Euler-Lagrange equation of motion:
has the solution:
, (432)
where the boundary conditions are
,
. For the classical action, it gives:
. (433)
Quantum-mechanical propagator has the form:
(434)
Equation (429), in a more explicit form, reads:
(435)
We note that the p-adic Gauss integral over the region of integration
is:
(436)
(437)
where
is defined as follows:
;
.
Using (437), for
, we obtain:
(438)
with condition
, i.e.
. For
, left hand side is equal to
if
holds. Applying also the (436) to (435), we have:
(439)
It becomes an equality if condition
take place.
Thence, we can rewrite Equation (435) as follow:
(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:
, (441)
where
and
is the scale factor. Instead of (441), we shall use:
. (442)
The corresponding v-adic action for this one-dimensional minisuperspace model is:
, (443)
where
. The classical equation of motion
, with the boundary conditions
and
gives:
. (444)
After substitution (444) into (443) and integration, one obtains that the classical action is:
. (445)
Since (445) is quadratic in
and
, quantum-mechanical propagator has the form:
. (446)
Equation (429) reads:
(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:
(448)
Thence, there is a mathematical connections between the slow-roll formula regarding the number of e-foldings
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:
; (449)
; (450)
; (451)
. (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:
(453)
; (454)
; (455)
. (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:
. (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:
;
. (458)
9. Conclusion
In this work, we have described mathematical connections between the slow-roll formula regarding the number of e-foldings
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.