Sub-Planckian scale and limits for $f(R)$ models

We study the Universe evolution starting from the sub-Planckian scale to present times. The requirement for an exponential expansion of the space with the observed metric as a final stage leads to significant restrictions on the parameter values of a $f(R)$-function. An initial metric of the Universe is supposed to be maximally symmetric with the positive curvature.


I. INTRODUCTION
It is generally believed that our Universe originated from Planck energies and evolved by expanding and cooling to its present state. The initial stage of quick expansion starting from the sub-Planckian energy density seems inevitable. We regard the sub-Planck scale as the highest energy scale in which classical behavior can dominate. The Planck scale is characterized by complete dominance of quantum fluctuations. The spontaneous creation of an inflationary universe is described in detail, for example, in Reference [1]. At the same time, the effects associated with the quantization of gravity may be responsible for model parameter alternation if the energy scale is large enough. Additionally, the gravity quantization leads to a nonlinear geometric extension of the Einstein-Hilbert action. The first and most successful formulation of the inflationary model, the Starobinsky model [2], considers nonlinear geometric terms belonging to the f (R) class of theories. Gravity with higher derivatives is widely used in modern research [3], despite the internal problems inherent in this approach [4,5]. Attempts were made to avoid Ostrogradsky instabilities [6], and f (R)-gravity was one of the simplest extensions of Einstein-Hilbert gravity free from Ostrogradsky instability. A necessary element of such models is the fitting of the model parameters to reconstruct the Einstein-Hilbert gravity at low energies [7]. For example, in Reference [8], the authors reconstructed the form of the function f (R) using the boundary conditions imposed on the scale factor so that it satisfied the observations in the early and late stages of the evolution of the universe. A variety of ways to study the nonlinear multidimensional gravity was discussed in Reference [9].
A wide variety of functions f (R) are presented in the literature. As explicit examples, it is worth citing a couple of functions that relate to a wide range of f (R) functions. The is considered in Reference [10]. This model unifies the early time inflationary era and the late time acceleration of the universe expansion. The authors investigated the viability of the model and obtained corresponding constraints on free parameters. The Tsujikawa model [11] is in agreement with the cosmological observations [12] but is slightly different from the ΛCDM model predictions. The function f (R) chosen there has the following form: Other attempts were undertaken to describe the whole period of evolution of the universe with the multiparametric f (R) function; see References [13][14][15][16][17][18][19].
The uncertainty in the parameter values is one of the common questions for such models.
In this article, we discuss the restrictions on the parameters of the following models based on the known behavior of the scale factor starting from the sub-Planck scale: In the last three cases, we neglect the cosmological constant compared to the energies we deal with and use a 0 = 0.
We assume that quantum fluctuations nucleate compact Planck-sized manifolds. Here, we rely on the quantum field theory, where a quantum transition is usually suppressed exponentially by a volume of nucleated systems. As the spatial part of the considered fourdimensional metric, we choose the metric of the three-dimensional sphere as the simplest representative: Other metrics are also nucleated on equal footing, and we plan to study some of them (compact hyperbolic and torus metrics) in the future.
Constraints on the parameters of the considered models of f (R) gravity, under which exponential growth of the scale factor is possible, are investigated. It is also necessary to determine the conditions under which the exponential growth of the scale factor is replaced by the observed stage of slow expansion. The parameters of the model are also limited by the condition that the current size of space must exceed the visible size of the universe.
During our study, we kept in mind the following issues: -the requirement of model stability, i.e., f (R) > 0 and f (R) > 0; -the quick growth of the space size. It must exceed the size of the visible Universe, ∼ 10 28 cm; -extremely small space expansion at the present time.
These requirements are in addition to those usually imposed on the models by the observations at low energies, in particular, inside the solar system.

