Integrability from Point Symmetries in a family of Cosmological Horndeski Lagrangians

For a family of Horndeski theories, formulated in terms of a generalized Galileon model, we study the integrability of the field equations in a Friedmann-Lema\^{\i}tre-Robertson-Walker spacetime. We are interested in point transformations which leave invariant the field equations. Noether's theorem is applied to determine the conservation laws for a family of models that belong to the same general class. The cosmological scenarios with or without an extra perfect fluid with constant equation of state parameter are the two important cases of our study. \ The De Sitter universe and ideal gas solutions are derived by using the invariant functions of the symmetry generators as a demonstration of our result. Furthermore, we discuss the connection of the different models under conformal transformations while we show that when the Horndeski theory reduces to a canonical field the same holds for the conformal equivalent theory. Finally we discuss how singular solutions provides nonsingular universes in a different frame and vice versa.


INTRODUCTION
The plethora of phenomena which have been discovered the last few years have led to the consideration of alternative/modified gravitational theories [1]. These extended theories of gravity provide additional terms in the field equations which -in conjunction to those of General Relativity -can explain the various phases of the universe. In physical science the existence of a fundamental axiom, such as Hamilton's principle, is of paramount importance. Among the theories that are generated by an action integral, those that are of second-order in respect to their equations of motion possess a distinguished role; from Newtonian Mechanics to General Relativity.
A particular family of theories which has drawn attention are the so-called Horndeski theories. Horndeski, in 1974 [2] derived the most general action for a scalar field in a four-dimensional Riemannian space in which the Euler-Lagrange equations are at most of second-order. All relevant scalar-tensor theories with this property, such as the Brans-Dicke [3][4][5], the Galileon [6,7] and the generalized Galileon [8] belong to the general family of Horndeski theories. Although the latter are of second-order, they provide non-canonical nonlinear field equations which -even for the simplest line element of the underlying space -might not be able to lead to solutions in terms of closed-form expressions. Due to the high level of nonlinearity, numerical methods are applied in order to approximate the evolution of the system. However, whether a solution actually exists is not always known. For that reason, in this work, we are motivated to study the integrability of the field equations for a class of families of Hordenski theories (or equivalently those of the generalized Galileon Lagrangian [9]).
There are various methods to study the integrability of a system of differential equations. Two of the most famous are: a) the existence of invariant transformations, i.e. symmetries and b) the singularity analysis; for a recent discussion and comparison of these two methods see [10]. Both of these have been applied widely in gravitational theories [11][12][13][14][15][16][17]. In several cases the gravitational field equations, after a specific ansatz for the line element is adopted, can be derived by point-like Lagrangians [18]. The application of Noether's theorem over the corresponding mini-superspace action has been utilized for the determination of conservation laws/analytical solutions in various models both at the classical as well as in the quantum level, for instance see [19][20][21][22][23][24][25][26][27][28][29][30] and references therein. Noether's Theorem is the main mathematical tool that we use in this study.
From the various different classes of generators of the invariant transformations we consider the most simple; the one corresponding to the so-called point symmetries. Point symmetries are the generators of the invariant transformations in the base manifold in which the dynamical system is defined. The application of Noether's theorem for point symmetries leads to conservation laws linear in the momentum. Well known conservation laws of this type from classical mechanics are those of the momentum and the angular momentum. The plan of the paper it follows.
In Section 2, we present our model which is a special consideration of the Hordenski Lagrangian and can be seen as a first generalization of the canonical scalar-tensor theories. Moreover we consider the cosmological scenario of an isotropic and homogeneous universe in which an extra fluid exists with constant equation of state parameter. For that model the field equations are calculated while the minisuperspace Lagrangian is discussed. Furthermore our model admits four unknown functions which define the specific form of the cosmological model. In Section 3, we apply Noether's theorem in the minisuperspace Lagrangian in order to specify the unknown form of the functions which define the model and derive the corresponding Noetherian conservation laws. For the completeness of our analysis we consider separately the cases with or without an ideal gas and the cases with zero or nonzero spatial curvature for the underlying spacetime. In order to demonstrate the usefulness of our results we present some closed-form (special) solutions in Section 4; which are derived from the invariant functions of the admitted symmetry vectors. Finally in Section 5 we discuss our results and specifically we discuss the relation between our models under conformal transformations and we show how a singular universe is mapped to a nonsingular universe under the change of the frame.

