Evolution of sub-spaces at high and low energies

TL;DR

This paper investigates how sub-spaces evolve at different energies within higher-derivative gravity, using numerical solutions to understand the stability and size of our 3D space versus extra dimensions.

Contribution

It provides numerical solutions to differential equations describing sub-space evolution and highlights the importance of initial conditions in space dynamics within higher-derivative gravity.

Findings

Initial conditions crucially influence space evolution.

Stable, large 3D space can coexist with small, stable extra dimensions.

Parameter values at high energies are highly uncertain due to quantum corrections.

Abstract

The evolution of sub-spaces in the framework of gravity with higher derivatives is studied. Numerical solutions to exact differential equations are found. It is shown that the initial conditions play crucial role in the space dynamic. Appropriate metrics describing an expanding and a stationary sub-space shed light on the well-known question: why our 3-dim space is large but extra space is small and stable (if exists)? It is assumed that the values of parameters at high energies strongly depend on uncontrolled quantum corrections and, hence, are not equal to their values at low energies. Therefore, there is no way to trace solutions throughout the energy range, and we restrict ourselves to the sub-Planckian and the inflationary energies.

Figures

Equations

$$\displaystyle S_{eff}=\frac{1}{2}m_{D}^{D-2}\int d^{D}x\sqrt{|g_{D}|}\Big{[}f(R)+c_{1}R_{AB}R^{AB}+c_{2}R_{ABCD}R^{ABCD}\Big{]}$$

$$\displaystyle{\cal L}_{GB}=k\sqrt{-g}\biggr{\{}R^{2}-4R_{AB}R^{AB}+R_{ABCD}R^{ABCD}\biggl{\}}$$

$$S_{gen}=\frac{1}{2}m_{D}^{D-2}\int d^{D}x\sqrt{|g_{D}|}\Big{[}\tilde{f}(R)+k\biggr{\{}R^{2}-4R_{AB}R^{AB}+R_{ABCD}R^{ABCD}\biggl{\}}\Big{]},$$

$$f(R)=a{R}^2+bR+c$$

$$d{s}^2=$$

$$\displaystyle-\frac{1}{2}\tilde{f}(R)\delta_{B}^{A}+(R_{B}^{A}+\nabla^{A}\nabla_{B}-\delta_{B}^{A}\square)\tilde{f}_{R}+k\Big{[}-8{R^{AC}}_{;BC}-12R^{AC}R_{CB}$$

$$+2{\delta }_B^A{R}_{CD}{R}^{CD}-\frac{{\delta }_B^A}{2}{R}^{CDEF}{R}_{CDEF}+2R^{ACDE}{R}_{BCDE}+4{R^{;A}}_{;B}$$

$$\displaystyle+4{R^{CAD}}_{B}^{\ \cdot}\ R_{CD}-\frac{\delta^{A}_{B}}{2}R^{2}+2RR^{A}_{B}\Big{]}=0,$$

$$-36k\left[{\dot{\alpha }}^3\dot{\beta }+3{\dot{\alpha }}^2{\dot{\beta }}^2+\dot{\alpha }{\dot{\beta }}^3+e^{-2\alpha }\dot{\beta }\left(\dot{\alpha }+\dot{\beta }\right)+e^{-2\beta }\dot{\alpha }\left(\dot{\alpha }+\dot{\beta }\right)+e^{-2\alpha }{e}^{-2\beta }\right]$$

$$-3{\tilde{f}}_{RR}\dot{R}\left(\dot{\alpha }+\dot{\beta }\right)+3{\tilde{f}}_R\left(\ddot{\alpha }+\ddot{\beta }+{\dot{\alpha }}^2+{\dot{\beta }}^2\right)-\frac{1}{2}\tilde{f}(R)=0,$$

$$-12k\{2\ddot{\alpha }\dot{\beta }\left(\dot{\alpha }+\dot{\beta }\right)+\ddot{\beta }\left({\dot{\alpha }}^2+4\dot{\alpha }\dot{\beta }+{\dot{\beta }}^2\right)+2{\dot{\alpha }}^3\dot{\beta }+6{\dot{\alpha }}^2{\dot{\beta }}^2+6\dot{\alpha }{\dot{\beta }}^3+{\dot{\beta }}^4$$

$$+e^{-2\alpha }\left[\ddot{\beta }+2{\dot{\beta }}^2\right]+e^{-2\beta }\left[2\ddot{\alpha }+\ddot{\beta }+3{\dot{\alpha }}^2+2\dot{\alpha }\dot{\beta }+{\dot{\beta }}^2\right]+e^{-2\alpha }{e}^{-2\beta }\}-{\tilde{f}}_{RRR}{\dot{R}}^2$$

