A simple $F({\cal R},\phi)$ deformation of Starobinsky inflationary model

TL;DR

This paper explores a modified inflationary model combining $F({ m R})$ gravity with a scalar field, resulting in two-field potentials that align with observational data, extending the Starobinsky model.

Contribution

It introduces a new $F({ m R}, \, \phi)$ deformation of the Starobinsky model with specific scalar potentials, demonstrating viability with current observations.

Findings

All three potential cases produce viable inflationary scenarios.

The models are conformally equivalent to Einstein-Hilbert theory with two scalar fields.

Results are consistent with current cosmological data.

Abstract

We study a model including a real scalar field $\phi$ non-minimally coupled to $F({\cal R})$ gravity, which is conformally equivalent to an Einstein-Hilbert theory, involving two real scalar fields. We consider three special cases of the potential of the field $\phi$ in the $F({\cal R})$-frame: a vanishing potential, a mass term and a Higgs potential. All these lead to non-trivial two-field potentials in the Einstein-frame which in particular directions resemble the well-known Starobinsky model. We find, that all these cases can yield viable inflationary models in complete agreement with current observational data.

Figures

Equations

$$n_s$$

$$r$$

$$f_{NL}$$

$$S=\int d^{4}x\sqrt{-\tilde{g}}\left[-\frac{1}{2}\tilde{\mathcal{R}}+\frac{1}{12M^2}{\tilde{\mathcal{R}}}^2\right],$$

$$S=$$

$$S=\int d^{4}x\sqrt{-\tilde{g}}\left[F(\tilde{\mathcal{R}},\varphi )+\frac{1}{2}{\tilde{g}}^{\mu \nu }{\partial }_{\mu }\varphi {\partial }_{\nu }\varphi -U(\varphi )\right].$$

$$S=$$

$$\frac{\partial \mathcal{L}}{\partial \psi }=0\Rightarrow \Phi =\tilde{\mathcal{R}}.$$

$$\frac{\partial \mathcal{L}}{\partial \Phi }=0\Rightarrow \frac{\partial F(\Phi ,\varphi )}{\partial \Phi }=-\psi \Rightarrow \Phi =\xi (\psi ,\varphi ).$$

$$S=$$

$$\tilde{\mathcal{R}}=2\psi \mathcal{R}+\frac{3}{\psi }{g}^{\mu \nu }{\partial }_{\mu }\psi {\partial }_{\nu }\psi +6{\psi }^2{\nabla }_{\mu }\left(\frac{{\partial }^{\mu }\psi }{{\psi }^2}\right).$$

$$S=$$

$$V(\varphi ,\rho )=$$

$$F(\mathcal{R},\varphi )=-\frac{1}{2}f(\varphi )\mathcal{R}+\frac{1}{12M^2(\varphi )}{\mathcal{R}}^2$$

$$V(\varphi ,\rho )=e^{-2\sqrt{\frac{2}{3}}\rho }\left[\frac{3}{4}{M}^2(\varphi ){\left(f(\varphi )-e^{\sqrt{\frac{2}{3}}\rho }\right)}^2+U(\varphi )\right],$$

$$d{s}^2=d{t}^2-{\alpha }^2(t)d{\vec{x}}^2$$

$$\ddot{\rho }+3H\dot{\rho }+\frac{1}{\sqrt{6}}{e}^{-\sqrt{\frac{2}{3}}\rho }(\dot{\varphi }{)}^2$$

$$\ddot{\varphi }+3H\dot{\varphi }-\sqrt{\frac{2}{3}}\dot{\rho }\dot{\varphi }$$

$$3H^2$$

$$\dot{H}$$

$$S=\int d^{4}x\sqrt{-g}\left[-\frac{R}{2}+\frac{1}{2}{\mathcal{G}}_{IJ}{g}^{\mu \nu }{\partial }_{\mu }{\varphi }^I{\partial }_{\nu }{\varphi }^J-V({\varphi }^I)\right].$$

$${\mathcal{D}}_J{A}_I=A_{I,J}-{\Gamma }_{IJ}^K{A}_K,$$

$${\Gamma }_{IJ}^K=\frac{1}{2}{\mathcal{G}}^{KL}({\partial }_I{\mathcal{G}}_{LJ}+{\partial }_J{\mathcal{G}}_{IL}-{\partial }_L{\mathcal{G}}_{IJ}).$$