THE MODEL
The gravitational action integral that we consider is of the form that falls into the class of a generalized Galileon (or Horndeski) model [9,31], where S m indicates any possible additional matter content. Furthermore, for the background geometry, we consider a four dimensional space of Lorentzian signature, isotropic and homogeneous with spatially flat three dimensional geometry; that is, a Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime line element The action integral (1) is the simplest generalization of those scalar-tensor theories where the kinetic energy is not coupled with functions entailing derivatives of φ. As we can see, when g (φ) → 0, an action integral of that type, including the Brans-Dicke case, is recovered. Our goal is to derive all admissible gravitation models in (1) that possess an integral of motion of the form previously described. In that respect, note that the existence of ω(φ) is not trivial since its absorption with a reparameterization. φ = f (ϕ) = ω(ϕ) −1/2 dϕ leads to the transformation of the Laplacian of φ as (primes denote differentiation with respect to the argument) introducing a new kinetic energy squared term in action (1). This is in contrast to what happens in the case of a scalar field whose action does not involve a coupling with derivatives of φ over its kinetic term, i.e. when in our case g (φ) = 0, where without loss of generality we can always select ω (φ) to be a constant. The mini-superspace Lagrangian that we can derive with the help of the gravitation plus scalar field part of (1) and (2) by integrating out the spatial degrees of freedom reads 1 where the dot stands for differentiation with respect to time t.
In what regards the matter part, we start our investigation by considering the S m = 2 L m d 4 x contribution in (1) to be that of a perfect fluid obeying the barotropic equation of state P = γρ, with P (t) and ρ(t) the pressure and energy density respectively. In order to perform the variation of L m = √ −gρ with respect to the metric we need to a priori define a continuity equation [32]. We choose the following well known relation for an ideal gas in a FLRW space-time:ρ where ρ 0 is a constant of integration. Thus, the relevant addition to the mini-superspace Lagrangian is As a result the total Lagrangian reads It can be easily verified that the three Euler-Lagrange equations are completely equivalent to the field equations of motion of (1) for the metric and the scalar field whenever the line-element (2) is substituted and the isometries of the space-time are inherited by the matter φ = φ(t).