$$-{\tilde{f}}_{RR}\left[\ddot{R}+\dot{R}\left(2\dot{\alpha }+3\dot{\beta }\right)\right]+{\tilde{f}}_R\left(\ddot{\alpha }+3{\dot{\alpha }}^2+3\dot{\alpha }\dot{\beta }+2e^{-2\alpha }\right)-\frac{1}{2}\tilde{f}(R)=0,$$

$$-12k\{\ddot{\alpha }\left({\dot{\alpha }}^2+4\dot{\alpha }\dot{\beta }+{\dot{\beta }}^2\right)+2\ddot{\beta }\dot{\alpha }\left(\dot{\alpha }+\dot{\beta }\right)+{\dot{\alpha }}^4+6{\dot{\alpha }}^3\dot{\beta }+6{\dot{\alpha }}^2{\dot{\beta }}^2+2\dot{\alpha }{\dot{\beta }}^3$$

$$+e^{-2\alpha }\left[\ddot{\alpha }+2\ddot{\beta }+{\dot{\alpha }}^2+2\dot{\alpha }\dot{\beta }+3{\dot{\beta }}^2\right]+e^{-2\beta }\left[\ddot{\alpha }+2{\dot{\alpha }}^2\right]+e^{-2\alpha }{e}^{-2\beta }\}-{\tilde{f}}_{RRR}{\dot{R}}^2$$

$$-{\tilde{f}}_{RR}\left[\ddot{R}+\dot{R}\left(3\dot{\alpha }+2\dot{\beta }\right)\right]+{\tilde{f}}_R\left(\ddot{\beta }+3\dot{\alpha }\dot{\beta }+3{\dot{\beta }}^2+2e^{-2\beta }\right)-\frac{1}{2}\tilde{f}(R)=0,$$

$$R=6\left(\ddot{\alpha }+\ddot{\beta }+2{\dot{\alpha }}^2+3\dot{\alpha }\dot{\beta }+2{\dot{\beta }}^2+e^{-2\alpha }+e^{-2\beta }\right).$$

$$-36k\left[{\dot{\alpha }}^3\dot{\beta }+3{\dot{\alpha }}^2{\dot{\beta }}^2+\dot{\alpha }{\dot{\beta }}^3+e^{-2\alpha }\dot{\beta }\left(\dot{\alpha }+\dot{\beta }\right)+e^{-2\beta }\dot{\alpha }\left(\dot{\alpha }+\dot{\beta }\right)+e^{-2\alpha }{e}^{-2\beta }\right]$$

$$-3\left(\dot{\alpha }+\dot{\beta }\right)\dot{R}{\tilde{f}}_{RR}+\left(-3{\dot{\alpha }}^2-9\dot{\alpha }\dot{\beta }-3\dot{\beta }-3e^{-2\alpha }-3e^{-2\beta }+\frac{R}{2}\right){\tilde{f}}_R-\frac{\tilde{f}}{2}=0,$$

$$\alpha (t)=H_{1}t,\beta (t)=H_{2}t,H_1>0,H_2>0,t\to \infty .$$

$$\displaystyle-H_{1}^{2}-3H_{1}H_{2}-H_{2}^{2}+\Big{(}24H_{1}^{3}H_{2}+18H_{1}^{2}H_{2}^{2}+24H_{1}H_{2}^{3}\Big{)}k$$

$$\displaystyle-\Big{(}36H_{1}^{3}H_{2}+54H_{1}^{2}H_{2}^{2}+36H_{1}H_{2}^{3}\Big{)}a-\frac{c}{6}=0,$$

$$\displaystyle-H_{1}^{2}-2H_{1}H_{2}-2H_{2}^{2}+\Big{(}4H_{1}^{3}H_{2}+18H_{1}^{2}H_{2}^{2}+24H_{1}H_{2}^{3}+20H_{2}^{4}\Big{)}k$$

$$\displaystyle-\Big{(}12H_{1}^{3}H_{2}+42H_{1}^{2}H_{2}^{2}+48H_{1}H_{2}^{3}+24H_{2}^{4}\Big{)}a-\frac{c}{6}=0,$$

$$\displaystyle-2H_{1}^{2}-2H_{1}H_{2}-H_{2}^{2}+\Big{(}20H_{1}^{4}+24H_{1}^{3}H_{2}+18H_{1}^{2}H_{2}^{2}+4H_{1}H_{2}^{3}\Big{)}k$$