$$R_{\mu \nu }-\frac{R}{2}{g}_{\mu \nu }=T_{\mu \nu },$$

$$T_{\mu \nu }=G_{IJ}{\partial }_{\mu }{\varphi }^I{\partial }_{\nu }{\varphi }^J-g_{\mu \nu }\left[\frac{G_{IJ}}{2}{g}^{\rho \lambda }{\partial }_{\rho }{\varphi }^I{\partial }_{\lambda }{\varphi }^J+V({\varphi }^I)\right].$$

$$g^{\mu \nu }{\varphi }_{;\mu ;\nu }^I+g^{\mu \nu }{\Gamma }_{JK}^I{\partial }_{\mu }{\varphi }^J{\partial }_{\nu }{\varphi }^K-G^{IK}{V}_{,K}=0$$

$$ϵ$$

$$\eta$$

$$ϵ\ll 1\text{and}\eta \ll 1.$$

$${\varphi }^I(x^{\mu })={\phi }^I(t)+\delta {\varphi }^I(x^{\mu }),$$

$$d{s}^2=(1+2A)d{t}^2-2\alpha (t){\partial }_{i}Bd{x}^{i}dt-{\alpha }^2(t)[(1-2\psi ){\delta }_{ij}+2{\partial }_i{\partial }_{j}E]d{x}^{i}d{x}^j.$$

$$Q^I=\delta {\varphi }^I+\frac{{\dot{\phi }}^I}{H}\psi .$$

Full text

1 Introduction

The mechanism of cosmological inflation [1, 2, 3, 4, 5] was first introduced in 1980’s in order to solve crucial problems of the Big Bang Cosmology, such as the horizon, the flatness and the Monopole problems. Inflation is a period of accelerated (quasi-de Sitter) expansion of the very early Universe, which elegantly allows for near large-scale homogeneity and spatial flatness of our Universe. An extra bonus of introducing inflation in Standard Cosmology is that it can explain the formation of large-scale structure, being the only known mechanism to do this. Quantum fluctuations during the inflationary epoch presumably seeded the perturbations which grew under gravitational instability into the structures we observe today [6]. Due to this, the inflationary mechanism has been intensively studied resulting to a better theoretical understanding of it. Also, in recent years, the interest in inflationary cosmology has grown considerably because of the great amount of data made available sourcing from various cosmological surveys. Despite its successful predictions the origin of inflation is not well understood, as yet. It is more like a phenomenological construction, whose origin should be sought in some fundamental theory, such as high-energy particle physics or gravity.

In particle physics, superstring theory and supergravity, multiple scalar fields are involved and some of them may play the role of inflaton. Furthermore, in curved spacetime and in the context of renormalization of scalar fields we have the arise of non-minimal couplings between scalar fields and the Ricci scalar [7, 8]. Thus it is reasonable to search for inflationary models which include many fields non-minimally coupled to gravity, whose potential energy dominates the energy-momentum tensor and drives inflation [9, 10, 11, 12, 13, 14, 15]. For a model of inflation to be viable it should be in agreement with the recent observational constraints for the spectral index, the tensor-to-scalar ratio and non-gaussianity (for multi-field inflation). These quantinties are defined in the context of cosmological perturbations [16, 17, 18, 19, 6, 20, 21, 22, 23, 24, 25, 26, 27] and their constraints obtained from observetions of cosmic microwave background (CMB), according to [28, 29], are:

[TABLE]

Among the single-field models of inflation, various classes of models can be in agreement with the aforementioned constraints. The first is the well studied non-minimally coupled Higgs inflation [30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49] which provides a particle origin to inflaton, but is also strongly connected with gravity. The second is the class of models of chaotic inflation [50, 51, 52] and its variants, and another is natural inflation models [53, 54]. Also a well-known motivated model is the Starobinsky model of inflation [1], which remarkably was proposed almost four decades ago and furnishes a gravitational origin to inflaton. This model can be seen as the simplest inflationary model within the context of $\cal F(R)$ theories of gravity [55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67], as the only extension from Einstein Gravity is the addition to the Hilbert - Einstein action of an extra ${\cal R}^{2}$ term

[TABLE]

where above $M$ is a parameter with dimensions of mass, $\tilde{g}_{\mu\nu}$ is the metric tensor and $\cal\tilde{R}$ is the Ricci scalar. In this paper, the reduced mass Planck is dimensionless and equal to 1, and we use the $(+,-,-,-)$ spacetime signature notation. This action is classically equivalent, through a conformal transformation, to the following scalar-tensor theory with a non-minimal coupling between the scalar field $\phi$ and gravity

