The method of generating functions in exact scalar field cosmology

Exact solutions construction in scalar fields cosmology is of growing interest. In this work we review the results which obtained with the help of one of the most effective method. Namely, the method of generating functions for exact solutions construction in the scalar field cosmology. We also included into debate the superpotential method which may be considered as the bridge to slow roll approximation equations. Basing on the review, we suggested classification for the generating functions and found connection for all of them with the superpotential.

The construction of exact solutions in scalar field inflationary cosmology is of growing interest. In this work, we review the results which have been obtained with the help of one of the most effective methods, viz., the method of generating functions for the construction of exact solutions in scalar field cosmology. We also include in the debate the superpotential method, which may be considered as the bridge to the slow roll approximation equations. Based on the review, we suggest a classification for the generating functions, and find a connection for all of them with the superpotential. PACS numbers: 98.80.Bp, 98.80.Cq, 98.80.Es, 98.80.Jk

I. EARLY INFLATION AND THE IMPLEMENTED SCALAR FIELD
Inflationary expansion of the universe during very early times, once the universe emerged from the quantum gravity (Planck) era, has been proposed in the late 1970's and, mainly in the beginning of the 1980's and is becoming more accepted as a necessary stage of the standard Big Bang theory model. In the work of Starobinsky (1978) [1], nonsingular isotropic cosmological models with a massive scalar field were investigated. Although this work was concerned more with a bouncing cosmology rather than inflation, the field equations and the corresponding slow-roll solutions were already derived in this paper. These were later employed in the chaotic inflation model of Linde [9]. The works by Starobinsky (1980) [2], Guth (1981) [3] Linde (1982) [4], Albrecht and Steinhardt (1982) [5] include the physical mechanism based on quantum corrections and phase transitions during the very early stage of the universe. Exponential (de Sitter) expansion is the feature of inflationary models which helped to solve the long standing problems of the standard Big Bang theory model: the horizon, flatness, homogeneity, isotropy and some other problems.
In Sato's work (1981) [6], the first-order phase transition model of the early Universe that leads to an exponential expansion which stretches domains much greater than the horizon scales was considered. Also, in [7] it was * chervon.sergey@gmail.com † ingvor@inbox.ru ‡ beeshama@unizulu.ac.za pointed out that fluctuations associated with the phase transition are exponentially stretched and then may play the role of seed fluctuations for large-scale structures. Further, the monopole problem was also shown to be resolved by exponential expansion by Einhorn and Sato (1981) [8]. The chaotic inflation scenario proposed by Linde (1983) [9] differs from other previous versions since it is not based on the theory of high-temperature phase transition in the very early universe, but contains the locally homogeneous scalar field which is slowly rolling down to the minimum of the scalar field potential. After that proposal, many investigations took place of the inflationary universe connected with a self-interacting scalar field as the source of gravitation in the Friedmann world. Let us briefly mention some interesting works concerning the study of a scalar field in inflationary cosmology.
Homogeneous isotropic cosmological models with a massive scalar field have been studied in the works [10,11]. It was shown that inflationary stages are a fairly general property of most solutions in the considered model. The general conditions for inflation were investigated in the work [12]. It was found that under the lower limit for the amplitude of a scalar field, the universe naturally enters into and exits out of an inflationary phase. What is important is that such behavior takes place under a large variety of scalar potentials which are polynomial, logarithmic or exponential. It was also stated that a scalar field is essential for inflation [12]; it is unlikely that a vector or other non-scalar field will lead to inflation. The difference between scalar potentials in particle physics and those in cosmology has been stressed in the work [13]. The author wrote: "... we do not really know which theory of particle physics best describes the very early universe. One should therefore keep an open mind as to the form of V (φ)." Halliwell chose the exponential potential and showed that it leads to a solution with power-law inflation, and that this solution is an attractor. Detailed investigations of power-law inflation have been carried out in the work [14]. The authors found the constraints on the model coming from the requirement of solving the horizon, flatness, reheating and perturbation-spectrum problems. It was stated also that these constraints can be suitably satisfied. An exact power-law inflationary solution possessing an exponential potential was given in the work [15]. The generic inhomogeneous generalization of this solution having both scalar and tensor superhorizon "hairs" was derived in the paper [16].
Let us mention also the investigation carried out by Ivanov (1981) [17] where he found exact solutions for a nonlinear scalar field in cosmology. The solutions he obtained included polynomial, trigonometric and exponential potentials. The method he used for searching for exact solutions was subsequently called the Hamilton-Jacobi-like approach.
From the observational point of view, most results which can be related to observational data have been obtained from the so-called slow roll approximation of the cosmological dynamical equations [4], [5]. Detailed investigations of various physical phenomena from particle physics and GUT theories for the period until the 1990's can be studied from the reviews [18][19][20]. Our attention will be concentrated on exact solutions of inflationary models, the study of which started about ten years later, after inflationary cosmology had been proposed.
Thus we are going to present a brief review of the construction of exact solutions in the inflationary universe, i.e., the solutions of self-consistent Einstein and scalar field equations in Friedmann cosmology. The direct connection between scalar field cosmology and cosmology based on the perfect fluid stress-energy tensor needs to be mentioned. This connection is always valid except in the case of dust matter. Therefore we included the case of exact solutions for perfect fluid as the source of gravitation.
The construction of exact solutions in inflationary cosmology started with the work by Muslimov (1990) [21]. The results presented in that article will be discussed in Sec. III. Here we would like to mention that the very method and many interesting exact solutions presented in [17] have been reproduced and generalized in [21]. New methods and new sets of exact solutions have been developed in the work [21] as well.
Barrow [22] found a simple way to solve exactly the cosmological dynamic equations in terms of a pressuredensity relationship. In this way he obtained the known power-law and de Sitter forms of inflation and new classes of behavior in which the expansion scale factor increases as the exponent of some power of the cosmic time coordinate. The double-exponential law solution was obtained as well.
The work by Ellis and Madsen (1991) [23] was the first where "the inverse problem" was considered in the framework of cosmology. Usually one suggests that we know the scalar potential in the very early universe from HEP, and our task is to find the scale factor and the scalar field as functions of time. However Ellis and Madsen (1991) [23] suggested starting from the given scale factor! Indeed, it is clear that the scale factor may be found from observational data. Then we may take into account this fact to find the potential and scalar field from the cosmological equations. This work was done and examples of exact solutions have been presented for the pure scalar field (without taking into account radiation which is also considered there). Further this approach was developed in the works [24,25].
Reconstruction of models with a scalar field (quintessence) and dust-like particles (baryons and dark matter) from observational data was further developed by Starobinsky (1998) [26], Huterer and Turner (1999) [27], Nakamura and Chiba (1999) [28] using the luminosity data and in Starobinsky [26] from the growth factor of inhomogeneities.
Our paper is organised as follows: In section II, we present the basic equations of scalar field cosmology. In section III, we discuss generating functions for finding solutions, and sections IV is devoted to the classification of the generating functions. The superpotential method is presented in section V and we conclude with section VI.