$$\displaystyle-\Big{(}24H_{1}^{4}+48H_{1}^{3}H_{2}+42H_{1}^{2}H_{2}^{2}+12H_{1}H_{2}^{3}\Big{)}a-\frac{c}{6}=0.$$

$$(66k-126a)H^4-5H^2-\frac{c}{6}=0,$$

$$H=H_1=H_2=\sqrt{\frac{5\pm \sqrt{25+44kc-84ac}}{12(11k-21a)}}.$$

$$\alpha (t)=Ht,\beta (t)={\beta }_{as},H>0,{\beta }_{as}>0,t\to \infty .$$

$$-\left(1+\frac{12a}{e^{2{\beta }_{as}}}\right)H^2+\frac{6(k-a)}{e^{4{\beta }_{as}}}-\frac{1}{e^{2{\beta }_{as}}}-\frac{c}{6}=0,$$

$$\left(60k-72a\right)H^4-\left(6+\frac{24a}{e^{2{\beta }_{as}}}\right)H^2-\frac{6(k-a)}{e^{4{\beta }_{as}}}-\frac{1}{e^{2{\beta }_{as}}}-\frac{c}{2}=0$$

$$R=R_4+R_3+P_k;P_k\equiv 6\ddot{\beta }+12{\dot{\beta }}^2+18\dot{\alpha }\dot{\beta }$$

$$P_k,R_4\ll R_3$$

$$S=$$

$$g_{ik}\rightarrow g_{ik}^{(E)}=|e^{3\beta}f^{\prime}(R_{3})|g_{ik};\quad R_{4}^{(E)}=\Big{|}e^{3\beta}f^{\prime}(R_{3})\Big{|}^{-1}R_{4}$$

$$S_E=\frac{v_3}{2}\int d^{4}x\sqrt{-g_4^{(E)}}\text{mbox}sign(f^{\prime })[R_4^{(E)}+K_E(\varphi ){\dot{\varphi }}^2-2V_E(\varphi )],$$

Full text

I Introduction

The origin of our Universe remains as an unresolved problem up to now. It is usually assumed that its nucleation is related to the quantum processes at high energies van Elst et al. (1994); Bousso and Linde (1998); Halliwell (2009); Yurov and Yurov (2005). The probability of its creation remains unclear in spite of wide discussion, see e.g. Kubyshin et al. (1989); Firouzjahi et al. (2004); Chopovsky et al. (2013). Here we are interested in the subsequent classical evolution of the metrics rather than a calculation of this probability. It is assumed that manifolds can be described by specific metrics after their nucleation. After nucleation, these manifolds evolve classically forming a set of asymptotic manifolds, one of which could be our Universe.

The complexity of the problem is greatly aggravated by two factors. First, the metric evolution should lead to the formation of our Universe with the strong fine-tuning of the observational parameters Donoghue (2007); Bauer et al. (2010). Second, the inclusion of extra dimensions is of particular interest because the idea of extra space is widely used in modern research. They shed light to such issues as the grand unification Arkani-Hamed et al. (1998); Dienes et al. (1999), neutrino mass Arkani-Hamed et al. (2002), the cosmological constant problem Starkman et al. (2001); Günther et al. (2003); Rubin (2018) and so on. In this regard, the immediately aroused question is: why specific number of dimensions are asymptotically compact and stable while others expand Greene and Levin (2007); Carroll et al. (2002); Nasri et al. (2002)? Which specific property of subspace leads to its quick growth? Sometimes one of the subspaces is assumed to be FRW space by definition Bronnikov and Rubin (2013). There are many attempts to clarify the problem, mostly related to introduction of fields other than gravity. It may be a scalar field Ketov and Nakada (2017); Greene and Levin (2007) (most used case) and gauge fields Kihara et al. (2009) for example. A static solution can be obtained using the Casimir effect Candelas and Weinberg (1984); Bolokhov and Bronnikov (2018) or form fields Freund and Rubin (1980). Another possibility was discussed in Yoshimura (1984); Bronnikov et al. (2013): it was shown that if the scale factor $a(t)$ of our 3D space is much larger than the growing scale factor $b(t)$ of the extra dimensions, a contradiction with observations can be avoided.

In our previous article Lyakhova et al. (2018) we studied evolution of manifolds after their creation on the basis of pure gravitational Lagrangian with higher derivatives. It was shown analytically and confirmed numerically that an asymptotic growth of the manifolds depends weakly on initial conditions. We have shown that the initial conditions can be a reason of nontrivial solutions (funnels) and studied their properties. A number of final states of metric describing our Universe is quite poor if we limit ourselves with a maximally symmetric extra space and the $f(R)$ gravity.

