Dynamically Generated Inflationary Lambda-CDM

Our primary objective is to construct a plausible unified model of inflation, dark energy and dark matter from a fundamental Lagrangian action first principle, where all fundamental ingredients are systematically dynamically generated starting from a very simple model of modified gravity interacting with a single scalar field employing the formalism of non-Riemannian spacetime volume-elements. The non-Riemannian volume element in the initial scalar field action leads to a hidden nonlinear Noether symmetry which produces energy-momentum tensor identified as a sum of a dynamically generated cosmological constant and a dust-like dark matter. The non-Riemannian volume-element in the initial Einstein-Hilbert action upon passage to the physical Einstein-frame creates dynamically a second scalar field with a non-trivial inflationary potential and with an additional interaction with the dynamically generated dark matter. The resulting Einstein-frame action describes a fully dynamically generated inflationary model coupled to dark matter. Numerical results for observables such as the scalar power spectral index and the tensor-to-scalar ratio conform to the latest 2018 PLANCK data.


Introduction
In the last decade or so a groundbreaking concept emerged about the intrinsic necessity to modify (extend) gravity theories beyond the framework of standard Einstein's general relativity. The main motivation for these developments is to overcome the limitations of the latter coming from: (i) Cosmology -for solving the problems of dark energy and dark matter and explaining the large scale structure of the Universe [1][2][3]; (ii) Quantum field theory in curved spacetime -because of the non-renormalizabilty of ultraviolet divergences in higher loops [4][5][6][7][8][9]; (iii) Modern string theory -because of the natural appearance of higher-order curvature invariants and scalar-tensor couplings in low-energy effective field theories [10][11][12][13][14].
Another parallel crucial development is the emergence of the theoretical framework based on the concept of "inflation", which is a necessary part of the standard model of cosmology, since it provides the solution to the fundamental puzzles of the old Big Bang theory, such as the horizon, the flatness, and the monopole problems [15][16][17][18][19][20][21][22]. It can be achieved through various mechanisms, for instance through the introduction of primordial scalar field(s) , or through correction terms into the modified gravitational action .
Additionally, inflation was proved crucial in providing a framework for the generation of primordial density perturbations [123,124]. Since these perturbations affect the Cosmic Background Radiation (CMB), the inflationary effect on observations can be investigated through the prediction for the scalar spectral index of the curvature perturbations and its running, for the tensor spectral index, and for the tensor-to-scalar ratio.
Various classes of modified gravity theories have been employed to construct viable inflationary models: f (R)-gravity; scalar-tensor gravity; Gauss-Bonnet gravity (see [125,126] for a detailed accounts); recently also based on non-local gravity ( [127] and references therein) or based on brane-world scenarios ( [128] and references therein). The first early successful cosmological model based on the extended f (R) = R + R 2 -gravity produces the classical Starobinsky inflationary scalar field potential [16].
Dynamically generated models of inflation from modified/extended gravity such as the Starobinsky model [18,126,129,130] still remain viable and produce some of the best fits to existing observational data compared to other inflationary models [131].
Unification of inflation with dark energy and dark matter have been widely discussed [79,81,[132][133][134][135][136][137][138][139][140][141][142]. It is indeed challenging to describe both phases of acceleration using a single scalar field minimally couple to gravity, without affecting the thermal history of the universe which has been verified to a good accuracy. In order to enable slow-roll behavior, the scalar field potential should exhibit shallow behaviour at early times followed by a steep region for most of the universe history turning shallow once again at late times. Although a simple exponential potential does not comply with the above picture, here we present a simple modified gravity model naturally providing a dynamically generated scalar potential, whose inflationary dynamics is compatible with the recent observational data. On the other hand, the task of describing particle creation will be discussed in our future work.
Another specific broad class of modified (extended) gravitational theories is based on the formalism of non-Riemannian spacetime volume-elements. It was originally proposed in [143][144][145][146][147], with a subsequent concise geometric formulation in [148][149][150]. This formalism was used as a basis for constructing a series of extended gravity-matter models describing unified dark energy and dark matter scenario [151,152], quintessential cosmological models with gravity-assisted and inflaton-assisted dynamical suppression (in the "early" universe) or dynamical generation (in the post-inflationary universe) of electroweak spontaneous symmetry breaking and charge confinement [153,154], as well as a novel mechanism for dynamical supersymmetric Brout-Englert-Higgs effect in supergravity [148].
In the present paper our principal aim is to construct a plausible unified model, i.e., describing (most of the) principal physical manifestations of a unification of inflation and dark energy interacting with dark matter, where the formalism of the non-Riemannian spacetime volume-elements will play a fundamental role. To this end we will consider a simple modified gravity interacting with a single scalar field where the Einstein-Hilbert part and the scalar field part of the action are constructed within the formalism of the non-Riemannian volume-elements -alternatives to the canonical Riemannian one √ −g. The non-Riemannian volume element in the initial scalar field action leads to a hidden nonlinear Noether symmetry which produces energy-momentum tensor identified as a sum of a dynamically generated cosmological constant and a dynamically generated dust-like dark matter. The non-Riemannian volume-element in the initial Einstein-Hilbert action upon passage to the physical Einstein-frame creates dynamically a second scalar field with a non-trivial inflationary potential and with an additional interaction with the dynamically generated dark matter. The resulting Einstein-frame action describes a fully dynamically generated unified model of inflation, dark energy and dark matter. Numerical results for observables such as the scalar power spectral index and the tensor-to-scalar ratio conform to the latest 2018 PLANCK data.
Let us briefly recall the essence of the non-Riemannian volume-form (volume-element) formalism. In integrals over differentiable manifolds (not necessarily Riemannian one, so no metric is needed) volume-forms are given by nonsingular maximal rank differential forms ω: where Our conventions for the alternating symbols ε µ 1 ,...,µ D and ε µ 1 ,...,µ D are: ε 01...D−1 = 1 and ε 01...D−1 = −1). The volume element Ω transforms as scalar density under general coordinate reparametrizations. In Riemannian D-dimensional spacetime manifolds a standard generally-covariant volume-form is defined through the "D-bein" (frame-bundle) canonical one-forms e A = e A µ dx µ (A = 0, . . . , D − 1): yields: To construct modified gravitational theories as alternatives to ordinary standard theories in Einstein's general relativity, instead of √ −g we can employ one or more alternative non-Riemannian volume element(s) as in (1) given by non-singular exact D-forms ω = dA where: so that the non-Riemannian volume element reads: Thus, a non-Riemannian volume element is defined in terms of the (scalar density of the) dual field-strength of an auxiliary rank D − 1 tensor gauge field A µ 1 ...µ D−1 . The modified gravity Lagrangain actions based on the non-Riemannian volume-elements formalism are of the following generic form (here and in what follows we will use units with 16πG Newton = 1): where Φ 0 (A) and Φ 1 (B) are of the form (6) (for D = 4), R is the scalar curvature, and the Lagrangian densities L 0,1,2 contain the matter fields (and possibly higher curvature terms, e.g. R 2 ). A basic property of the class of actions (7) is that the equations of motion w.r.t. auxiliary gauge fields, defining the non-Riemannian volume-elements Φ 0 (A) and Φ 1 (B) as in (6), produce dynamically generated free integration constants M 1 , M 0 : (cf. Eqs. (15) and (28) below) whose appearance will play an instrumental role in the sequel.
Further, let us stress on the following important characteric feature of the modified gravity-matter actions (7). When considering the gravity part in the first order (Palatini) framework (i.e., R = g µν R µν (Γ) with a priori indepedent metric g µν and affine connection Γ λ µν ), then the auxiliary rank 3 tensor gauge fields defining the non-Riemannian volume-elements in (7) are almost pure-gauge degrees of freedom, i.e. they do not introduce any additional propagating gravitational degrees of freedom when passing to the physical Einstein-frame except for few discrete degrees of freedom with conserved canonical momenta appearing as arbitrary integration constants. This has been explicitly shown within the canonical Hamiltonian treatment [149,153].
On the other hand, when we treat (7) in the second order (metric) formalism (the affine connection Γ λ µν is the canonical Levi-Civitta connection in terms of g µν ), while passing to the physical Einstein-frame via conformal transformation (cf. Eq.(30) below): the first non-Riemannian volume element Φ 1 (A) in (7) is not any more (almost) "pure gauge", but creates a new dynamical canonical scalar field u via , which will play the role both of an inflaton field at early times, as well as driving late-time de Sitter expansion (see Section 3 below).
In Section 2 we briefly review our construction in [152] of a simple gravity-scalar-field model -specific member of the class of modified gravitational models (7) of the form (10) below, which yields an explicit dynamical generation of independent (non-interacting among each other) dark energy and dark matter components in an unified description as a manifestation of a single material entity ("darkon" scalar field)a simplest realization of a ΛCDM model.
In section 3 we extend the previous construction to dynamically generate, apart from dark matter, also early-time inflation and late-time de Sitter expansion -via dynamical creation of an additional canonical scalar field u ("inflaton") out of a non-Riemannian volume-element with the following properties: (i) u aquires dynamically a non-trivial inflationary type scalar field potential driving inflation at early times of the universe' evolution; (ii) At late times the same evolving u flows towards a stable critical point of the pertinent dynamical system describing the cosmological evolution, driving a late-time de Sitter expansion in a dark energy dominated epoch; (iii) In this case the field u, induces a specific interaction between the dark energy and dark matter.
In Section 4 we study the cosmological implications of the latter dynamically generated inflationary model with interacting dark energy and dark matter. In Section 5 several plots of the numerical solutions for the evolution of the dynamical inflationary field and for the behavior of the relevant inflationary slow-roll parameters and the corresponding observables are presented. Section 6 contains our conclusions and outlook.

