Classical evolution of subspaces

We study evolution of manifolds after their creation at high energies. Several kinds of gravitational Lagrangians with higher derivatives are considered. It is shown analytically and confirmed numerically that an asymptotic growth of the maximally symmetric manifolds depends strongly on their dimensionality. A number of final metrics describing our Universe is quite poor if we limit ourselves with a maximally symmetric extra space. We show that the initial conditions can be a reason of nontrivial solutions (funnels) and study their properties.


Introduction
The compact extra spaces is widely used idea. Their inclusion into physical theories helps to move forward on such issues as the grand unification [1,2], neutrino mass [3], the cosmological constant problem [4,5] and so on. Any multi-dimensional model has to lead to the effective 4-dim theory at low energies. This would imply relations between the observable fourdimensional physics and a metric of the higher dimensions.
One of the question remaining not clarified yet is: why specific number of dimensions are compactified and stable while others expand [6][7][8]? Which specific property of subspace leads to its quick growth? There are many attempts to clarify the problem, mostly related to introduction of fields other than gravity. It may be a scalar field [6,9] (most used case), gauge fields [10]. A static solutions can be obtained using the Casimir effect [11] or form fields [12,13]. Sometimes one of the subspace is assumed to be FRW space by definition [14]. Another possibility was discussed in [15,16]: 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.
The origin of our Universe is usually related to its quantum creation from the space-time foam at high energies [17,18]. a e-mail: sergeirubin@list.ru The probability of its creation is widely discussed, see e.g. [19]. Here we are interested in the subsequent classical evolution of the metrics rather than a calculation of this probability. Manifolds are nucleated having specific metrics. The set of such metrics is assumed to be very rich. After nucleation, these manifolds evolve classically forming a set of asymptotic manifolds, one of which could be our Universe. In this paper the asymptotic set of the maximally symmetric manifolds with positive curvature is studied in the framework of pure gravity with higher derivatives. We consider models of the f (R) gravity and a more general model acting in 5 and 6 dimensions. No other fields are attracted to stabilize an extra space. We have found out that a number of asymptotic solutions is quite limited. This conclusion was confirmed both analytically and numerically. There is a set of initial conditions that lead to a common asymptote of classical solutions. In Sect. 3 we have elaborated a method for prediction the asymptotic behavior of metric judging on the specific form of the initial metric. We also study the funnel solution [20] as the result of an inhomogeneity of initial metric.
The gravity with higher derivatives is widely used in modern research despite the internal problems inherent in this approach [21]. Attempts to avoid the Ostrogradsky instabilities are made [22] and extensions of the Einstein-Hilbert action attract much attention. Promising branch of such models is based on the Gauss-Bonnet Lagrangian and its generalization to the Lovelock gravity. These models were adjusted to obtain differential equations of the second order so that the Ostrogradsky theorem is not dangerous for such models.
A lot of papers devoted to the f (R)-gravity -the simplest extension of the Einstein-Hilbert gravity. Reviews [23,24] contain description of the f (R)-theories including extension to the Gauss-Bonnet gravity. Examples of research with specific form of the function f (R) can be found in [25,26]. Most of the research assume positive curvature of extra space metric, but as was shown in [27], hyperbolic manifolds can also be attracted to explain the observable acceleration of the Universe.
In the framework of the gravity with higher derivatives, a variety of regimes with expanding three and contracting extra dimensions has been found in [28][29][30][31][32]. The power-law and the exponential analytical behavior of scale-factors are studied in [33,34]. Stability of specific extra space metrics is discussed in [5,[34][35][36]. The conclusion is that stable metrics do exist but their fraction is quite small. Here we consider a wider class of metrics depending on the initial conditions. It was found out that different initial metrics can lead to one and the same asymptotic solution.
Throughout this paper we use the conventions for the curvature tensor 2 Destiny of subspaces: exact results

Setup and classical equations
In this section we analyze the classical behavior of the extra space metrics. This can be done on the basis of two well known frames -the Jordan frame and the Einstein one, which are connected by the conformal transformation [37]. There are intensive debates on the selection of the frame that should be used for appropriate description of the Nature [38]. Our analysis is mainly based on the Jordan frame.
has been nucleated due to some quantum processes at high energies. The probability of this process is a subtle point and we do not discuss it in this paper. The entropy growth leads to evolution of subspaces, to those which are the maximally symmetric [39]. In this paper we study the classical evolution of subspaces M d 1 and M d 2 whose metric is assumed to be maximally symmetric with a positive curvature. It is supposed that manifolds are born with accidental shape. Subsequently they acquire symmetries due to the entropy growth [39]. We start our study after the process of symmetrization is finished. In this section, we consider the following action where R is the scalar curvature of a D-dimensional spacetime. This action appears to be an appropriate tool to study the behavior of the system just after its nucleation.
Einstein's equations of this theory are Using the results given in Appendix A, we can write the nontrivial equations of this system as in terms of metric (1). Here we have kept in mind According to (A13), the Ricci scalar is For calculations it is convenient to consider the Ricci scalar R(t) as additional unknown function and the definition (7) as the fourth equation. Three equations of this system [for example, (4), (5), (7)] can be solved with respect to the higher derivativesβ 1 ,β 2 ,R. Then substitutingβ 1 andβ 2 into equation (6), we obtain the equation which plays the role of restricting the solutions of the coupled second order differential equations. This can be checked, for example, by writing the set of four Eqs. (4), (5), (7), (8) as an equivalent set of (six) coupled first-order equations plus one algebraic equation. The Eq. (8) reduces to the algebraic transcendental equation, i.e., it is a constraint. The complete set of initial conditions therefore requires specifying six pieces of information, namely β 1 (t 0 ), β 2 (t 0 ), R(t 0 ),β 1 (t 0 ),β 2 (t 0 ), andṘ(t 0 ). These initial conditions are not independent due to Eq. (8). The latter will be used to derive an exact relation between these initial data.

Analysis of numerical solutions and their asymptotes
The system of differential equations is highly nonlinear so that one could expect a rich set of its solutions. In this section a set of solutions is discussed. It will be also shown that solutions and their asymptotes depend on the dimensionality of the extra spaces. To perform numerical simulation one needs to specify the form of the function f (R) and the initial values of the functions β 1 (t), β 2 (t),β 1 (t),β 2 (t),Ṙ(t) at t = 0. The initial value R(0) is found from auxiliary condition (8). Let with the parameter values Our choice being quite arbitrary doesn't contain small parameters. Now we are ready to fix the subspace dimensions and perform numerical calculations. Variation of the boundary conditions does not alter the asymptotes -see Fig. 2. This observation can be proved analytically. Suppose that the asymptotic behavior is as follows which is usual for the De Sitter metric. In this case we can strongly simplify the equations of motion. At t → ∞ Eq. gives Equations (4), (5), (6) at t → ∞ are transformed into the system of algebraic equations According to these equations, the subspaces are expanded with equal speed, If we substitute these expressions in (12), then we obtain the equation We note that the results at t → ∞ do not depend on the initial conditions and H 1 , H 2 , R are consistent with the numerical results. Nevertheless, the conclusion that the asymptotic behavior is independent of the initial conditions in their whole range is hasty. Indeed, if we change the parameteṙ β 2 (0) in the set displayed in Fig. 1 to the valueβ 2 (0) = 0.1 or smaller, the behavior of solution changes drastically. As shown in Fig. 3 the stable solution is absent.

The case d
Numerical simulation reveals the fact that both subspaces expand equally (see the middle and the right plots in Fig. 4). Moreover, the rate of expansion at t → ∞ does not depend on the initial conditions. However, in contrast to the previous case, the expansion rate is proportional to exp (Bt 2 ).
Let us check analytically that the asymptotes have the form In this case Eq. (7) gives and Eqs. (4), (5), (6) look as follows Here we analyze theories with f (R) such that which is true, for example, for theories with a polynomial function though another class of models that does not satisfy this condition exists, see e.g. [40]. For our case, expressions are true if we take into account (18). In turn, the Eqs. (19), (20), (21) acquire the simple form It gives and the final expression for the Ricci scalar (18) is as follows According to (23) The latter being substituted into (28) gives the relation between the degree of the polynomial "n" in (23) and the dimensionality d 1 and d 2 The choice made in this subsection (n = 3, d 1 = 3, d 2 = 2) satisfies this condition. The value of the parameters B 1 = B 2 = 0.00256 was obtained numerically. We conclude that interplay between the dimensionality of extra spaces and the form of the function f (R) influence the asymptotic behavior of the extra space metrics. This conclusion is true for those initial conditions that lead to stable solutions. As was shown numerically, if the derivativeṡ β 1 (0),β 2 (0) are small there are no stable solutions.

Extra space with constant volume, β 2 (t) = const
The discussion above indicates that time dependence of extra space metric is determined by the initial conditions, though asymptotes have more universal character. The latter leads to growing volumes of both extra spaces. This means that such solutions are hardly applicable to description of our Universe. A more realistic case relates to the situation when one of the subspaces (say M d 2 ) has a stationary radius e β c , while the other one, M d 1 , expands. In this case, the system of Eqs. (4)- (7) admits an analytic solution. More definitely, the combination For a given function f (R) this equation determines R(t) Substituting expressions (31) and (33) into Eqs. (4)- (7) leads to the following system of 4 equations Subtracting Eq. (35) from (37), we obtain that gives the connection for The solution of Eqs. (40) and (39) with respect to β 1 (t) is where H, R 0 and β c can be found from the last relations (40). Analytical solution (41) at t → ∞ describes the de Sitter space with nonzero Hubble parameter H . The result looks promising, but our numerical simulation indicates that this solution is instable.
In this section we have obtained the set of asymptotic metrics in the framework of f (R) gravity. Dependence on the initial conditions appears to be nontrivial. Nevertheless, two issues should be clarified. Namely, the question on the stability of these solutions remains. Connection between the initial conditions and the asymptotes of the metric is also not clear.

Destiny of subspaces and initial conditions
This section is devoted to the problems mentioned above on more realistic basis. The fact is that f (R) theory with compact maximally symmetric extra dimensions can not reproduce the Minkowski space at low energies [41]. This means that modern de Sitter stage (which is extremely close to the Minkowski metric) can be achieved by the price of ultra-fine tuning of Lagrangian parameters. Therefore we have to start with a more general form of action.
The Gauss-Bonnet Lagrangian is the appropriate starting point because of the absence of higher derivatives in the equations of motion. Lagrangian containing the Gauss-Bonnet term plus f (R) term (a function of the Ricci scalar) was used to describe the dark energy phenomenon [42]. Nevertheless we will use more general form of the pure gravitational action [43] S gen = m D−2 It is assumed that such Lagrangian is the basis of an effective theory [44]. In the following, the parameter b = 1 without the lost of generality. The Gauss-Bonnet term can be restored if c 1 = −4, c 2 = 1.
In this section we develop the approach to study simultaneous evolution of two subspaces analytically. To make nota-tions more familiar, we will denote

for our Universe) and consider metric of the form
where the indexes i, j, k, . . . refer to the d 0 -dim part of the metric, and a, b, . . . to its extra-dimensional part; (x) and (y) mark the dependence on x i and y a , respectively. The subspaces M d 1 and M d 2 are assumed to be maximally symmetric.
A substantial simplification of the field equations is achieved if we consider only those events for which the radius r 1 ≡ e 2β 1 (t) of the subspace M d 1 is much larger size than the radius r 2 ≡ e 2β 2 (t) of the space M d 2 , or in terms of the Ricci scalars We also suppose that the metrics of both spaces M d 1 and M d 2 vary slowly along d 0 coordinates as compared to the extra coordinates y a . More specifically, that is, each derivative ∂ k contains a small parameter . This choice breaks the equivalence of these two subspaces and we can follow the method elaborated in [43]. It will be shown that the larger subspace grows constantly while the destiny of the smaller subspace depends on initial conditions. Then metric (1) leads to the Ricci scalar in the form Additional inequality followed from agreement (46) means that the functions β 1 (t), β 2 (t) varies slowly. According to this inequality and relation (47) at hand we can perform the Tailor decomposition of the function f (R) in the units m D = 1. Here The Planck mass may be determined as a multiplier of the Ricci scalar R d 0 in expression (49). It is assumed that the field β 2 has been settled in a potential minimum in the modern epoch of the Universe evolution and hence β 2 = β c = const for the second subspace. In this case the observable Planck mass is Here the asymptotes β 2 (t = ∞) = β c and φ c = d 2 (d 2 − 1)e −2β c are written in the Jordan frame. It is more familiar to work in the Einstein frame. To this end we have to perform conformal transformation of the metric describing the subspace M d 0 . That leads to the action in the Einstein frame in the form [43] We have got Einstein gravity with the uniformly distributed scalar field where our physical intuition works properly. The extra space volume v d 2 plays the role of the Planck mass square in the Einstein frame in m D units. Remind for the future that the observable Planck mass (50) is written in terms of the Jordan frame. The action (52) leads to the Einstein-Hilbert equations with a scalar field. Solutions to these equations, β 1 (t) and φ 2 (t), evidently depends on initial conditions. Due to relation (55), knowledge of time dependence of φ 2 (t) means knowledge of the metric β 2 (t). Our nearest aim is to study asymptotic behavior of functions β 1 (t → ∞), β 2 (t → ∞).
The potential (54) with specific parameter values is represented in Fig. 5. We remind that the volume of the subspace M d 1 is much greater than the same for M d 2 from the beginning, see inequality (45).
There are two different types of the metrics evolution depending on initial conditions. The first case is realized when the space M d 2 is nucleated with the Ricci scalar φ 2,in from the right of the potential maximum, see Fig. 5. The final stationary state is characterized by the field in the potential minimum at φ 2 = φ c = const = 0.04. It means that the Ricci scalar R d 2 ≡ φ c = const and hence the size of extra space M d 2 is also constant, On the other hand, the constant term V E (φ c ) in the action (52) plays the role of the cosmological constant. Therefore, another subspace V d 0 is described by the De Sitter metric at t → ∞ provided that V E (φ c ) > 0. Its scale factor grows exponentially.
The second case is realized when the space M d 2 is nucleated with the Ricci scalar φ 2,in from the left of the potential maximum, see Fig. 5 and the field φ 2 starts its classical motion to zero value. That means that the size of the extra space is increasing with time according to relation (56).
We conclude that a destiny of nucleated subspaces drastically depends on initial conditions, see also discussion in [38]. There are only two ways of metric evolution if one subspace is much greater than the other from the beginning.