II. BASIC EQUATIONS OF SCALAR FIELD COSMOLOGY
We consider the model of a self-gravitating scalar field φ with the potential of self-interaction V (φ). The action of such a model is where R is the curvature scalar, φ the scalar field, φ µ = ∂ µ φ the short representation of the partial derivative dφ/dx µ , κ is Einstein's gravitational constant, and Λ is the cosmological constant, which will mainly be included in the scalar field potential V (φ) as the constant part of it.
In the standard way one can obtain the energymomentum tensor (EMT) and the Einstein equation may be represented through the trace of the EMT in the form Varying the action (1) with the scalar field φ, we obtain the dynamic equation of the scalar field We consider the homogeneous and isotropic Universe as the spacetime with the Friedmann-Robertson-Walker (FRW) metric where ǫ = 0, ǫ = 1, ǫ = −1 for the spatially-flat, closed and open universe, respectively.
The Einstein Eq. (3) and the equation of the scalar field dynamics (5) in the FRW metric (6) lead to the system of equations Eqs. (7) and (8) can, in an equivalent way, be replaced by thier sum and the linear combination 3(7)+ (8). Including the Hubble parameter H =ȧ/a, the system (7)-(9) can be rewritten in the form We will refer to the system (10)-(12) as the Scalar Cosmology Equations (SCEs). The representation above, Eqs. (10)-(12) has some advantages for the derivation of any of the three equations (10)-(12) from the other two, and differential consequences of them.
Another representation of the SCEs was first proposed by G. Ivanov [17]. Suggesting the dependence of the Hubble parameter H on the scalar field φ, the transformation of the equations (10)- (12) for the spatially-flat universe (ǫ = 0) to the form, which was called later the Hamilton-Jacobi-like form, was made. Eq. (11) is transformed to Squaring the above equation and making the substitutioṅ φ 2 /2, expressed in term of H ′2 , and substituting into (10), one can obtain 2 3κ It is worthwhile to note that this procedure and the very equation (14) have been obtained by G. Ivanov in 1981 [17]. Unfortunately this result was published in limited editions (in Russian) and it was not familiar outside the USSR. Fortunately, in 1990, in the work of A. Muslimov [21], the Ivanov result was reproduced (Muslimov referenced the Ivanov article), and some of the solutions were generalized as well. In the same year, within the detailed investigation of long-wavelength metric fluctuations in inflationary models, D. Salopek and J. Bond [29] obtained the "separated Hamilton-Jacobi equation that also governs the semiclassical phase of the wave functional". The obtained equation contains a couple of scalar fields and definitely can be applied to a single scalar field. Therefore we suggest that Eq. (14) in the cosmological context should be called the Ivanov-Salopek-Bond (ISB) equation.

