Analytic Solution of the Starobinsky Model for Inflation

We prove that the field equations of the Starobinsky model for inflation in a Friedmann-Lema\^{\i}tre-Robertson-Walker constitute an integrable system as the field equations pass the singularity test. The analytical solution in terms of a Painlev\'{e} Series for the Starobinsky model is presented for the case of zero and nonzero spatial curvature. In both cases the leading-order term describes the radiation era provided by the corresponding higher-order theory.


INTRODUCTION
In the so-called modified/extended theories of gravity [1] new dynamical terms, of geometric origin, are introduced which force the evolution of the gravitational field equations in order to explain various phenomena which were raised by recent observations [2,3]. However, in the modified gravitational theories the new terms increase the complexity of the field equations and even in the simplest models, such as that of an isotropic and homogeneous universe, the existence of an analytical solution is not obvious. Although numerical methods can be applied to approximate the evolution of the field equations that is not sufficient for the complete study of a theory; while the analysis of the critical points it is not sufficient to provide us with information for the evolution of a system far from the critical points. Consequently the existence of analytical solutions for the field equations has lead to the application of various techniques from the analysis of dynamical systems for the study of the integrability 1 .
One of the simplest modifications of the Einstein-Hilbert Action which consider quantum corrections is the Starobinsky model of inflation [9] with Action Integral where R 2 describes the quantum-gravitational effects in the early universe and L m is the Lagrangian of the matter source. The latter Action Integral corresponds to the family of the so-called quadratic theories instance [10][11][12].
The gravitational field equations are of fourth order and in the case of a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) universe with line element 2 are calculated to be and where ρ m , p m are the energy density and the pressure of the matter source. In the case of the vacuum the field equations (3)-(4) admit an unstable (special) nonsingular solution [9]. Moreover it is important to mention that, when the equation of state parameter of the matter source, ρ m , p m is that of an ideal gas, i.e., p m = (γ − 1) ρ m , then equation (4) can be integrated to (3), while in general the conservation lawρ m + 3H (ρ m + p m ) = 0 exists. The Action Integral (1) corresponds to the f (R)-theories of gravity [13], where f (R) = R + qR 2 , while a more general consideration of the Starobinsky model is the 3 f (R) = R + qR n theory 4 [28] or [29]. Exact solutions of power-law f(R) theories can be found in [30,31]. The importance of the Starobinsky model (1) is that it provides an inflationary scenario which is favored by the observations [32]. Furthermore it has been shown that various models of inflation are identical with the Starobinsky model when the inflationary phase takes place [33] whereas the Mixmaster universe provides nonchaotic trajectories [34].
The introduction of a Lagrange Multiplier in f (R)-theories [35] can be used to reduce the order of the theory from a fourth-order to a second-order theory by increasing at the same time the number of degrees of freedom [14]. In particular a new field is introduced which is a equivalent to that of a Brans-Dicke scalar field with zero Brans-Dicke parameter [15,16] the so-called O'Hanlon theory [36]. Therefore the field equations in a FLRW background form a two-dimensional canonical Hamiltonian system which describes a particle moving in a flat space while the potential which forces the evolution of the system is related with the form of the f (R) function. Because the scalar field description is that of Brans-Dicke field the theory is defined in the Jordan frame. Hence under a conformal transformation a minimally coupled scalar field is defined and the theory is defined now in the Einstein frame 5 . Therefore the Starobinsky model can be seen as a mechanics to provide a minimally scalar field [39] to drives the inflationary phase of the universe, for a review see [40] By using the property that the field equations describe a canonical Hamiltonian system various functions f (R) have been determined in which the field equations admit conservation laws which are linear or quadratic in the momentum [8], while recently in [41] it was found that the cosmological model f (R) = R + qR n passes the singularity test and is integrable for some values of the power n. However, the case n = 2, which is that of the Starobinsky model has been ruled out and the main reason is that for n = 2 the field equations admit singular special solutions following from the the R n term; in contrast to the Starobinsky model in which the R 2 term provides a nonsingular solution as mentioned above.
A specific f (R) theory provides a de Sitter universe if there exists R = R 0 such that the Barrow-Ottowill [42] condition holds It is straightforward to see that for arbitrary R 0 , that is, R 0 → R, the latter condition can be seen as a first-order differential equation with solution the quadratic function f (R) = f 0 R 2 , where f 0 is a constant of integration.