Fixing parameters
The action (43) contains several parametersc 1 , c 2 and those containing in the function f (R). Some connections between them should be imposed if we intend to consider this model as the ordinary scalar field acting in the Minkowski space. The first conditions supply the energy density of the Universe be zero. The inequalities are needed for the stability reasons. We also assume that the curvature of extra space is positive and the Planck mass is a real number -see (50). This leads to the additional inequalities The algebraic equations (57) together with definition (54) give position of the potential minimum and connection between the Lagrangian parameters valid for the Minkowski metric of the space M d 0 . Here we assume b = 1 for generality. Consider inequalities (58) and (59) in more detail. First inequality in (58) gives Second inequality in (58) leads to the following expression: are the result of expressions (59). The action parameters must satisfy conditions (60)-(64) for the (43) model, so that it is realistic. The numerical values of the parameters used in Fig. 5 do satisfy these conditions.

Analysis of the metric dynamic
As shown above, there are different ways of the metric evolution depending on the initial conditions. In this Section, we analyze the space expansion before the modern horizon has been appeared. Numerical values of the physical parameters are chosen such that they satisfy conditions (60)-(64). Suppose that the second subspace M d 2 was nucleated with the Ricci scalar φ 2 = φ 2,in = 0.02, see Fig. 5 while the first extra space M d 0 was nucleated having much larger the Ricci scalar. The dynamic of the field φ 2 is evident due to mechanical analogy: it will move to the potential minimum at φ 2 = φ c = 0.04. Our interest is to evaluate the character of this motion as well as the dynamic of the larger subspace M d 0 .
According to (52), the Lagrangian has the form For our estimations, it is enough to approximate the kinetic term by the constant, K E (φ 2,in ) = K E (0.02) 6 · 10 −6 and the potential by a linear function . These approximations are valid in a region around φ 2,in . In terms of new field the Lagrangian acquires the form in the vicinity of φ 2,in . Equation of motion for the scalar field is as follows Following the spirit of the inflationary ideas we use here the slow roll approximation (the second derivative in (68) is omitted) and the gravitational part is reduced to equation for d 1 = 3 with denotation H =β 1 (t). Therefore, the size of the first extra space grows as while the size of the second extra space evolves along the slope of the potential in Fig. 5. Its motion is described by the solution A remark is necessary. Up to now we study the rate of expansion in the Einstein picture. In the Jordan picture the expansion rate is different and can be found by application of formula (51)-(70)