A Simple Model of Unification of Dark Energy and Dark Matter
In [152] we started with the following non-conventional gravity-scalar-field action -a simple particular case of the class (7) -containing one metric-independent non-Riemannian volume-element alongside with the standard Riemannian one: with the following notations: • The first term in (10) is the standard Einstein-Hilbert action with R(g) denoting the scalar curvature w.r.t. metric g µν in the second order (metric) formalism; • Φ 0 (A) is particular representative of a D = 4 non-Riemannian volume-element density (6): • L(ϕ, X) is general-coordinate invariant Lagrangian of a single scalar field ϕ(x): Varying (10) w.r.t. g µν , ϕ and A µνλ yield the following equations of motion, respectively: where M 0 is arbitrary integration constant (the factor 2 is for later convenience). As stressed in [152], the scalar field dynamics is determined entirely by the first-order differential equation -the dynamical constraint Eq. (15). The usual second order differential equation (14) for ϕ is in fact a consequence of (15) together with the energy-momentum conservation: Also, as exhibited in [152], the specific form of the scalar field potential V(ϕ) does not affect the dynamics of the system (10), see the remark below following (18). The same phenomenon occurs in the extension of (10) to the model (24) in Section 3 and 4 below.
The canonical Hamiltonian analysis in [152] of the action (10) reveals that the auxiliary gauge field A µνλ is in fact an almost pure-gauge, i.e., it is a non-propagating field-theoretic degree of freedom with the integration constant (−2M 0 ) identified with the conserved Dirac-constrained canonical momentum conjugated to the "pure gauge" "magnetic" component of A µνλ . For a general canonical Hamiltonian treatment of Lagrangian action with one or more non-Riemannian volume-elements, we refer to [155].
A crucial property of the model (10) is the existence of a hidden nonlinear Noether symmetry revealed in [152]. Indeed, both Eqs.(14)- (15) can be equivalently rewritten in the following current-conservation law form: The covariantly conserved current J µ (17) is the Noether current corresponding to the invariance (modulo total derivative) of the action (10) w.r.t following hidden nonlinear symmetry transformations: dual" components of the auxiliary gauge field A µνλ (11).
Remark. We notice that the existence of the hidden nonlinear symmetry (18) of the action (10) does not depend on the specific form of the scalar field potential V(ϕ).
The next important step is to rewrite T µν (13) and J µ (17) in the relativistic hydrodynamical form (again taking into account (15)): Here the integration constant M 0 appears as dynamically generated cosmological constant and: We now find that the covariant conservation laws for the energy-momentum tensor (19) ∇ µ T µν = 0 and the J-current (17) acquire the form: Eqs. (21) imply in turn the geodesic equation for the "fluid" 4-velocity u µ : Therefore, comparing (19) with the standard expression for a perfect fluid stress-energy tensor T µν = ρ + p)u µ u ν + pg µν , we see that T µν (19) consists of two additive parts which have the following interpretation according to the standard Λ-CDM model [156][157][158][159][160][161][162][163] (using notations p = p DM + p DE and ρ = ρ DM + ρ DE ): • Dynamically generated dark energy part given by the second cosmological constant term in T µν (19) due to (15), where p DE = −2M , ρ DE = 2M; • Dynamically generated dark matter part given by the first term in (19), where p DM = 0 , ρ DM = ρ 0 with ρ 0 as in (20), which in fact according to (21) and (22) describes a dust fluid with fluid density ρ 0 flowing along geodesics. Thus, we will call the ϕ scalar field by the alias "darkon".
The conservation laws (21) due to the hidden nonlinear Noether symmetry (18) imply that in the model (10) there is no interaction between dark energy and dark matter -they are separately conserved.