III. GENERATING FUNCTION FOR SOLVING THE IVANOV-SALOPEK-BOND EQUATION
The structure of Eq. (14) prompts us to find the form of the potential which can give the exact solution. Indeed, let the potential V (φ) be of the form where F = F (φ) is a C 1 function on φ and F * = const. We will call the function F (φ) the generating function. Comparing Eq. (14) with (15), it is easy to find the solution of (14) Thus, one can directly from (14) obtain the potential if the Hubble parameter is given. And vice versa, if one sets the potential in the form (15), then the solution of (14) will be defined by the Eq. (16).

A. Potential in polynomial form
As an example, let us choose the generating function F (φ) as the finite series on degrees of the field φ Under this circumstance, the potential V (φ) takes the following form (18) Let us consider the simple case when F * = 0, k = 0, p = 1. Then the potential becomes The generating function F (φ) and the Hubble parameter are If we additionally set λ 0 = 0, then we obtain the solution for the massive scalar field as in [17], with λ 2 1 = m 2 /2 (for the sake of simplicity we also set c = = 1). Thus the potential takes the form Solving Eq.(13) one can obtain the evolution of the scalar field The index "s" ("singularity") here is related to the values at the initial time t = 0, i.e., for a singularity in accordance with big bang theory. The Hubble parameter has a dependence φ on time (22), an this gives us a possibility to perform integration and obtain the dependence of the scale factor on time Thus we obtained the exact solution for the potential (21), which is represented by the dependence of φ on time (22), and the scale factor on time (24).
Also, for the first time, solution (24) was obtained in [9] by using the slow-roll approximation for the potential V (φ) = m 2 φ 2 /2 in contrast to the shifted potential (21).
It is interesting to note that the same solution and its application for the calculation of the number of e-folds and scalar spectral parameter were found and developed later by Wang [30].
When λ 0 = 0 the solution for the scale factor will differ by the factor a s in front of the exponent Hereã s = a s exp λ 0 κ/3 . The potential V (φ) then takes the form (26) Let us note that the linear transformation of the field without changing of the mass leads to the potential (21) for the fieldφ. Let us consider the case when k = 1, 2. The generating function F (φ) takes the form The potential V (φ) is We can make a simplification by considering λ 1 = 0. The potential then takes the Higgs form Using the relation (16), we may find H and H ′ expressed through φ Eq. (13) takes the form Performing the integration we find the dependence of φ on time Substituting this result into (31) and performing the integration, we will find the scale factor a(t) This is the double-exponental law solution [17], [22].
To make a comparison with Ivanov's results (with subscript "I" in his notation) [17], let us display the relations between the parameters of the model The case when λ 1 = 0 leads to the scalar field The scale factor then takes the following form (36) It is useful to note the role of the addition of the constant F * to the function F (φ). In our presentation for H(φ) (16), we can extract the constant part of the Hubble parameter The presence of the constant H * will be exhibited as the additional factor for a(t) In the considered examples above (25) and (36), such factors can be extracted explicitly.