INTEGRABILITY OF THE FIELD EQUATIONS
In the case of the vacuum the field equations in f (R) = R 2 theory are and The two equations are not independent and derivation of (6) gives the fourth-order equation (7). There are various ways in which equation (6) can be written as a first-order ordinary differential equation 6 . If we select the new dependent variable w = 1 u du dv and independent variable v, where u =ȧ, v = a then equation (6) becomes the following Riccati equation where w 0 is a constant of integration. Therefore it follows that H(t) where H (t) =ȧ a , and for initial conditions such that v 0 = 0 provides the closed-form solution a (t) ≃ t −1 . That is not the unique case. In order to see that consider now the new variables {x, y} = H (t) , d dt (H (t)) equation (6) equation (6) takes the form of the linear equation , that is, a (t) = a 0 (t − t 0 ), which is an ideal gas solution which mimics radiation solution, while it is a singular (special) solution. This singular solution is used below in order to prove the integrability of the Starobinsky model. The existence of the radiation solution it is not a surprise in the sense that f (R)-gravity can provides always a radiation epoch in the evolution of the universe [43]. However the radiation solution have been investigated before in a higher-order theory which include the Starobinsky term as also other terms follows from the Gauss-Bonnet invariant in [54][55][56]. Moreover the radiation solution in quadratic theories has been found that can describes a past isotropic singularity for the Bianchi I universe [57].
The method that we apply is that of the singularity analysis and specifically we follow the ARS algorithm [44][45][46]. Singularity analysis is a powerful method which has been applied in cosmological studies for the reconstruction of the analytical solution of various models [41,[47][48][49][50]. We omit the properties of the singularity analysis and we refer the reader to the extended review [51].
We continue by firstly applying the method for the quadratic theory f (R) = R 2 and consider now equation (7). We find that the leading-order behavior is the power-law solution a (t) = a 0 τ 1/2 , where τ = t − t 0 and t 0 denotes the position of the singularity. The application of the ARS algorithm provides the resonances to be s 1 = −1 , s 2 = 0 , s 3 = 3 2 and s 4 = 5 2 ,which means that the analytic solution is expressed by the Right Painlevé Series [52] a (t) = a 0 τ where the constants of integration are the a 0 , a 3 , a 5 and the position of the singularity t 0 . However, with the use of (6) we find that a 5 = 0, while the calculation of the first coefficient constants gives the solution to be a (t) = a 0 τ 1 2 + a 3 τ 2 + 19 32 For the field equations of the Starobinsky model we apply the same algorithm and we find the same resonances as those of the quadratic model, which means that the analytic solution is given by expression (10) or specifically by calculation the first nine coefficient constants the solution is where the constants of integration are again the coefficients a 0 , a 3 , a 5 and the position of the singularity t 0 , while the constraint equation (3) gives that a 5 = 0 or, if we assume the existence of a dust fluid, that is, p m = 0 and ρ m = ρ m0 a −3 , it follows that ρ m0 = 315 2 qa 5 (a 0 ) 2 . While in the latter scenario it is important to mention that the term a 0 t 1/2 describes the leading-order behaviour.
From the values of the resonances it is easy to see that the radiation solution is an unstable solution 7 , while the field equations of the Starobinsky model for inflation in a spatially flat FLRW spacetime pass the singularity test and are integrable.

DISCUSSION
Singularity analysis is a powerful method to study the integrability of dynamical systems. However, it has a basic disadvantage in that it is coordinate dependent. That is the reason that the Starobinsky model did not pass the singularity analysis in the consideration of [41]. The reason is that in the space of variables {a, R}, in which usually the f (R)-gravity is referred, the leading-order behaviour a (t) = a 0 t 1 2 , provides a singular behaviour for only one of the dynamical variables while for the Ricci Scalar it is a constant. However, we overpassed that problem by working directly on the fourth-order differential equation and without using the Lagrange Multiplier.
We now consider the case of nonzero spatially curved spacetime. Hence for the Action Integral (1) the field equations are derived to be where for the matter source we assumed that of a dust fluid. We apply the ARS algorithm and we find that the solution is expressed again by the Right Painlevé Series (11) where now the coefficient constants depend also upon the curvature k. For instance the first terms of the solution are where from (13) follows ρ m0 = 315 2 qa 5 (a 0 ) 2 + 30a 3 qk. We conclude that the Starobinsky model for inflation in a FLRW spacetime with or without spatial curvature it is an integrable system. Last but not least from the singularity analysis we found that the radiation era is described by a unstable point which is in agreement with the dynamical analysis for a higher-order theory [54,55].