In this article, we continue to study the Universe evolution at the sub-Planckian scale. The space $V_{D}$ is assumed to be the direct product $T\times V_{3}\times W_{3}$ of the time and the two maximally symmetric manifolds with positive curvature. Both sub-spaces are born with the size of the order of the Planck scale or more. Most of the resulting sub-spaces are characterized by initial metrics, which lead to the growth of both 3-dimensional sub-spaces, which clearly contradicts the observations. We have found a set of metrics that could lead to the observable space-time metric.

The action used should not contradict the observations at the low and intermediate energies. More definitely, the Lagrangian parameters at low energies should be chosen in such a way to supply (almost) Minkowski space for the modern Universe, the stationarity of the extra space metric at the modern epoch, and reproduce the inflationary stage with the Hubble parameter of the order of $H\sim 10^{13}$GeV. On the contrary, the parameter values at the sub-Planckian energies are free from such restrictions. Indeed, the quantum corrections to the parameter values cause their dependence on the energy scale. There is a lot of literature devoted to this subject, see e.g. Babic et al. (2002). The quantum corrections to the physical parameters are large at the very high energies where interactions of all kinds of fields must be taken into account. Our knowledge of the physical parameters at the low energies is blurring at sub-Planckian scale.

The gravity with higher derivatives is widely used in modern research despite the internal problems inherent in this approach Woodard (2015). Most promising model of inflation is the Starobinsky model based on pure gravitational action. Attempts to avoid the Ostrogradsky instabilities are made Paul (2017) and extensions of the Einstein-Hilbert action attract much attention. Promising branch of such models is based on the Gauss-Bonnet Lagrangian Ivashchuk (2010); Pavluchenko and Toporensky (2018) and its generalization to the Horndeski gravity Starobinsky et al. (2016). These models were adjusted to obtain differential equations of the second order so that the Ostrogradsky theorem is not dangerous for them. Our model contains $f(R)$ gravity with the Gauss-Bonnet term acting in $D=7$ dimensions. The similar research was performed in Shafi and Wetterich (1987). The authors considered the action containing all scalars made from the Riemann tensor up to the second term multiplied by an arbitrary factor. They discuss asymptotic evolution of the two subspaces with one of them tending to a stationary state. The equations were simplified by the slow motion approximation. Our analysis is based on exact solutions to the nonlinear equations and hence can be applied to the whole variety of initial metrics. This ability is important for our study, which aims to study the role of the initial conditions.

Throughout this paper we use the conventions for the curvature tensor $R_{ABC}^{D}=\partial_{C}\Gamma_{AB}^{D}-\partial_{B}\Gamma_{AC}^{D}+\Gamma_{EC}^{D}\Gamma_{BA}^{E}-\Gamma_{EB}^{D}\Gamma_{AC}^{E}$ and for the Ricci tensor $R_{MN}=R^{F}_{MFN}$. The units $\hbar=c=1$ are also used.

II The model

The action

[TABLE]

can be considered as a basis of an effective theory Burgess (2007). Here $f(R)$ is an arbitrary function of the Ricci scalar $R$, $m_{D}$ is the D-dimensional Planck mass and $c_{1},c_{2}$ are parameters of the Lagrangian. Such action was used in Shafi and Wetterich (1987) to study the extra space dynamic in the slow motion limit. Its particular case, the Gauss-Bonnet Lagrangian

[TABLE]

is the appropriate starting point because of the absence of higher derivatives in the equations of motion. $f(R)$ extensions of the gravitational action is also applied in this content Troisi (2017).

The action

[TABLE]

containing the Gauss-Bonnet term plus $\tilde{f}(R)$ term (a function of the Ricci scalar) will be used throughout the paper. Such action was used to describe the dark energy phenomenon Nojiri et al. (2007) for example. The action (3) is the particular case of the action (1) provided that $\tilde{f}(R)=f(R)-kR^{2}$ and $c_{1}=-4k,c_{2}=k$. In what follows we restrict ourselves to the quadratic function

[TABLE]

($b=1$ without the loss of generality).

A separate problem is the values of the physical parameters. It is well known that values of Lagrangian parameters depend on the energy scale. There are at least two approaches to calculate this dependence - the Renormalization group analysis Babic et al. (2002) and the Effective field theory Burgess (2007). In any case, the physical parameters are functions $a(M),k(M),c(M)$ of an energy scale $M$. The next Section relates to the highest energy scale of the D-dim Planck mass, $M\simeq m_{D}$ and hence we may find the specific values $a(m_{D}),k(m_{D}),c(m_{D})$ but not functions $a(M),k(M),c(M)$. Unfortunately, there are no way to connect accurately these functions at the inflationary ($M_{infl}$) scale and at the sub-Planckian scale. Some restrictions to the parameters $a(M_{infl}),k(M_{infl}),c(M_{infl})$ are found in Section IV.