Inflation and Unified Dark Energy and Dark Matter
Now we will extend the simple model (10) of unified dark energy and dark matter by introducing another metric-independent non-Riemannian volume-element: inside the gravity (Einstein-Hilbert) part of the action (using again units with 16πG Newton = 1): Here Λ 0 is a dimensionful parameter to be identified later on as energy scale of the inflationary universe' epoch.
The specific form of the action (24) may be justified by the requirement about global Weyl-scale invariance under the transformations: and provided we choose V(ϕ) = ϕ 4 . Concerning global Weyl-scale invariance let us note that it played an important role already since the first original papers on the non-canonical volume-form formalism [146].
In particular, models with spontaneously broken dilatation symmetry have been constructed along these lines, which are free of the Fifth Force Problem [147].
The equations of motion of the action (24) w.r.t. ϕ and A µνλ are the same as in (14)- (15), therefore once again (24) is invariant under the hidden nonlinear Noether symmetry (18) with the associated Noether conserved current (17), which we rewrite here for later convenience taking into account (15): On the other hand, the equations of motion w.r.t. g µν and B µνλ now read: where T µν is the same energy-momentum tensor as in (13) or (19), which taking into account (15) and using short-hand notation χ 0 in (26)) reads T µν = −2M 0 g µν + 1 + χ 0 ∂ µ ϕ∂ ν ϕ , and M 1 is another free integration constant similar to M 0 in (15). Taking trace of (27) together with (28) imply a dynamical equation for χ 1 (χ 0 and χ 1 as defined in (26) and (28), respectively): The passage to the Einstein-frame is accomplished via the conformal transformation: on Eqs. (27) and (29), and upon using the known formulae for conformal transformations of Ricci curvature tensor and covariant Dalambertian (see e.g. [164]; bars indicate magnitudes in theḡ µν -frame): In the process we introduce the field redefinition χ 1 → u: so that u appears as a canonical scalar field in the Einstein-frame transformed equations (27), (29) and (15): and most importantly u acquires a non-trivial dynamically denerated potential: due to the appearance of the free integration constants from the equations of motion of the original-frame non-Riemannian spacetime volume-elements. The hidden nonlinear Noether symmetry current conservation (17), equivalent to the ϕ-equation of motion, becomes in the Einstein-frame: Thus, the Einstein-frame Lagrangian action producing the Einstein-frame equations of motion (34)-(38) reads: with U eff (u) as in (37) and where now χ 0 (from (26)) becomes a simple Lagrange multiplier. The upper line in S EF (39) represents an inflationary Lagrangian action with dynamically generated inflationary potential U eff (u) (37) obtained in [107] from a pure gravity initial action (without any matter fields) in terms of non-Riemannian volume-elements: which is graphically depicted on Fig.1, is a generalization of the classic Starobinsky inflationary potential [16]. In fact, the latter is a special case of (37) for the particular values of the parameters: Λ 0 = M 0 = 1 4 M 1 .
Inflation -Slow roll U eff (u) (37) possesses two main features relevant for cosmological applications: • (i) U eff (u) (37) has an almost flat region for large positive u: U eff (u) 2Λ 0 for large u. This almost flat region correspond to "early" universe' inflationary evolution with energy scale 2Λ 0 as it will be evident from the autonomous dynamical system analysis of the cosmological dynamics in Section 4; • (ii) U eff (u) (37) has has a stable minimum for a small finite value u = u * : ∂U eff ∂u = 0 for u ≡ u * , where: • (iii) As it will be explicitly exhibited in the dynamical system analysis in Section 4, the region of u around the stable minimum at u = u * (41) correspond to late-time de Sitter expansion of the universe with slightly varied late-time Hubble parameter (dark energy dominated epoch), where the minimum value of the potential: is the asymptotic value at t → ∞ of the dynamical dark energy density [165,166].
The lower line in S EF (39) represents the interaction between the dynamical inflaton field u and the "darkon" field ϕ, in other words here we have unification of inflation, dark energy and dark-matter. This is reflected in the structure of the Einstein-frame energy-momentum tensorT µν (34) -the first two terms being the stress-energy tensor of u and the last term being the "darkon" stress-energy tensor coupled to u.