[TABLE]

Using the slow-roll inflation formalism and the agreement with the observables the value of $M$ is restricted to $M\simeq 1.3\cdot 10^{-5}$.

The great success of Starobinsky inflation model and its elegant physical interpretation of inflaton has led to an intensive study of inflationary models that are extensions or modifications of this model in the framework of $F({\cal R})$ theories or $F({\cal R},\phi)$ theories studied in the metric formalism [68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95] (see for instance [96] for a beautiful discussion of inflation in the Jordan frame). The latter have the meaning of being $F({\cal R})$ theories in the presence of scalar fields, which are in general non-minimally coupled to gravity. An alternative variational principle leading to the equations of motion of General Relativity, which has been intensively studied for cosmological purposes [97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112], is the Palatini formalism [57, 61, 113, 114, 115], in which the metric tensor $g_{\mu\nu}$ and the connection $\Gamma_{\mu\nu}^{\lambda}$ are treated as independent variables. Motivated by the multi-field scope of particle physics and the viability of Starobinsky inflation, in this paper we study the robustness of Starobinsky inflation in the presence of a scalar field non-minimally coupled to gravity both in the $\cal R$ and the ${\cal R}^{2}$ term.

This paper is organized as follows. In Section 2 we present the general theoretical setup of the model at hand. In Section 3 we specialize our study and obtain numerical results for the observables for the case where the pre-existing scalar field is massless and its potential is zero. In Section 4 we do the same work with the addition of a mass term for the pre-existing scalar field. In Section 5 we identify $\phi$ with the SM Higgs boson at the electroweak scale. We conclude and discuss potential extensions and future study of this model in Section 6.

2 Theoretical setup

Our starting point is an inflationary model including a real scalar field $\phi$ non-minimally coupled to $F({\cal R})$ gravity. This model, is described by the action, in the $F({\cal R})$-frame,

[TABLE]

For this theory we can obtain a better physical understanding by working in the Einstein frame. It can easily be seen that this theory is classically equivalent to

[TABLE]

as the equations of motion for the field $\psi$, which plays the role of a Lagrange multiplier, yield

[TABLE]

The corresponding equations of motion for $\Phi$ are then

[TABLE]

Thus the action (2.2) can be written in the following form in the Jordan frame

[TABLE]

In order to pass to the Einstein frame we perform a conformal transformation of the metric tensor $\tilde{g}_{\mu\nu}=g_{\mu\nu}/2\psi$ [7, 8]. Under this transformation the Ricci scalar transforms as

[TABLE]

The action, after eliminating a total derivative, is given by

[TABLE]

Finally, using the field redefinition $2\psi=e^{\sqrt{\frac{2}{3}}\rho}$, which leads to a canonical kinetic term for the field $\rho$, we obtain the following form for the action in the Einstein frame:

[TABLE]

In this the potential $V(\phi,\rho)$ is given by

[TABLE]

An interesting model belonging to this class of theories is that proposed in [79] where

[TABLE]

In this $f(\phi)$ and $M^{2}(\phi)$ are two generic functions, which should be positive defined in order to avoid ghosts. In this model one can easily find, using equation (2.9), that the potential $V(\phi,\rho)$ is given by the following expression

[TABLE]

which however is in disagreement with the result derived in [79] 111The potential found therein is of the form $V(\phi,\rho)=\frac{3}{4}M^{2}(\phi)\left[f(\phi)-e^{\sqrt{\frac{2}{3}}\rho}\right]^{2}+e^{2\sqrt{\frac{2}{3}}\rho}U(\phi)$.. In the following we shall study the cosmological prediction of this model with the potential as given by (2.11).

The equations of motion for the two scalar fields $(\rho,\phi)$ in a spatially flat Friedman-Robertson-Walker spacetime

[TABLE]

assuming that the fields are only time-dependent, based on the observed homogeneity and isotropy of our Universe at large scales, following from the action (2.8) are given by

[TABLE]

where we denote $\dot{}\equiv d/dt$, $V_{,I}=\partial V/\partial\phi^{I}$ with $I=\phi,\rho$ and $H=\dot{\alpha}/\alpha$ is the Hubble rate. The Einstein equations for the action (2.8) lead to the following Friedmann equations

[TABLE]