III Subspaces evolution at the Planck energies

In this section, the evolution of extra spaces at highest energies $m_{D}$ is analyzed on the basis of action (3). Both extra spaces are described by 3-dim maximally symmetrical metrics with positive curvature

[TABLE]

Below, we show that one of the sub-spaces can evolve into a space with a large volume while the second one - into a small static extra space.

Einstein’s equations for this model are

[TABLE]

Our aim is to study the behaviour of the two scale factors at highest energies $E\lesssim m_{D}$ by solving the main system (III), (III), (III), (10). The main goal is to find a solution - $\alpha(t)=Ht,\ H=const$ and $\beta(t)=const$. The Hubble parameter $H$ at high energies is not related to those at the inflationary stage.

The nontrivial system of equations (III) are

[TABLE]

[TABLE]

[TABLE]

where we have kept in mind denotations $\partial_{t}\tilde{f}_{R}=\tilde{f}_{RR}\dot{R}$ and $\partial^{2}_{t}\tilde{f}_{R}=\tilde{f}_{RRR}\dot{R}^{2}+\tilde{f}_{RR}\ddot{R}$. The Ricci scalar is

[TABLE]

Here and in the following the units $m_{D}=1$ are assumed.

For calculations, it is convenient to consider the Ricci scalar $R(t)$ as the additional unknown function and interpret definition (10) as the fourth equation. Three equations of this system (for example, (III), (III), (10)) can be solved with respect to the higher derivatives $\ddot{\alpha},\ddot{\beta},\ddot{R}$. Then, substitution $\ddot{\alpha}$ and $\ddot{\beta}$ into equation (III) gives equation

[TABLE]

which plays the role of restriction to the solutions of the coupled second order differential equations. This can be checked, for example, by writing the set of four equations (III), (III), (10), (III) as an equivalent set of (six) coupled first-order equations plus one algebraic equation. The equation (III) reduces to the algebraic transcendental equation, i.e., it is a constraint. The complete set of initial conditions may be chosen in the form $\alpha(t_{0}),\beta(t_{0}),R(t_{0}),\dot{\alpha}(t_{0}),\dot{\beta}(t_{0})$, and $\dot{R}(t_{0})$. These initial conditions are not independent due to equation (III). The latter will be used to derive an exact relation between these initial data.

Natural values of the parameters are assumed to be of the order of the Planck scale: $a\sim k\sim m_{D}^{-2},c\sim m_{D}^{2}$. They are not related to the observational values because of strong and uncontrolled contribution of the quantum corrections at sub-planckian energies. Hence, they are considered as free parameters.

Our analysis revealed a complex dynamic of the sub-spaces depending on the values of the parameters $k,a,c$ and the initial metric. There are several variants for the metric evolution:

i) both sub-spaces expand at equal rates;

ii) both sub-spaces expand at different rates;

iii) one of the sub-spaces expands while the other remains constant.

iv) one or both sub-spaces shrink. Looking ahead, we note that such solutions do exist (see solid line in the left panel of Fig. 4). This means that there exists such a set of initial metrics for which manifolds just nucleated, come back to the space-time foam. The destiny of manifold depends on its initial metric. Our main aim is to find a set of those metrics that could correspond to the evolution of our Universe and hence validate the issue iii)

Let us find the asymptotic solution of the expanding sub-spaces in the form

[TABLE]

In this case we can strongly simplify the equations of motion (III), (III), (III) which are transformed into the system

[TABLE]

[TABLE]

[TABLE]

Only two of the three equations (III, III, III) are independent (the combination $H_{1}\cdot$ Eq.(III) +$H_{2}\cdot$ Eq.(III) $-(H_{1}+H_{2})\cdot$ Eq.(III) = 0 is an identity). These equations have two different solutions. The simplest one is characterized by equal asymptotes. In this case, equations (III), (III) and (III) are reduced to a single equation

[TABLE]

with the solution

[TABLE]

In this paper, we are interested in the dependence of solutions to system (III), (III), (III) on the initial conditions. In this particular case, the initial conditions leading to asymptotic behaviour (17) are found. An example of solution with such asymptotes is represented in FIG. 1.