SYMMETRIES AND CONSERVATION LAWS
By considering the mini-superspace action to be form invariant under point transformations generated by 2 one is naturally led to the well known infinitesimal criterion where pr (1) Y is the first prolongation of Y , i.e. the extension of the generator in the first jet space spanned by (t, a, φ, N,ȧ,φ,Ṅ ). Clearly, sinceṄ does not appear in (4), we can disregard the relevant term and just write where Φ i = dξi dt −q i dχ dt , q = (a, φ), i = 1, 2. Application of (13) leads to an overdetermined system of partial differential equations to be solved for χ(t, a, φ, N ), ξ i (t, a, φ, N ), i = 1, ...3, and F (t, a, φ, N ). The former is formed by gathering the coefficients of all terms involving the derivatives of a, φ and N , to which the functions entering the generator (12) have no dependence. By gradually integrating the equations, we obtain the following results (the basic steps to the solution can be found in appendix A): • An infinite dimensional symmetry group generated by with χ(t) an arbitrary function of time. Its appearance reflects the fact that the mini-superspace action for (4) is form invariant under arbitrary time transformations. Its existence leads, trough Noether's second theorem, to a differential identity between the Euler-Lagrange equations of motion (8). The latter implies the well known fact that not all of them are independent, i.e. L tot is a constrained (or singular) Lagrangian. It is well known that when an infinite dimensional symmetry group is present (i.e. Noether's second theorem is applicable) the system is necessarily singular, however the inverse is not true [34]. In our case this particular symmetry is a remnant of a more general group, the four dimensional diffeomorphism, χ(x) µ ∂ ∂x µ , under which the full gravitational action (1) is form invariant.
• Additionally to the previous group -which always exists for time reparameterization invariant Lagrangians [35] -we obtain a symmetry generator if h(φ), ω(φ), g(φ) and V (φ) satisfy certain criteria. In particular, we see satisfies (13) whenever in which λ is a constant. Note that the corresponding gauge function appearing in (13) is trivial in the calculation, i.e. F (t, a, φ, N ) =const. and this remains that way for every case that we examine later on in the analysis. Note also, that (16) and (17) describe a particular solution of (13) in the special case when γ = −1. When the fluid contribution plays the role of a cosmological constant, it can be absorbed inside the potential V (φ) (as can be seen by the form of L tot ) and the general solution of the γ = −1 case is the same as the one that we get if we completely omit the fluid.
On the other hand, for the special case where the perfect fluid describes radiation, i.e. γ = 1/3 the situation changes and the symmetry generator assumes the general form while the functions entering the action need to be whereg(φ) is an arbitrary non-constant function. As a result, for γ = 1/3, there exists an infinite family of models belonging to the general class of actions of the form (1) that possesses an integral of motion of the type we are investigating.
3.1. Particular case: ρ0 = 0 Let us see how the situation alters if we remove the ideal gas from our considerations. In other words let us use the Lagrangian given by (4), in the criterion (13) instead of the Lagrangian L tot . Then, in addition to the diffeomorphism group characterized by (15) -which always exists for time reparameterization invariant Lagrangians [35] -we obtain a symmetry generator that satisfies (13) when with λ being a constant andg(φ) again an arbitrary (non-constant) function of φ. Once more, we have an infinite set of physically different models that are characterized by a functiong(φ). We can see that, in comparison to the γ = 1/3 case, result (21) is identical to (19) in the special case when λ = −3.
It is useful to study how the inclusion of spatial curvature (either positive or negative) may alter the conditions under which a symmetry generator appears. If we consider the line element where r 2 = x 2 + y 2 + z 2 , Lagrangian (4) is modified by the addition of an extra term in the potential and reads Thus, for the same class of point transformations we considered earlier and with the application of the infinitesimal symmetry criterion (13), it is straightforward to derive the result given by (20) and (21) under the condition that λ = −3 or equivalently result (18) and (19) corresponding to the radiation fluid case where γ = 1/3. In other words when k = 0 we have the same situation as in the γ = 1/3 case. The functiong(φ) once more remains arbitrary, under the restriction of course of not being constant. As we can see, even though the non vanishing of k imposes a restriction (λ = −3) in comparison to the k = 0 case, there still exists an infinite number of models admitting an integral of motion of the type we consider in this work. Moreover, for the combined case where we have both a non-zero spatial curvature and a perfect fluid, in other words when the Lagrangian is being given by L k + 2L m , it comes as no surprise that the existence of an integral of motion implies that γ = 1/3. In fact the result is once more exactly the same with what we see in (18) and (19).
Having derived the previous results, we can use Noether's second theorem and derive the integrals of motion corresponding to each case. It can be easily verified that where the conserved charge is while A i E i denotes a linear combination of the two spatial equations of motion. Of course, in place of L there can be either L tot of (7), L of (4) or L k of (23), depending on the generator that we use and the case that we examine. In every situation, we have on mass shell the corresponding conserved quantity (25) for each symmetry vector field Y . The fact that in every case but the generic fluid with γ = 1/3 an arbitrary functiong(φ) is involved in the generator as well as in the action itself, implies that we possess an infinite collection of models -corresponding to different sets of functions h(φ), g(φ), V (φ) and ω(φ) -which admit at least one integral of motion of this type; hence being integrable. Note here, that the fact that we are led to an infinite number of models through the arbitrariness ofg(φ) is owed to the use of the reparameterization invariant Lagrangian (4). Had we chosen to adopt the gauge N = 1 at the Lagrangian level, then only a particular case would have emerged whereg(φ) = φ κ , with the corresponding χ(t) for this symmetry being χ(t) ∝ t. Of course, due to the system being autonomous, ∂ t also exists; the latter being the only remnant of (15) after fixing the gauge. The difference appearing here in the arbitrariness ofg(φ) lies in the consideration of Lagrangian (4), which naturally generates the constraint equation E 0 = ∂L ∂N = 0 in (13) and (24), allowing the derivation of a larger class of symmetries.
The existence of an integral of motion of the form (25) is of great interest in the search of solutions. The system at hand possesses two degrees of freedom that are bound by a constraint equation. Thus, the existence of I implies that, in principle, we need only to solve the constraint equation together with I =const. in order to fully integrate the system. As a result, our problem is immediately reduced to one that involves only first order differential equations. The usefulness of I =const. also lies in the fact that it is linear inȧ, an advantage which the constraint equation does not have since it is non-linear in both velocitiesȧ andφ. Nevertheless, even for the simplest of models defined byg(φ), the situation can be highly complicated. In most of cases it may be more useful if the result is used in a way to distinguish any existing solutions that are invariant under the action of the generator, i.e. derive the characteristic of Y for a given modelg(φ) so as to deduce a possible relation between a and φ and use the latter together with the integral of motion and the constraint equation.