In order to find the time-evolution of the fields $\rho$, $\phi$ and the scale factor we need just to solve the equations of motion (2.13), (2.14) and (2.15). The equation (2.16) is not independent but it is related to the other three equations of motion.

It is worth noting that the action (2.8) can be seen as a special case of the well-studied generalized non-linear sigma model of multifield inflation

[TABLE]

In the above expression, the Latin indices account for the number of the fields of the theory and ${\cal G}_{IJ}$ is the metric tensor of the curved field space manifold. Then, in correspondence with the theory of Gravity, we have the definition for the covariant derivative:

[TABLE]

with $A_{I}$ being a vector in the field space and $\Gamma^{K}_{IJ}$ being the corresponding Christoffel symbols in the curved field space, calculated by the expression

[TABLE]

The Einstein equations for the action (2.17) are of the usual form

[TABLE]

where the stress tensor is given by the expression

[TABLE]

The equations of motion for the fields $\phi^{I}$ are given by the expression

[TABLE]

In this generalized model the slow-roll parameters are defined [89] as222In this work we assume slow-roll, slow-turn solutions. For highly curved field spaces the dominant behaviour can be of different nature (see e.g. [116]).

[TABLE]

and should obey the ordinary (from one-field inflation models) slow-roll conditions

[TABLE]

Also for this case of inflationary models, there have been derived some well known expressions for the spectral index and the tensor-to-scalar ratio that we measure today from the CMB spectrum. The definition of $n_{s}$ and $r$ happens in the context of cosmological perturbations. In the framework of this theory we have quantum fluctuations of the fields sourcing perturbations in the metric and vice versa. Let us briefly mention, using the description used in [10, 11, 12], the basic concepts for the understanding of the quantities needed for the calculation of $n_{s}$ and $r$ in the concept of multi-field inflation. For a more complete view of this theory we recommend the reader to see these articles and the references therein [19, 6, 21, 9, 10, 11, 12, 13, 14, 15, 20, 22, 23, 24, 25, 26].

The results that we will quote in the following are obtained by keeping only first order terms in the perturbative expansion of the spacetime-dependent fields, $\varphi^{I}(x^{\mu})$, around the time-dependent background fields, $\phi^{I}(t)$

[TABLE]

and of the spacetime metric around the spatially flat FRW background

[TABLE]

In order to have a gauge-invariant formalism and due to the fact that the fluctuations $\delta\phi^{I}$ do not transform covariantly under spacetime gauge transformations, it is necessary to introduce the gauge-invariant Mukhanov-Sasaki variables

[TABLE]

Using these variables we can rewrite the equations of motion (2.22) separating to background and first-order expressions as

[TABLE]

and

[TABLE]

where $k$ is the comoving wave number and the derivative ${\cal D}_{t}$ is defined as

[TABLE]

while the mass-squared matrix ${\cal M}_{J}^{I}$ is defined as

[TABLE]

A useful quantity for the simplification of the analysis of the multi-field cosmological perturbations is the length of the velocity vector for the background fields given by

[TABLE]

from which one can define the adiabatic (curvature) and entropy (isocurvature) directions in the curved field space via the unit vectors

[TABLE]

where the turn-rate vector is given by $\omega^{I}(t)\equiv{\cal D}_{t}\hat{\sigma}^{I}$ and $\omega=|\omega^{I}|$. In a two-field model any vector in field space can be decomposed into components along these two unit vectors

[TABLE]

Using this we can decompose the Mukhanov-Sasaki variables along adiabatic and entropy directions

[TABLE]

Using these new definitions of the Mukhanov-Sasaki variables and making a Fourier transformation of the form $\partial_{i}\partial^{i}F(t,x^{i})=-(k^{2}/\alpha^{2})F_{k}(t)$ we obtain from equation (2.30) the following two equations

[TABLE]

and

[TABLE]

with $\Psi$ the gauge-invariant Bardeen potential

[TABLE]

Above we used the following definitions for the mass-squared matrix

[TABLE]

The effective mass of the isocurvature perturbations is defined by

[TABLE]

If $\mu_{s}^{2}/H^{2}\gg 1$, the isocurvature perturbations are heavy during slow-roll, so the growth of $Q_{s}$ is strongly suppressed. On the other hand a tachyonic mass, $\mu_{s}^{2}<0\,,$ leads to the rapid amplification of the isocurvature modes.

The gauge-invariant curvature perturbation and the renormalized entropy perturbation are expressed in terms of the $Q_{s}$ and $Q_{\sigma}$ via the relations