The choice of another set of physical parameters can change the picture. Sub-spaces grow exponentially at different rates ($H_{1}\neq H_{2}$) as is represented in FIG. 2. In this case, as well as in the previous one, both sub-spaces appear to be large at the present time and hence such growth is inadmissible.

According to the observations, the most promising case is iii) - one of the sub-spaces expands while the other remains constant:

[TABLE]

Asymptotic regime (18) can be obtained by substituting these expressions into equations (III), (III), (III). This leads to the algebraic system

[TABLE]

which is used for the determination of the Hubble parameter $H$ and the asymptotic size $e^{\beta_{as}}$.

The results of numerical calculation of equations (III), (III), (III) are represented in FIG. 3. One of the sub-spaces grows exponentially while the second sub-space remains constant. The solution does not vary for a long time and looks stable. The results of numerical calculations with slightly different initial conditions $\beta(0)$ are shown in FIG. 4 (the first sub-space is growing permanently and is not represented in the figure). There are three types of metric behaviour depending on the initial value of $\beta(0)$ - after some time, the space is growing; the subspace tends to stay constant; the subspace shrinks and finally returns to the space-time foam (left panel).

The stability region exists (left panel FIG. 4) but is very small, $\delta\beta\simeq 10^{-11}$ for the chosen set of parameters (right panel). The quantum fluctuations could easily push out the solution from the stability region and the second sub-space passes to the expanding or shrinking regime. Therefore, it is worth seeking for another set of physical parameters with more appropriate stability conditions. The successful result is shown in FIG. 5. The size of one sub-space is expanded while the metric function of the second sub-space remains constant. Such solutions are realized in a broad interval of initial conditions.

In this section we have shown that the behavior of the extra spaces depends on initial conditions. We have found out promising conditions and the Lagrangian parameters for which one of the extra spaces remains small while the other expands quickly. We also have found initial conditions that lead to the expansion of both extra spaces. This case, unlike the previous one, contradicts observations. The following sections are devoted to studying the model at the electro-weak and inflationary energy scales. In this case, we face the problem of choice of the Lagrangian parameters $a,k,c$ at the low (intermediate) energy scale, $M\lesssim 10^{13}$ GeV. Indeed, the quantum corrections to the parameter values are known to be extremely important at the Planck scale. Parameter values at high energies could differ drastically from those at low energies. It means that we can not use the values obtained earlier. Instead, we are able to obtain information from the observations. The metric of larger subspace is flat with great accuracy in the present time due to exceedingly small value of the dark energy density. Also, the Hubble parameter during the slow rolling inflation is more or less firmly established, $H\sim 10^{13}$GeV. This information is used below to impose restrictions to the Lagrangian parameters.

It is assumed that the parameters $a(M),k(M),c(M)$ vary slowly for the scale $M$ smaller than the scale of inflation and we will neglect this dependence below.

IV Subspaces evolution at intermediate and low energies

Previous study has shown that one can choose such parameter values and initial conditions for which one extra space exponentially grows and the other one remains constant starting from the sub-Planckian energies ($\sim 10^{18}$GeV). In this case, the beginning of the inflation is characterized by the two 3-dim sub-spaces - one of large and the other of extremely small size. The latter is treated as the extra space in the following.

The main goal of this section is to find those parameter values for which the extra space remains stable and its size is small. Restrictions to the parameters included in the Lagrangian (3) can be imposed if we suppose that the sub-space $V_{3}$ is (almost) flat at present time due to extremely small value of the cosmological constant. Also, the knowledge of the Hubble parameter during the inflationary stage gives additional restriction which is discussed below.

IV.1 Parameter values at low energies

There are several model parameters in action (3) and some of them can be fixed by the low energy physics. It was shown in the previous Section that two subspaces $V_{3}$ and $W_{3}$ can evolve differently starting from the sub-Planckian scale of energy. It is assumed here that subspace $V_{3}$ acquires large volume and the subspace $W_{3}$ remains small up to the beginning of the inflation. This choice breaks the equivalence of these two sub-spaces and we may follow the method elaborated in Bronnikov and Rubin (2006) to facilitate analysis.

Metric (5) of our 7-dim space leads to the Ricci scalar in the form

[TABLE]

We assume that inequalities

[TABLE]

hold. This follows from the facts that the Ricci scalar $R_{4}$ of the sub-space $T\times V_{3}$ is small as compared to the Ricci scalar $R_{3}$ of the compact sub-space $W_{3}$ and the function $\beta(t)$ varies slowly during the inflation.

Expressions (21) and (22) validate the Tailor decomposition of the function $f(R)$ in action (1)

[TABLE]