a(t) e Ht
If the field φ 2 moves slowly to its minimum as it is in most inflationary models, exponential grows of the main space remains the same qualitatively while the size of the second extra space tends to constant. We have considered the first type of the metric evolution. The second type (the subspace M d 2 is also increasing) is realized if this subspace is nucleated with a small value of the Ricci scalar (from the left of the potential maximum in Fig. 5). To simplify our analysis, consider motion near φ 2 = 0. In this case the functions K E and V E may be approximated by the first term in the Tailor decomposition, Substitution of (55) leads to classical equation for β 2 (t) which can be solved analytically. The result is where C i is a constant depending on initial conditions. Therefore the size of subspace M 2 grows exponentially. As the result, the first subspace M d 1 exponentially grows while the second space M d 2 has two variants of evolution. If it was born with the Ricci scalar R d 2 from the left of the potential maximum in Fig. 5, this subspace is expanded (its radius is proportional to √ 2/φ 2 ). If it was born with the Ricci scalar R d 2 = φ 2 from the right of the potential maximum, this subspace volume tends to the finite value.

Funnel as a connection of two subspaces
Discussion in this section is devoted to consequences of the most promising case we have found -one extra space expands while the volume of the second one tends to constant. The particular coordinate could experience expansion or contraction depending on the initial conditions. The question is: what happens if the initial conditions are different in two areas of the manifold M?
Consider a manifold M with topology T × M 1 × M 2 × M 3 , where M 1 is 1-dimensional infinite flat space, and M 2 , M 3 are 2-dimensional spheres. We study the pure gravitational field action in the form (43) with the metric Here A(u), β 1 (u) and β 2 (u) are functions of the Schwarzschild radial coordinate u, −∞ < u < ∞. The Ricci scalar has the form with prime denoting differentiation with respect to u. For the Ricci tensor squared we have where whereas the Kretschmann scalar is as follows where We are going to describe the structure of the funnel, which is the transition between the domain with (large M 2 /small M 1 ) subspaces and the domain with subspaces (large M 1 /small M 2 ) [20]. It is implied that the Minkowski metric is realized asymptotically which imposes the following restrictions Here r 0 = e β c is the radius of extra space at u → ∞. The value of r 0 = 1/ √ −c is connected with the physical parameter c according to formulas (55), (60), (61).
Numerical simulations were performed using the Ritz method which means that metric functions A(u), β 1 (u), β 2 (u) in (74) are approximated by trial functions. Let us choose the trial functions in the form keeping in mind conditions (80), (81). Parameters ξ 1 , ξ 2 , ξ 3 are not independent. Indeed, the integrand in (43) has to decrease faster than 1/u when u tends to infinity in order to avoid divergences. The determinant grows as g ∼ u 4 so the content of square brackets in (43) should tend to zero faster than 1/u 3 . It can be obtained by setting all the coefficients of terms u n , n ≥ −3 equal zero in the Tailor series of the integrand. For terms u n , n ≥ −2 it is achieved by the very choice of trial functions (82). For the term 1/u 3 , the coefficient equals zero if the relation holds. This relation between ξ 1 , ξ 2 and ξ 3 makes them be dependent. In what follows we consider ξ 1 , ξ 3 as independent values, and ξ 2 as their function. Parameters ξ 1 , ξ 2 , ξ 3 are to be defined by the action minimization. Qualitative behavior of the extra space radii is represented in (Fig. 6). Now we proceed the search of approximate metric in the following manner. The action S(ξ 1 , ξ 3 ) appears to be a function of two variables after the substitution of trial functions (82) into action (43). Classical solutions should satisfy equations δS δg ab = 0 ( 8 4 ) or, according to the Ritz method, [45] One of the ways to solve this system is to find a minimum of the auxiliary function Note that if the trial functions are chosen absolutely correctly by accident, a minimum of the function min = 0. Probability of such accurate choice is almost zero. So we just search for the values of parameters ξ 1 , ξ 2 (ξ 1 , ξ 3 ), ξ 3 which give the minimum of , and consider the resulting trial functions as approximate solutions of (85). The described minimization procedure for (ξ 1 , ξ 2 (ξ 1 , ξ 3 ), ξ 3 ) gives the minimum at ξ * 1 = 9.96, ξ * 2 = 9.57, ξ * 3 = 1.85 (see Figs 7,8). One can also notice that the minimum is very deep. It means that the trial functions are chosen properly. Note that parameters ξ 2 , ξ 3 have concrete physical sense. The asymptote of g 00 is ∼ 1 − ξ 2 /u, whereas for weak gravitational field we have g 00 ∼ 1 − 2m/u. It means that the distant observer can "feel" the parameter ξ 2 as the mass of the object (Funnel): m = ξ 2 /2. The second parameter, ξ 3 is related to the size of the funnel transition as d = r 0 + ξ 3 . Thus it is of interest to present values of these parameters in dimensional form. Using (50) and the set of values from the capture of Fig. 5 we obtain For an external observer the funnel looks like a microscopic object with the mass of the order of the Planck scale. It does not interact with matter other than gravitationally and may serve as a dark matter candidate with nontrivial internal structure (Fig. 9).