[TABLE]

In the long-wavelength and slow-roll limits these two functions can be found (from the differential equations of the Mukhanov-Sasaki variables) to obey the differential equations

[TABLE]

with the function $\beta(t)$ given by

[TABLE]

In the above equation $n_{ss}$ and $n_{\sigma\sigma}$ are the slow-roll parameters defined by

[TABLE]

The solutions of the equations (2.43) and (2.44) can be written in the form

[TABLE]

where

[TABLE]

are the transfer functions which relate the gauge-invariant perturbations at the time $t_{*}$, when the perturbations of pivot scale during the inflationary era crossed outside the Hubble volume for the first time, to their value at some later time $t$.

The dimensionless primordial spectra ${\cal P}_{R}(k)$, ${\cal P}_{S}(k)$ and $C_{RS}(k)$ are defined at some time $t$ after inflation during the radiation-dominated era when all cosmological scales are still outside the horizon

[TABLE]

Due to the fact that we are interested in the calculation of the power spectrum ${\cal P}_{R}(k)$ to first-order in the slow-roll parameters we keep only term to zeroth-order in ${\cal P}_{R}(k,t_{*})$, ${\cal C}_{RS_{*}}(k,t_{*})$ and ${\cal P}_{S}(k,t_{*})$, where we have

[TABLE]

Thus replacing back to (2.53) we obtain the equation which relates the power spectrum of curvature perturbations at time $t_{*}$ to its value a latter time $t$

[TABLE]

The spectral index can be determined via the power spectrum of the curvature perturbations from the expression

[TABLE]

Using the fact that the derivative with respect to $\ln k$ can be converted to a time derivative via the relation $d\ln k/dt=(1-\epsilon)H$ and that for the second and the third term of (2.57) holds true

[TABLE]

we find that (2.57) takes the form

[TABLE]

where we have used the definitions

[TABLE]

As it regards the tensor-to-scalar ratio, $r$, it can be determined from the ratio between the power spectrum of the tensor perturbations and that of the curvature perturbations, meaning

[TABLE]

Due to the vanish of $T^{i}_{j}$ terms for $i\neq j$ to first-order the tensor perturbations evolve in the same way as in the single-field models and thus ${\cal P}_{T}$ is given by the expression

[TABLE]

Thus for the tensor-to-scalar ratio in two-field models from equations (2.56) and (2.63) we obtain

[TABLE]

From the above expression it is obvious that if the maximum value $r_{\text{max}}=16\epsilon(t_{*})$ is in agreement with the observed constraint (1.1), there is no need of calculation of the transfer functions, as far as $r$ is concerned.

For the system at hand one can easily infer, by comparing (2.8) and (2.17), identifying $\phi^{1}$ by $\rho$ and $\phi^{2}$ by $\phi$, that the metric tensor ${\cal G}_{IJ}$ has the following form

[TABLE]

Then, by the use of (2.19) we find that the only non-zero components of $\Gamma^{K}_{IJ}$ are

[TABLE]

Thus, using the definitions (2.23), (2.24) and (2.18) for our case, we are leaded to the following expressions for the slow-roll parameters

[TABLE]

where the unit vectors $\hat{\sigma}^{\phi}$, $\hat{\sigma}^{\rho}$ appearing in the equations above are

[TABLE]

So far we have not specified the functions $U(\phi)$, $f(\phi)$ and $M(\phi)$ defining the function $F(\tilde{{\cal R}},\phi)$, (2.10), which defines the model given in (2.1). In order to proceed to predictions of the inflationary observables within the context of this type of models in the following sections, and for different forms of the potential $U(\phi)$, we choose the functions $M^{2}(\phi)$ and $f(\phi)$ to be the same as those employed in [79],

[TABLE]

In these $\alpha$ and $\beta$ are some constants signaling deviations from Starobinsky model. We mainly focus on the parameter regime where $\alpha,\beta>0.$ Negative values for the parameters $\alpha$ and $\beta$ are also acceptable provided that the functions $M^{2}(\phi)$ and $f(\phi)$ remain positive to avoid ghosts.