B. Trigonometric potential
The solution for the potential which leads to the Sine-Gordon type equation was obtained in [17]. Such a setting used the special choice of the additional parameter. Let us consider this point in detail.
We choose the corresponding generating function F (φ) as Then the potential is To obtain the potential suggested in [17], it is enough to choose the parameter λ in the following way: λ 2 = 3κ 2 . Such a choice leads to the potential Equating the parameter A 2 = µ, we find the correspondence of the potential function with that presented in [17]. The Hubble parameter in terms of the scalar field can be defined from (16) Integrating equation (13), it is necessary to consider the integral dx cos x which has various functional representations In [17] the first representation was shown. We take the second representation in which the scalar field is defined from the relation Applying this result to (42) and performing the integration over time t, we obtain the scale factor C. Exponential potential The exponential potential in [17] is given as If we set then the potential takes the form Comparing this result to the original potential (46), we can find the relations In accordance with the general procedure explained at the beginning of Sect. III, one can obtain Then the dependence of the scalar field on time t has a logarithmic character The scale factor is evaluated via a power law An addition of the constant F * to F (φ) leads to the generalization of the solution (52) This is the exponential power law solution. Then the potential acquires the additional terms Muslimov (1990) [21] found the generalization of Ivanov's solution for the exponential potential. Let us represent this, which contains both solutions.
If we take the generating function F (φ) in the form (47) with the potential (54) then We can find that Integrating (13), we obtain where v * is a constant of integration. Finally we find To obtain Ivanov's solution [17], we set v * = 0 and take into account the relations (49).
To obtain Muslimov's solution [21], we set v * = 1 and take into account the relations below The solution (58) without restrictions on the parameter v * gives some generalization.

D. The solution with an inverse potential
The potential in [21] was presented in the following way It is not difficult to check that the same potential can be obtained from the generating function Therefore it is clear that the solution can be obtained by a general scheme. The Hubble parameter is As we know the influence of H * on the result, we may take H * = 0 for the sake of simplicity.
Integrating (13) we can find the dependence of the scalar field on time This result leads to the time dependence of the Hubble parameter Then the scale factor is (64) The solution of such a type can be confronted with both a very early and late time universe. Similar solutions have been obtained in [22].

E. The solution with an intermediate (hyperbolic) function
In the range of the results described above, Muslimov [21] suggested a new original approach for solving the scalar field cosmology equation. To simplify calculations, let us, following [21], introduce a new variable and the potential function For this notation, the Ivanov-Salopek-Bond (ISB) equation (14) reduces to The upper case corresponds to the positive sign of the potential. It is interesting to mention that an equation of this type was studied by Mitrinovitch in 1937 [31].

Let us search for a solution in the form
The other choice is for the lower sign Here we included inverse hyperbolical tangents instead of inverse hyperbolical cotangents as in Muslimov's work with the aim of avoiding plus-minus signs in the final equation. Using the formulae above, for transition to the function y(x), we obtain Here, for the sake of briefness, we introduce the function u in the following way We can shift the second term on the left hand side of (70) to the right hand side and, using the property of the hyperbolic function, take the root on the left and right hand side of the equation (considering all values as positive). As the result, we obtain the equation The transition to the function y(x) is performed by inverse substitution of (71) and by introducing the derivative Finally we arrive at the following relation After simple algebraic transformations, we acquire the Abel equation Repeating the same procedure for the lower case, we once again arrive at the equation (75).

F. Kim's exact solutions
H.-C.Kim in the article [32] proposed the generating function G(φ), which we can compare with the Ivanov generating function F (φ) (15)- (16) with F * = 0. The Hubble parameter in the work [32] is represented as To obtain the relation between the generating functions F (φ) and G(φ) we may simply equalise F (φ) ≡ G(φ). With the representation (76) the solution of the SCEs can be obtained with the following formulae In the work [32], two generating functions were considered. First of them is for the constant potential V (φ) = Λ > 0 : This gives De Sitter space-time with the scale factor expanding exponentially.
The second generating function is given by The scalar field and the scale factor behave as The third generating function leads to the power-like potential where µ, n are constants. The scalar field and the scale factor are where we set φ 0 = φ(0) = 0 at the time t = t 0 .

G. Exact solutions for constant-roll inflation
Now, we consider inflationary models with the additional conditionφ where α is arbitrary constant parameter. Standard slowroll inflation occurs when α ≃ −3 while for the other ultra-slow-roll regime, one has α = 0. Such models correspond to the constant-roll inflation models which interpolate between these two regimes, and were considered in the papers [33,34]. Also, this approach was used for analysis of cosmological models in f (R) gravity on the basis of conformal transformations from the Jordan frame to the Einstein frame which lead to Ivanov-Salopek-Bond (Hamilton-Jacobi-like) equations (13)-(14) [35]. In the works [35,36], it was shown that these types of models satisfy the latest observational constraints. Eq. (13), in the system of units with κ = 1, and Eq. (87) gives the condition on the Hubble parameter [34] with general solution From Eq. (14) one has the potential On the basis of the general solutions (89)-(90), one can find the specific ones [34]: • The solutions for α > −3 only: Power-law inflation 1 , which was already considered in subsection (III C) where M is an integration constant.
• The solutions for both cases α > −3 and α > −3: The inflation model with hyperbolic scale factor in the 1 In the original work [34], the time t was omitted in the expression for the Hubble parameter (94), but the scale factor was correct.
case of α > −3: This solution is equivalent to a solution found in [48] but in a different context. The inflation model with hyperbolic scale factor in the case of α < −3: This is the a particular case of hilltop inflation which is considered in paper [37]. The model with oscillating scale factor: From the scale factor (110) one hasä(t) < 0; nevertheless, this solution corresponds to the oscillatory dynamics of the early universe.

IV. THE CLASSIFICATION OF GENERATING FUNCTIONS
We propose a classification of generating functions as they appeared in the literature in chronological order. As the first class, we name the Ivanov generating function F (φ) (15)- (16). It was not represented in direct form in [17], but in the present publication for the first time, using this generating function, we recover and generalise all solutions of Ivanov's work [17]. Now we continue the classification of the generating functions that occur in the literature.
A. The second class of generating functions L. P. Chimento, A. E. Cossarini, and A. S. Jakubi (1993)(1994)(1995) in the papers [38,39] represented the potential of the scalar field in the form where F (a) is a new type of generating function. The scalar field equation can be integrated and it yields where C is an arbitrary integration constant. Thus, the problem of generation of exact solutions has reduced to the quadratures: where ∆t ≡ t − t 0 , ∆φ ≡ φ − φ 0 and t 0 , φ 0 are arbitrary integration constants. For the generating function where B(> 0), b(> 0), s and n are constants and s(n + 1) = 6. If we take C = 0, the potential is simplified to This potential has a non-vanishing minimum at ∆φ = 0 for s > 0, which is equivalent to an effective cosmological constant. When s < 0 , the origin becomes a maximum, and the potential vanishes exponentially for large φ.
In [38] the equation (113) was evaluated for some values of s: For s > 0, the evolution begins from a singularity as ∆t 1/3 and is asymptotically De Sitter with ∆φ → 0 for t → ∞. On the other hand, for s < 0 the evolution has a deflationary behaviour from a De Sitter era in the far past to a Friedmann behavior ∆t 1/3 when t → ∞.

B. The third class of generating functions
In the paper F. E. Schunck and E. W. Mielke (1994), [40] the equations (7)-(9) are written in the following formḢ where The potential as a function of the Hubble parameter V = V (H) is defined as By choice of the graceful exit function g(H), the exact solutions of the SCES can be generated. The choice for power-law and intermediate inflation are where n is real and A is a positive constant. For n = 0, the following solutions were found: For n = 1 For n = 2 For n = 0, 1, 2 (138)

C. The fourth class of generating functions
A. T. Kruger and J. W. Norbury (2000) [41] proposed the generating function F = F (φ) viȧ Then the scalar field dynamic equation transforms tö Equations (139) and (140) form a set of coupled equations which are equivalent to equation (12). They can be uncoupled by differentiating (139) with respect to time and substitutingφ in (140). Using the evident relation φ = d(φ 2 )/(2dt), by taking time derivative of the square of Eq. (139), we havë where F ′ ≡ dF/dφ. Then, inserting the expression forφ into (140) we obtain It can be seen that if one chooses F ≡ F (φ), then equation (142) is always separable and the potential is given by where F ≡ F (φ) and β is a constant. The choice of the generating function in the work [41], is F = cosh(λφ). The the solution can be obtained from (143) and (139): These equations are only consistent if all upper or all lower signs are taken, i.e., one should not mix upper and lower signs. For the sake of simplicity, we chose λ = √ 6. The upper " − " sign in equation (144), corresponding to the upper " + " sign in the definition (139), gives The lower " + " sign in (144), corresponding to the lower " − " sign in the definition (139), leads to Using this result (146) in equation (139), we obtaiṅ The solution for φ(t) yields where D is a constant. The form of the potential suggests that φ(t) be a function that decreases from an initial maximum value similar to the chaotic inflation model. We can choose D = 0 and so (148) can be written as This solution can be used to determine the evolution of the scale factor with time for the expanding (" + " square root) solution D. The fifth class of generating functions T. Charters and J. P. Mimoso (2010) in the paper [42] proposed a generating function in the following form The scalar field's dynamical equations are transformed to the following Further, in [42] the known and new exact solutions were obtained by the special choice of the generating function x(φ). For example, the following type of x(φ) = λφ corresponds to the generalization of the Easther solution [43]: where Ei is the exponential integral function.
where the function G can be defined via the scalar field from the equation (159) as In [44] the authors considered the case when the scalar field potential can be represented as the function of G in the form where α 1 is an arbitrary constant. With this choice, the evolution equation takes the simple form Performing the integration one can obtain the general solution where φ 0 is an arbitrary constant of integration. Using the form (163) of G, one can obtain the self-interaction potential of the scalar field and the scale factor A simple solution of the gravitational field equations for a power-law type scalar field potential can be obtained by assuming for the function G the following form With this choice of G, Eq. (158) immediately provides the scalar field potential given by where V 0 is an arbitrary constant of integration. The time dependence of the scalar field is given by a simple power law, The scale factor can be obtained from da/dφ = (1/ √ 6) coth G a = (φ/2α 2 ) a, and it has exponential dependence on the scalar field and time were a 0 is an arbitrary constant of integration.

V. THE SUPERPOTENTIAL METHOD
The superpotential method for the standard inflationary model (1) was successfully applied for solving the SCEs (7)- (9). The main idea was [45] to represent the SCEs in the form of the slow roll approximation. i.e., in the equations for a spatially-flat universe (we obtain them from (7)-(9) by setting ǫ = 0) we should omitφ,φ 2 . After that, the SCEs in the slow roll regime take the following form To obtain the desired form of the equations, the potential of total energy [46] (or the superpotential [25], [47]) was introduced, as the sum of kinetic energy (in terms of scalar field argument U (φ) =φ) and the potential energy After this substitution, the SCEs take the following form: As we can see the system of equations above exactly reproduce the SCEs in slow roll form (173)-(175) if we substitute W (φ) instead of V (φ) and U (φ) instead ofφ.
As we know, any one from the presented exact SCEs in the slow roll form (177)-(179) can be derived as a consequence (or a differential consequence) of the remaining two. In our approach, we exclude from consideration equation (178) in the first step.
Excluding from (177) the Hubble parameter H, and inserting it into (179), taking into account the superpotential definition (176), we have the consequence of (177) and (179) in the form Integrating (180) with respect to W , we obtain the relation: which leads to the new method of exact solution construction in cosmology, viz., by suggesting that evolution of the scalar field is given, one can determine the superpotential by solving the integral on the right hand side of (181) φ2 dt. Knowing W , one can find H from the Friedmann equation (177) with the following relation Then, by integration, one can find the scale factor a(t). Thus, the proposed method presents some combination of the two methods: slow roll-like presentation of the exact equations [46] and obtaining cosmological solutions for the given scalar field evolution [48].
The advantage of the proposed method lies in the essential simplification of the integration procedure: one needs to calculate only one integral for obtaining a superpotential and Hubble parameter. Then the potential V (176), as well as the scale factor a(t), are calculated from related definitions. The exhibition of the simplicity of the procedure, and its effectiveness can be found in the applications of the method to cosmology on the brane, in phantom and tachyon fields [49], [50], [51]. The two last have a very restricted number of exact solutions which can be essentially extended by virtue of the superpotential method.
We can have a look at the superpotential method from another position. The system (170) -(172) has three unknowns, φ(t), V (φ) and H(t) (or a(t)). To solve it, one of these variables has to be given a priori. It is customary to look for the solution for a given V (φ), but as it is known, it is very difficult to solve the SCEs for a given potential exactly.
In the superpotential approach, it is proposed to reduce the equations to a simpler form which helps to solve them exactly. In order to do this, let us consider the superpotential function W (φ) defined in (176) Now, with the change of variable dt = dφ/φ and using the inverse transformation from U (φ) toφ, one can obtain Hence, the SCEs (170) -(172) can be rewritten in the slow roll form (173) -(175). Therefore we can try to solve SCEs in the superpotential presentation because the slow roll approximation was intensively studied. The superpotential W (φ) (> 0) shows up as the main part of the potential function, driving the dynamics of the Hubble parameter H or the scale factor. To solve them, note that Eq. (177) definesȧ/a as a function of φ, H(φ), which when inserted into Eq.(179), gives the scalar field φ(t) as a function of t, at least in quadratures Finally, inserting φ(t) into eqs. (176) and (179) gives V (φ) and a(t), respectively, and the solution is completed.
Obviously, one could simply have begun by giving H=H(φ), but it is usually desirable to have some description of the potential instead, and for this reason it is preferable to give W (φ). One could also use H(t) to determine φ(t), since implies that where ∆φ(t) = φ−φ * , φ * is the integration constant. So since W = 3H 2 (t), a complete knowledge of H(t) fully determines the solution to the problem. Also, we can obtain potential V from the superpotential W . (187)

A. Examples of exact solutions
Firstly, we consider the superpotential W = λφ 2n . For this type of superpotential, we have V (φ) = λφ 2n − λ 2n 2 3 φ 2(n−1) , (190) Another potential used in the literature is that of a hyperbolic cosine W = V 0 (cosh(βφ) − 1) [52]. The solution found there, with φ 0 =0, is V (φ) = V 0 (cosh(βφ) − 1)−V 0 β 2 6 (cosh(βφ) + 1) , (192) The expressions for generating functions by means of superpotential W are represented in the Table 1. The space of exact solutions in scalar field cosmology is very huge. In this work, we review exact solutions in inflationary cosmology and the method of finding such solutions. There are basically three methods of constructing such solutions. For any scale factor a(t), we can find the potential and kinetic energy to be satisfied in the selfconsistent system of Einstein and scalar field equations. Secondly, from the evolution of the scalar field φ = φ(t), we can find the scale factor a(t) and the potential V (t), thus defining the exact solution. Thirdly, choosing the generation function, represented in this given review, one can, once again, obtain a great number of exact solutions.
Thus, as the next step of the investigation, we suggest analysing the method of confrontation of theoretical predictions from exact solutions with observational data.