Here $\epsilon^{2}=R_{4}/R_{3},v_{3}=2\pi^{2}$ is the volume of $3$-dimensional sphere of the unit radius.

It is more familiar to work in the Einstein frame. To this end we have to perform conformal transformation

[TABLE]

of the metric describing the subspace $M_{4}=T\times V_{3}$. That leads to the action in the Einstein frame in the form Bronnikov and Rubin (2006)

[TABLE]

[TABLE]

where our physical intuition works properly. In the modern epoch, the field $\beta$ is settled in a potential minimum.

The observable Planck mass is

[TABLE]

in the units $m_{D}=1$. The $D$-dim Planck mass is slightly smaller than the 4-dim Planck mass which doesn’t contradict the limit $m_{D}>10^{13}$GeV obtained in Nikulin and Rubin (2019).

The average metric of our Universe is the de Sitter metric with the extremely small Hubble parameter. This permits us to approximate the metric by the Minkowski one that strongly facilitate the analysis. Some relations to the model parameters can be imposed in this case. The first condition

[TABLE]

supplies the energy density of the Universe be zero. The inequalities

[TABLE]

are needed for the stability reasons. We also assume that the curvature of extra space is positive

[TABLE]

In the following, we put $b=1$ in $f(R)=aR^{2}+bR+c$ without the loss of generality. The algebraic equations (30) together with definition (27) give a position of the potential minimum

[TABLE]

and connection between the Lagrangian parameters valid for the Minkowski metric of the space $M_{4}=T\times V_{3}$

[TABLE]

The first inequality in (31) gives

[TABLE]

The second inequality in (31) leads to the following expression:

[TABLE]

and hence is true for any values of the parameters.

Finally, we have two inequalities

[TABLE]

that lead to $f^{\prime}(\phi_{0})>0$ and formulas (33), (34) which should be taken into account below.

We will continue to use the Gauss-Bonnet model, for which $c_{1}+c_{2}=-3k$, that is,

[TABLE]

The parameter "$c$" is fixed by (34), the Ricci scalar of the compact static extra space is known, see (33) and (28).

IV.2 Moderate energies. Inflation

We have obtained relations (37) between the parameters, connection (33) and expression (34) at the energy scale $M_{low}=M_{ew}\sim 10^{2}$GeV. Another restrictions can be obtained at the energy scale $M_{infl}\sim 10^{13}$GeV where the inflation takes place. We will assume that the parameters of our model varies slowly within the energy range $0-M_{infl}$ which is the widespread approach for inflationary models. This means that the restrictions obtained in the previous section hold at the inflation stage.

The idea is to find those parameter values that reduce potential (27) to the form suitable for the inflationary scenario. The Ricci scalar $R_{3}$ plays the role of the inflaton $\phi$. The potential should have at least one minimum which is responsible for the (re-)heating after 50-60 e-folds during the inflation.

It is known that the inflation is finished when equality

[TABLE]

and

[TABLE]

are true Morris (2001). The range of the Lagrangian parameters can be limited from approximate equalities

[TABLE]

Here, the expression (29) is used and $\phi_{end}$ is the inflaton value at the end of the inflation. We remind that the volume of the subspace $V_{3}$ is much greater than of the subspace $W_{3}$ at the beginning of the inflation.

The preliminary simulations indicate that appropriate form of the potential is realized for the parameter values $k(M_{low})=200,a(M_{low})=150$, see FIG. 6. The parameter $c$ is fixed by condition (34).

The inflation is finished at $\phi_{end}\simeq 0.0095$. This value was found by the solution of equation (42). The present horizon arises when the field $\phi$ slightly differs from $\phi_{end}$. At the present time, the inflaton is placed in the potential minimum

[TABLE]

One can conclude that the field moves near the potential minimum during the inflationary stage and after it. Keeping in mind the connection (28) between the inflaton and the Ricci scalar of the compact extra space it can be concluded that this extra space is stable.

What about inequality (22)? It is certainly true nowadays when the curvature is as small as $R_{4}\sim 12H_{0}^{2}\sim 10^{-60}M_{P}^{2}$. During the inflation, inequality (22) also remains true. Indeed, for the chosen parameter set $R_{4}=12H^{2}\sim 10^{-11}M_{P}^{2}\sim 10^{-10}$ is much smaller than $R_{3}=\phi_{0}=0.01$ during inflation.

Our analysis was performed on the basis of the Einstein frame. Meantime, there is no firm opinion which frame is realized in the Nature the Einstein frame or the Jordan one. In the latter case, we have to turn to equation (23) after the evolution of the extra space is finished and the Ricci scalar $R_{3}=\phi_{0}$. The expression for the 4-dim Planck mass can be obtained by equating the value $M_{P}^{2}/2$ to the multiplier to $R_{4}$