In the next sections we will be concentrated on the calculation of the spectral index $n_{s}$ and the tensor-to-scalar ratio $r$ in the pivot scale. The main difference of our two-field inflationary model against the Starobinsky model is the presence of isocurvature perturbations. Those perturbations turn out to be very small and undetectable so far. So, we did not examine primordial non-Gaussianities in our deformation of Starobinsky inflationary model, because we expect it to be negligible, like in the pure Starobinsky model of inflation and thus in agreement with the constraints [28, 29]. Also, for all the cases which will be studied in the following sections we found that the effective mass, $\mu_{s}$, of the isocurvature perturbations is small and negative and thus the formalism we presented in this section is applicable.

3 The Case $U(\phi)=0$

In this section we review the special case, already studied in [79], when $\phi$ is a free massless field, in the $F(\cal{R})$-frame. We have $U(\phi)=0$ and thus the scalar potential $V(\phi,\rho)$ becomes

[TABLE]

which is semi-positive definite. From the above expression we observe that for $\alpha=\beta=0$ we retain the Starobinsky potential, as we expected. For non-zero values of $\alpha$ and $\beta$ we have a deformation of the potential in the $\phi$ direction, leading to the stabilization of the field $\phi$. In the case where only $\alpha=0\,,$ the parameter $\beta$ cannot grow well beyond 50. Actually, the marginal case where $\beta=50\,,$ is marginally out due to the non agreement with the $n_{s}$ data. The increase of the parameter $\beta$ does not have a significant effect on the value of $r$. On the other hand, if $\beta=0\,,$ then $\alpha$ cannot be bigger than $\sim 10^{-3}\,.$ Also, according to our numerical calculations if both $\alpha$ and $\beta$ grow well beyond 10 the inflation becomes again non-viable for the same reason, namely our model fails to get the observed value of $n_{s}$. From now on and in the following, for the parameters $\alpha,\beta$ and $M$ we choose for definiteness the values $\alpha=0.01,$ $\beta=0.001$ and $M=1.3\cdot 10^{-5}$, like in the pure Starobinsky model. In Figure 1 we show the profile of the potential for this choice of parameters. It is worth observing that in this potential we do not have a unique Minkowski vacuum but rather a valley of Minkowski vacua on the contour determined by the points $(\phi_{\text{min}},\rho_{\text{min}})=(c,\sqrt{\frac{3}{2}}\ln(1+\alpha c^{2})),$ where $c\in{\rm I\!R}.$ As we shall see below, the choice of the initial values of the fields $\phi$ and $\rho$ lead us to the vacuum $(\phi_{\text{min}},\rho_{\text{min}})=(0,0)$. The choice of different initial conditions will have as a result an altered vacuum which can lead us to successful inflation.

In order to calculate the scalar spectral index and the tensor-to-scalar ratio in the pivot scale, we solve numerically the equations of motion (2.13) , (2.14) and (2.15) to find the time evolution of the fields $\phi$ and $\rho$. Our results are shown in Figure 2. In all graphs displayed the time $t$ is normalized by $M$. From Figure 2 we can see that at the end of inflation the field $\rho$ oscillates near in its approach to the minimum of the potential, whereas the field $\phi$ does not. In fact $\phi$ drops off rapidly. The oscillation of the field $\rho$ can be seen after the time $t=605$. Using the numerical results for $\phi$ and $\rho$, we then solve the equation (2.15), find the time dependence of the scale factor $a(t)$ and plot the time evolution of the Hubble parameter and the number of e-foldings, as shown in Figure 4. From the bottom diagram of Figure 4 we can find the pivot scale used for the calculation of the spectral index and the tensor-to-scalar ratio. Furthermore, using the relations (2.67) and (2.68) we calculate the slow-roll parameters $\epsilon$ and $\eta$, respectively, and we present their behaviour as a function of time in Figure 5. From that figure we can see that the slow-roll parameters $\epsilon$ and $\eta$ indeed obey the slow-roll conditions (2.25). Finally, using relations (2.61) and (2.64) we calculate the spectral index $n_{s}$ and the max value of the tensor-to-scalar ratio in the pivot scale. For $N\simeq 50-60$ we find that333The calculation of $n_{s}$ and $r$ has been made between 50 and 60 e-foldings, so the $\pm$ or $\mp$ signs correspond to 50 (down sign) and 60 (upper sign) e-foldings respectively.

[TABLE]

These values for $n_{s}$ and $r$ are in agreement with the constraints (1.1) and thus lead to viable inflation. Scales of cosmological interest first crossed the Hubble radius between $50$ and $60$ e-foldings before the end of inflation. As indicated in Figure 3, $T_{RS}$ remains small between $N=50\,\,\text{and}\,\,60,$ so corrections to $n_{s}$ and $r$ from $T_{RS}$ remain negligible.