II. BASIC EQUATIONS
Consider the theory described by action: The corresponding extended field equations are as follows: This system of equations coincides with Einstein's field equations for f (R) = R. Throughout this paper, we use the conventions for the curvature tensor σν Γ σ µα and the Ricci tensor is defined as R µν = R α µαν . Let us suppose that the action and metric have the forms (4) and (3) consequently. In this case, the nontrivial Equation (5) acquirse the following form: where Equation (6) correspond to the (tt)-component and Equation (7) corresponds to the coinciding components (xx) = (yy)=(zz) of system (5). The definition of the Ricci scalar for metric (3) is Substitutingα from (8) into the Equation (6), we obtain an equation that does not contain the second derivatives of the functions α and R: There are three Equations (6)-(8) with respect to the unknown functions α(t) and R(t), but only two of them are independent. It is technically easier to solve Equations (7) and (8).
Equation (9) plays the role of a restriction to the solutions of second-order differential Equations (7) and (8). This equation was used twofold. Firstly, this equation should be the identity when the solution of systems (7) and (8) are substituted. Secondly, applied at t = 0, it was used to fix one of the initial variables.
We look for those solutions to this system of equations that have "correct" asymptotic behavior. The latter are those that could describe our universe at present time. The space size should be not smaller than the size of the universe. Therefore, space should expand extremely quickly, at least in the beginning. The asymptotic value of the Hubble parameter should not be bigger than the observable one. Due to its smallness, compared to the sub-Planckian energies, we use H Several solutions of this equation could take place for specific values of the physical parameters of function f (R). The Ricci scalar averaged over large scale is negligibly small at present time. Therefore, our aim is the asymptotic solution R c = R U niverse 0.
Restriction (9) is used to fix the initial value of the curvature R(0) = R 0 .
We are interested in the dynamics of the maximally symmetric manifold starting from the sub-Planck scale. Therefore, the natural choice of the initial conditions is Further, we work in the Planck units, m P l = 1.
The sections below describe the rate of space growth for several forms of the f (R) function depending on the initial data and physical parameters.

A. R Gravity
In the case f (R) = R − 2Λ, we have a well-known solution: The observations indicate [29] that the parameter Λ ∼ 10 −122 in the Planck units. Therefore, the initial size of the manifold nucleated is of the order 10 61 . The nucleation probability of such a huge volume due to the quantum effects is negligible. Therefore, this model does not satisfy the considered assumptions.

B. R 2 Gravity
For a well-studied model the asymptote of the curvature is zero (R c = 0), which is the solution of Equation (10).  The model predicts the substantial growth of the space size that looks quite evident. Nevertheless, the model parameter a 2 is too large to be natural. This means that there are some processes that occur above the inflationary scale that strongly influence the parameter value. It is the nontrivial subject of future research.

C. R 3 Gravity
Our next choice is the function There are three types of asymptotes following from algebraic Equation (10) The first one is realized in our Universe. Our immediate task is to find the solutions to Equations (6) and (7)  are not acceptable since the other conditions are not satisfied (the solution R c = 0 is not stable and/or the space growth is too slow). The boundaries are smooth due to a possible variation in the initial conditions in Equation (12). It is assumed that the manifolds are nucleated due to the quantum effects at the sub-Planckian scale so that the probability of large sized manifolds is negligible. The part of the acceptable parameter region derived in Reference [31] is marked by a black dashed line (the right panel). The common area belonging to both restrictions is much less than each of them. Appropriate results can be obtained not only for the trivial initial conditionṘ(0) = R 1 = 0; see as an example Figure 3 withṘ(0) = R 1 = 0.01.  and a 2 = 10 9 (right side); and the initial conditions α 0 = 2.3, α 1 = 0.1, R 0 = 0.29, and R 1 = 0.01.
As a final example, consider the function with the most realistic estimation of the parameter a 2 ∼ 10 9 according to the discussion in Section III B. A typical behavior of the metric is shown in Figure 4. The phase diagram for parameters a 4 and a 3 at fixed a 2 ∼ 10 9 is presented in Figure 5.

IV. CONCLUSIONS
In this paper, we discuss new restrictions imposed on the parameters of some f (R) models of the gravity. These restrictions are the result of studying the Universe evolution at high energies. We suppose that our Universe was nucleated with the size of the Planck scale order. It must expand rapidly to reach a size no smaller than that of our Universe at present time. We also choose the 3 dimensional spherical metric from the very beginning as the additional assumption.
These suppositions being quite natural lead to new restrictions compared to limits based on the observations in the Solar system. For example, the parameter range of R 3 gravity is severely tightened if we apply both our restriction and those in the paper of [31]. In all models discussed here, the parameter ranges depend on the initial conditions that lead to their slight uncertainties. Nevertheless, these restrictions should be taken into account in the considered models based on gravity with higher derivatives. It is worth mentioning that the pure Einstein-Hilbert gravity with the Λ term is not realized in the framework of our approach.