A Real Scalar Field Unifying the Early Inflation and the Late Accelerating Expansion of the Universe through a Quadratic Equation of State: The Vacuumon

In a previous paper we introduced a cosmological model describing the early inflation, the intermediate decelerated expansion, and the late accelerating expansion of the universe in terms of a single barotropic fluid characterized by a quadratic equation of state. We obtained a scalar field representation of this fluid and determined the potential V(ϕ) connecting the inflaton potential in the early universe to the quintessence potential in the late universe. This scalar field has later been called the `vacuumon’ by other authors, in the context of the Running Vacuum model. In this paper, we study how the scalar field potential is modified by the presence of other cosmic components such as stiff matter, black-body radiation, baryonic matter, and dark matter. We also determine the mass m and the self-interaction constant λ of the scalar field given by the second and fourth derivatives of the potential at its extrema. We find that its mass is imaginary in the early universe with a modulus of the order of the Planck mass MP=(ℏc/G)1/2=1.22×1019GeV/c2 and real in the late universe with a value of the order of the cosmon mass mΛ=(Λℏ2/c4)1/2=2.08×10−33eV/c2 predicted by string theory. Although our model is able to describe the evolution of the homogeneous background for all times, it cannot account for the spectrum of fluctuations in the early universe. Indeed, by applying the Hamilton–Jacobi formalism to our model of early inflation, we find that the Hubble hierarchy parameters and the spectral indices lead to severe discrepancies with the observations. This suggests that the vacuumon potential is just an effective classical potential that cannot be directly used to compute the fluctuations in the early universe. A fully quantum field theory may be required to achieve that goal. Finally, we discuss the connection between our model based on a quadratic equation of state and the Running Vacuum model which assumes a variation of the cosmological constant with the Hubble parameter.


Introduction
The universe displays three main periods of evolution: an early phase of inflation during which the scale factor increases exponentially rapidly with time, an intermediate phase of decelerated expansion during which the scale factor increases algebraically, and a late phase of accelerating expansion during which the scale factor increases exponentially rapidly again. 1 The idea that a period of accelerated expansion (early inflation) may have occured in the early universe was introduced by Guth [1] in 1981 to explain the observed isotropy, homogeneity, and flatness of the universe in a natural way. The early inflation also explains the near scale-invariant spectrum of cosmological perturbations. The present acceleration of the universe was discovered at the end of the twentieth century [2][3][4][5] and was a surprise.
The intermediate phase of decelerated expansion is relatively well-understood. It corresponds to a relativistic radiation era followed by a nonrelativistic matter era. These two periods are described by a linear equation of state where P is the pressure and ρc 2 is the energy density. The coefficient is equal to α = 1/3 for radiation (corresponding to a gas of photons or other ultrarelativistic particles such as neutrinos) and to α = 0 for pressureless matter (baryonic matter and dark matter). One can also consider a small value of α = k B T/mc 2 ∼ 10 −7 in the matter era in order to take into account thermal effects (to make this estimate we have used v c ∼ (k B T/m) 1/2 ∼ 100 km/s obtained from the rotation curves of the galaxies). The early inflation and the late accelerating expansion of the universe are less well understood. The fundamental constant that describes quantum mechanics is the Planck constanth = 1.05 × 10 −31 m 2 g s −1 [6]. From the Planck constant, and from the other fundamental constants of physics (the speed of light c and the gravitational constant G), one can construct a density ρ P = c 5 hG 2 = 5.16 × 10 99 g m −3 , called the Planck density. This density is extremely high. It is expected to play a fundamental role in the early universe (which is very dense) and be responsible for the first phase of inflation. It is usually believed that inflation corresponds to a de Sitter stage during which the density of the universe is constant and equal to the Planck density. 2 On the other hand, Einstein [7] introduced a cosmological constant Λ in the equations of general relativity in order to obtain a static universe. Although Einstein rejected this constant after the discovery of the expansion of the universe, calling it his "biggest blunder" [8], we now know that a nonzero value of the cosmological constant Λ = 1.00 × 10 −35 s −2 is favored by current observations. The effect of the cosmological constant is equivalent to the effect of a constant density called the cosmological density. This density is extremely low. It is expected to play a fundamental role in the late universe (which is very dilute) and be responsible for its observed present-day acceleration. The origin of this acceleration is generally called dark energy. It is usually believed that the late accelerating expansion of the universe corresponds to a second de Sitter stage during which the density of the universe is constant and equal to the cosmological density. The ratio between the Planck density (early universe) and the cosmological density (late universe) is They differ by 123 orders of magnitude. The cosmological constant is usually interpreted in terms of the vacuum energy. The vacuum is described by an equation of state P = −ρc 2 [9][10][11] which implies a constant energy density leading to a de Sitter exponential expansion. However, particle physics predicts that the vacuum energy should be of the order of the Planck scale. Therefore, this interpretation leads to a discrepancy of 123 orders of magnitude with the measured value of Λ. This is the so-called "cosmological constant problem" [12,13]. Various models have been proposed to describe the primordial inflation. The model of Starobinsky [14] takes into account quantum gravitational effects. It has a purely geometric nature and consists in a generalization of the Einstein-Hilbert action to contain an R 2 contribution, where R is the Ricci scalar curvature. It leads to a nonsingular de Sitter stage instead of the initial big bang singularity predicted by the classical Einstein equations. The phase of inflation in the very early universe can also be described by a hypothetical scalar field φ, called inflaton, running down a potential. This scalar field has its origin in the quantum fluctuations of the vacuum and is usually associated with a nonequilibrium phase transition [15]. Interestingly, the Starobinsky model can be mapped on a scalar field model with a special form of potential (see, e.g., [16] and references therein).
Various models of dark energy have also been proposed to describe the late acceleration of the universe. The simplest model is the cold dark matter model with a cosmological constant (ΛCDM model), which describes dark matter as a pressureless fluid and dark energy as a cosmological constant [17,18]. The ΛCDM model provides a very good description of the large scale structure of the universe and can account for the observations of the Planck mission [19,20]. However, the ΛCDM model suffers from the cosmological constant problem mentioned above and from the "cosmic coincidence problem" [21][22][23], namely why the fractions of dark matter and dark energy turn out to be of the same order of magnitude at the present epoch although they scale differently with the universe's expansion. 3 The CDM model also suffers from small-scale problems (at the galactic scale) such as the core-cusp problem of dark matter halos [24], the missing satellite problem [25][26][27], and the "too big to fail" problem [28]. This leads to the so-called small-scale crisis of CDM [29]. In order to solve or alleviate these difficulties, other models of dark energy have been introduced. Inspired by the models of inflation, some authors have proposed to describe the dark energy in terms of a self-interacting scalar field called quintessence [30][31][32][33] which can be interpreted as a dynamical vacuum energy. Other authors [34] have invoked an exotic fluid with a negative pressure called the Chaplygin gas [35]. This model provides a unification of dark matter and dark energy in terms of a single dark fluid. Additional models of unified dark matter and dark energy (called UDM or quartessence models [36]), such as the generalized Chaplygin gas [34,37], the modified Chaplygin gas [38], the polytropic gas [39][40][41][42][43][44][45][46][47], or the logotropic dark fluid [48][49][50][51][52], have been introduced subsequently. These models are consistent with observation data only if they are extremely close to the ΛCDM model [53] which is equivalent (in its UDM interpretation) to a single dark fluid with a constant negative pressure P = −ρ Λ c 2 [42,53,54].
In a series of papers [39][40][41][42][43][44][45][46][47], we have proposed to describe the whole evolution of the universe, from its early inflation to its late accelerating expansion, by a quadratic equation of state of the form in which the coefficients are the Planck density ρ P (see footnote 2 ) and the cosmological density ρ Λ . The left term describes the early inflation, the middle term describes the decelerated expansion (corresponding to a fluid with a linear equation of state P = αρc 2 ), and the right term describes the late accelerating expansion. This equation of state can be simplified in some limits. In the early universe (ρ ρ Λ ), the last term in Equation (5) is negligible and the quadratic equation of state P = −(α + 1) ρ 2 ρ P c 2 + αρc 2 (6) describes the transition between the early inflation and the phase of decelerated expansion.
On the other hand, in the late universe (ρ ρ P ), the first term in Equation (5) is negligible and the affine equation of state P = αρc 2 − (α + 1)ρ Λ c 2 (7) describes the transition between the phase of decelerated expansion and the late accelerating expansion (the ΛCDM model, corresponding to a constant equation of state P = −ρ Λ c 2 , is recovered for α = 0). Interestingly, the equations of state (6) and (7) can be viewed as generalized polytropic equations of state of the form P/c 2 = αρ + kρ 1+1/n involving a linear term P/c 2 = αρ and a polytropic term P/c 2 = kρ 1+1/n with a negative pressure (k < 0) and an index n = +1 (resp. n = −1) in the early (resp. late) universe. These generalized polytropic equations of state have been studied at a general level in [41][42][43]. They can be viewed as generalized (or modified) Chaplygin gas models. Our main results can be summarized as follows: (i) In Ref. [41], we studied in detail the equation of state (6) with α = 1/3 describing the smooth transition between the early inflation and the radiation era. This equation of state provides a "graceful exit" to the de Sitter era. We considered more general models with an arbitrary value of α (instead of α = 1/3) and an arbitrary positive polytropic index n e (instead of n e = +1). We showed that the results remain qualitatively the same in these more general situations. (ii) In Ref. [42], we studied in detail the equation of state (7) with α = 0 describing the smooth transition between the matter era and the late inflation. We showed that this equation of state returns the ΛCDM model. We considered more general models with an arbitrary value of α (instead of α = 0) and an arbitrary negative polytropic index n l (instead of n l = −1). We showed that the results remain qualitatively the same in these more general situations. (iii) In Refs. [39,40,42,47], we described the complete history of the universe. It involves two de Sitter eras (early and late inflation) bridged by an intermediate decelerated era. These results are reported in Figure 14 of [42]. They have been obtained by connecting the results valid in the early universe (see Equation (6)) to the results valid in the late universe (see Equation (7)). This approach describes the successive phases of early inflation, radiation, matter, and late inflation. In this manner, our model smoothly connects the primordial inflation to the ΛCDM model. In our model, the early and late evolution of the universe is remarkably symmetric. It is described by two polytropic equations of state of index n e = +1 and n l = −1 respectively. In addition, the cosmological density ρ Λ in the late universe is the counterpart of the Planck density ρ P in the early universe. As the universe expands, the density decreases from the Planck density ρ P = 5.16 × 10 99 g m −3 to the cosmological density ρ Λ = 5.96 × 10 −24 g m −3 , spanning 123 orders of magnitude (see Figure 15 of [42]). The resulting model of universe is non-singular and non-phantom. There is no big bang singularity in the past, nor big rip singularity in the future. The early and late behaviors of the universe are described by two de Sitter eras with density ρ P and ρ Λ , respectively. The universe exists eternally in the past and in the future. There is no question such as "What happens for t < 0 before the big bang?". We called this nonsingular and fully symmetric model of universe, exhibiting two extreme de Sitter eras bridged by a period of decelerated expansion, the "aioniotic" universe (see Section 7.4 of [42]). (iv) In Refs. [39][40][41][42]47], we studied the thermal history of the universe. As the Friedmann equations are dissipationless, the total entropy of the universe (including all kinds of matter and energy) is constant. It has the very large value S/k B = 5.04 × 10 87 [41]. We obtained a generalized Stefan-Boltzmann law valid in the early universe (see Equation (84a) of [41]). In our model, the temperature T increases exponentially rapidly during the inflation up to the Planck temperature T P = 1.42 × 10 32 K, then decreases algebraically during the radiation and matter eras (see Figure 16 of [42]). This is very different from other models of inflation where the temperature drops drastically during the exponential inflation and one has to invoke a phase of re-heating by various high energy processes (that are not very well-understood) in order to restore the initial temperature. (v) In Refs. [39,40,42,47], we developed a scalar field representation of our model. We determined the "inflaton" potential (see Equation (123) of [42]) associated with the equation of state (6) which describes the smooth transition between the early inflation and the radiation era, and we determined the "quintessence" potential (see Equation (125) of [42]) associated with the equation of state (7) which describes the smooth transition between the matter era and the late inflation.
(vi) In Refs. [45][46][47], we considered the possibility that the cosmic history of the universe involves an additional stiff matter era after the inflation and prior to the radiation era. This stiff matter era is described by an equation of state of the form P = ρc 2 where the speed of sound, given by c 2 s = P (ρ), is equal to the speed of light (c s = c) [45,55,56]. We proposed to describe the transition between the inflation and the stiff matter era in the primordial universe by an equation of state of the form of Equation (6) with α = 1, and we derived the corresponding scalar field potential (see Equation (140) in [45] and Equation (F.42) in [46]). 4 (vii) In Refs. [39,40,42,44], we solved the general model described by the quadratic equation of state (5). We first derived the explicit relation between the energy density and the scale factor (see Equation (86) of [44]). We then obtained an exact analytical solution of the Friedmann equations giving the complete temporal evolution of the scale factor a(t) and energy density ρ(t) from t = −∞ to t = +∞ (see Equation (106) of [44]). This solution describes the early inflation, the intermediate decelerated expansion, and the late accelerating expansion of the universe. The quadratic equation of state (5) therefore provides a unification of the early and late inflation of the universe. We determined the general scalar field potential associated with this equation of state.
We obtained its exact analytical expression in terms of Jacobian Elliptic functions (see Equation (121) of [44]) and proposed a simple approximate expression obtained by using matched asymptotic expansions (see Equation (131) of [44]). Interestingly, our scalar field theory describes, with a unique potential, the whole evolution of the universe, from its early inflation to its late accelerating expansion, passing through a phase of algebraically decelerating expansion. In this sense, the scalar field potential V(φ) unifies the inflaton potential in the early universe and the quintessence potential in the late universe.
Very similar results have been obtained in parallel by J. Solà and his collaborators in the context of the Running Vacuum Model (RVM). The basic idea behind this model (see the review [57] for an exhaustive list of references) is that the cosmological constant actually depends on time. Using results of particle physics and the renormalization group approach, they argued that the cosmological constant Λ(H) is related to the Hubble constant H by a quartic equation. When combined with radiation and matter, this model produces a phase of early inflation, a phase of decelerated expansion, and a phase of late inflation which are very similar to the ones obtained in our model (see Appendix C for a comparison between the two approaches). Recently, Basilakos et al. [58] developed a scalar field representation of the RVM. They proposed to call the scalar field that accounts for the temporal evolution of the vacuum energy the "vacuumon". As they obtained exactly the same potentials in the early and late universe as the ones obtained previously in our papers (compare Equations (4.18) and (4.42) of [58] with Equations (123) and (125) of [42]), we will also call "vacuumon" the scalar field associated with the quadratic equation of state (5) of our model [39][40][41][42][43][44][45][46][47]. 5 Following this terminology, what we referred to as "inflaton" in the early universe (associated with the equation of state (6)) will be called "early vacuumon" and what we referred to as "quintessence" in the late universe (associated with the equation of state (7)) will be called "late vacuumon". Even more recently, these authors developed a new model that they called String-inspired Running Vacuum [59]. They managed to make a connection between the phenomenological RVM and string theory. Interestingly, their new model includes a stiff matter era in the primordial universe similar to the one introduced heuristically in our papers [45,46]. They described the transition between the inflation and the stiff matter era in a manner similar to the one described in Section XI of [45] or in Appendix F of [46]. They also derived a scalar field potential coinciding with the one previously obtained in our papers [45,46] (compare Equation (40) of [59] with Equation (140) of [45] or Equation (F.42) of [46]). 6 Therefore, the RVM [57][58][59] and our model [39][40][41][42][43][44][45][46][47] are consistent, valuable, and complementary to each other (see Appendix C).
We have argued above that the quadratic equation of state (5) describes the whole evolution of the universe from its early inflation to its late accelerating expansion. There is, however, a difficulty with this description which is connected to the value that one should ascribe to the parameter α. If we take α = 1/3, the first two terms in Equation (5) describe the transition between the inflation and the radiation and the last term describes the late accelerating expansion of the universe. However, this equation of state does not account for the matter era. Indeed, it describes a universe undergoing early inflation, radiation era, and late inflation. Alternatively, if we take α = 0, the first term in Equation (5) describes the early inflation and the last two terms describe the transition between the matter era and the late accelerating expansion of the universe. However, this equation of state does not account for the radiation era. Indeed, it describes a universe undergoing early inflation, matter era, and late inflation. This implies that the quadratic equation of state (5) is not able to describe the whole content of the universe. Indeed, it cannot describe simultaneously the radiation era and the matter era. At that point, we have two possibilities: (A) The first possibility is to assume that the coefficient α that appears in the equation of state (5) depends on the density in such a way that α → 1/3 at high densities and α → 0 at low densities. In the early universe (high densities), we can neglect matter and dark energy and use Equation (6) with the coefficient α = 1/3 (radiation).
In the late universe (low densities), we can neglect inflation and radiation and use Equation (7) with the coefficient α = 0 (matter). We then have to match these two asymptotic limits. This is the point of view adopted in [39,40,42,47]. This point of view is also consistent with the RVM (see Appendix C) provided that the density ρ is interpreted as the total density of the universe, i.e., the sum of the running vacuum energy density plus the energy density of radiation in the early universe or the matter energy density in the late universe (the same comment applies to the pressure P). (B) Another possibility is to assume that the quadratic equation of state (5) with a fixed coefficient α describes only one cosmic fluid. This exotic fluid could correspond to a scalar field in its hydrodynamic representation. Then, we must consider, in addition, the contributions of other species treated as independent noninteracting fluids. These additional species correspond to standard fluids (stiff matter, radiation, and baryonic or dark matter) described by a linear equation of state. Different choices are possible. A first choice is to take α = 1/3. In that case, the quadratic equation of state (5) characterizes a scalar field which is responsible for a phase of early inflation, a phase of radiation (it could be the standard radiation corresponding to photons or relativistic particles or a "dark radiation" different from the standard radiation), and a phase of late inflation. This scalar field provides a unification of inflation, radiation, and dark energy. Then, we have to add standard radiation (α r = 1/3), baryonic matter (α b = 0), and dark matter (α dm = 0) as additional species. Another choice is to take α = 0. In that case, the quadratic equation of state (5) characterizes a scalar field which is responsible for a phase of early inflation, a phase of pressureless dark matter, and a phase of late inflation. This scalar field provides a unification of inflation, dark matter, and dark energy. Then, we have to add standard radiation (α r = 1/3) and baryonic matter (α b = 0) as additional species. We can also choose another value of α, different from 1/3 or 0, such as α = 1 corresponding to stiff matter [45,46]. In that case, the quadratic equation of state (5) characterizes a scalar field which is responsible for a phase of early inflation, a stiff matter era, and a phase of late inflation. This scalar field provides a unification of inflation, stiff matter and dark energy. Then, we have to add standard radiation (α r = 1/3), baryonic matter (α b = 0) and dark matter (α dm = 0) as additional species.
In conclusion, we are led to a considering a model of universe involving a scalar field described by a quadratic equation of state given by Equation (5) plus zero, one, or several X-fluids described by a linear equation of state P = α X ρc 2 . We assume that these different species are independent from each other and noninteracting. The case of a scalar field alone described by the quadratic equation of state (5) has been treated in our previous papers [39,40,42,44]. We found that this scalar field has a potential given by Equation (108) in the general case, and by Equations (115) and (132) in the early and late universe, respectively. This is the potential of the "bare" vacuumon. If we now consider a scalar field described by the quadratic equation of state (5) in the presence of other species, such as X-fluids described by a linear equation of state P = α X ρc 2 , the potential of the scalar field will change. 7 In this paper, we explain how this potential can be calculated and we give its general expression under the form of an integral. This is the potential of the "dressed" vacuumon due to the presence of other species. Unfortunately, its general expression cannot be obtained analytically, except in particular limits.
Once the potential V(φ) has been obtained, we can determine the main characteristics of the scalar field such as its mass m and self-interaction constant λ which are given by the second and fourth derivatives of the potential at its extrema [47]. The potential of the vacuumon presents a maximum V = ρ P c 2 in the early universe at φ = 0 and a minimum V = ρ Λ c 2 in the late universe at φ = φ max . Correspondingly, the early vacuumon (inflaton) has an imaginary mass of the order of the Planck mass M P = (hc/G) 1/2 = 1.22 × 10 19 GeV/c 2 and the late vacuumon (quintessence) has a real mass of the order of the cosmon mass m Λ =h √ Λ/c 2 = 2.08 × 10 −33 eV/c 2 predicted by string theory. We find that the mass of the vacuumon in the early universe satisfies a fundamental quantization rule while its mass in the late universe does not. Finally, we apply the Hamilton-Jacobi formalism [60] to our model of early inflation in order to obtain the Hubble hierarchy parameters and the spectral indices. We show that it leads to severe discrepancies with the observations. This suggests that the scalar field potential of the vacuumon is just an effective classical potential that cannot be directly used to compute the spectrum of fluctuations in the early universe. A fully quantum field theory may be required to achieve that goal.
The paper is organized as follows. In Section 2, we recall the basic equations that describe the cosmic evolution of a fluid with a linear equation of state and the basic equations that describe the cosmic evolution of a canonical self-interacting real scalar field. In Section 3, we briefly review our model of universe [39][40][41][42][43][44][45][46][47] based on the quadratic equation of state (5). In Section 4, we determine the potential of a scalar field associated with an arbitrary barotropic equation of state P(ρ) in the presence of X-fluids and apply these general results to the quadratic equation of state (5). In Section 5, we show how these results can be simplified in the case of a scalar field alone in the universe and we recover the results of our previous papers [39][40][41][42][43][44][45][46][47]. In Section 6, we determine the parameters of the scalar field (vacuumon) such as its mass and its self-interaction constant in the early and late universe. In Section 7, we determine the potential of the scalar field in the presence of an additional fluid in the intermediate regime between the early and late inflation. In that case, the scalar field potential can be obtained analytically. In Section 8, we consider the evolution of the scalar field in the presence of one fluid in the late universe. In Section 9, we determine the Hubble hierarchy parameters corresponding to our model of early inflation (early vacuumon) and conclude that a fully quantum field theory is required to achieve a form of agreement with observations.

Basic Equations
In this section, we recall the basic equations determining the cosmological evolution of a fluid described by a linear equation of state (called here X-fluid) and the basic equations determining the cosmological evolution of real canonical self-interacting scalar field φ described by a potential V(φ).

Friedmann Equations
If we consider an expanding homogeneous background and adopt the Friedmann-Lemaître-Robertson-Walker (FLRW) metric, the Einstein field equations reduce to the Friedmann equations where H =ȧ/a is the Hubble parameter, a(t) is the scale factor, Λ is the cosmological constant, and k determines the curvature of space. The universe is closed if k > 0, flat if k = 0, or open if k < 0. Equation (9) can also be written as Combining Equations (8) and (9), we obtain the energy conservation equation This equation can be directly obtained from the conservation of the energy-momentum tensor D µ T µν = 0 which is included in the Einstein equations through the contracted Bianchi identities. It can be rewritten as Introducing the volume V ∝ a 3 and the energy E = ρc 2 V, Equation (12) takes the form dE = −PdV. It can be interpreted as the first principle of thermodynamics for an adiabatic evolution of the universe dS = 0 [41]. The equation of state parameter is defined by According to Equation (11), the energy density decreases with the scale factor if w > −1 and increases with the scale factor if w < −1 (it is constant if w = −1). The case where the energy density increases with the scale factor corresponds to a phantom universe. The deceleration parameter is defined by The universe is decelerating if q > 0 and accelerating if q < 0. In this paper, we consider a flat universe (k = 0) in agreement with the inflation paradigm [1] and the observations of the cosmic microwave background (CMB) [19,20]. On the other hand, we set the cosmological constant to zero (Λ = 0) because, in our model, the acceleration of the expansion of the universe will be taken into account in the equation of state of the scalar field. The Friedmann equations then reduce to the form In that case, the deceleration parameter is related to the equation of state parameter by The universe is decelerating if w > −1/3 and accelerating if w < −1/3. When w = −1/3 the scale factor increases linearly with time.
We assume that the universe contains different X-fluids of density ρ X and a scalar field of density ρ φ . The total density and the total pressure are ρ = ∑ X ρ X + ρ φ and P = ∑ X P X + P φ . The Friedmann Equation (16) can then be written as

X-Fluids
We assume that the X-fluids are described by a linear equation of state of the form where α X is constant. We also assume that these fluids are independent from each other (and from the scalar field) so that they individually satisfy the energy conservation equation Solving this equation with the equation of state (20), we find where ρ 0 c 2 denotes the present value of the energy density of the universe and Ω X,0 denotes the proportion of the X-fluid at the present time (a = 1). Assuming that the X-fluid is alone in the universe (or dominates the other species) and solving the Friedmann Equation (16) with the density from Equation (22), we find that the scale factor and the density of the X-fluid evolve with time as if α X > −1, and as if α X = −1.
For the radiation (α r = 1/3): For baryonic matter (α b = 0): For dark matter (α dm = 0): For stiff matter (α s = 1): For the dark energy (α de = −1): The equation of state of the dark energy (vacuum energy) leads to a constant energy density implying a phase of exponential expansion (de Sitter). However, we will not need the equation of state of the dark energy in our model because the acceleration of the expansion of the universe will be taken into account in the equation of state (or in the potential) of the scalar field.

Canonical Scalar Field
A canonical scalar field minimally coupled to gravity evolves according to the Klein-Gordon (KG) equation 8φ where V(φ) is the potential of the scalar field. The scalar field tends to run down the potential towards lower energies. The density and the pressure of the scalar field are given by We can check that these equations imply the energy conservation Equation (see Appendix A) When the kinetic term dominates the potential term (φ 2 /2 V(φ)), we obtain the equation of state P φ ∼ ρ φ c 2 of stiff matter. This is the kination regime [63]. When the potential term dominates the kinetic term (φ 2 /2 V(φ)), we obtain the equation of state P φ ∼ −ρ φ c 2 of vacuum or dark energy. This is the inflation (de Sitter) regime [15]. Inversely, if the scalar field is described by the stiff equation of state P φ = ρ φ c 2 , we find V(φ) = 0 and φ = (c 2 /12πG) 1/2 ln t + cst corresponding to the kination regime. If the scalar field is described by the vacuum or dark energy equation of state P φ = −ρ φ c 2 , we findφ = 0, i.e., φ = cst and V(φ) = cst corresponding to the de Sitter regime (see Section 5.6).

General Equations
We assume that the scalar field is described by a quadratic equation of state of the form [39][40][41][42][43][44][45][46][47] where ρ P = 5.16 × 10 99 g m −3 is the Planck density (see footnote 2 ) and ρ Λ = 5.96 × 10 −24 g m −3 is the cosmological density. The linear term, characterized by a constant α, may represent radiation (α = 1/3), dark matter (α = 0 or α 0), stiff matter (α = 1), or even be a new species. 10 The equation of state parameter w φ = P φ /(ρ φ c 2 ) and the squared speed of sound Solving the energy conservation Equation (32) with the equation of state (33), we find that the density of the scalar field is given in excellent approximation by (see Appendix B) where a 1 is a constant of integration. To obtain Equation (36), we have used the fact that ρ Λ /ρ P 1. This equation is essentially exact as ρ Λ /ρ P ∼ 10 −123 is extremely small. It can be rewritten as where ρ 0 c 2 is the present value of the energy density and we have introduced the notations Ω P,0 ≡ ρ P /ρ 0 and Ω Λ,0 ≡ ρ Λ /ρ 0 (Ω Λ,0 represents the present proportion of dark energy in the universe while Ω P,0 is just a convenient notation). After the period of inflation (see below), we can make the approximation In this manner, we see that Ω P,0 a 3(α+1) 1 represents the present proportion Ω α,0 of the α-fluid in the universe. Therefore, the constant a 1 is given by We can then rewrite Equation (37) under the equivalent form Combining Equations (33)-(36) and using again the fact that ρ Λ /ρ P 1, we obtain in excellent approximation We note that the denominator in Equation (43) can be replaced by ρ Λ (a/a 1 ) 6(α+1) + ρ P (a/a 1 ) 3(α+1) + ρ P with the same degree of approximation.
If the universe contains only the scalar field, the deceleration parameter q is given by Equation (18). Using Equations (34) and (43), we obtain and The complete analytical solution a(t) of the Friedmann Equation (16) with Equation (36) is given by [42,44] with t P = 1/ Gρ P = (hG/c 5 ) 1/2 = 5.39 × 10 −44 s (Planck time). This solution describes the evolution of the universe from an early de Sitter era to a late de Sitter era bridged by a decelerating algebraic expansion (α-era). For the particular values α = 0, 1/3, 1, the density evolves with the scale factor as To better understand the physical meaning of the different terms involved in the foregoing equations, we successively study how these equations can be simplified in the early and in the late universe. We only give the equations that will be needed in the subsequent analysis, and we refer to our previous papers [39][40][41][42][43][44][45][46][47] for a more thorough discussion.

Early Universe
In the early universe, where the density is high, the equation of state (33) reduces to This amounts to neglecting the contribution of dark energy (ρ Λ = 0). The equation of state parameter and the squared speed of sound are given by Solving the energy conservation Equation (32) with the equation of state (49), we find that the density of the scalar field is It can be rewritten as Then, we find P φ If the universe contains only the scalar field, the deceleration parameter is given by When a a 1 , the density of the scalar field given by Equation (52) tends to a constant ρ φ ρ P equal to the Planck density. This leads to a phase of early inflation where the scale factor increases exponentially rapidly with time (early de Sitter era). When a a 1 , the density of the scalar field decreases algebraically as ρ φ /ρ 0 = Ω α,0 /a 3(α+1) . In that case, it behaves as an α-fluid with a linear equation of state P φ = αρ φ c 2 . The scale factor increases algebraically as t 2/[3(1+α)] (the expansion is decelerating if α > −1/3). The equation of state (49) thus describes the smooth transition between a phase of inflation and an α-era in the early universe. The transition takes place at a a 1 . This equation of state is studied in detail in [41,44]. The temporal evolution of the scale factor a(t) is given analytically by For the particular values α = 0, 1/3, 1, the density evolves with the scale factor as If we follow approach (A) of the Introduction, it is relevant to take α = 1/3 (radiation) or possibly α = 1 (stiff fluid) in the early universe.

Late Universe
In the late universe, where the density is low, the equation of state (33) reduces to This amount to neglecting quantum effects (ρ P → +∞). The equation of state parameter and the squared speed of sound are given by Solving the energy conservation Equation (32) with the equation of state (61), we find that the density of the scalar field is We note that a 1 a 2 Equation (64) can be rewritten as Then, we find P φ If the universe contains only the scalar field, the deceleration parameter is given by When a a 2 , the density of the scalar field given by Equation (64) tends to a constant ρ φ ρ Λ equal to the cosmological density. This leads to a phase of late accelerating expansion (or late inflation) where the scale factor increases exponentially rapidly with time (late de Sitter era). When a a 2 , the density of the scalar field decreases algebraically as ρ φ /ρ 0 = Ω α,0 /a 3(α+1) . In that case, it behaves as an α-fluid with a linear equation of state P φ = αρ φ c 2 . The scale factor increases algebraically as t 2/[3(1+α)] (the expansion is decelerating if α > −1/3). The equation of state (61) thus describes the transition between an α-era and a phase of accelerating expansion (dark energy) in the late universe. The transition takes place at a a 2 . This equation of state is studied in detail in [42,44]. The temporal evolution of the scale factor a(t) and density ρ φ (t) is given analytically by with t Λ = 1/ Gρ Λ = (8π/Λ) 1/2 = 1.46 × 10 18 s (cosmological time). For the particular values α = 0, 1/3, 1, the density evolves with the scale factor as If we follow approach (A) of the Introduction, it is relevant to take α = 0 (pressureless matter) in the late universe. In that case, the equation of state is constant P φ = −ρ Λ c 2 . This UDM model is equivalent to the ΛCDM model not only to 0-th order in perturbation theory (background) but to all orders, even in the nonlinear clustering regime [42,53,54].

Intermediate Regime
In the intermediate regime, valid after the early inflation and before the late accelerating expansion of the universe, the equation of state of the scalar field reduces to This amounts to neglecting quantum effects (ρ P → +∞) and dark energy (ρ Λ = 0). The equation of state parameter and the squared speed of sound are given by Solving the energy conservation Equation (32) with the equation of state (73), we find that the density of the scalar field is If the universe contains only the scalar field, the deceleration parameter is given by The temporal evolution of the scale factor a(t) and density ρ φ (t) is given analytically by The expansion is decelerating for α > −1/3. For the particular values α = 0, 1/3, 1, the density evolves with the scale factor as

Complete Evolution of the Universe
Regrouping the foregoing results, we see that the quadratic equation of state (33) describes in a unified manner an early inflation followed by a phase of decelerated expansion (α-era) and, finally, a late accelerating expansion. This equation of state is studied in detail in [44]. As discussed in the introduction, the equation of state has a drawback in the sense that it cannot describe both the radiation and matter eras. To solve this problem, we have proposed two approaches.
If we follow approach (A), we have to consider that α changes in the course of time (i.e., with the density of the universe). In that case, we are led to considering the equation of state (49) in the early universe with α = 1/3 (or α = 1) and the equation of state (61) in the late universe with α = 0. We can then solve the Friedmann Equation (16) with Equation (53) or (66) in these two epochs and then match the results to obtain the whole evolution of the universe. This is equivalent to solving the Friedmann equation with the density 11 This is the procedure adopted in [39,40,42,47]. Equations (80) and (81) determine the evolution of the scale factor a(t). They describe successively the phases of inflation, radiation, pressureless matter, and dark energy. The ΛCDM model is recovered by ignoring the early inflation, i.e., by taking Ω P,0 → +∞ (i.e., ρ P → +∞ orh → 0) in Equation (81). If we follow approach (B), we can consider that α is fixed in the quadratic equation of state (33) of the scalar field (the value α = 1/3 or α = 1 may be the most relevant). Then, we have to account for the presence of other species such as stiff matter, radiation, baryonic matter, and dark matter. Finally, we have to solve the Friedmann Equation (19) with Equation (40) for the scalar field and Equations (25)- (28) for the other species (except the one that has been incorporated in the equation of state of the scalar field). In the simplest model, ignoring stiff matter, taking α = 1/3 in the equation of state of the scalar field and treating dark matter as an independent species with α X = 0, this leads again to Equations (80) and (81). Alternatively, taking α = 1 in the equation of state of the scalar field and treating radiation and dark matter as independent species with α X = 1/3 and α X = 0 respectively, we find Equations (80) and (82) determine the evolution of the scale factor a(t). They describe successively the phases of inflation, stiff matter, radiation, pressureless matter, and dark energy. As in footnote 11 , in order to avoid a spurious divergence of the energy density at a = 0, the radiation and matter component terms Ω r,0 /a 4 and Ω m,0 /a 3 have to be introduced at a sufficiently late time, i.e., after the inflation era when ρ ρ P . The possibility that a radiation (or a stiff matter) era occurs before the inflation era, leading to a big bang singularity, is considered in Appendix E. See also Ref. [47] for more general models generalizing Equations (81) and (82).
The complete history of the universe determined by Equations (80) and (81) has been described in [39,40,42,44,47] (see in particular Figure 14 of [42]). Although the evolution of the background density is the same in approaches (A) and (B), the scalar field potential is different in these two approaches, especially in the late universe, as discussed in Sections 4 and 5 below.

Scalar Field in the Presence of X-Fluids
In this section, we derive the potential V(φ) of a real scalar field characterized by an equation of state P φ (ρ φ ) in the presence of other species. We provide general results using standard techniques [64,65] and apply them to the quadratic equation of state (33) of our model in the presence of X-fluids (radiation, matter, stiff matter. . .).

General Results
Let us consider a canonical scalar field defined by Equations (30)- (32). From Equation (31), we findφ where we have introduced the equation of state parameter w φ = P φ /ρ φ c 2 . Usingφ = (dφ/da)Ha and the Friedmann Equation (19), we find that the relation between the scalar field and the scale factor is given by 12 On the other hand, according to Equation (31), the potential of the scalar field is given by Therefore, we find that the potential of the canonical scalar field in the presence of other species is determined in parametric form by the equations where we have taken the origin of the scalar field (φ = 0) at a = 0.

Quadratic Equation of State
For the quadratic equation of state (33), using the results of Section 3, we have in excellent approximation (using the fact that ρ Λ /ρ P 1): On the other hand, for the X-fluids we have As a result, we find that the scalar field potential of the vacuumon unifying the early inflation and the late accelerating expansion of the universe in the presence of X-fluids is determined by the equations In the early universe (Ω Λ = 0), we obtain φ(a) = 3c 2 8πG In the late universe (Ω P → +∞), we obtain 13 φ(a) = 3c 2 8πG These equations determine the potential of the vacuumon in the presence of X-fluids. Although it is not possible to obtain the scalar field potential V(φ) analytically in the general case, it can always be obtained numerically. 14 The scalar field potential V(φ) can be obtained analytically in the absence of other species (see Section 5). In the presence of Xfluids, it can be obtained analytically in particular situations of physical interest depending on the approach considered.
In approach (A), the value of α changes with the epoch during the evolution of the universe. It is equal to α = 1 (stiff matter) or α = 1/3 (radiation) in the early universe and to α = 0 (matter) in the late universe. In this approach, we do not have to take into account the presence of other species, as their effect has been incorporated in the equation of state of the scalar field. Therefore, we can consider that the scalar field is alone in the universe but that α has a different value in the early and in the late universe. The scalar field potential is determined by Equations (96) and (97) with Ω X,0 = 0 and α = 1 or α = 1/3 in the early universe (leading to Equations (119) and (120)) and by Equations (98) and (99) with Ω X,0 = 0 and α = 0 in the late universe (leading to Equation (137)). These cases are treated in Sections 5.2 and 5.3.
In approach (B), the value of α is fixed. We then have to consider two situations: (i) We first consider a universe without stiff matter. In that case, it is relevant take α = 1/3 in the quadratic equation of state (33) so that the scalar field accounts for the inflation era, the radiation era and the late acceleration of the universe. Then, we have to add baryonic matter and dark matter as independent species (they can be treated as a single species with α X = 0). In the early universe, we can approximate the equation of state of the scalar field by Equation (49). As matter is negligible in the early universe, we can consider that the scalar field is alone in the universe at that epoch. The scalar field potential is then determined by Equations (96) and (97) with Ω X,0 = 0. This situation, which describes the transition between the inflation era and the radiation era, is treated in Section 5.2. It leads to the same potential (120) as in approach (A).
In the late universe, we can approximate the equation of state of the scalar field by Equation (61). On other hand, we have to take into account the presence of matter as an independent species. In that case, the scalar field potential is determined by Equations (98) and (99) with α X = 0. Unfortunately, the integral in Equation (98) cannot be performed analytically (see Section 8). However, at very late time, the contribution of matter becomes negligible and the scalar field is alone in the universe.
In that case, its potential is given by Equation (136). We note that the scalar field potential in the late universe in approach (B) is very different from its expression in approach (A) as it corresponds to Equation (132) with α = 1/3 instead of α = 0. Finally, in the intermediate era between the early inflation and the late accelerating expansion of the universe, we can approximate the equation of state of the scalar field by Equation (73). We also have to take into account the presence of matter as an independent species. The scalar field potential is then determined by Equations (96) and (97) with Ω P,0 → +∞ and α X = 0 or by Equations (98) and (99) with Ω Λ,0 = 0 and α X = 0. This situation, which describes the period where the universe contains radiation and matter, is treated in Section 7 leading to Equation (188) with α = 1/3 and α X = 0. (ii) We now consider a universe with stiff matter. In that case, it is relevant to take α = 1 so that the scalar field accounts for the inflation era, the stiff matter era and the late acceleration of the universe. Then, we have to add radiation (α X = 1/3) and matter (α X = 0) as independent species. In the early universe, we can approximate the equation of state of the scalar field by Equation (49). As radiation and matter are negligible in the (very) early universe, we can consider that the scalar field is alone in the universe at that epoch. The scalar field potential is then determined by Equations (96) and (97) with Ω X,0 = 0. This situation, which describes the transition between the inflation era and the stiff matter era, is treated in Section 5.2. It leads to the same potential (119) as in approach (A). 15 In the late universe, we can approximate the equation of state of the scalar field by Equation (61). On other hand, we have to take into account the presence of matter (α X = 0) as an independent species. In that case, the scalar field potential is determined by Equations (98) and (99) with α X = 0. The integral in Equation (98) can be performed analytically (see Section 8). However, the potential is independent of φ and given by Equation (135) as when the scalar field is alone in the universe (which is the case at very late time when the contribution of matter becomes negligible). We note that the scalar field potential in the late universe in approach (B) is very different from its expression in approach (A) as it corresponds to Equation (132) with α = 1 instead of α = 0. Finally, in the period between the early inflation and the late accelerating expansion, we can approximate the equation of state of the scalar field by Equation (73). We also have to take into account the presence of radiation and matter as independent species. The scalar field potential is then determined by Equations (96) and (97) with Ω P,0 → +∞ and α X = {1/3, 0} or by Equations (98) and (99) with Ω Λ,0 = 0 and α X = {1/3, 0}.
This situation describes the period where the universe contains stiff matter, radiation, and matter. In this general situation involving a scalar field (stiff matter) and two external fluids (radiation and matter), the scalar field potential cannot be obtained analytically. However, if we neglect the matter contribution (at sufficiently early times) or the radiation contribution (at sufficiently late times), we have just one external fluid (radiation or matter) and the scalar field potential can be obtained analytically. This situation, which describes the period where the universe contains stiff matter and radiation or stiff matter and matter, is treated in Section 7, leading to Equation (188) with α = 1 and α X = 1/3 or α X = 0. Note that even if the stiff matter (played by the scalar field) is subdominant in that period, it remains important because it will ultimately lead to the late accelerating expansion of the universe.
Noting that x = (a/a 1 ) 3(α+1)/2 , where a 1 is defined by Equation (39), and using Equation (106), the relation between the scalar field and the scale factor is We can then express all the parameters of Section 3.1 as a function of ψ instead of a. In particular, the energy density and the pressure of the scalar field are given by The temporal evolution of the scale factor and of the scalar field is discussed in detail in our previous papers [39][40][41][42][43][44][45][46][47].

Early Vacuumon (Inflaton)
In the early universe (Ω Λ = 0), Equations (102) and (103) reduce to φ(a) = 3c 2 8πG Introducing the notations from Equation (104), we find which corresponds to Equations (106) and (107) with κ = 0. This returns Equations (122) and (123) of [44]. Eliminating x between these equations we obtain [44] V This is the potential of the vacuumon in the early universe (playing the role of the inflaton). It is associated to the equation of state (49) valid in the early universe (see Section 8.1 of [42]). For ψ → 0, it can be expanded in Taylor series as For ψ → +∞, we obtain the exponential asymptotic behaviors We note that the coefficient α = 1 (stiff matter) plays a special role as it leads to a faster decay of the potential. In that case [45,46] V(ψ) = ρ P c 2 cosh 4 ψ (α = 1).
We can then express all the parameters of Section 3.2 as a function of ψ instead of a. In particular, the energy density and the pressure of the scalar field are given by The temporal evolution of the scalar field is discussed in detail in our previous papers [39][40][41][42][43][44][45][46][47]. According to Equations (59) and (122) we have Remark 1. For a → 0, we can simplify Equation (111) into yielding φ(a) = 2 3 3c 2 8πG Substituting Equation (126) into Equation (112) we obtain the quadratic potential This corresponds to the first term in the expansion of the potential from Equation (116).

Late Vacuumon (Quintessence)
In the late universe (Ω P → +∞), Equations (102) and (103) reduce to φ(a) = 3c 2 8πG These expressions are valid for sufficiently large values of a. Introducing the notations from Equations (104) and (105), we find which correspond to Equations (106) and (107) with x 1. Here, ψ max appears as a constant of integration that can be obtained by matched asymptotics (see below). We recover Equations (125) and (126) of [44]. Eliminating x between these equations we obtain [44] This is the potential of the vacuumon in the late universe (playing the role of the quintessence). It is associated to the equation of state (61) valid in the late universe (see Section 8.1 of [42]). For ψ → ψ max , it can be expanded in Taylor series as For ψ → −∞, we obtain the exponential asymptotic behavior For α = 1 (stiff matter), the scalar field potential is constant For α = 1/3 (radiation), we have For α = 0 (pressureless matter), we have The potential corresponding to α = 0 (relevant in the late universe) is plotted in Figure 2. This is the scalar field potential associated with the ΛCDM model viewed as a UDM model described by the constant equation of state P φ = −ρ Λ c 2 .  Noting that x √ κ = (a/a 2 ) 3(α+1)/2 , where a 2 is defined by Equation (64), and using Equation (130), the relation between the scalar field and the scale factor is (a 2 /a) 3(α+1)/2 = sinh(ψ max − ψ).
We can then express all the parameters of Section 3.3 as a function of ψ instead of a. In particular, the energy density and the pressure of the scalar field are given by The temporal evolution of the scalar field is discussed in detail in our previous papers [39][40][41][42][43][44][45][46][47]. According to Equations (71) and (138) we have Remark 2. For a → +∞, we can simplify Equation (128) into Substituting Equation (142) into Equation (129) we obtain the quadratic potential This corresponds to the first term in the expansion of the potential from Equation (133).

Remark 3.
For α = 1, the approximate potential from Equation (144) takes the particularly simple form

Intermediate Regime
In the intermediate regime valid after the early inflation and before the late accelerating expansion of the universe (Ω P → +∞ and Ω Λ = 0), the scalar field is described by the linear equation of state (73). In that case, Equations (102) and (103) reduce to φ(a) = 3c 2 8πG These equations can also be obtained from Equations (111) and (112) with Ω P → +∞ or from Equations (128) and (129) with Ω Λ = 0. Integrating Equation (146), we obtain where φ * is a constant of integration. Combining Equations (147) and (148), we obtain the exponential potential This result is equivalent to Equations (117) and (134) obtained long after the primordial inflation or long before the late accelerating expansion of the universe. This determines the constant of integration ψ * = (1/2) ln(4Ω P,0 /Ω α,0 ). For α = 1 (stiff matter), we find corresponding to the regime of kination where the potential vanishes. Comparing Equation (150) with Equations (118) and (135) we see that the limits α → 1 and φ → ±∞ do not commute. For α = −1 (vacuum energy), we find corresponding to the de Sitter regime whereφ = 0. Combining Equations (73), (76), and (148), we can express the density and the pressure as a function of the scalar field as According to Equations (78) and (148), the temporal evolution of the scalar field is

Constant Scalar Field Potential
We have found different solutions corresponding to a constant scalar field potential (see Equations (135), (150) and (151)): (i) For the equation of state P φ = ρ φ c 2 of stiff matter (α = 1), we have (see Sections 3.4 and 5.5) This solution describes a pure stiff matter era.
This solution describes a pure de Sitter era. (iii) For the affine equation of state P φ = ρ φ c 2 − 2ρ Λ c 2 , we have (see Sections 3.3 and 5.3) This solution describes a stiff matter era followed by a de Sitter era (late inflation). 16 Therefore, several laws of evolution may correspond to the same potential.

Parameters of the Scalar Field
In this section, we determine the parameters (mass, self-interaction constant, and scattering length) of the scalar field in the early and late universe by comparing the form of the potential close to its maximum and minimum with the normal form of a quartic potential given by Equation (A54).

Early Vacuumon (Inflaton)
In the early universe, the potential of the vacuumon is given by Equation (115). For ψ → 0, it can be expanded up to fourth order yielding Equation (116). Comparing Equation (116) with Equation (A54), we obtain the following results (see Appendix D): (i) The maximum value of the potential (corresponding to ψ = 0) is equal to the Planck energy density where ρ P = c 5 /hG 2 = 5.16 × 10 99 g m −3 is the Planck density. (ii) The squared mass of the scalar field is given by with f e (α) = − 9 4 (α + 1)(α + 3).
We find g e = 9 for α = 0, g e = 512/27 for α = 1/3, and g e = 56 for α = 1. The self-interaction constant of the early vacuumon (inflaton) is positive (g e > 0), corresponding to a repulsive self-interaction of order 1. (iv) The dimensional self-interaction constant of the scalar field is given by We find g e /| f e | = 4/3 for α = 0, g e /| f e | = 256/135 for α = 1/3, and g e /| f e | = 28/9 for α = 1. The dimensional self-interaction constant of the early vacuumon (inflaton) is of order Gh 2 /c 2 = 5.15 × 10 −71 eV cm 3 . (v) The scattering length of the bosons associated with the scalar field is given by We find g 2 e /| f e | = 12 for α = 0, g 2 e /| f e | = 131072/3645 for α = 1/3, and g 2 e /| f e | = 1568/9 for α = 1. The scattering length a s of the bosons associated with the early vacuumon (inflaton) is of the order of the Planck length l P = GM P /c 2 = (Gh/c 3 ) 1/2 = 1.62 × 10 −35 m which corresponds to the semi Schwarzschild radius associated with the Planck mass M P .
In theories of extended supergravity, the mass of a scalar field is quantized according to the rule [66] where n * is an integer and H is the Hubble factor in the de Sitter era. As H 2 = 8πGρ P /3 during the early inflation, we can rewrite Equation (166) as The quantum of mass (n * = 1) in the early universe is the rescaled Planck mass M * P = √ 8π/3M P . Comparing Equation (167) with Equation (160), the quantization rule implies that f e (α) should be an integer n * . We see that this quantization rule is realized for the integer value n * = −10 when α = 1/3 (radiation) and for the integer value n * = −18 when α = 1 (stiff) [47,67]. By contrast, f e (α) is not an integer when α = 0 (matter). This may be connected to the fact that the index α = 0 is not justified in the early universe as there is no matter there in the usual sense.

Late Vacuumon (Quintessence)
In the late universe, the potential of the vacuumon is given by Equation (132). For ψ → ψ max , it can be expanded up to fourth order yielding Equation (133). 17 Comparing Equation (133) with Equation (A54), with the substitution φ → φ max − φ, we obtain the following results (see Appendix D): (i) The minimum value of the potential (corresponding to ψ = ψ max ) is equal to the cosmological energy density where ρ Λ = Λ/8πG = 5.96 × 10 −24 g m −3 is the cosmological density. (ii) The squared mass of the scalar field is given by We find f l = 9/4 for α = 0, f l = 2 for α = 1/3, and f l = 0 for α = 1. The mass of the late vacuumon (quintessence) is real ( f l > 0) and of the order of the cosmon mass m Λ =h √ Λ/c 2 = 2.08 × 10 −33 eV/c 2 , except for α = 1 (stiff matter) where m = 0. Our approach provides therefore a physical interpretation of the cosmon mass as being the mass of the scalar field responsible for the dark energy in the late universe. To the best of our knowledge, this interpretation has not been given before. (iii) The dimensionless self-interaction constant of the scalar field is given by We find g l = 1 for α = 0, g l = 32/27 for α = 1/3, and g l = 0 for α = 1. The self-interaction constant of the late vacuumon (quintessence) is positive (g l > 0), corresponding to a repulsive self-interaction of order ρ Λ /ρ P ∼ 10 −123 , except for α = 1 (stiff matter) where λ = 0. (iv) The dimensional self-interaction constant of the scalar field is given by We find g l / f l = 4/9 for α = 0, g l / f l = 16/27 for α = 1/3, and g l / f l → 8/9 for α → 1. The dimensional self-interaction constant of the late vacuumon (quintessence) is of order Gh 2 /c 2 = 5.15 × 10 −71 eV cm 3 , except for α = 1 (stiff matter) where λ s = 0 (note that the result obtained in the limit α → 1 is different from the result obtained for α = 1). (v) The scattering length of the bosons associated with the scalar field is given by We find g l / f l = 2/3 for α = 0, g l / f l = 16 √ 2/27 for α = 1/3, and g l / f l → 0 for α → 1. The scattering length a s of the bosons associated with the late vacuumon (quintessence) is of the order of the cosmon radius r Λ = Gm Λ /c 2 = Gh √ Λ/c 4 = 2.75 × 10 −96 m, which corresponds to the semi Schwarzschild radius associated with the cosmon mass m Λ , except for α = 1 (stiff matter) where a s = 0.
As H 2 = 8πGρ Λ /3 during the late inflation, we can rewrite the quantization rule (166) as The quantum of mass (n * = 1) in the late universe is the rescaled cosmon mass m * Λ = 1 √ 3 m Λ . 18 Comparing Equation (169) with Equation (175), the quantization rule implies that f l (α) should be an integer n * . We see that this quantization rule is realized for the integer value n * = 2 when α = 1/3 (radiation) and for the integer value n * = 0 when α = 1 (stiff matter) [47,67]. By contrast, f l (α) is not an integer when α = 0 (matter). We recall that our model of late universe with α = 0 is equivalent to the ΛCDM model interpreted as a UDM model. Therefore, the quintessence field associated with the standard ΛCDM model (α = 0) does not satisfy the quantization rule (166) [47,67]. This suggests (provided that this quantization rule must hold, which is not firmly established) that the relevant value of α to take is α = 1/3 or α = 1, and that dark matter should be treated as an independent species, as in scenario (B).

Remark 4.
In this paper, we have described the late universe by a polytrope of index n = −1 [42]. For α = 0 and n = −1, our model is equivalent to the ΛCDM model. If we consider a more general polytropic equation of state with α = 0 and an arbitrary index n < 0 (see Equation (A45) and Ref. [42]), we find that the mass of the scalar field in the late universe is given by (see Appendix D and [47,67]) We note that this mass presents a maximum as a function of n for n = −1. Interestingly, this maximum selects the ΛCDM model among all the polytropic models of the form of Equation (A45) with α = 0 and n < 0 [47,67].

General Results
In this section, we consider the complete model of Section 4 in the intermediate regime between the two de Sitter eras. In this regime, the scalar field has a linear equation of state P φ = αρ φ c 2 . We use approach (B) in which the value of α is fixed. As we have seen, relevant values of α are α = 1 (stiff matter) or α = 1/3 (radiation). On the other hand, we have to take into account the presence of X-fluids which also have a linear equation of state P X = α X ρ X c 2 . If the scalar field evolves in the presence of radiation, one has α r = 1/3. If it evolves in the presence of matter, one has α m = 0. For the sake of generality, we consider arbitrary values of α and α X but we assume that α ≥ −1 and α X ≥ −1 so that the universe has not a phantom behavior.
The energy density of the scalar field and of the X-fluids evolve with the scale factor as The total energy density is therefore The time evolution a(t) of the scale factor is determined by the Friedmann equation where the right hand side is given by Equation (178). 19 On the other hand, the scalar field potential is determined by the equations 20 φ(a) = 3c 2 8πG In the following, in order to have analytical results, we assume that there is just one X-fluid in addition to the scalar field. It can represent for example radiation (α X = 1/3) or matter (α X = 0). The case of two X-fluids (radiation and matter) in addition to the scalar field could be treated numerically. We obtain the form of the potential in the limiting case where one perfect fluid dominates over the other. In our model, the foregoing equations are valid only in the intermediate regime between the two de Sitter eras, where the equation of state of the scalar field can be approximated by a linear law, so they have a restricted domain of validity. However, in our mathematical analysis below, we shall study Equations (177)-(181) for all times because they can be relevant in other contexts, for example in the case where the fluid describes dark matter and the scalar field describes dark energy (see Section 7.3).

Hyperbolic Potential
We first assume α < α X . In that case, the X-fluid dominates at early times and the scalar field dominates at late times. Equation (180) can be rewritten as φ(a) = 3c 2 8πG . (182) Making the change of variables we obtain φ(a) = 3c 2 8πG Using the identity we find φ(a) = 3c 2 8πG Inversely, Combining Equation (187) with Equation (181), we obtain the hyperbolic potential We now assume α > α X . In that case, the scalar field dominates at early times and the X-fluid dominates at late times. We note that the integral from Equation (182) diverges at x = 0. We can easily circumvent this problem by redefining the additive constant in the scalar field so that . (189) With this convention, we can check that the previous equations (186)-(188) remain valid. Finally, when α = α X , Equation (180) reduces to φ(a) = 3c 2 8πG where we have chosen the origin of the scalar field at a = 1 in order to avoid spurious divergences. This yields φ(a) = 3c 2 8πG Inversely, Combining Equation (192) with Equation (181) we obtain the exponential potential 1+α Ω X,0 +Ω α,0 Ω α,0 For a scalar field alone (Ω X,0 = 0) we recover Equation (149). Below, we consider particular limits of the hyperbolic potential (188).

Power-Law Potential
When the terms in brackets in Equations (186)-(188) go to zero, we find φ(a) = 3c 2 8πG This corresponds to the situation where the X-fluid dominates the scalar field. In that case, we obtain a power-law potential. When α < α X , this approximation is valid at sufficiently early times (R → 0, φ → 0 + and V → +∞). When α > α X , it is valid at sufficiently late times (R → +∞, φ → 0 − and V → 0). Pure (inverse) power-law potentials were introduced by Peebles and Ratra [30,31]. They were used in relation to tracker solutions [22,23,33].

Exponential Potential
When the terms in brackets in Equations (186)-(188) go to +∞, we find φ(a) = 3c 2 8πG This corresponds to the situation where the scalar field dominates the X-fluid. In that case, we obtain an exponential potential as when the scalar field with a linear equation of state is alone in the universe (see Section 5.5). 21 When α < α X , this approximation is valid at sufficiently late times (R → +∞, φ → +∞ and V → 0). When α > α X , it is valid at sufficiently early times (R → 0, φ → −∞ and V → +∞). Pure exponential potentials were introduced by Halliwell [69] (see also [31,[70][71][72][73][74][75]). Exponential potentials arise very naturally in all models of unification with gravity such as Kaluza-Klein theories, supergravity theory or string theory. Most theories undergoing dimensional reduction to an effective four-dimensional theory yield exponential potentials.

Comparison with Other Works
The hyperbolic potential (188) has been obtained by several authors [18,[76][77][78][79][80][81] using different methods. In these models, the fluid characterized by an equation of state P X = α X ρ X c 2 with 0 ≤ α X ≤ 1 describes radiation (α X = 1/3) or dark matter (α X = 0) and the scalar field characterized by an equation of state P φ = αρ φ c 2 with −1 < α ≤ −1/3 (e.g., α = −0.65 [79] or α = −2/3 [80]) describes DE. As α < α X these models describe a radiation or dark matter era followed by a dark energy era. During the radiation or dark matter era, the universe is decelerating and the scale factor increases algebraically as t 1/2 or t 2/3 (see Equation (23)). During the dark energy era, the universe is accelerating and the scale factor increases algebraically as t 2/[3(1+α)] (power-law inflation). In the radiation and dark matter eras, the scalar field is subdominant and the potential V(φ) has an inverse power-law behavior (tracking). In the dark energy era, the scalar field dominates the other species (X-fluids) and the potential V(φ) has an exponential behavior as when a scalar field with a linear equation of state is alone in the universe. We thus achieve a tracker solution that can drive the universe into its current inflationary state. The behavior of the potential as an inverse power-law potential at early times avoids the fine-tuning problem and the cosmic coincidence problem [22]. Its exponential behavior at late time drives the universe into a power-law inflationary stage in agreement with the observations. In this model, the scalar field describes only dark energy, and we have to add radiation and (dark) matter as additional species.
In the present work, the hyperbolic potential (188) is obtained in a particular limit of a more general model based on a scalar field with a quadratic equation of state in the presence of X-fluids. It is valid in the intermediate regime between the early and late inflation. The scalar field models dark radiation (α = 1/3) or stiff matter (α = 1) and the fluid models radiation (α X = 1/3) or matter (α X = 0). As α ≥ α X , we are always in the situation where the scalar field dominates at early times and the X-fluid dominates at late times (note that the scalar field dominates again in the late universe due to the constant term −(α + 1)ρ Λ c 2 in its equation of state which has been ignored in the intermediate regime). Therefore, in this intermediate period, the scalar field potential is initially exponential and then becomes algebraic.

Remark 5.
The ΛCDM model corresponds to α X = 0 and α = −1. In that case, φ = cst and V(φ) = cst. The ΛCDM model can also be obtained by taking α X = −1 and α = 0. In that case, the scalar field has a potential

Scalar Field in the Presence of One Fluid in the Late Universe
In this section, we consider the complete model of Section 4 in the late universe (we consider just one X-fluid in addition to the scalar field). In that case, the scalar field has an affine equation of state P φ = αρ φ c 2 − (α + 1)ρ Λ c 2 . The scalar field potential is determined by Equations (98) and (99) which can be rewritten as φ(a) = 3c 2 8πG The scalar field is described by the index α = 1/3 (radiation) or α = 1 (stiff) and the X-fluid is described by the index α X = 0 (matter). Unfortunately, the integral cannot be calculated analytically in general except in the degenerate case α = α X (see below) which is not relevant in our situation. However, in the very late universe, we can ignore the term in Ω X,0 in Equation (201). In that case, we are led back to the equations (128) and (129) valid for a scalar field alone in the universe leading to the potential from Equation (132). Note that for α = 1 (stiff matter), the potential from Equation (202) is constant: whether or not an X-fluid is present. In addition, for α = 1 and α X = 0 the integral in Equation For a 1, it reduces to This expression provides the correction to Equation (128) due to the presence of the X-fluid.

Remark 6.
It is interesting at a general level to consider the case α = α X in Equation (201) where the integral can be calculated analytically. This corresponds to a situation where a scalar field described by the affine equation of state P φ = αρ φ c 2 − (α + 1)ρ Λ c 2 evolves in the presence of an X-fluid described by the linear equation of state P X = α X ρ X c 2 with α X = α. 22 In that case, Equation (201) reduces to φ(a) = − 3c 2 8πG where we have chosen the origin of the scalar field at a → +∞ in order to avoid spurious divergences at x = 0. 23 Making the change of variables we obtain Using the identity from Equation (185), we find φ(a) = − 3c 2 8πG argsinh Inversely, a = Ω X,0 + Ω α,0 Ω Λ,0 Finally, combining Equation (210) with Equation (202), we obtain the hyperbolic potential For a scalar field alone (Ω X,0 = 0), we recover Equation (132).

Spectrum of Fluctuations in the Primordial Universe
In this work and in previous ones [39][40][41][42][43][44][45][46][47], we have developed a model of early inflation based on the quadratic equation of state (49). This model is able to describe the evolution of the homogeneous background and to account for a primordial phase of exponential expansion (de Sitter era) followed by a stiff matter era (α = 1) or a radiation era (α = 1/3). It also provides a graceful exit to the de Sitter era. However, explaining the evolution of the homogeneous background is not sufficient. A relevant model of inflation must also reproduce the observed spectrum of fluctuations in the primordial universe. In this section, we apply the Hamilton-Jacobi formalism of inflation [60] to a scalar field described by the equation of state (49) and we determine the Hubble hierarchy parameters and the spectral indices. We show that the obtained results are in severe disagreement with the observations. These negative results suggest that the vacuumon potential V(φ) is just an effective classical potential that cannot be used to compute the fluctuations in the primordial universe. A fully quantum mechanics approach may be required.

Hamilton-Jacobi Formalism
Let us first briefly review the Hamilton-Jacobi formalism of inflation following Ref. [60]. This formalism is very general. In particular, it does not rely on the slow roll approximation. This is important because, as we shall see, our model of inflation is never in the slow roll regime. 24 In this section, we assume that the scalar field is alone in the primordial universe. The basic equations of the problem are the KG equation and the Friedmann equation where is the energy density of the scalar field. In the following, we assume that the function H(φ) is known. It is called the generating function. Specifying H(φ) is equivalent to specifying the equation of state P φ (ρ φ ) or the potential V(φ) of the scalar field. Taking the time derivative of Equations (213) and (214), combining the results with Equation (212), where denotes derivative with respect to φ. Using aH =ȧ = a φ and Equation (215), we obtain Integrating this equation, we find that This equation determines the scale factor as a function of φ. On the other hand, combining Equations (213)-(215), we find This equation determines the scalar field potential. In the slow-roll approximation wherė φ 2 V(φ) and |φ| V (φ), the scalar field potential is given by The first Hubble hierarchy parameter is defined by where we have used Equation (216) to obtain the second equality. Usingä = (d/dt)(Ha) = aḢ + Hȧ, implyingä/a = H 2 +Ḣ = H 2 (1 +Ḣ/H 2 ) = H 2 (1 +Ḣa/Hȧ), we find thaẗ Therefore, H is related to the deceleration parameter from Equation (14) by The parameter H gives information about the acceleration of the universe. During inflation (ä > 0), we have H < 1. Inflation ends whenä = 0 yielding H (t end ) = 1. The number of e-folding before inflation ends is defined by It can be written in terms of the scalar field as where we have used Equation (215) to obtain the second equality and Equation (220) to obtain the last equality. Here, φ end is the value of the scalar field at the end of inflation. It is determined by the condition H (φ end ) = 1. The pressure of the scalar field is given by Summing Equations (214) and (225), we find Substituting Equations (213) and (215) into Equation (226), we obtain This equation determines the pressure as a function of φ. Finally, the energy density ρ φ (φ) of the scalar field is given by Eliminating φ between Equations (227) and (228), we obtain the equation of state P φ (ρ φ ).

Application to Our Model of Inflation
We now apply the Hamilton-Jacobi formalism to our model of early inflation. The generating function of our model is 25 Using Equations (227) and (228), we obtain the pressure and the energy density This returns the results from Equation (123). Eliminating ψ between these equations, we obtain the equation of state This returns the result from Equation (49). Using Equation (218), we obtain the scalar field potential This returns the result from Equation (115). Using Equation (217), we obtain the scale factor This returns the result from Equation (122). Using Equations (223) and (234), we find that the number of e-folding before inflation ends is given by In our model, we find N i ∼ ln(a 1 /l P ) ∼ 67 [44]. Using Equations (220) and (222), we obtain the first Hubble hierarchy parameter and the deceleration parameter through the relation The condition H = 1 (i.e., q = 0) gives tanh 2 (ψ end ) = 2 3(α + 1) .
This equation determines the value ψ end of the scalar field at the end of the inflation. Therefore, inflation takes place between ψ = 0 and ψ end . The second and third Hubble hierarchy parameters are defined by and where we have used Equation (216) to obtain the second equality in Equation (238). For the generating function from Equation (229), we find and We note that η H → − 3 2 (α + 1) = 0 when ψ → 0, so we are never in the slow roll regime unless α −1 (which is a very peculiar situation [41]). The scalar spectral index n s is given by [60] n Using Equations (236) and (238), we obtain In our model, the scalar spectral index n s turns out to be independent of φ. The running scalar spectral index n run is given by [60] Using Equations (236), (238), and (241), we find In our model, the running scalar spectral index vanishes whatever the values of φ and α.
Comparison with Planck data: The values of the parameters of inflation obtained by the Planck collaboration [82] are n s = 0.9603 ± 0.0073, n run = −0.0134 ± 0.0090, r < 0.11 and N = 60 − 70. A value of n s close to 1 means that we have nearly scale-dependent density perturbations. For α = 1 (stiff matter), α = 1/3 (radiation) or α = 0 (matter), the value of n s predicted by our model is very far from unity. This implies that we are not in the slow roll regime where H and η H are small. For α = 1, we find n s = −5; for α = 1/3, we find n s = −3; and for α = 0, we find n s = −2. It is only for α −1 that n s becomes close to unity. However, in that case, we cannot describe the graceful exit to inflation (i.e., the smooth transition from inflation to a stiff matter or a radiation era). Therefore, for relevant values of α, our inflationary model cannot account for the fluctuations in the primordial universe. This is probably because the vacuumon potential V(φ) (see Equation (233)) is just an effective classical field that cannot be used to compute the spectrum of fluctuations in the primordial universe. It may be therefore necessary to develop a fully quantum model of inflation.

Conclusions
We have developed the model of universe introduced in our previous papers [39][40][41][42][43][44][45][46][47]. This model assumes that the universe is filled with an exotic fluid described by a quadratic equation of state. The presence of this fluid accounts for a phase of early inflation, followed by a decelerated expansion, and finally a late accelerating expansion. Therefore, our quadratic equation of state unifies inflation, relativistic or nonrelativistic matter, and dark energy in the spirit of a generalized Chaplygin gas. Stiff matter, radiation, and (dark or baryonic) matter may be incorporated in this exotic fluid at different periods of its evolution or treated as independent species. We have given the scalar field representation of this exotic fluid and determined its potential in the absence or presence of other species. For a pure scalar field the potential can be expressed in terms of the Jacobian Elliptic function. We have shown that the mass of the early vacuumon (inflaton) is equal to the Planck mass M P while the mass of the late vacuumon (quintessence) is equal to the cosmon mass m Λ . Our model is able to describe the complete history of the universe from an early de Sitter era with density ρ P to a late de Sitter era with density ρ Λ , bridged by a phase of decelerated expansion. In our model, the periods of early and late inflation are described by two polytropic equations of state with index n e = +1 and n l = −1, respectively. This makes our model of universe very symmetric, the Planck density ρ P in the early universe playing a role similar to the cosmological density ρ Λ in the late universe. They represent fundamental upper and lower density bounds differing by 123 orders of magnitude. Our model is also fully consistent with the ΛCDM model at late times and completes it by incorporating in a natural manner a phase of early inflation that avoids the primordial (big bang) singularity. Therefore, in our model, the universe exists eternally in the past and in the future and does not present singularities. This has been called the aioniotic universe [42]. We have also made the connection between our model [39][40][41][42][43][44][45][46][47] based on a quadratic equation of state and the RVM [57][58][59] based on a quartic dependence of the cosmological constant on the Hubble parameter. In this connection, our exotic fluid can be viewed as a mixture of a barotropic fluid (representing stiff matter, radiation, or dark matter) and vacuum energy. Furthermore, our scalar field [42,44] corresponds to what Basilakos et al. [58] have called later the vacuumon. Despite its interest, our model cannot account for the spectrum of fluctuations in the early universe. This suggests that the vacuumon potential is just an effective classical potential that cannot be directly used to compute the fluctuations in the early universe. A fully quantum field theory may be required to achieve that goal.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Energy Conservation Equation for a Real Scalar Field
Starting from the KG Equation (30), we can derive the energy conservation Equation (32) for a real scalar field as follows. Taking the time derivative of the energy density from Equation (31), we obtain Combining Equation (A1) with Equation (30), we find that Summing the density and the pressure from Equation (31), we finḋ Substituting Equation (A3) into Equation (A2) we obtain the energy conservation Equation (32). Inversely, from Equations (31) and (32) we can derive the KG Equation (30).

Appendix B. Expression of the Energy Density of the Scalar Field Described by our Quadratic Equation of State
Let us consider a scalar field described by the quadratic equation of state (33). From this equation, we find The two roots of the term in brackets are As κ = ρ Λ /ρ P ∼ 10 −123 1, we can make the approximations These approximations are essentially exact as κ is extremely small. Therefore, Equation (A4) can be written in excellent approximation as Substituting this expression into the energy conservation Equation (32), we obtain where a 1 is a constant of integration. This equation can be rewritten as where we have used ρ P ρ Λ . We then find giving Using ρ P ρ Λ again, this relation can be rewritten as The constant of integration a 1 can be obtained by applying Equation (A12) at a = 1 where ρ = ρ 0 . This yields where we have used ρ P ρ Λ . We note that ρ 0 − ρ Λ = ρ α,0 can be interpreted as the present density of the α-fluid. Equation (A12) can finally be written as Remark A1. Actually, it is possible to integrate Equation (32) with Equation (A4) without making any approximation as detailed in Appendix C of [44]. However, in view of the extremely small value of κ = ρ Λ /ρ P ∼ 10 −123 1, this exact expression is not particularly useful for our purposes.

Appendix C. Connection between the RVM and our Quadratic Equation of State
In this appendix, we discuss the connection between the RVM [57][58][59] and our model [39][40][41][42][43][44][45][46][47]. The RVM assumes that the cosmological constant Λ(H) depends on the Hubble parameter and that this dependence is specified by a quartic function of H (actually, a quadratic function of H 2 ). This relationship is motivated by particle physics and the renormalization group approach. On the other hand, our phenomenological model assumes that the universe is filled with an exotic fluid described by the quadratic equation of state (33). We show that the two models are consistent and complementary to each other.

Appendix C.1. Approach Based on a Quartic Λ(H)
Let us first recall the main lines of the RVM [57][58][59]. For an N-species system, the conservation law reads where ρ N c 2 and P N denote the energy density and the pressure of each species, respectively. This equation expresses the conservation of the total energy (summed over all species present in the universe). If we consider a system composed of just one X-fluid + vacuum, we finḋ In the general case, the cosmological constant Λ(t) and, correspondingly, the vacuum energy density may depend on time. 26 The equation of state of vacuum is given by Therefore, Equation (A16) reduces tȯ We also recall the Friedmann equation where the right hand side includes the energy density of the X-fluid + the energy density of vacuum. We note that, according to Equation (A19), the X-fluid and the vacuum energy are coupled to each other and are not treated as independent species. Introducing the total density and the total pressure 27 , we can rewrite Equations (A19) and (A20) aṡ and We now assume that the X-fluid has a linear equation of state where w X is constant. Combining Equations (A22)-(A24), we obtaiṅ Using Equation (A21), this equation can be rewritten aṡ Finally, using Equation (A23), we finḋ and Equation (A31) can be solved by making the same approximations as in Appendix B 28 , yielding . (A32) Appendix C.2. Approach Based on a Quadratic P(ρ) Let us now recall the main lines of our approach [39][40][41][42][43][44][45][46][47]. We assume that the universe is filled with an exotic fluid characterized by a barotropic equation of state P(ρ). This fluid may correspond to a scalar field 29 , but this is not compulsory. The energy conservation equation is given byρ As H =ȧ/a, we can rewrite this equation as where ρ = dρ/da. We also recall the Friedmann equation Based on heuristic considerations [39][40][41][42][43][44][45][46][47], we assume that the pressure P(ρ) is a quadratic function of the energy density given by Substituting P(ρ) from Equation (A36) into Equation (A34), we obtain a differential equation for ρ (see below). The solution of this differential equation gives ρ(a). We can then obtain a(t) by solving the Friedmann Equation (A35). Specifically, for the quadratic equation of state P(ρ) from Equation (A36), the differential Equation (A34) for ρ reads which is solved in Appendix B to give This solution was first obtained in Refs. [39,40,42,44,47]. The complete analytical solution a(t) of the Friedmann Equation (A35) with Equation (A38) is given in [42,44]. This solution (see Equation (47)) describes a universe going from an early de Sitter era to a late de Sitter era passing by an intermediate phase of decelerated expansion.
3. Comparing Equations (A36) and (A42), we see that the coefficient of the linear term is given by In the matter era, we have w X = 0 implying α = −ν. Therefore, the RVM suggests that α = 0 in the matter era even though w X = 0. A nonvanishing value of α could be accounted for in our model by assuming that dark matter has an effective temperature so that P = αρc 2 with α = k B T eff /mc 2 . Taking v c ∼ (k B T/m) 1/2 ∼ 100 km/s from the rotation curves of the galaxies, we find α ∼ 10 −7 . 30 This estimate is consistent with the RVM provided that ν < 0 and |ν| ∼ 10 −7 . It is argued in [57,58] that ν can be positive or negative and that |ν| 1, probably in the range |ν| = 10 −6 − 10 −3 , and with a typical value |ν| = O(10 −3 ). 4. Except for the slight differences mentioned above, the RVM [57][58][59] and our model [39][40][41][42][43][44][45][46][47] are consistent and complementary to each other and both represent an interesting modelling of the evolution of the universe.
The scalar field potential corresponding to the equation of state (A45) is [42] where we have defined For ψ → 0, it can be expanded into For ψ → ±∞, we obtain the asymptotic behaviors We note that the coefficient α = 1 of stiff matter plays a special role as it leads to a faster decay of the potential. In that case, Equation (A47) becomes The relation between the scalar field and the scale factor is [42] (a/a * ) 3(α+1)/2n = sinh ψ.
This relation allows us to express ρ φ , P φ , etc. as a function of ψ instead of a.

Appendix D.3. Normal form of the Potential
The approximate expression (A49) of the scalar field potential for ψ → 0 can be compared with the normal form of a quartic potential where V 0 is the value of the potential at φ = 0, m is the mass of the scalar field, and λ is the dimensionless self-interaction constant. When the scalar field describes the wave function of a Bose-Einstein condensate (BEC), we have the relation 31 where a s is the scattering length of the bosons and λ C =h/(|m|c) is their Compton wavelength. One can also introduce the dimensional self-interaction constant Comparing Equation (A49) with Equation (A54), and recalling Equation (A48), we obtain the following results: is the cosmon radius (the semi Schwarzschild radius of a particle of mass m Λ ). Finally, Remark A3. We note that

Appendix E. Scalar Field in the Presence of One Fluid in the Early Universe
In the main text, we have assumed in approach (B) that the early universe is dominated by the scalar field and that the X-fluids appear later in the evolution of the universe. In that case, there is first a phase of inflation followed by an α-era both due to the scalar field. As the scalar field is alone in the early universe, its potential is given by Equation (115). The α-era may correspond to a stiff matter era (α = 1) or to a radiation era (α = 1/3). If the α-era corresponds to a stiff matter era (α = 1), we have to introduce an X-fluid with α X = 1/3 to represent the subsequent phase of radiation while this is not necessary if the α-era corresponds to a radiation era (α = 1/3). 33 Then, there is a matter era due to another X-fluid with α X = 0, and finally a late inflation due to the constant pressure of the scalar field.
In this appendix, we consider the possibility that an X-fluid coexists with the scalar field since the beginning of the universe. In that case, the X-fluid dominates the scalar field for a → 0 and leads to a big bang singularity where ρ X ∝ a −3(1+α X ) . This X-era is followed by an inflation era due to the scalar field when ρ X ρ P . Then, the evolution of the universe is dominated by the X-fluid or by the α-fluid (different possibilities can arise depending on the values of the parameters). 34 In any case, we exit this period through a radiation era due to the scalar field (if α = 1/3) or to the X-fluid (if α X = 1/3). Then, there is a matter era due to another X-fluid with α X = 0, and finally a late inflation era due to the constant negative pressure of the scalar field.
Substituting Equation (A85) into Equation (A83), we find that the scalar field potential is given in this regime by V(φ) = ρ P c 2 − 1 2 ρ P c 2 (α + 3) Ω P,0 Ω α,0 8πG 3c 2 1/2 3 2 α + α X + 2 √ 1 + α Ω α,0 Ω X,0 Ω P,0 φ 2(α+1) α+α X +2 . (A86) Notes 1 It is expected, but not observationally established, that the periods of acceleration are exponential, corresponding to a de Sitter stage. 2 The Planck density corresponds to a mass scale ∼ 10 19 GeV/c 2 . Actually, the scale of primordial inflation could be three orders of magnitude smaller, corresponding to the grand unification theory (GUT) scale ∼ 10 16 GeV/c 2 . In the following, for convenience, we shall identify the scale of primordial inflation to the Planck scale. If another scale turns out to be more relevant for our problem, we just have to replace the Planck density ρ P by the corresponding density in the equations. 3 This can be viewed as an "initial conditions problem" or as a "fine tuning problem". Indeed, as dark matter and dark energy evolve at different rates with the universe expansion, conditions in the early universe must be set very carefully in order for them to be comparable to the ones existing today. 4 We stress that there is no initial singularity in our model [45,46] as the stiff matter era starts after the inflation era (de Sitter) during which the density is constant (finite). 5 Note that Basilakos et al. [58] did not derive the complete energy density evolution nor the complete scalar field potential obtained in Equations (86), (106), and (121) of [44], but only their asymptotic expressions in the early and late universe. 6 Apparently, their formula (40) in [59] contains a mistake. The potential associated with the inflation + stiff matter era should read like Equation (F.42) of [46]. 7 This is because the other species enter into the Friedmann equation as additional components of the energy density and therefore alter the evolution of the scale factor and of the Hubble constant with respect to the free scalar field. 8 In this paper, we assume that the scalar field is real. The cosmological evolution of a complex scalar field is considered in [52,61,62].
On the other hand, in order to simplify the equations, we make the change of notation φ → cφ. 9 We will see in Sections 4 and 5 how to relate the potential V(φ) of the scalar field to its equation of state P φ (ρ). 10 As discussed in approach (A) of the Introduction, α may change with the density of the universe. Therefore, its value may depend on the epoch under consideration. 11 To avoid a spurious divergence of the energy density at a = 0, the matter component term Ω m,0 /a 3 has to be introduced at a sufficiently late time, i.e., after the inflation era when ρ ρ P . 12 We assume a non-phantom universe w φ > −1. We also assume that the scalar field φ increases with the scale factor a so that dφ/da ≥ 0. 13 We have left the lower limit of integration undetermined as the expression of the integrand is only valid for sufficiently large values of a. The lower limit of integration has to be obtained by matching the solutions in the early and late universe. 14 For example, we can take α = 1 (stiff matter) in the equation of state of the scalar field and add radiation and matter as additional species. 15 The case where a radiation era occurs before the inflation era, leading to a big-bang singularity, is considered in Appendix E. 16 This is a particular case of the general solution given in [42]. 17 For α = 1 (stiff matter), the scalar field potential is constant V = ρ Λ c 2 (see Equation (135)). 18 A similar quantization rule was introduced by Wesson [68]. By using the dimensional reduction from higher dimensional relativity and by assuming that the Compton wavelength of a particle cannot take any value, he proposed that the mass is quantized according to the rule m = (n * h /c 2 ) √ Λ/3 = n * m * Λ , where n * is an integer (this differs from Equation (175) in that it involves n * instead of √ n * ). Hence, m * Λ is the minimum mass corresponding to the ground state n * = 1. In our model, the mass of the scalar field associated with the ΛCDM model (α = 0) is (3/2)m * Λ . 19 Some analytical solutions of the Friedmann equation involving two or more fluids with a linear equation of state (e.g., stiff matter, radiation, matter, or dark energy) are given in [45]. 20 These equations can be obtained from Equations (86) and (87) by using the fact that both the scalar field and the X-fluids have a linear equation of state. They can also be recovered from Equations (94) and (95) by taking Ω Λ,0 = 0 and Ω P,0 → +∞. 21 The two potentials (149) and (199) are equivalent provided that the terms depending on X in Equation (199) are included in the constant of integration φ * appearing in Equation (149). 22 In the framework of our model, this corresponds to the situation where the scalar field describes dark radiation (α = 1/3) and the X-fluid describes normal radiation (α X = 1/3). This also corresponds to the situation where the scalar field describes dark matter (α = 0) and the X-fluid describes baryonic matter (α X = 0). 23 This assumes ρ Λ = 0. The case ρ Λ = 0 leads to Equations (190)-(193). 24 This is basically why we do not find a good agreement with the observations. 25 The function H(φ) can be obtained from Equations (104), (123), and (213). However, in line with the Hamilton-Jacobi formalism, we shall proceed here the other way round and take Equation (229) as a starting point from which all the results of our model can be derived. 26 In this section, the time-dependent vacuum energy density ρ Λ (t) should not be confused with the constant cosmological density ρ Λ = 5.96 × 10 −24 g m −3 appearing in Equation (33). 27 In the scalar field representation of the RVM [58,59], the scalar field is associated with the total density and pressure. This implies that the X-fluid is part of the scalar field. In other words, the vacuumon is not associated to the vacuum alone, but to the vacuum + X-fluid. We also recall that the X-fluid changes with the epoch considered (early or late universe). 28 We show in Appendix C.3 that the calculations in the RVM are equivalent to those of Appendix B. 29 In that case, ρ and P represent the energy density and the pressure of the scalar field denoted ρ φ and P φ in the main text. 30 This value is consistent with the condition α ≤ 10 −7 necessary to avoid the presence of oscillations in the matter power spectrum (see Appendix C in [52]). 31 The factor 2/3 arises because we consider a real scalar field (see [83] for more details). In principle, Equation (A55) makes sense only if m is a positive real number. However, in order to treat all possible situations, we formally extend this formula to the case where m is imaginary (m 2 < 0) by taking its modulus |m|. 32 This mass scale is often interpreted as the smallest mass of the elementary particles predicted by string theory [84] or as the upper bound on the mass of the graviton [85]. The mass m Λ also represents the quantum of mass in theories of extended supergravity [66]. The mass scale m Λ is simply obtained by equating the Compton wavelength of the particle λ C =h/mc with the Hubble radius R Λ = c/H 0 (the typical size of the visible universe) giving m Λ =hH 0 /c 2 ∼h √ Λ/c 2 (as H 2 0 ∼ Gρ Λ ∼ Λ). The mass m Λ corresponds to Wesson's [68] minimum mass interpreted as a quantum of dark energy (Wesson's maximum mass M Λ = (4/3)πρ 0 R 3 Λ = c 3 /2GH 0 = 9.20 × 10 55 g is of the order of the mass of the universe). The mass scales M Λ and m Λ also appear in Refs. [51,86] and represent the mass of the visible universe and the minimum mass of the bosons. Böhmer and Harko [87] proposed to call the elementary particle of dark energy having the mass m Λ the "cosmon". Cosmons were originally introduced by Peccei et al. [88] to name scalar fields that could dynamically adjust the cosmological constant to zero (see also [89][90][91]). The name cosmon was also used in a different context [92] to designate a very light scalar particle (dilaton) of mass ∼ 10 −3 eV/c 2 which could mediate new macroscopic forces in the submillimeter range. 33 We could, however, introduce an X-fluid with α X = 1/3 to distinguish the radiation due to the scalar field (dark radiation) from the ordinary radiation due to photons or other relativistic particles.