It is worth pointing out, that trying different values for the initial conditions of the fields we observed that for bigger values of $\phi$ we have smaller duration of the inflationary period (number of e-foldings) but always in agreement with the observational constraints (1.1). In fact we have a wide range of allowed initial values for the field $\phi$, for every initial value of $\rho$, to obtain at least $50-60$ e-foldings of inflation. Therefore we do not have fine tuning for the initial conditions since a wide range of them leads to realistic inflation.

4 The Case $U(\phi)=\frac{1}{2}m^{2}\phi^{2}$

Now we make a simple but physical extension of the calculation we made in the previous section, by adding a mass term for the field $\phi$ in our action. The scalar potential $V(\phi,\rho)$ in this case takes the form

[TABLE]

This model for the case of $\alpha=\beta=0$ was studied in references [80, 89] where it was found that it yields viable inflation without the need of fine-tunning. Specifically, in [89], using the $\delta N$ formalism, the non-gaussianity was calculated and was found to be in agreement with the observed constraints [28, 29]. Thus our expectation of feasible non-gaussianity for the case $\alpha=0.01$ and $\beta=0.001$ is strengthened. In this model, values for $\alpha$ and $\beta$ greater than 10 but less than $10^{2}$ are also acceptable as we have checked numerically.

In Figure 6 we see the profile of the potential for the choice of parameters $\alpha=0.01$, $\beta=0.001$ and $m=M=1.3\cdot 10^{-5}$. From this graph we observe, in contrast to the potential considered in the previous section, that in this potential we have a unique Minkowski vacuum, at the point $(\phi_{\text{min}},\rho_{\text{min}})=(0,0).$ For the mass $m$ we study the cases between $m/M=1$ and $m/M=10^{3}$. Our results for the fields, the Hubble parameter and the number of e-foldings, as functions of time for three typical values of the ratio $m/M$ are presented in Figure 7 and Figure 8. In all the cases we use the same initial conditions for the fields $\rho=7.4$ and $\phi=6$.

From these figures and our results we find that with increasing of the ratio $m/M$ we have new effects that are absent in the massless case. First, we observe that we have the damped oscillation of the field $\phi$, whose frequency is increased with increasing the ratio $m/M$. Thus we may have a possible contribution of the field $\phi$ to the reheating as $m$ increases. This effect can be also observed for the massless case we studied in the previous section for larger values of the constant $\alpha$ (larger than 1), but for these values the viability of the inflationary model breaks down. Second, we see that in the beginning we have a valid increase of the field $\rho$ which first reaches a maximum value and then starts the slow-roll inflationary process. This maximum value is increased with increasing the ratio $m/M$. Specifically, in the case $m/M=1$ the field $\rho$ does not increase, in the case $m/M=10^{2}$ increases till the maximum value $7.6$ and in the case $m/M=10^{3}$ till $8.4$. Last, we notice that with increasing the ratio $m/M$ the duration of the inflationary period becomes longer.

From our results we find that the values for the observables in the case of including the mass term in the potential agree with the constraints (1.1), for values of $m/M$ ranging from $1$ to $10^{3}$, and thus also this case leads to viable inflation. For example for the case of $m/M=1$ the spectral index and the tensor-to-scalar ratio are $n_{s}=0.963\pm 0.004$ and $r=0.0037\mp 0.0007$, while in the case of $m/M=10^{3}$, they are $n_{s}=0.963\pm 0.003$ and $r=0.0037\mp 0.0007$.

5 The Case $U(\phi)=\frac{1}{4}\lambda(\phi^{2}-\upsilon^{2})^{2}$

In the following, we shortly study the case where $\phi$ is identified as the Standard Model Higgs boson. The scalar potential $V(\phi,\rho)$ in this case takes the form

[TABLE]

We focus on the case where $\lambda=m_{h}^{2}/2\upsilon^{2}\simeq 0.13$ is fixed by the measured Higgs vacuum expectation value $\upsilon\simeq 246\,GeV$ and Higgs boson mass $m_{h}\simeq 125\,GeV$ at the electroweak scale. The case $\beta=0$ was exhaustively studied in [84] where it is also being studied the case of smaller values of $\lambda$ that are suggested by the Standard Model RG flow, which drives the running coupling $\lambda$ to very small values at high energy scales. As discussed in [84] and [86] in the case where $\beta=0$ the constraints in the parameter $\alpha$ somehow relaxed, compared to the previous models, in order to be in agreement with all data. That is $\alpha\lesssim 10^{3}-10^{4}\,.$ Such values for the parameter $\alpha$ are also acceptable for $\beta>0$.