INVARIANT SOLUTIONS
In what follows, and in order to demonstrate the importance of our results, we present a few illustrative applications of invariant solutions that we derive in the gauge N = 1 with and without a perfect fluid. For a treatment of how the conserved quantity can be utilized to lead to more general solutions expressed in an arbitrary gauge we refer the reader to appendix B. Note that the solutions which we present here are derived with a mentality of keeping all terms inside the action, i.e. we do not express any additional solutions which may have a vanishing V 0 , g 0 or ω 0 .
The on mass shell value of the conserved quantity I for solution (26) is It is interesting to note what occurs when the theory assumes a Brans-Dicke-like form, i.e. when h(φ) = φ and ω(φ) = ω 0 φ −1 . This happens when λ = 3γ−1 3(γ−1) and it leads to a particular solution with γ = −1 (cosmological constant contribution in the action) that reads where now ω 0 = −5/3, V 0 = −2ρ 0 , while g 0 is kept arbitrary. We can observe that solution (27) satisfies the same condition as the vacuum Brans-Dicke solution (for a spatially flat universe), namely that a(t) 3 φ(t) ∝ t (in the gauge N = 1) [36]. This is explained by the fact that for γ = −1 the potential in (17) becomes a constant with the value V 0 = −2ρ 0 so that it effectively cancels the cosmological constant contribution of the fluid in (7).
Nevertheless, for an arbitrary γ, we can see how the symmetry that is present in the Brans-Dicke model -satisfying a 6 φ=constant -is generalized in the class of models we are considering here. In particular we encounter the more general relation a 1−3γ λ φ =constant, that is implied by the existing symmetry generator. Of course in our case, (26) is not a general solution of the equations of motion but a particular one, which however, as a power-law, is of special cosmological interest.

2.
In the absence of a fluid,g(φ) = e −µφ case. For a spatially flat FLRW line-element, we choose the function appearing in (21) to beg(φ) = e −µφ . After an appropriate redefinition of the constants ω 0 , g 0 and V 0 , we derive the corresponding model to be A special solution in terms of a power-law for the scale factor exists when the constants ω 0 and V 0 are given by while the rest parameters remain free. It can be seen that on this solution the conserved quantity I (the Noetherian conservation law) becomes zero.
3. Caseg(φ) = φ −µ , without a perfect fluid. Again in (21), we make the choiceg(φ) = φ −µ and obtain (once more with an appropriate redefinition of the constants V 0 , g 0 , ω 0 in the action and λ → λ/µ) the subsequent set of functions entering the action Two solutions for which the conserved quantity is again zero can be easily derived for this model.
The first solution is again a power-law where λ = 2 κ − 2µ + 1, are the constants appearing in (30). Once more we can also distinguish a Brans-Dicke sub-case that appears when µ = 0. It is interesting to observe that condition a(t) 3 φ(t) ∝ t of the vacuum Brans-Dicke model does not apply here, although it was still relevant in the corresponding model of case 1. As we can see, the power σ in the scale factor expression is not connected to the power of time in the φ(t). Only if we choose to turn the action into the pure Brans-Dicke form does (31) with µ = 0 satisfy the aforementioned condition: By setting g 0 = 0 and V 0 = 0, which immediately results in κ = 1 − 3σ or κ = −2σ. Another solution we can derive for theg(φ) = φ −µ model leads to a space-time of a constant scalar Ricci curvature, a de Sitter universe that is: when the following relations hold 4. Case k = 1,g(φ) = e −µφ , without a perfect fluid. Let us furthermore consider the case with nonzero spatial curvature k = 0. Then, from the two particular types of solutions we examined above in cases 2 and 3, only for g(φ) = e −µφ arises the following configuration: which, as long as the scale factor is concerned, it corresponds to the k = −1 solution of Einstein's gravity in vacuum, that is the Milne solution. Here, it is given for any sign of k with the difference being carried in the coupling (through ω 0 and V 0 ) that is necessary for its existence. The integral of motion -calculated with the help of (25), (23) and (18) -assumes on mass shell the value I = 2 k As is obvious, we are able -just by choosing a particular functiong(φ) -to find through (21) or (19) the corresponding class of integrable models that admit a conserved quantity of the aforementioned form. As we already stated, it is due to the existence of the constraint equation (8a) that one can in principle obtain the general solution by considering only the first order system E 0 = ∂L ∂N = 0 and I =const. However, because of the complexity of the equations this is a highly nontrivial task.