Conclusion
General aim of this research is to analyze the role of initial conditions in asymptotic behavior of subspace metrics. On the basis of pure f (R) gravity we have shown that asymptotic behavior is the same for a wide set of the initial conditions. Exact asymptotes of the subspaces where obtained in an ana- In this case the extra space sizes grow as e Bt 2 at t → ∞. We also obtain analytic result for asymptote of one extra space -formula (81) -provided that the size of another one is constant. This result is of purely academic interest. Indeed, as was shown in [41], 4dim Minkowski space-time is incompatible with maximally symmetric compact extra spaces of nonzero curvature in the framework of pure f (R) gravity.
The results based on more complicated model (43) are more promising. The 4-dim Minkowski space-time and maximally symmetric compact extra spaces with positive curvature could coexist. The price is connection (61) between the parameters of the Lagrangian which may be considered as the strong fine tuning of the model. Three kinds of the solutions are found. One of them is characterized by expansion of both subspaces while the second solution describes expansion of one extra space and stabilization of the other. The third solution takes place when the initial conditions are different in separate points of a manifold. In this case, the nontrivial solution (funnel) could be formed. For observers, it looks like a point-like object of the order of the Planck mass and size. They can be considered as the dark matter candidates.
One can conclude that a number of final states of metric describing our Universe endowed by extra dimensions is quite poor if we limit ourselves with a maximally symmetric extra space. Inclusion of inhomogeneous extra spaces [41,46,47] improves the situation significantly. tion and Science of the Russian Federation, MEPhI Academic Excellence Project (Contract N 02.a03.21.0005, 27.08.2013). The work of S.G.R. was also supported by the Ministry of Education and Science of the Russian Federation, Project N 3.4970.2017/BY. The work of S.G.R. and A.A.P. is performed according to the Russian Government Program of Competitive Growth of Kazan Federal University.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .