Dynamically Generated Inflation from Non-Riemannian Volume Forms

TL;DR

This paper introduces a modified gravity model using non-Riemannian volume forms that naturally generates an inflationary potential, leading to cosmological predictions consistent with current observations.

Contribution

It presents a novel gravity model with non-Riemannian volume elements that produce inflationary dynamics without initial matter fields.

Findings

Predicts spectral index n_s ≈ 0.96

Predicts tensor-to-scalar ratio r ≈ 0.002

Compatible with observational data for 60 e-folds

Abstract

We propose a simple modified gravity model without any initial matter fields in terms of several alternative non-Riemannian spacetime volume elements within the metric (second order) formalism. We show how the non-Riemannian volume-elements, when passing to the physical Einstein frame, create a canonical scalar field and produce dynamically a non-trivial inflationary-type potential for the latter with a large flat region and a stable low-lying minimum. We study the evolution of the cosmological solutions from the point of view of theory of dynamical systems. The theory predicts the spectral index $n_s \approx 0.96$ and the tensor-to-scalar ratio $r \approx 0.002$ for 60 e-folds, which is in accordance with the observational data. In the future Euclid and SPHEREx missions or the BICEP3 experiment are expected to provide experimental evidence to test those predictions.

Figures

Equations

$$\sqrt{-g}\equiv \sqrt{-\det \Vert g_{\mu \nu }\Vert }$$

$$\displaystyle\int_{{\mathcal{M}}}\omega\bigl{(}\ldots\bigr{)}=\int_{{\mathcal{M}}}dx^{D}\,\Omega\bigl{(}\ldots\bigr{)}\;\;,$$

$$\omega =\frac{1}{D!}{\omega }_{{\mu }_1\dots {\mu }_D}d{x}^{{\mu }_1}\wedge \dots \wedge d{x}^{{\mu }_D},$$

$${\omega }_{{\mu }_1\dots {\mu }_D}=-{\epsilon }_{{\mu }_1\dots {\mu }_D}\Omega ,$$

$$\omega =e^0\wedge \dots \wedge e^{D-1}=\det \Vert e_{\mu }^A\Vert d{x}^{{\mu }_1}\wedge \dots \wedge d{x}^{{\mu }_D}$$

$$⟶\Omega =\det \Vert e_{\mu }^A\Vert =\sqrt{-\det \Vert g_{\mu \nu }\Vert }.$$

$$A=\frac{1}{(D-1)!}{A}_{{\mu }_1\dots {\mu }_{D-1}}d{x}^{{\mu }_1}\wedge \dots \wedge d{x}^{{\mu }_{-1}}$$

$$⟶\Omega \equiv \Phi (A)=\frac{1}{(D-1)!}{\epsilon }^{{\mu }_1\dots {\mu }_D}{\partial }_{{\mu }_1}{A}_{{\mu }_2\dots {\mu }_D}.$$

$$S=\int d^{4}x\Bigl{\{}\Phi_{1}(A)\Bigl{[}R+L^{(1)}\Bigr{]}+\Phi_{2}(B)\Bigl{[}L^{(2)}+\frac{\Phi_{0}(C)}{\sqrt{-g}}\Bigr{]}+\ldots\Bigr{\}}\;.$$

$$\int d^{4}x\,\sqrt{-g}\Bigl{[}R+\ldots\Bigr{]}$$

$$g_{\mu \nu }\to {\bar{g}}_{\mu \nu }={\chi }_1{g}_{\mu \nu },{\chi }_1\equiv \frac{{\Phi }_1(A)}{\sqrt{-g}},$$

$$S=\int d^{4}x\Bigl{\{}\Phi_{1}(A)\Bigl{[}R-2\Lambda_{0}\frac{\Phi_{1}(A)}{\sqrt{-g}}\Bigr{]}+\Phi_{2}(B)\frac{\Phi_{0}(C)}{\sqrt{-g}}\Bigr{\}}\;,$$

$${\Phi }_1(A)\equiv \frac{1}{3!}{\epsilon }^{\mu \nu \kappa \lambda }{\partial }_{\mu }{A}_{\nu \kappa \lambda },{\Phi }_2(B)\equiv \frac{1}{3!}{\epsilon }^{\mu \nu \kappa \lambda }{\partial }_{\mu }{B}_{\nu \kappa \lambda },$$

$${\Phi }_0(C)\equiv \frac{1}{3!}{\epsilon }^{\mu \nu \kappa \lambda }{\partial }_{\mu }{C}_{\nu \kappa \lambda }.$$

$$g_{\mu \nu }\to \lambda g_{\mu \nu },$$

$$A_{\mu \nu \kappa }\to \lambda A_{\mu \nu \kappa },B_{\mu \nu \kappa }\to {\lambda }^2{B}_{\mu \nu \kappa },C_{\mu \nu \kappa }\to C_{\mu \nu \kappa }.$$

$$R-4{\Lambda }_0\frac{{\Phi }_1(A)}{\sqrt{-g}}=-M_1\equiv \mathrm{const},$$

$$\frac{{\Phi }_0(C)}{\sqrt{-g}}=-M_2\equiv \mathrm{const},\frac{{\Phi }_2(B)}{\sqrt{-g}}={\chi }_2\equiv \mathrm{const}.$$

$$R_{\mu\nu}-\Lambda_{0}\chi_{1}\,g_{\mu\nu}+\frac{1}{\chi_{1}}\bigl{(}g_{\mu\nu}\Box{\chi_{1}}-\nabla_{\mu}\nabla_{\nu}\chi_{1}\Bigr{)}-\frac{\chi_{2}M_{2}}{\chi_{1}}g_{\mu\nu}=0\;,$$

$$3\frac{\square {\chi }_1}{{\chi }_1}-\frac{4{\chi }_2{M}_2}{{\chi }_1}-M_1=0.$$

$$R_{\mu \nu }(g)=R_{\mu \nu }(\bar{g})-3\frac{{\bar{g}}_{\mu \nu }}{{\chi }_1}{\bar{g}}^{\kappa \lambda }{\partial }_{\kappa }{\chi }_1^{1/2}{\partial }_{\lambda }{\chi }_1^{1/2}$$

$$\quad+\chi_{1}^{-1/2}\bigl{(}{\bar{\nabla}}_{\mu}{\bar{\nabla}}_{\nu}\chi_{1}^{1/2}+{\bar{g}}_{\mu\nu}{\bar{\Box}}\chi_{1}^{1/2}\bigr{)}\;,$$

$$\Box\chi_{1}=\chi_{1}\Bigl{(}{\bar{\Box}}\chi_{1}-2{\bar{g}}^{\mu\nu}\frac{\partial_{\mu}\chi_{1}^{1/2}\partial_{\nu}\chi_{1}}{\chi_{1}^{1/2}}\Bigr{)}\;,$$

$$R_{\mu \nu }(\bar{g})-\frac{1}{2}{\bar{g}}_{\mu \nu }R(\bar{g})=$$

$$\displaystyle\frac{1}{2}\Bigl{[}\partial_{\mu}u\partial_{\nu}u-{\bar{g}}_{\mu\nu}\bigl{(}\frac{1}{2}{\bar{g}}^{\kappa\lambda}\partial_{\kappa}u\partial_{\lambda}u+U_{\rm eff}(u)\bigr{)}\Bigr{]}\;,$$

$$\Phi_{1}(A)/\sqrt{-g}\equiv\chi_{1}=\exp{\bigl{(}u/\sqrt{3}\bigr{)}}$$

$$U_{\rm eff}(u)=2\Lambda_{0}-M_{1}\exp{\bigl{(}-\frac{u}{\sqrt{3}}\bigr{)}}+\chi_{2}M_{2}\exp{\bigl{(}-2\frac{u}{\sqrt{3}}\bigr{)}}\;.$$

$$\bar{\square }u+\frac{\partial U_{\mathrm{eff}}}{\partial u}=0$$

$$S_{\rm EF}=\int d^{4}x\sqrt{-{\bar{g}}}\Bigl{[}R({\bar{g}})-\frac{1}{2}{\bar{g}}^{\mu\nu}\partial_{\mu}u\partial_{\nu}u-U_{\rm eff}(u)\Bigr{]}\;.$$

$$\exp\bigl{(}-\frac{u_{*}}{\sqrt{3}}\bigr{)}=\frac{M_{1}}{2\chi_{2}M_{2}}\quad,\quad,\frac{\partial^{2}U_{\rm eff}}{\partial u^{2}}\bigg{|}_{u=u_{*}}=\frac{M_{1}^{2}}{6\chi_{2}M_{2}}>0\;.$$

$$U_{\mathrm{eff}}(u_{\ast })=2{\Lambda }_0-\frac{M_1^2}{4{\chi }_2{M}_2}\equiv 2{\Lambda }_{\mathrm{DE}}$$

$$\displaystyle S_{\rm FLRW}=\int d^{4}x\Bigl{[}-6\frac{a\stackrel{{\scriptstyle.}}{{a}}^{2}}{N}+Na^{3}\Bigl{(}\frac{1}{2}\frac{\stackrel{{\scriptstyle.}}{{u}}^{2}}{N^{2}}$$

$$\displaystyle+M_{1}e^{-u/\sqrt{3}}-M_{2}\chi_{2}e^{-2u/\sqrt{3}}-2\Lambda_{0}\Bigr{)}\Bigr{]}\;.$$

$$H^2=\frac{1}{6}\rho ,\rho =\frac{1}{2}{\stackrel{.}{u}}^2+U_{\mathrm{eff}}(u),$$

$$\stackrel{.}{H}=-\frac{1}{4}(\rho +p),p=\frac{1}{2}{\stackrel{.}{u}}^2-U_{\mathrm{eff}}(u),$$

Full text

1 Introduction

Developments in cosmology have been influenced to a great extent by the idea of inflation Starobinsky:1979ty ; Starobinsky:1980te ; Linde:1981mu ; Guth:1980zm ; Albrecht:1982wi , which provides an attractive solution of the fundamental puzzles for the standard Big Bang model, as the horizon and the flatness problems. In addition, providing a framework for sensible calculations of primordial density perturbations were discussed in Mukhanov:1981xt ; Guth:1982ec . However, it has been recognized that a successful implementation requires some very special restrictions on the dynamics that drives inflation. In particular, in New Inflation Linde:1981mu , a potential with a large flat region, which then drops to zero (or almost zero) in order to reproduce the vacuum with almost zero (in Planck units) cosmological constant of the present universe, is required.

In a parallel development, extended (modified) gravity theories as alternatives/generalizations of the standard Einstein General Relativity are being extensively studied in the last decade or so. The main motivation for this development comes from:

• •

(a) Cosmology – modified gravity may solve the problems of dark energy and dark matter and explain the large scale structure and the accelerated expansion of the universe Perlmutter:1998np ; Copeland:2006wr );

• •

(b) Quantum field theory in curved spacetime – because of non-renormalizability of standard general relativity in higher loops it fails to describe the universe at quantum scales weinberg-79 ;

• •

(c) Modern String theory – because of the natural appearance of scalar-tensor couplings and higher-order curvature invariants in low-energy effective field theories aimed at phenomenologically realistic description of particle physics GSW-1 .

The principal approaches to construct modified gravity theories include $f(R)$-gravity, scalar-tensor theories, Gauss-Bonnet gravity models. For detailed accounts, see the book extended-grav-book and the extensive reviews extended-grav ; odintsov-1 ; berti-etal ; odintsov-2 , as well as for further details Refs.Dvali:1998pa -Benisty:2019jqz .

One broad class of actively developed modified/extended gravitational theories is based on employing alternative non-Riemannian spacetime volume-forms, i.e., metric- independent generally covariant volume elements in the pertinent Lagrangian actions on spacetime manifolds with an ordinary Riemannian geometry, instead of the canonical Riemannian volume element $\sqrt{-g}\,d^{4}x$ whose density is given by the square-root of the determinant of the Riemannian metric:

[TABLE]

originally proposed in TMT-orig-0 ; Hehl ; TMT-orig-1 ; TMT-orig-2 ; 5thforce . For a concise geometric formulation, see susyssb-1 ; grav-bags .

This formalism was used as a basis for constructing a series of extended gravity-matter models describing unified dark energy and dark matter scenario dusty ; dusty-2 , quintessential cosmological models with gravity-assisted and inflaton-assisted dynamical suppression (in the “early” universe) or generation (in the post-inflationary universe) of electroweak spontaneous symmetry breaking and charge confinement grf-essay ; varna-17 ; bpu-10 , and a novel mechanism for the supersymmetric Brout-Englert-Higgs effect in supergravity susyssb-1 .

In the present paper we propose a very simple gravity model without any initial matter fields involving several non-Riemannian volume-forms instead of the standard Riemannian volume element $\sqrt{-g}\,d^{4}x$. We show how the non-Riemannian volume-elements, when passing to the physical Einstein frame, generate a canonical scalar field $u$ and manage to create dynamically a non-trivial inflationary-type potential for $u$ with a large flat region for large positive $u$ and a stable low-lying minimum, i.e., $u$ will play the role of a dynamically created “inflaton”. This dynamically generated inflationary potential turns out to be a generalization of the well-known Starobinsky potential Starobinsky:1979ty .

We study the evolution of the cosmological solutions from the point of view of the theory of dynamical systems and calculate the spectral index $n_{s}$ and the tensor-to-scalar ratio $r$ in our model whose values are in accordance with the observational data.

In Section 2 below we briefly review the general notion of volume forms on arbitrary differential manifolds. Section 3 briefly presents the general construction of Lagrangian actions on Riemannian manifolds employing metric-independent (non-Riemannian) volume forms (volume elements). Our main results are contained in Sections 4, 5 and 6. In Section 4 we propose our simple modified gravity model in terms of several non-Riemannian volume elements without any matter fields and derive the corresponding Einstein-frame description with the associated dynamical creation of a canonical scalar field with a non-trivial effective inflationary potential. In Section 5 we study the cosmological evolutionary solutions within the Friedmann-Lemaitre-Robertson-Walker framework. In Section 6 we derive the explicit expressions for the Hubble slow-roll parameters and use them to obtain analytic results for the scalar power spectral index and the tensor-to-scalar ratio which we compare with the available observational data.The last Section 7 contains our conclusions.

2 Non-Riemannian Volume-Forms Formalism

Let us first recall the general notion of volume-forms (volume elements) in integrals over arbitrary differentiable manifolds – not necessarily Riemannian one, so no metric is needed. Volume forms are given by nonsingular maximal rank differential forms $\omega$ (see e.g. Ref.spivak ):

[TABLE]

(our conventions for the alternating symbols $\varepsilon^{\mu_{1},\ldots,\mu_{D}}$ and $\varepsilon_{\mu_{1},\ldots,\mu_{D}}$ are: $\varepsilon^{01\ldots D-1}=1$ and $\varepsilon_{01\ldots D-1}=-1$).

The volume element density $\Omega$, as it is evident from its definition in (2), transforms as scalar density under general coordinate reparametrizations on the manifold.

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}_{\mu}dx^{\mu}$ ($A=0,\ldots,D-1$):

[TABLE]

To construct modified gravitational theories as alternatives to ordinary standard theories in Einstein’s general relativity, instead of $\sqrt{-g}$ we can employ one or more alternative non-Riemannian volume element densities as in (2) given by non-singular exact $D$-forms $\omega=dA$ where:

[TABLE]

Thus, the non-Riemannian volume element density $\Phi(A)$ is defined in terms of the (scalar density of the) dual field-strength of an auxiliary rank $D-1$ tensor gauge field $A_{\mu_{1}\ldots\mu_{D-1}}$ and it transforms as scalar density under general coordinate transformations, which is evident from its definition (4). Accordingly, the integration element $\int d^{4}x\Phi(A)$ is manifestly invariant under general coordinate transformations.

Let us stress that the term “non-Riemannian” relates only to the nature of the volume element density (4), whose definition does not involve the metric. Otherwise the geometry of the spacetime is a regular Riemannian one – scalar products of vector fields are given as usual by the Riemannian metric $g_{\mu\nu}$, the connection $\Gamma_{\mu\nu}^{\lambda}$ is the usual Levi-Civita one in terms of $g_{\mu\nu}$, there is no torsion, etc..

3 The Action

In general, modified gravity Lagrangian actions based on the non-Riemannian volume-form formalism have the following generic form (here and below we are using units with $16\pi G_{\rm Newton}=1$):

[TABLE]

Here $\Phi_{1}(A),\Phi_{2}(B),\Phi_{0}(C)$ are several different non-Riemannian volume element densities of the form (4), i.e., defined by auxiliary rank 3 tensor gauge fields $A_{\mu\nu\lambda},B_{\mu\nu\lambda},C_{\mu\nu\lambda}$; $R$ denotes the scalar curvature in either first-order (Palatini) or second order (metric) formalism; $L^{(1)}$ and $L^{(2)}$ are some matter field Lagrangians; the dots indicate possible additional terms containing higher powers of the non-Riemannian volume element densities e.g., \bigl{(}\Phi_{1}(A)\bigr{)}^{2}/\sqrt{-g}. The specific forms of $L^{(1)}$ and $L^{(2)}$ can be uniquely fixed via the requirement for invariance of (5) under global Weyl-scale invariance (see (10) below).

Let us stress that the modified gravity action (5), in complete analogy with

[TABLE]

which is the standard Einstein-Hilbert action, is explicitly invariant under general coordinate reparametrizations since, as mentioned above, non-Riemannian volume element densities transform as scalar densities similarly to $\sqrt{-g}$.

A characteristic feature of the modified gravitational theories (5) is that when starting in the first-order (Palatini) formalism all non-Riemannian volume-forms are almost pure-gauge degrees of freedom, i.e. they do not introduce any additional physical (field-propagating) gravitational degrees of freedom except for few discrete degrees of freedom with conserved canonical momenta appearing as arbitrary integration constants. The reason is that the modified gravity action in Palatini formalism is linear w.r.t. the velocities of some of the auxiliary gauge field components defining the non-Riemannian volume element densities, and does not depend on the velocities of the rest of auxiliary gauge field components. The (almost) pure-gauge nature of the latter is explicitly shown in Refs.grav-bags ; grf-essay (appendices A) employing the standard canonical Hamiltonian treatment of systems with gauge symmetries, i.e., systems with first-class Hamiltonian constraints a’la Dirac (e.g., henneaux-teitelboim ; rothe ).

Unlike Palatini formalism, the above situation changes significantly when we treat (5) in the second order (metric) formalism. In the latter case the “Einstein-Hilbert” part $\int d^{4}x\,\Phi_{1}(A)R$ of the modified gravity action (5) contains second order time derivative terms of the metric in $R$, which is in sharp contrast with the case of ordinary Riemannian volume element $\int d^{4}x\,\sqrt{-g}R$ where the corresponding second-order time derivatives amount to a total derivative. According to the general canonical Hamiltonian treatment of systems with higher-order time derivatives on the canonical variables (see e.g., gitman-tyutin – modern version of the classical Ostrogradsky formalism ostrogradski ) the presence of the latter implies the appearance of some of the corresponding velocities as additional physical degrees of freedom. In the present case this is reflected in the fact that (as we will see below, Eqs.(19)-(21)) upon passing to the physical Einstein frame via conformal transformation:

[TABLE]

the first non-Riemannian volume element density $\Phi_{1}(A)$ in (5) is not any more a “pure gauge”, but creates a new dynamical canonical scalar field $u$ via $\chi_{1}=\exp{\frac{u}{\sqrt{3}}}$. In the following Section we will see how a non-trivial inflationary potential for $u$ is dynamically generated.

4 Einstein Frame - the Effective Scalar Potential

Let us now consider the simplest member in the class of modified gravitational models (5) with no original matter fields. i.e., $L^{(1)}=0$ and $L^{(2)}=0$, and where we only add a quadratic term w.r.t. non-Riemannian volume element density $\Phi_{1}(A)$:

[TABLE]

Here $R$ is the scalar curvature in the second order (metric) formalism and:

[TABLE]

The specific form of the action (8) is dictated by the requirement about global Weyl-scale invariance under:

[TABLE]

where $\lambda=\text{const}$.

Scale invariance has always played an important role since the original papers on the non-canonical volume-form formalism TMT-orig-1 .

In a more general context let us recall, that scale invariance is a symmetry which relates small scales to large scales. As such (together with conformal symmetry) it plays fundamental role in quantum field theory and modern string theory in particle physics at (ultra)high energies as it dynamically generates (via spontaneous breakdown) mass scales hierarchies. On the other hand it plays an important role in cosmology as well, where it leads naturally to flat inflationary potentials (in the present context this is because it introduces a shift symmetry of the scalar field(s)) and produces candidates for dark matter (see the lectures at CERN’s Workshop on “Scale invariance in particle physics and cosmology”, Ref.CERN-workshop ). Also let us note another specific application of spontaneously broken dilatation symmetry in combination with application of the non-canonical volume form formalism: elimination of the Fifth Force Problem in a quintessential inflationary scenario 5thforce .

The equations of motion resulting from (8) upon variation w.r.t. the auxiliary gauge fields $A_{\mu\nu\lambda},\,B_{\mu\nu\lambda}\,,C_{\mu\nu\lambda}$ yield, respectively:

[TABLE]

Here $M_{1},M_{2}$ and $\chi_{2}$ are (dimensionful and dimensionless, respectively) integration constants. The appearance of $M_{1},M_{2}$ indicate spontaneous breaking of global Weyl symmetry (10).

The equations of motion w.r.t. $g_{\mu\nu}$ from (8) read:

[TABLE]

with $\chi_{1}$ as in (7). On the other hand, taking the trace of (13) and using Eq.(11) we obtain the equation of motion for $\chi_{1}$:

[TABLE]

We now transform Eqs.(13) and (14) via the conformal transformation (7) and show that the transformed equations acquire the standard form of Einstein equations w.r.t. the new “Einstein-frame” metric ${\bar{g}}_{\mu\nu}$. To this end we are using the known formulas for the conformal transformations of $R_{\mu\nu}$ and $\Box\Psi$, the latter being an arbitrary scalar field, in particular $\Psi\equiv\chi_{1}$ (see e.g. Ref.dabrowski ; bars indicate magnitudes in the ${\bar{g}}_{\mu\nu}$-frame):

[TABLE]

[TABLE]

and

[TABLE]

Following the analogous derivation in Ref.benisty-1 , upon using (15)-(16) we rewite Eqs.(13) as:

[TABLE]

where we have redefined:

[TABLE]

in order to obtain a canonically normalized kinetic term for the scalar field $u$, and where:

[TABLE]

On the other hand, using (16) we rewrite Eq.(14) in terms of the canonical scalar field $u$:

[TABLE]

with $U_{\rm eff}$ as in (19).

Accordingly, the corresponding Einstein-frame action reads:

[TABLE]

We now observe an important result – in (21) we have a dynamically created scalar field $u$ with a non-trivial effective scalar potential $U_{\rm eff}(u)$ (19) entirely dynamically generated by the initial non-Riemannian volume elements in (8) because of the appearance of the free integration constants $M_{1},M_{2},\chi_{2}$ in their respective equations of motion (11)-(12).

The qualitative shape of (19) is depicted on Fig.1. The effective potential $U_{\rm eff}(u)$ has two main features relevant for cosmological applications. First, $U_{\rm eff}(u)$ (19) possesses a flat region for large positive $u$ and, second, it has a stable minimum for a small finite value $u=u_{*}$:

(i) $U_{\rm eff}(u)\simeq 2\Lambda_{0}$ for large $u$;

(ii) $\frac{{\partial{U_{\rm eff}}}}{{\partial{u}}}=0\;$ for $u\equiv u_{*}$ where:

[TABLE]

The flat region of $U_{\rm eff}(u)$ for large positive $u$ correspond to “early” universe’ inflationary evolution with energy scale $2\Lambda_{0}$. On the other hand, the region around the stable minimum at $u=u_{*}$ (22) correspond to “late” universe’ evolution where the minimum value of the potential:

[TABLE]

is the dark energy density value Angus:2018tko ; Zhang:2018gbq .

Let us note that the effective potential $U_{\rm eff}$ (19) generalizes the well-known Starobinsky inflationary potential Starobinsky:1979ty ((19) reduces to Starobinsky potential upon taking the following special values for the parameters: $\Lambda_{0}=\frac{1}{4}M_{1}=\frac{1}{2}\chi_{2}M_{2}$).

5 Evolution of the homogeneous solution

We now consider reduction of the Einstein-frame action (21) to the Friedmann-Lemaitre-Robertson-Walker (FLRW) setting with metric $ds^{2}=-N^{2}dt^{2}+a(t)^{2}d{\vec{x}}^{2}$, and with $u=u(t)$. In order to study the evolution of the scalar field $u=u(t)$ and the Friedmann scale factor $a=a(t)$, it is useful to use the method of autonomous dynamical systems.

The FLRW action describes a minimally coupled canonical scalar field $u$ with specific potential $U_{\rm eff}(u)$ (19) (using again units with $16\pi G_{\rm Newton}=1$):

[TABLE]

Variations w.r.t. $N,a,u$ (and subsequently using the gauge $N=1$ for the lapse function) yield the pertinent Friedmann and field equations ($H=\stackrel{{\scriptstyle.}}{{a}}/a$ being the Hubble parameter):

[TABLE]

In the treatment of Eqs.(25)-(27) it is instructive to rewite them in terms of a set of dimensionless parameters (following the approach in Ref.Bahamonde:2017ize ):

[TABLE]

with $L_{\rm DE}$ as in (23). In these coordinates the system defines a closed orbit:

[TABLE]

which is equivalent to the first Friedmann equation (25).

Employing the variables $(x,y,z)$ in Eqs.(25)-(27) and taking into account the constraint (29) we obtain the autonomous dynamical system w.r.t. $(x,z)$:

[TABLE]

where the primes denote derivative w.r.t. the parameter ${\mathcal{N}}=\log a$, and the function $\xi(x,z)$ is defined as:

[TABLE]

The phase portrait of the system (30) is depicted on Fig.2. There are two critical points in the system. The stable point $A\left(x=0,z=1\right)$ corresponds to the “late” universe de Sitter solution with the asymptotic cosmological constant $\Lambda_{\rm DE}$ (23).

The second point $B\left(x=0,z=\sqrt{\Lambda_{\rm DE}/\Lambda_{0}}\right)$ is unstable corresponding to the beginning of the universe’ evolution in the “early” universe at large $u$. If the evolution starts at any point close to $B$, initially the evolution is of de Sitter type with effective cosmological constant $\approx\Lambda_{0}$. Then the dynamics drives the system away from $B$ all the way towards the stable point $A$ at late times.

Numerical solutions are demonstrated in Fig. 3. One can see that the Hubble parameter begins and ends with two different values. The first one is related to the inflationary epoch and the other related to the dark energy in the late universe. The scalar field $u$ oscillates around the minimum point $u_{*}$ (22) of $U_{\rm eff}$ (19), which corresponds to particle creation in the reheating epoch.

6 Perturbations

In order to check the viability of the model we investigate the perturbations of the above background evolution, in particular focusing on the inflationary observables such as the scalar spectral index $n_{s}$ and the tensor-to-scalar ratio $r$. As usual, we introduce the Hubble slow-roll parameters Martin:2013tda ; Bamba:2016wjm , which in our case using the potential (19) read:

[TABLE]

where:

[TABLE]

with $\Lambda_{\rm DE}$ – the dark energy density (23), and therefore $\delta$ very small.

Inflation ends when $\epsilon(u_{f})=1$ for some $u=u_{f}$ where ($\zeta_{f}\equiv\frac{M_{2}\chi_{2}}{M_{1}}e^{-u_{f}/\sqrt{3}}$):

[TABLE]

For the number of $e$-foldings ${\mathcal{N}}=\frac{1}{2}\int_{u_{i}}^{u_{f}}du\;U_{\rm eff}/U_{\rm eff}^{\prime}$ we obtain:

[TABLE]

where $\zeta_{i}\equiv\frac{M_{2}\chi_{2}}{M_{1}}e^{-u_{i}/\sqrt{3}}$ and $u=u_{i}$ is very large corresponding to the start of the inflation. Ignoring $\delta$ and using the last equality (35) we have approximately:

[TABLE]

Using the slow-roll parameters, one can calculate the values of the scalar spectral index and the tensor-to-scalar ratio respectively as Nojiri:2019kkp ; Dalianis:2018frf :

[TABLE]

Taking into account Eqs.(32),(33) (ignoring $\delta$) and (37) we find:

[TABLE]

and

[TABLE]

where $u_{i}({\mathcal{N}})$ is the solution of the transcedental Eq.(37) for $u_{i}$ as a function of ${\mathcal{N}}$.

One viable example in our model is to take ${\mathcal{N}}=60$ $e$-folds. Eq.(37) yields ${\mathcal{N}}=60$ provided we choose $\frac{M_{1}}{M_{2}\chi_{2}}\simeq 10^{-2}$, which yields $u_{i}\simeq 17$. In such a way the observables are predicted to be:

[TABLE]

which are well inside the PLANCK observed constraints Akrami:2018odb :

[TABLE]

Fig. 4 demonstrates the relation between the number of $e$-folds and the dimensionless parameters. One can see that all those values fit the latest PLANCK collaboration constraints.

7 Conclusions

We propose a very simple gravity model without any initial matter fields in terms of several alternative non-Riemannian spacetime volume elements within the second order (metric) formalism. We show how the non-Riemannian volume-elements, when passing to the physical Einstein frame, create a canonical scalar field and produce dynamically a non-trivial inflationary-type potential for the latter with a large flat region and a stable low-lying minimum. We study the evolution of the cosmological solutions from the point of view of the theory of dynamical systems. Our model predicts scalar spectral index $n_{s}\approx 0.96$ and tensor-to-scalar ratio $r\approx 0.002$ for 60 $e$-folds, which is in accordance with the observational data.

A natural next step is to consider two-field inflation (Refs. Bahamonde:2017ize ; 2-field-inflation-1 ; 2-field-inflation-1a ; 2-field-inflation-2 ; Ben-Dayan:2018mhe ; Artymowski:2019vfy ; Ben-Dayan:2014isa ; Benisty:2019bmi ; Aragam:2019khr ; Kubo:2018kho , for a geometric treatment see Refs. lili ; MadrizAguilar:2019wgd ; Barvinsky:2019qzx ; Ahmad:2019jbm ; Granda:2019jqy , and references therein) by adding a new scalar field $\varphi$ with non-trivial potentials in the starting modified gravity action (8) built in terms of several non-Riemannian volume elements and subject to preserving the requirement of global Weyl-scale invariance (10).In this case the non-Riemannian volume elements will again generate a second scalar field $u$ and create dynamically a non-trivial two-field scalar potential with a very specific geometry of the field space of $\varphi,u$. This is studied in more detail in our subsequent work Benisty:2019bmi , where it is shown that the latter dynamically generated two-field inflationary model similarly conforms to the observational data.

Acknowledgements.

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 the action CA18108. D.B., E.N. and S.P. acknowledge illuminating discussions with Lilia Anguelova. Finally we would like to thank the referees whose comments contributed to significant improvement of the presentation.