Cosmological Implications
Let us now consider reduction of the Einstein-frame action (39) to the Friedmann-Lemaitre-Robertson-Walker (FLRW) framework with metric ds 2 = −N 2 dt 2 + a(t) 2 d x 2 , where u = u(t) and ϕ = ϕ(t): The equations of motion w.r.t. χ 0 and ϕ from (43) are equivalent to the FLRW reduction of the dynamical constraint (36) and of the Noether current conservation (38), respectively: .
which imply the relation: with c 0 a free integration constant. Taking into account (45), the FLRW reduction of the Einstein-frame energy-momentum tensor (34) becomes:T Relations (47) explicitly show that the last term in ρ: represents the "dust" dark matter part of the total energy denisty -it is "dust" because of absence ot corresponding contribution for the pressure p in (47). The equation of motion from (43) a /a being the Hubble parameter): ..
and, finally, the two Friedmann equations (varying (43) w.r.t lapse N and a) read: .
Remark. We observe that due to the hidden nonlinear Noether symmetry current conservation (45), the FLRW dynamics given by (49)-(51) does not depend on the explicit form the "darkon" part of the FLRW action (43) -the only trace of the "darkon" is embodied in the integration constant c 0 .
It is instructive to analyze the system FLRW equations (49)-(51) as an automous dynamical system. To this end it is useful to rewrite the system (49)-(51) in terms of a set of dimensionless coordinates (following the approach in [167]): with L DE as in (42) and ρ DM as in (48). In these coordinates the system defines a closed orbit: which is equivalent to the first Friedmann equation (50). Then Eqs. (49) and (51) can be represented as a 3-dimensional autonomous dynamical system for the (x, y, H variables (cf. (52)): where the primes indicate derivatives w.r.t. number of e-folds N = log(a) (meaning d dN = 1 H d dt ). The dynamical system (54)-(56) possesses two critical points: • (A) Stable critical point: where all three eigenvalues of the stability matrix are negative or with negative real parts (λ 1 = . The stable critical point (57) corresponds to the late-time asymptotics of the universe' evolution where according to the definitions (52) u(t) → u * -the stable minimum of the effective potential U eff (u) (37), so that U eff (u) → 2Λ DE , the dark matter energy density (48) ρ DM → 0, and .
H→ 0 accordin to (56), i.e., late-time accelerated expansion with where one of the three eigenvalues of the stability matrix is zero ( 16M 0 Λ 0 ). According to the definitions (52), in the vicinity of the unstable critical point (58) u(t) is very large positive (u → ∞), so that U eff (u) 2Λ 0 , ρ DM is vanishing ρ DM ≈ 0, and we have there a slow-roll inflationary evolution with inflationary scale Λ 0 where the standard slow-roll parameters are very small:

Numerical Solutions
Going back to the system of equations (49)-(51) we can use (50) to replace the term ρ DM ≡ 2 c 0 a 3 e −u/2 √ 3 in (49) and (51) so that we will obtain a closed system of two coupled nonlinear diffential equations for u(t), H(t) of second and first order, respectively: .. .
where U eff (u) is given by (37): Below we present several plots qualitatively illustrating the evolutionary behavior of the numerical solution of the system (61)  To obtain plausible values for the observables -the scalar power spectral index n s and the tensor to scalar ratio r [53,105,171]. we need the functional dependence of the slow-roll parameters and η w.r.t. N = log(a) -number of e-folds. More specifically, in N f is the number of e-folds at the end of inflation defined as (N f ) ≈ 1, then we need the values (N i ) and η(N i ) at N i -e-folds at the start of inflation, where it is assumed that N f − N i ∼ 60. Then, according to [53,105]: where the corresponding slow roll parameter read: and where H = H(N ) is the functional dependence of Hubble parameter w.r.t. the number of e-folds. To this end we employ numerical simulation of the autonomous dynamical system equations (54)- (56).
From the inflationary scenario we know that the observed value of the inflationary scale Λ 0 ∼ 10 −8 M 4 Pl is way larger than the current value (∼ 10 −122 M 4 Pl ) of the cosmological constant Λ DE (42). So, as in the numerical example above for the numerical solution of the system for u(t), H(t) (61)-(62), we will take again the values for the parameters according to (63) meaning that we set the initial condition for the Hubble parameter to be according to (58) With those numerical values we obtain for the observables (64) to be: which are well inside the last PLANCK observed constraints [131]: 0.95 < n s < 0.97 , r < 0.064 . In order to see the pattern of the general behavior depending on the initial conditions, we employ here Monte Carlo simulation with 10 4 samples for the initial conditions using a normal distribution: Λ 0 = 50 ± 10, M 1 = 20 ± 10, while the error bar is taken to be 1 σ. Fig.7 shows how different values of initial conditions yield different number of e-folds until end of inflation (where = 1) and, accordingly, different values for the observables r and n s , whereas Fig.8 depicts the corresponding relation between r and n s . Nevertheless, all the values of the latter fall within the constraint (67).  Figure 7. The scalar to tensor ratio r and the scalar spectral index n s vs. the number of e-folds for different values of the initial conditions. The sampling of the latter is done with a normal distribution Λ 0 = 50 ± 10, M 1 = 20 ± 10.

Conclusions and Outlook
In the present paper, starting from the basic first principle of Lagrangian field-theoretic actions combined with a non-canonical modification of gravity via employing non-Riemannian spacetime volume forms as alternatives to the standard Riemannian one given by √ −g, we have constructed a unified model of dynamically generated inflation with dark energy and dark matter coupled among themselves. Upon In the present model dark matter in the form of a dust-like fluid is created already in the early universe inflationary epoch without a significant impact on the inflationary dynamics. After end of inflation the dust-like dark matter apart from a short period of matter domination still does not exert a sufficient impact, which means that one has to further extend the present formulation in order to take properly into account the full dark matter contribution to the evolution.
One subject that has to be addressed is the "reheating of the universe", since of course we need temperature in the early universe to account for processes like Bing Bang Nucleosynthesis. There are many way to achieve this, due to the oscillating nature of inflaton solutions near the minimum of the inflaton potential, which leads in general to particle creation. For example, one possible way to complement the modified gravity-scalar field model (24) in order to incorporate the effect of radiation after end of inflation is to include a coupling to the "topological" density of a electromagnetic field A µ with field strength F µν = ∂ µ A ν − ∂ ν A µ in the following way: Upon passage to the Einstein-frame via the conformal transformation (30) the action (68) becomes (cf. (39)): The coupling term e −u/ √ 3 ε µνκλ F µν F κλ is suppressed in the inflationary stage where the derivative of u is small (because of the slow roll regime), whereas after end of inflation it may produce pairs of photons out of u due to the appreciable time-derivative of u resulting from the oscillations near the minimum of the effective potential. Of course, many other possible interaction terms can be introduced. Finally, in the reheating stage many particles can be produced, some of them could be no standard-model particles. If those are stable, they could provide additional "dark matter" apart from the "darkon" dust-like dark matter discussed here. Of course, if all created particles beyond those of the standard models turn out to be unstable, then we will be left with the "darkon" as the unique source of dark matter.

Acknowledgments
We gratefully acknowledge support of our collaboration through the Exchange Agreement between Ben-Gurion University, Beer-Sheva, Israel and Bulgarian Academy of Sciences, Sofia, Bulgaria. E.N. and S.P. are thankful for support by Contract DN 18/1 from Bulgarian National Science Fund. D.B., E.G. and E.N. are also partially supported by COST Actions CA15117, CA16104 and CA18108.