DISCUSSION
We remark that for the cases without the extra matter fluid where we have obtained an infinite set of models owed to the arbitrariness ofg(φ), there exists an interesting connection among them. These models can be mapped to each other by conformal transformations. Take for example two different models whose action is characterized by two different scalar field functionsg 1 (φ 1 ) andg 2 (φ 2 ). Then, the transformation where g (1) µν and g (2) µν are the two corresponding spacetimes, maps the one action to the other. What is more, if we try to map the relevant form of the action to the Einstein frame we can see that an interesting "degeneration" occurs 3 : The different models that we get for the variousg(φ) all collapse to a single action. Take for example the case where the initial action is characterized by the functions of φ given by (21). It can be seen to be mapped into an action of the formS by a conformal transformation in the metric and a re-parametrization of φ(Φ). The new to the old variables are related by the expressionsḡ with V 1 a constant and σ, ω 1 and ω 2 in (35) being related to the initial constants g 0 , V 0 , ω 0 and λ. We can see that the difference of this action lies in the existence of a term that is quadratic to the kinetic energy of the scalar field, while in (1) we considered a theory that has at most linear expressions in X = − 1 2 φ ,µ φ µ . However, what is also interesting here is that when V (Φ) ≃ e σΦ dominates then the contribution of the terms e −σΦ ω 1 Φ ,µ Φ ,µ¯ Φ + ω 2 (Φ ,µ Φ ,µ ) 2 ≃ 0, which means that (35) reduces to the Action Integral of a minimally coupled scalar field. The latter is exactly the limit which relates a canonical scalar field between the Jordan and the Einstein frames.
In general, great care needs to be taken when using conformal transformations. The fact that all the different models for the various choices ofg(φ) can be mapped to (35) does not make them equivalent to the latter, or even to each other (due to (34)) for that matter. Two actions in order to be physically equivalent, they need to be mapped by gauge transformations of the theories under consideration. In the case of gravitational actions, those are the four dimensional diffeomorphism of space-time. As we know, conformal transformations cannot always be attributed to coordinate changes. As a result, the gravitational space-time arising in each situation is generally different.
In order to see that consider the nonsingular solution (32) corresponding to a de Sitter universe, which came out of a model characterized byg(φ) = φ −µ . If we choose to map this model to (35) then, by virtue of (36) we obtain whereN ,ā are the new lapse and scale factor respectively, while c 1 is a constant. It is easy to see that if you go from this metric to the gauge whereN = 1 (so as to compare with what we have in (32)) by performing the time transformation we derive, with the appropriate scalings in x, y and z, the line element where ψ is a constant. What was an exponential solution with a constant Ricci scalar in the theory we are investigating, has now become a power law with R(τ ) ∝ τ −2 and a space-time with curvature singularity at τ = 0 in the system described by action (35). Hence we can see that, whenever the conformal transformation does not correspond to a general coordinate transformation, the gravitational properties of the system are bound to change and the solutions represent different geometries. For discussions on the relation between analytical solutions and physical quantities between the Jordan and the Einstein frames see [37][38][39] and references therein.
Additionally to the previous geometrical considerations, extra care needs to be taken when a fluid is used as a matter source. This is owed to the fact that you need to pre-define a continuity equation in order to derive a rule for the variation of the energy density ρ with respect to the metric g µν . In the beginning of our analysis we considered a perfect fluid that is characterized by continuity equation (5). It is a well known fact that the latter is not conformally invariant, while its solution is necessarily utilized at the level of the minisuperspace Lagrangian. Henceforth, after a conformal transformation is being made, one needs to take into account a different fluid, which is now interacting with the scalar field, in order to make such a correspondence possible. This also results in a change of the physical behaviour of the system, since the properties of the matter source need to be altered.
In a future work we plan to study the integrability of other Horndeski theories by including more terms in the action integral and for other kind of transformations which leaves the field equations invariants, such as the generalized symmetries as also to investigate the effects of the conformal transformations in Hordenski theories.
Appendix A: Calculation of the symmetry generator As indicated in the main text, application of (13) leads to an overdetermined system of partial differential equations. The latter is formed by gathering and demanding that are zero the coefficients multiplying terms involving derivatives of a, φ and N . For Lagrangian L tot as written in (7) we infer from the coefficients ofȧ 3 ,ȧ 2Ṅ andȧ 2φ that we need to set respectively The fourth order coefficientsȧ 2φ2 andȧφ 2Ṅ each imply which means that ξ 2 = ξ 2 (t, φ). However, with the help of χ = χ(t), we can get a further restriction from the coefficient ofȧφ 2 that leads to Due to (A1) and (A3) , the terms involvingφ 3Ṅ andφ 3 bring about the conditions Thanks to the above restrictions (A1), (A3) and (A4) we also get fromṄ ,φ andȧ respectively. The equation emanating from the coefficient ofȧ 2 can be solved algebraically with respect to ξ 3 (t, a, φ, N ) the latter being a consequence of the fact that no derivative of N enters in the Lagrangian. After this step one can see that the equation extracted from the coefficient ofȧφ 3 is and it can be immediately integrated to yield With the dependence with respect to a being now completely specified, the relation produced by the coefficient oḟ aφ involves only unknown functions of φ. Hence, we can now start gathering coefficients with respect to powers of a that appear inside it. With the help of a useful reparametrization g(φ) = g ′ 1 (φ) 3 one can straightforwardly obtain: where the c i 's indicate constants of integration. The only appearance of time in the remaining coefficient is inside the zero-th order component (that multiplies no derivative ofȧ,φ orṄ ) and it implies that Thus, the gauge function becomes trivial.
With the results up to here all dependence of the functions of the generator with respect to a, t and N is specified. Inside each of the remaining equations we can gather coefficients with respect to a, since now the unknown functions involve only φ. However, due to the existence of γ the way that the coefficients are to be gathered depends on its value. We can distinguish two main cases: • Case γ = 1/3. The equation produced by the zero-th order coefficient leads to We note that, since we want the most general result, we avoid any special solution that leads to vanishing of any of the functions involved in our starting action. Theφ 2 component leads to while theφ 4 gives rise to the third order equation which under a transformation g 1 = exp g 2 (φ)dφ becomes with general solution g 2 = c8 φ+c9 . Thus, the final function needed so that all equations are satisfied is At this point the system of equations is completely satisfied. By absorbing the trivial constant c 9 inside φ with a translation and with an appropriate parametrization of the rest of the constants we obtain result (16), (17). Additionally, we can observe that the function χ(t) remained arbitrary through the calculation and this explains the existence of (15).
• Case γ = 1/3. With this choice of γ the zero-th order equation implies The relations from theφ 2 andφ 4 coefficients respectively lead to It is a matter of reparametrizing the constants of integration and an introduction of a new functiong(φ) as g 1 (φ) = 1 g(φ) to obtain result (18), (19) and of course the same comments as in the previous case, for the arbitrariness of χ(t), hold.
For the sake of completeness we have to point out that apart from the case γ = 1/3 there are also other values of γ that can lead to a different gathering of terms in the coefficients of powers of a; namely γ = −1 and γ = −7/5. The first gives the same result as the case without fluid (see eqs. (20), (21)) with the sole difference of a potential that is V (φ) = V 0g λ+3 g ′ − 2ρ 0 in place of the one appearing in (21). That is one that cancels the cosmological constant role of the fluid. The second case, γ = −7/5, after appropriate reparametrizations leads to the exact same result as the generic γ = 1/3 result and that is why we do not make a separate presentation of it here.

Appendix B: More general solutions
In section 4 we presented some invariant solutions where the functional dependence between a and φ can be extracted with the help of the infinitesimal symmetry generator. Although these are particular solutions which can be derived in a rather simple manner, they are the most cosmologically interesting. Even though we have proven the existence of an integral of motion I, thus reducing the problem to solving two first order differential equations: the constraint and I =const. The acquirement of the general solution is still a very difficult task, especially due to the non-linearity of the constraint equation in both derivatives involvedȧ andφ.
Here, we want to exhibit a method with which the integral of motion I can be used to derive the solution for the spatially flat case in the absence of a perfect fluid, when I = 0. Some of the solution we derived in section 4 lead to I = 0 but still they are not the general solution of this case, but rather particular solutions of it. Although this method does not lead to the full solution (where I =const.) of the system, it can be applied for an arbitrary functioñ g(φ). Thus, giving in closed form the full solution with the property I = 0 for any model of the integrable type we are considering.