[TABLE]

Here expression (28) at the potential minimum $\phi=\phi_{0}$ is taken into account and the parameter $m_{D}$ is restored. For chosen parameter values $M_{P}\simeq 700m_{D}$ and previous estimations based on the Einsten frame remains the same.

We conclude that the compact space volume could remain small enough during the whole period of its evolution - from its nucleation at the (sub-)Planckian energies up to the modern epoch.

V Quantum fluctuations and stability of extra space

In this section, we shortly discuss the role of quantum fluctuations on the stability of the extra space metrics discussed above. We start with the low energy scale $M_{low}$ where the present horizon is formed. The quantum fluctuations of the scalar field $\phi$ have been intensively studied Kofman and Linde (1987); Starobinskij (1992). The common conclusion is that in spite of their smallness they are the reason of the large scale structure formation in our Universe. The amplitude of the fluctuations during the slow roll inflation is of the order of the $\delta\phi\sim H/2\pi\sim 10^{-7}M_{P}$. They are also very small in the $m_{D}$ units: $\delta\phi\sim 10^{-6}$.

In this research, the Ricci scalar $R_{3}(\beta)$ plays the role of the scalar field $\phi$ - the inflaton - with the potential (27). Estimations made above indicate that the inflationary stage is performed near the bottom of the potential far from the potential extrema, see FIG. (6). Therefore we may not worry about the influence of the fluctuations on the final state of inflaton. If the initial value of the inflaton is to the right of the maximum, it will inevitably move to nonzero minimum. This indicates the stability of the metric evolution $\beta(t)=-\frac{1}{2}\ln[R_{3}(t)/6]=-\frac{1}{2}\ln[\phi(t)/6]\rightarrow-\frac{1}{2}\ln[\phi_{end}/6]$ with respect to the radial quantum fluctuations at the inflationary scale of energy.

The problem arises when we shift the scale up to the sub-Planckian energies where the quantum fluctuations $\delta\beta$ are of the order of the unity. Their role is twofold. If the region of stability is small, the quantum fluctuations could easily break the stationary behaviour of the extra space. The latter starts expanding or shrinking so that the probability of staying in a stable region is very small. On the other side, suppose that there exists the set of parameters $a,k,c$ leading to quasi stable classical solutions with slowly expanding sub-spaces $W_{3}$. Quantum fluctuations are able to turn its metric back to a stationary regime in some causally connected domains of the large sub-space $V_{3}$. Such domains could survive up to the beginning of the inflation.

VI Conclusion

The appearance of manifolds with different metrics as a result of quantum effects at high energies is a well known paradigm. After their creation, some manifolds evolve classically. The originated metrics serve as the initial conditions for their subsequent classical evolution. The measure of any metric originated from space-time foam is assumed nonzero though uncertain due to the absence of the Theory of Quantum Gravity.

The analysis performed in this paper indicates that there are several regimes of sub-spaces evolution at sub-Plankian energies. There are regimes characterized by the expansion of both sub-spaces at equal rates as well as at different rates depending on the Lagrangian parameters. We also have shown that some sub-spaces come back to the space-time foam.

The observable fact is that the only one 3-dim sub-space is large. It is those space where the modern physical processes are performed. Therefore, the regime characterized by only one growing sub-space is of most interest. We have found that such a regime is realized at the highest energies for specific values of the Lagrangian parameters $a(m_{D}),k(m_{D}),c(m_{D})$ and the specific initial metrics.

The parameter values of the Lagrangian $a(M_{infl}),k(M_{infl}),c(M_{infl})$ at the inflationary scale $M_{infl}$ were also discussed on the basis of chaotic inflation. They do not coincide with those at high energies $m_{D}$ due to uncontrolled quantum corrections at the sub-Planckian energy scale.

Shortly, the general picture is as follows. Sub-spaces are nucleated with different initial metrics. There is a class of multidimensional models with the higher derivatives for which some of the sub-spaces form pairs evolving classically in proper manner - one of the sub-spaces expands while the other remains constant.

Acknowledgments

The work is performed according to the Russian Government Program of Competitive Growth of Kazan Federal University and MEPhI Academic Excellence Project (contract № 02.a03.21.0005, 27.08.2013). The work of A.P. was also supported by the Russian Foundation for Basic Research Grant No 19-02-00496. The work of S.G.R. was also supported by the Ministry of Education and Science of the Russian Federation, Project № 3.4970.2017/BY.