The potential (5.1) is positive defined and it contains a two-fold degenerated Minkowski minimum and a saddle-point, being a minimum for $\rho$ and a maximum for $\phi$, which respectively are given by

[TABLE]

At the extrema, the potential acquires the values

[TABLE]

From the above equations we observe that the places of the extrema of the potential and the values of the potential calculated there do not depend on the constant $\beta$ and thus are the same with the places found for the case $\beta=0$ studied within [84]. Furthermore, solving the equation $\frac{dV}{d\phi}=0\,,$ it is trivial to find the trajectory that the inflaton takes place in the potential valley, but it is too complicated to be presented. Nevertheless in the limit where $\beta\rightarrow 0$ and for negligibly small $\upsilon$ the known trajectory, $\phi^{2}=\frac{e^{\sqrt{\frac{2}{3}}\rho}-1}{\alpha+\lambda/3\alpha M^{2}}$, is recovered.

We study this model for the case where $\alpha=0.01$ and $\beta=0.001$ as in the previous sections and for different initial values for $\phi$, $\phi=6,\,1,\,0.6\,\,\text{and}\,\,10^{-3}$ and $\rho=7.4.$ The behaviour of the field $\rho$ does not change in this range of initial values of $\phi.$

In all cases the spectral index and the tensor-to-scalar ratio are $n_{s}=0.964\pm 0.004$ and $r=0.0037\mp 0.0007$, leading thus to viable inflation.

6 Conclusion

Inspired by the success of the Starobinsky model and the multifield nature of theories of particle physics we study the inflationary model proposed in [79]. In the begining, we review the feauture of a general $F({\cal R},\phi)$ theory of gravity, which is conformally equivalent to an Einstein-Hilbert theory including two scalar fields. Specializing to the case $F({\cal R},\phi)=-\frac{1}{2}f(\phi){\cal R}+\frac{1}{12M^{2}(\phi)}{\cal R}^{2}$ and $U(\phi)=0$, as considered in [79], we found that it leads to a two-field potential differing from the one given in [79]. However this also yields a viable inflationary model, as concluded therein.

Moreover, we considered other cases, as well, by adding a mass term for the field $\phi$ in the action. We found that this case also yields a viable inflationary model, while the evolution of the fields depends mostly on the ratio $m/M$. More specifically, we found that with increasing the ratio $m/M$ we have new effects that are absent in the massless case, such as the observation of damped oscillations of the field $\phi$, whose frequency is increased with increasing the ratio $m/M$. The ratio $m/M$ affects the duration of the inflationary period, too. Furthermore, in a Higgs-like potential, identifying the field $\phi$ with the Standard Model Higgs boson, we also found a viable inflation.

Finally, in the constraints we have provided for $\alpha$ and $\beta$ in the studied models we expect the non-gaussianities to be relatively small because $\phi$ is almost stabilized to its vacuum value before the last $50-60$ e-folds such that in these last e-folds the Starobinsky trajectory of $\rho$ drives the inflation. This observation also depends on the initial conditions for $\phi$ and $\rho$, but we checked that for initial conditions such that we have sufficient e-folds for inflation and the field $\phi$ contributes to the last $50-60$ of them, the predictions for $n_{s}$ were not in the range predicted by the observations.

Within the goals of future work, in the context of this model, and in order to check its viability, would be the study of reheating and preheating after the end of inflation, but this lies beyond the scope of this paper. Regarding possible extensions of this model these mainly include considerations of other physical scalar potentials, in the place of the $U(\phi)$ studied here and/or the insertion of more complicated non-minimal couplings between gravity and the pre-existing scalar field which are motivated by higher dimensional theories, such as non-minimal kinetic terms and d’ Alembertian of Ricci scalar terms, or Supergravity theories.

Acknowledgments

We are grateful to thank A.B. Lahanas, V.C. Spanos and N. Tetradis for useful discussions, comments and suggestions on the draft. D.C. thanks P. Christodoulidis for illuminating discussions. The research work was supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the “First Call for H.F.R.I. Research Projects to support Faculty members and Researchers and the procurement of high-cost research equipment grant” (Project Number: 824).