This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Reconstructing G-inflation: From the attractors nS(N)n_{S}(N) and r(N)r(N)

Ramón Herrera ramon.herrera@pucv.cl Instituto de Física, Pontificia Universidad Católica de Valparaíso, Avenida Brasil 2950, Casilla 4059, Valparaíso, Chile.
(September 3, 2025)
Abstract

The reconstruction of an inflationary universe in the context of the Galileon model or G-model, considering as attractors the scalar spectral index nS(N)n_{S}(N) and the tensor to scalar ratio r(N)r(N) as a function of the number of e-folding NN is studied. By assuming a coupling of the form G(ϕ,X)=g(ϕ)XG(\phi,X)=g(\phi)\,X, we obtain the effective potential VV and the coupling parameter gg in terms of the cosmological parameters nSn_{S} and rr under the slow roll approximation. From some examples for nS(N)n_{S}(N) and r(N)r(N), different results for the effective potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi) are found.

pacs:
98.80.Cq

I Introduction

It is well known that during the dynamic evolution of the early universe, it presented a period of rapid growth called inflationary stage or merely inflation R1 ; R102 ; R103 ; Rin . The inflationary universe gives an elegant solutions to long standing cosmological problems present in the standard hot big bang model. Nevertheless, inflation not only solves the problems of the hot big bang, but also gives account of the Large-Scale Structure (LSS) R2 ; R203 , together with a causal description of the anisotropies observed in the Cosmic Microwave Background (CMB) radiation of the early universeastro ; Hinshaw:2012aka ; Ade:2013zuv ; Planck2015 ; Ob2 ; DiValentino:2016foa .

In the context of the different models that give account of the dynamical evolution of an inflationary universe, we can stand out one general class of inflationary models where the inflation is driven by a minimally coupled scalar field. In the literature these models are called the Galilean inflationary models or simply G-inflation and its generalization, known as G2{}^{2}-inflation which also corresponds to a subclass of the Horndeski theoryHo , was developed in Ref.Kobayashi:2011nu . In this context, the inclusion of the canonical and non-canonical scalar field in the model of G-inflation, is known as kinetic gravity braiding modelG1 ; G2 . From the observational point the view, the detection of gravitational waves by GW170817 and the γ\gamma-ray burst TheLIGOScientific:2017qsa ; Monitor:2017mdv ; GBM:2017lvd give a strong constrain on the speed gravitational waves. In this context, the G-model is consistent with the GW170817, since the speed gravitational waves is equal to the speed of light. Thus, the Galilean action includes an extra term of the form G(ϕ,X)ϕG(\phi,X)\square\phi to the standard action, but this term on the action does not modify the speed gravitational waves and it is equivalent to the speed of light G1 ; G2 . We should also mention that the field equations still have derivatives only up to second order, see Ref.Nic . In this way, different cosmological models have been developed in the framework of G-inflation. In particular, assuming the slow roll approximation and considering some effective potentials G-inflation was studied in Ref.288 . In relation to the Higgs field, the model of Higgs G-inflation viewed as a modification of the standard Higgs inflation, in which the function G(ϕ,X)ϕXG(\phi,X)\propto\phi X, was developed in Ref.289a , see also289 . In this modified gravity, the specific case in which the effective potential V(ϕ)=V(\phi)= constant together with the slow roll approximation, is known as ultra slow roll G-inflation and it was studied in Ut . For the case in which the scalar potential is of the power-law type was studied in Ref.Pw . The model of warm inflation and its thermal fluctuations in the context of G-inflation was developed in Ref.Herrera:2017qux . The reheating mechanism in this model was studied in Ref.Br and from string gas cosmology in Ref.agr , see alsoReh2 ; DE1 ; DE2 ; Ka .

On the other hand, the idea of the reconstruction of the physical variables present on the background dynamics of inflationary models, from observational parameters such as the scalar spectrum, scalar spectral index nSn_{S} and the tensor to scalar ratio rr, have been studied by several authorsH1 ; H2 ; H3 ; H4 ; M ; Chiba:2015zpa ; H5 . In this context, a reconstruction mechanism to obtain the physical variables in the inflationary stage considering the slow roll approximation, is though of the parametrization of the cosmological observables nS(N)n_{S}(N) and r(N)r(N) or commonly called attractors, in which the parameter NN denotes the number of e-folds.

From the observational point of view, the scalar spectral index nSn_{S}, is well supported by the parametrization in terms of the number of e-fonds NN, given by the attractor nS12/Nn_{S}\sim 1-2/N for large NN, in which the number N5070N\simeq 50-70 at the end of the inflationary scenario from the results of Planck and BICEP2-Keck Array Collaborations Planck2015 ; Ob2 . In the framework of the General Relativity (GR), different models can be reconstructed considering the parametrization or attractor nS(N)12/Nn_{S}(N)\sim 1-2/N (for large N𝒪(102)N\sim\,\mathcal{O}(10^{2})), to name of few we have; the hyperbolic tangent model or simply the T-model T , E-modelE , R2R^{2}-modelR102 , the chaotic inflationary modelR103 , the Higgs inflation Higgs ; Higgs2 , etc. In the context of warm inflation and its reconstruction was necessary to introduce the attractors nS(N)n_{S}(N) and r(N)r(N) (unlike cold inflation), in order to build the scalar potential and the dissipation coefficient in terms of the scalar field Herrera:2018cgi .

Another methodology used in the literature in order to reconstruct the scalar potential, scalar spectral index and the tensor to scalar ratio is by means of the slow-roll parameter ϵ(N)\epsilon(N), as function of the number of e-folds NN Huang:2007qz ; M ; Gao:2017owg ; M . Also, the use of two slow roll parameters ϵ(N)\epsilon(N) and η(N)\eta(N), for the reconstruction of the scalar potential and spectral index was assumed in Ref.Roest:2013fha , see also Refs.N1 ; N2 .

The goal of this investigation is to reconstruct the Galilean inflationary model or G-inflation, given the parametrization of the scalar spectral index and the tensor to scalar ratio, in terms of the number of e-folds. In this sense, we investigate how the Galilean inflationary model in which the function G(ϕ,X)G(\phi,X), is given by G(ϕ,X)=g(ϕ)XG(\phi,X)=g(\phi)\,X, modifies the reconstructions of the scalar potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi), from the attractor point nS(N)n_{S}(N) and r(N)r(N). Thus, we will determine the structure of the function g(ϕ)g(\phi) and in particular of V(ϕ)V(\phi), in order to in account of the observables nSn_{S} and the ratio rr given by the observations.

By considering the domination of the Galilean effect, we develop a general formalism in order to obtain the effective potential VV and coupling parameter gg, from the parametrization of the cosmological attractors nS(N)n_{S}(N) and r(N)r(N), under the slow roll approximation.

For the application of the developed formalism, we will study different examples for the attractor point nS(N)r(N)n_{S}(N)-r(N). From these attractors, we will reconstruct the effective potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi) in terms of the scalar field ϕ\phi. Also, we will obtain different constraints on the parameters present in the reconstruction.

The outline of the paper is as follows: The next section we give a brief description of the model of G-inflation. Here, the background dynamics and cosmological perturbations are presented. In the section III, we develop a general formalism in order to reconstruct the scalar potential and coupling parameter in function of the attractors nS(N)n_{S}(N) and r(N)r(N), respectively. In section IV we apply the method for different examples of nS(N)n_{S}(N) and r(N)r(N) so as to construct the effective potential V(ϕ)V(\phi) and the coupling g(ϕ)g(\phi) in terms of the scalar field ϕ\phi. Finally, in section V we give our conclusions. We chose units so that c==8π=1c=\hbar=8\pi=1.

II The model of G-inflation

In this section we give a brief description of the model of G-inflation. We start with the 4-dimensional action for the Galilean model given by

S=g4d4x(MP22R+K(ϕ,X)G(ϕ,X)ϕ),S=\int\sqrt{-g_{4}}d^{4}x\,\left(\frac{M^{2}_{P}}{2}R+K(\phi,X)-G(\phi,X)\square\phi\right), (1)

where g4g_{4} denotes the determinant of the space-time metric gμνg_{\mu\nu}, RR is the Ricci scalar and the quantity XX corresponds to X=gμνμϕνϕ/2X=-g^{\mu\nu}\partial_{\mu}\phi\partial_{\nu}\phi/2, where ϕ\phi denotes the scalar field. The quantities KK and GG are arbitrary functions of XX and the scalar field ϕ\phi, respectively. Here, the quantity MpM_{p} corresponds to the Planck mass.

By assuming a spatially flat Friedmann Robertson Walker (FRW) metric, along with a scalar field homogeneous in which ϕ=ϕ(t)\phi=\phi(t), the Friedmann equation can be written as

3H2=κρ.3H^{2}=\kappa\,\rho. (2)

Here, the parameter H=a˙aH=\frac{\dot{a}}{a} denotes the Hubble rate, aa corresponds to the scale factor and ρ\rho is the energy density. In the following, the dots denote differentiation with respect to the time and the quantity κ=1/Mp2\kappa=1/M_{p}^{2}.

From the action (1), we can identify that the energy density and the pressure associated to the scalar field ϕ\phi are given by G1 ; G2

ρ=2KXXK+3GXHϕ˙32GϕX,\rho=2K_{X}\,X-K+3G_{X}H\dot{\phi}^{3}-2G_{\phi}X, (3)

and

p=K2(Gϕ+GXϕ¨)X,\displaystyle p=K-2(G_{\phi}+G_{X}\ddot{\phi})X, (4)

respectively. In the following, we will assume that the notation KXK_{X} denotes KX=K/XK_{X}=\partial K/\partial X, GϕG_{\phi} corresponds to Gϕ=G/ϕG_{\phi}=\partial G/\partial\phi, KXX=2K/X2K_{XX}=\partial^{2}K/\partial X^{2}, etc.

In this framework, the continuity equation for the energy density ρ\rho can be written as ρ˙+3H(ρ+p)=0,\dot{\rho}+3\,H\,(\rho+p)=0, or equivalently

KXϕ+2KXXXϕ¨+2KXϕXKϕ2(GϕGXϕX)ϕ+6GX(H˙X+X˙H+3H2X)K_{X}\square\phi+2K_{XX}X\ddot{\phi}+2K_{X\phi}X-K_{\phi}-2(G_{\phi}-G_{X\phi}X)\square\phi+6G_{X}(\dot{H}X+\dot{X}H+3H^{2}X)
+6HGXXXX˙2GϕϕX4GXϕXϕ¨=0,whereϕ=ϕ¨+3Hϕ˙.+6HG_{XX}X\dot{X}-2G_{\phi\phi}X-4G_{X\phi}X\ddot{\phi}=0,\;\;\;\;\mbox{where}\,\,\,\,\;\;\square\phi=\ddot{\phi}+3H\dot{\phi}. (5)

Here, we have used Eqs.(3), and (4).

In particular, for the special cases in which K=XV(ϕ)K=X-V(\phi) and G=0G=0, in which V(ϕ)V(\phi) corresponds to the effective potential, we recovered the standard General Relativity (GR).

In order to study the reconstruction for the G-model, we will analyze the special case in which the functions K(ϕ,X)K(\phi,X) and G(ϕ,X)G(\phi,X) are given by

K(ϕ,X)=XV(ϕ),andG(ϕ,X)=g(ϕ)X,K(\phi,X)=X-V(\phi),\;\;\;\;\mbox{and}\;\;\;G(\phi,X)=-g(\phi)\,X, (6)

respectively. Here, g(ϕ)g(\phi) is a function that depends exclusively on the scalar field ϕ\phi.

From Refs.G1 ; G2 , we will assume the slow roll approximation. In this context, the potential V(ϕ)V(\phi) dominates over the quantities XX, |GXHϕ˙3||G_{X}H\dot{\phi}^{3}| and |GϕX||G_{\phi}X|, wherewith the energy density ρϕV(ϕ)\rho_{\phi}\sim V(\phi) and then the Friedmann equation (2) becomes

3H2κV(ϕ).3H^{2}\approx\,\kappa\,V(\phi). (7)

Introducing the slow-roll parameters ϵ1\epsilon_{1}, ϵ2\epsilon_{2}, ϵ3\epsilon_{3} and ϵ4\epsilon_{4} given by G2

ϵ1=H˙H2,ϵ2=ϕ¨Hϕ˙,ϵ3=gϕϕ˙gH,andϵ4=gϕϕX2Vϕ,\epsilon_{1}=-\frac{\dot{H}}{H^{2}},\,\,\,\epsilon_{2}=-\frac{\ddot{\phi}}{H\dot{\phi}},\,\,\,\epsilon_{3}=\frac{g_{\phi}\dot{\phi}}{gH},\,\,\mbox{and}\,\,\,\epsilon_{4}=\frac{g_{\phi\phi}X^{2}}{V_{\phi}}, (8)

then the Eq.(5) can be rewritten as

3Hϕ˙(1ϵ2/3gHϕ˙[3ϵ12ϵ22ϵ2ϵ3/3])=(12ϵ4)Vϕ.3H\dot{\phi}(1-\epsilon_{2}/3-gH\dot{\phi}[3-\epsilon_{1}-2\epsilon_{2}-2\epsilon_{2}\epsilon_{3}/3])=-(1-2\epsilon_{4})V_{\phi}. (9)

Here, we have used the functions given by Eq.(6) and the slow roll parameters given by Eq.(8).

By considering that the slow-roll parameters ϵ1\epsilon_{1}, |ϵ2||\epsilon_{2}|, |ϵ3||\epsilon_{3}|, |ϵ4|1|\epsilon_{4}|\ll 1, then the Eq.(9) is reduced to

3Hϕ˙(1+𝒜)Vϕ,where𝒜=3gHϕ˙.3H\dot{\phi}(1+{\cal{A}})\simeq-V_{\phi},\;\;\;\;\mbox{where}\;\;\;{\cal{A}}=-3gH\dot{\phi}. (10)

Here, we mention that in relation to the slow roll equation (10), we have two limiting instances. The situation in which |𝒜|1|{\cal{A}}|\ll 1 corresponds to the standard equations of GR in the framework of slow roll inflation. Instead, the inverse case in which |𝒜|1|{\cal{A}}|\gg 1, the Galileon effect changes the dynamic equation of the scalar field ϕ\phi and hence the dynamics of inflationary model. In this context, we are interested in the latter situation in which the Galileon effect modifies the dynamics of the G-model and its reconstruction. Thus, 3Hϕ˙𝒜Vϕ3H\dot{\phi}{\cal{A}}\simeq-V_{\phi} and then 9H2ϕ˙2(Vϕ/g)9H^{2}\dot{\phi}^{2}\simeq(V_{\phi}/g) suggesting that the ratio (Vϕ/g)>0(V_{\phi}/g)>0. Therefore, in the case in which Vϕ>0V_{\phi}>0 then the quantity g>0g>0 and vice versa. In the following we shall take Vϕ>0V_{\phi}>0 and g>0g>0.

Typically, if the scalar field roll down potential, then the velocity of the scalar field can be written as

ϕ˙Vϕ3κgV.\dot{\phi}\simeq-\sqrt{\frac{V_{\phi}}{3\kappa\,g\,V}}. (11)

Here, we have considered Eq.(7). Also, we note that the parameter 𝒜>0{\cal{A}}>0, since we have assumed that ϕ˙<0\dot{\phi}<0.

In relation to the expansion, we define the number of e-folding NN in order to give a measure of the inflationary expansion. In this way, we assume two different values of cosmological times tt and tet_{e}. Here, the time tet_{e} denotes the end of inflationary epoch. Thus, the number of e-folds NN in the limit in which 1|𝒜|1\ll|{\cal{A}}| is given by

N=tteH𝑑t=ϕϕeHdϕϕ˙κϕeϕV[gVϕ]1/2𝑑ϕ.N=\int_{t}^{t_{e}}H\,dt^{\prime}=\int_{\phi}^{\phi_{e}}H\,\frac{d\phi^{\prime}}{\dot{\phi}}\simeq\kappa\,\int_{\phi_{e}}^{\phi}\,V\,\left[\frac{g}{V_{\phi^{\prime}}}\right]^{1/2}\,d\phi^{\prime}. (12)

On the other hand, the cosmological perturbations together with the scalar and tensor spectrums were obtained in Refs.Kobayashi:2011nu ; GG1 ; GG2 ; Ohashi:2012wf for the model of G-inflation. In this sense, from the action (1), the amplitude of scalar perturbations 𝒫𝒮{\cal{P_{S}}} generated during the inflationary epoch for a flat space and assuming the slow roll approximation we have

𝒫𝒮κ3V312π2Vϕ2(1+𝒜)2(1+2𝒜)1/2(1+4𝒜/3)3/2.{\cal{P_{S}}}\simeq\frac{\kappa^{3}\,V^{3}}{12\pi^{2}\,V_{\phi}^{2}}\,\frac{(1+{\cal{A}})^{2}(1+2{\cal{A}})^{1/2}}{(1+4{\cal{A}}/3)^{3/2}}. (13)

Because the scalar spectral index nSn_{S} is defined as nS=dln𝒫𝒮/dlnkn_{S}=d\ln{\cal{P_{S}}}/d\ln k, then from Eq.(13) the index nSn_{S} in terms of the standard slow roll parameters can be written as

nS16ϵ1+𝒜+2η1+4𝒜/3[1𝒜6(1+2𝒜)2],n_{S}-1\simeq\,-\frac{6\epsilon}{1+{\cal{A}}}+\frac{2\eta}{1+4{\cal{A}}/3}\left[1-\frac{{\cal{A}}}{6(1+2{\cal{A}})^{2}}\right], (14)

where the standard parameters ϵ\epsilon and η\eta are defined as

ϵ=12κ(VϕV)2,andη=VϕϕκV.\epsilon=\frac{1}{2\kappa}\left(\frac{V_{\phi}}{V}\right)^{2},\;\,\,\,\,\,\mbox{and}\,\,\,\eta=\frac{V_{\phi\phi}}{\kappa\,V}. (15)

Note that in the limit g0g\rightarrow 0 (or equivalently 𝒜0{\cal{A}}\rightarrow 0), the spectral index given by Eq.(14) coincides with the GR in which nS16ϵ+2ηn_{S}-1\simeq-6\epsilon+2\eta. By assuming the limit in which |𝒜|1|{\cal{A}}|\gg 1, the scalar spectral index reduces to

nS16ϵ𝒜+3η2𝒜.n_{S}-1\simeq-\frac{6\epsilon}{{\cal{A}}}+\frac{3\eta}{2\cal{A}}. (16)

On the other hand, in relation to the tensorial perturbations, the amplitude of the tensor mode was determined in refs.Kobayashi:2011nu ; GG1 ; GG2 ; Ohashi:2012wf , and then the tensor spectrum 𝒫𝒯{\cal{P_{T}}} is defined as 𝒫𝒯(2κ2V/3π2){\cal{P_{T}}}\simeq(2\kappa^{2}V/3\pi^{2}). In this form, the tensor to scalar ratio rr in the framework of G-inflation can be written as

r=𝒫𝒯𝒫𝒮 16ϵ[(1+4𝒜/3)3/2(1+𝒜)2(1+2𝒜)1/2].r=\frac{{\cal{P_{T}}}}{{\cal{P_{S}}}}\simeq\,16\epsilon\,\left[\frac{(1+4{\cal{A}}/3)^{3/2}}{(1+{\cal{A}})^{2}(1+2{\cal{A}})^{1/2}}\right]. (17)

One again, note that in the limit 𝒜0{\cal{A}}\rightarrow 0 ( or equivalently g0g\rightarrow 0), the ratio rr coincides with that corresponding to GR in which r=16ϵr=16\epsilon.

Taking the limit |𝒜|1|{\cal{A}}|\gg 1, the tensor to scalar ratio rr results

r4233/216ϵ𝒜.r\simeq\frac{4\sqrt{2}}{3^{3/2}}\,\frac{16\epsilon}{\cal{A}}. (18)

In the following we will analyze the reconstruction of the effective potential V(ϕ)V(\phi) and the coefficient g(ϕ)g(\phi) in the framework of G- inflation. In order to realize the reconstruction we will assume the limit 1|𝒜|1\ll|{\cal{A}}|, together with an attractor point from the index nS(N)n_{S}(N) and the ratio r(N)r(N) on the rnSr-n_{S} plane.

III Reconstructing G-inflation

In this section we develop the method to follow in order to reconstruct, the scalar potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi), assuming the scalar spectral index nS(N)n_{S}(N) and the tensor to scalar ratio r(N)r(N) as attractors. In order to reconstruct analytically the potential VV and the coupling parameter gg, we shall take the limit |𝒜|1|{\cal{A}}|\gg 1. Following Refs.Chiba:2015zpa ; Herrera:2018cgi , we rewrite the spectral index and the tensor to scalar ratio given by Eqs.(16) and (18), in terms of the number of e-folds NN and its derivatives. Thus, from these relations and giving nS=nS(N)n_{S}=n_{S}(N) and r=r(N)r=r(N), we should find the potential VV and the coupling parameter gg as a function of the number NN. Later, with the help of Eq.(12) we should obtain the e-folding NN in terms of the scalar field ϕ\phi in order to reconstruct finally, the scalar potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi), respectively.

In this sense, we begin by rewriting the standard slow roll parameters as a function of the number of e-folds, taking into account that

Vϕ=dVdϕ=κVgVϕVN,V_{\phi}=\frac{dV}{d\phi}=\kappa V\sqrt{\frac{g}{V_{\phi}}}\,V_{\,N}\,,

then we get

Vϕ=(κVgVN)2/3,wherewithVN=dVdN>0.V_{\phi}=(\kappa V\sqrt{g}\,\,V_{\,N})^{2/3}\,,\,\,\,\,\mbox{wherewith}\,\,\,\,V_{\,N}=\frac{dV}{dN}>0. (19)

In the following, we will assume the subscript VN=dV/dNV_{\,N}=dV/dN, VNNV_{NN} to VNN=d2V/dN2V_{NN}=d^{2}V/dN^{2}, gN=dg/dNg_{{}_{N}}=dg/dN etc.

Similarly for VϕϕV_{\phi\phi} we get

Vϕϕ=23(κ4V4g2VN)1/3[VNV+(g)Ng+VNNVN].V_{\phi\phi}=\frac{2}{3}\,(\kappa^{4}\,V^{4}\,g^{2}\,V_{N})^{1/3}\,\left[\frac{V_{N}}{V}+\frac{(\sqrt{g})_{N}}{\sqrt{g}}+\frac{V_{\,NN}}{V_{\,N}}\right]. (20)

In this form, the standard slow roll parameters ϵ\epsilon and η\eta are rewritten as

ϵ=12κV2(κVgVN)4/3,\epsilon=\frac{1}{2\kappa V^{2}}\,(\kappa V\sqrt{g}\,V_{\,N})^{4/3}, (21)

and

η=23κV(κ4V4g2VN)1/3[VNV+(g)Ng+VNNVN],\eta=\frac{2}{3\kappa V}\,(\kappa^{4}\,V^{4}\,g^{2}\,V_{\,N})^{1/3}\,\left[\frac{V_{\,N}}{V}+\frac{(\sqrt{g})_{N}}{\sqrt{g}}+\frac{V_{\,NN}}{V_{\,N}}\right], (22)

respectively.

Now, the relationship between the e-folds NN and the scalar field ϕ\phi, from Eq.(12) can be written as

[VNκ2gV2]1/3𝑑N=𝑑ϕ.\int\left[\frac{V_{\,N}}{\kappa^{2}\,g\,V^{2}}\right]^{1/3}\,dN=\int\,d\phi. (23)

In this context, by considering Eqs.(16), (21) and (22) we find that the scalar spectral index nSn_{S} can be rewritten in terms of the e-folding NN, such that

nS1=2VNV+[(g)Ng+VNNVN],n_{S}-1=-2\frac{V_{\,N}}{V}+\left[\frac{(\sqrt{g})_{N}}{\sqrt{g}}+\frac{V_{\,NN}}{V_{\,N}}\right],

or equivalently

nS1=[ln(gVNV2)]N.n_{S}-1=\left[\ln\left(\frac{\sqrt{g}\,\,V_{\,N}}{V^{2}}\right)\right]_{\,N}. (24)

From Eq.(18) we obtain that the tensor to scalar ratio becomes

r=𝒫𝒯𝒫𝒮=CVNV,in whichC=3269 8.71.r=\frac{{\cal{P_{T}}}}{{\cal{P_{S}}}}=C\,\frac{V_{\,N}}{V},\,\,\,\,\mbox{in which}\,\,\,\;\;\,C=\frac{32\sqrt{6}}{9}\simeq\,8.71. (25)

Here, we have considered that the function 𝒜\cal{A} can be rewritten as 𝒜=3gHϕ˙=g(κVgVN)1/3.{\cal{A}}=-3gH\dot{\phi}=\sqrt{g}\,(\kappa\,V\,\sqrt{g}\,V_{\,N})^{1/3}. Also, we note that in the context of the reconstruction and assuming 𝒜1{\cal A}\gg 1, the ratio rr given by Eq.(25) does not depend of the parameter gg. Here, we noted that this ratio rr is similar to the one obtained in the standard GR, where r=8VN/Vr=8V_{N}/V, see Ref.Chiba:2015zpa .

Now, from Eqs.(24) and (25) we obtain that the effective potential in terms of the e-folding NN results

V(N)=exp[1Cr𝑑N],V(N)=\,\exp\left[\frac{1}{C}\int\,r\,dN\right], (26)

and the coupling parameter g(N)g(N) becomes

g(N)=C2r2exp[2(rC+[nS1])𝑑N].g(N)=\frac{C^{2}}{r^{2}}\,\exp\left[2\int\,\left(\frac{r}{C}+[n_{S}-1]\right)dN\right]. (27)

Here, we emphasize that the potential V(N)V(N) depends only of the tensor to scalar ratio r(N)r(N).

In fact, the Eqs.(23), (26) and (27) are the fundamental equations in order to reconstruct of the effective potential V(ϕ)V(\phi) and the parameter g(ϕ)g(\phi), giving the attractors nS(N)n_{S}(N) and r(N)r(N), respectively.

In the following we will study some specific examples in order to reconstruct the scalar potential V(ϕ)V(\phi) and coupling parameter g(ϕ)g(\phi), from the cosmological parameters or attractors nS(N)n_{S}(N) and r(N)r(N), respectively.

IV Some examples of reconstruction

In order to apply the formalism of above, we shall first consider the simplest example for the attractors nS(N)n_{S}(N) and r(N)r(N), so as to reconstruct analytically the effective potential V(ϕ)V(\phi) and coupling parameter g(ϕ)g(\phi). In this sense and following Refs.T ; Chiba:2015zpa , we assume that the spectral index in terms of the number of e-folds is given by

nS1=2N,n_{S}-1=-\frac{2}{N}, (28)

and the tensor to scalar ratio as

r=αN,r=\frac{\alpha}{N}, (29)

where α>0\alpha>0 corresponds to a dimensionless constant. We mention that in the framework of GR, for the chaotic model (in which V(ϕ)ϕ2V(\phi)\propto\phi^{2}) R103 , the parametrization in terms of NN of the scalar index is given by Eq.(28), where the value of the parameter α=8\alpha=8, i.e., r(N)=8/Nr(N)=8/N, see Ref.Chiba:2015zpa . In particular if we consider that the number NN before the end of inflationary scenario occurs at N60N\simeq 60, then the tensor to scalar ratio given by Eq.(29) is well corroborated by observational data when the constant α<4.2\alpha<4.2. Here, we have used that the ratio r<0.07r<0.07 from BICEP2 and Keck Array CollaborationsOb2 .

From the attractor given by Eq.(29), we find that the effective potential (26) results

V(N)=V0NαC,V(N)=V_{0}\,N^{\frac{\alpha}{C}}, (30)

where V0>0V_{0}>0 corresponds to the integration constant (with units of Mp4M_{p}^{4}). By utilizing the spectral index (28), we obtain that the integral expression exp[(ns1)]dN=β/N2\exp[\int(n_{s}-1)]dN=\beta/N^{2}, in which β\beta denotes an integration constant (with units of Mp11/2M_{p}^{-11/2} or κ11/4\kappa^{11/4}). In this form, we obtain that the coupling parameter g(N)g(N) in terms of e-folds from Eq.(27) can be written as

g(N)=g0N2(Cα)C,whereg0=(βCV0α)2.g(N)=g_{0}\,N^{\frac{-2(C-\alpha)}{C}},\,\,\,\,\,\mbox{where}\,\,\,\,\,g_{0}=\left(\frac{\beta CV_{0}}{\alpha}\right)^{2}. (31)

By using Eqs.(30) and (31), we find that the parameter 𝒜\cal{A} in terms of the number of e-folds results

𝒜(N)=𝒜0N(5C6α)3C,in which𝒜0=[CV02κ1/3β4/3α].{\cal{A}}(N)={\cal{A}}_{0}\,N^{\frac{-(5C-6\alpha)}{3C}}\,,\;\;\;\;\mbox{in which}\;\;\;\;\;{\cal{A}}_{0}=\left[\frac{C\,V_{0}^{2}\,\kappa^{1/3}\beta^{4/3}}{\alpha}\right]. (32)

From the condition in which predominate the Galileon effect, such that 𝒜1{\cal{A}}\gg 1, then from Eq.(32) we find a lower limit for the parameter β\beta given by βα3/4N5C6α4C(CV02κ1/3)3/4\beta\gg\frac{\alpha^{3/4}N^{\frac{5C-6\alpha}{4C}}}{(CV_{0}^{2}\kappa^{1/3})^{3/4}}. In particular for N=60N=60, α=4\alpha=4 and assuming that the potential at the end of inflation V(N=60)1011Mp4V(N=60)\simeq 10^{-11}M_{p}^{4} (in which V01.5×1012Mp4V_{0}\simeq 1.5\times 10^{-12}M_{p}^{4}), we find that the lower bound for the parameter β\beta is given by βMp11/22.9×1018𝒪(1018)\beta\,M_{p}^{11/2}\simeq 2.9\times 10^{18}\sim\,\mathcal{O}(10^{18}).

Now, combing Eqs.(23), (30) and (31), we find that the relationship between the e-folding NN and the scalar field ϕ\phi results

N=N0(ϕϕ0)γ1,N=N_{0}\,(\phi-\phi_{0})^{\gamma_{1}}, (33)

where ϕ0\phi_{0} corresponds to an integration constant and the quantities N0N_{0} and γ1\gamma_{1} are given by

N0=((4C3α)3C[κ2Cg0V0α]1/3)γ1,andγ1=3C4C3α,N_{0}=\left(\frac{(4C-3\alpha)}{3C}\,\left[\frac{\kappa^{2}Cg_{0}V_{0}}{\alpha}\right]^{1/3}\right)^{\gamma_{1}},\;\;\;\;\;\mbox{and}\,\,\,\,\,\gamma_{1}=\frac{3C}{4C-3\alpha}\,,

respectively.

In this form, we find that the reconstruction of the effective potential V(ϕ)V(\phi) can be written as

V(ϕ)=V¯0(ϕϕ0)αγ1C,V(\phi)=\bar{V}_{0}\,(\phi-\phi_{0})^{\frac{\alpha\,\gamma_{1}}{C}}, (34)

where the quantity V¯0\bar{V}_{0} is defined as V¯0=V0N0\bar{V}_{0}=V_{0}\,N_{0}. Here, we note that the effective potential V(ϕ)V(\phi) corresponds to a power law potential in which the exponent αγ1C>0\frac{\alpha\,\gamma_{1}}{C}>0. In particular for case in which the parameter α=4\alpha=4, we find that the effective potential V(ϕ)ϕ0.52V(\phi)\propto\phi^{0.52}.

Analogously, from Eqs.(31) and (33), we obtain that the coupling parameter gg in terms of the scalar field results

g(ϕ)=g¯0(ϕϕ0)γ2=g¯01(ϕϕ0)γ2,g(\phi)=\bar{g}_{0}\,(\phi-\phi_{0})^{-\,\gamma_{2}}=\bar{g}_{0}\,\frac{1}{(\phi-\phi_{0})^{\gamma_{2}}}, (35)

where the constant γ2=6(Cα)4C3α>0\gamma_{2}=\frac{6(C-\alpha)}{4C-3\alpha}>0 and the quantity g¯0=g0N0\bar{g}_{0}=g_{0}\,N_{0}. We noted that in particular if the parameter α=4\alpha=4, then the power γ21.2\gamma_{2}\simeq 1.2 and the coupling parameter gg decays as g(ϕ)ϕ1.2g(\phi)\propto\phi^{-1.2}. Also, we noted that for values of α0\alpha\simeq 0 (or equivalently r0r\simeq 0) and the power γ23/2\gamma_{2}\simeq 3/2, then the coupling parameter decays as g(ϕ)ϕ3/2g(\phi)\propto\phi^{-3/2}. Also, we note that for the value γ23/2\gamma_{2}\simeq 3/2 (or equivalently α0\alpha\simeq 0), the scalar potential V(ϕ)V(\phi)\sim constant, leading to an exponential expansion R1 . Thus, we find that the range for the parameter γ2\gamma_{2} is given by 1.1<γ23/21.1<\gamma_{2}\lesssim 3/2.

In order to give an account of the end of inflationary stage in this reconstruction, we consider that inflation ends when the slow-roll parameter ϵ1(ϕ=ϕe)=H˙/H2=1\epsilon_{1}(\phi=\phi_{e})=-\dot{H}/H^{2}=1 or equivalently a¨=0\ddot{a}=0. Under slow roll approximation and considering the limit 𝒜1{\cal{A}}\gg 1, we write ϵ1\epsilon_{1} in terms of standard slow roll parameter ϵ\epsilon as ϵ1ϵ(gVϕ)1/2\epsilon_{1}\simeq\frac{\epsilon}{(g\,V_{\phi})^{1/2}}. In this way, combining Eqs.(15), (34) and (35), we obtain that the scalar field at the end of inflation becomes

ϕe=ϕ0+ϕ0~,whereϕ0~=[12κ(g0¯V0¯)1/2(αγ1C)3/2]23+αγ1/Cγ2>0.\phi_{e}=\phi_{0}+\tilde{\phi_{0}},\,\,\,\,\,\mbox{where}\,\;\;\;\,\,\tilde{\phi_{0}}=\left[\frac{1}{2\kappa\,(\bar{g_{0}}\bar{V_{0}})^{1/2}}\,\left(\frac{\alpha\gamma_{1}}{C}\right)^{3/2}\right]^{\frac{2}{3+\alpha\gamma_{1}/C-\gamma_{2}}}>0.

Also, the condition for that inflation takes place is ϵ1<1\epsilon_{1}<1 (or equivalently a¨>0\ddot{a}>0). Thus, during inflation the scalar field is such that ϕ>ϕ0+ϕ0~\phi>\phi_{0}+\tilde{\phi_{0}}. This results suggest that the coupling parameter g(ϕ)g(\phi) given by Eq.(35) does not present a singularity at ϕ=ϕ0\phi=\phi_{0}, during the inflationary scenario. Thus, at the end of inflation g(ϕe)=g0¯/ϕ0~γ2g(\phi_{e})=\bar{g_{0}}/\tilde{\phi_{0}}^{\gamma_{2}} and the scalar potential V(ϕe)=V0¯ϕ0~αγ1/CV(\phi_{e})=\bar{V_{0}}\,\tilde{\phi_{0}}^{\alpha\gamma_{1}/C}.

Other type of the attractor for the ratio rr studied in the literature is given by r=1/N(1+ξN)r=1/N(1+\xi N), where ξ\xi is a free parameterChiba:2015zpa ; Herrera:2018cgi , such that ξ>4/315\xi>-4/315 in order to obtain from BICEP2 and Keck Array Collaborations r<0.07r<0.07, see Ref.Ob2 . By considering this attract together with nsn_{s} given by Eq.(28), we find a transcendental equation from Eq.(23) to express the number NN in terms of the field ϕ\phi, wherewith the reconstruction does not work.

In order to find a simple relationship between the number of e-folding NN and the scalar field, we consider Eq.(28) together with the attractor r(N)r(N) given by

r=CN(3+ξN1/3),r=\frac{C}{N(3+\xi N^{1/3})}, (36)

where ξ\xi corresponds to a constant (dimensionless) and it satisfies the lower bound ξ>0.24\xi>-0.24, such that in particular r(N=60)<0.07r(N=60)<0.07. Note that the relation between the ratio rr and the scalar spectral index nSn_{S}, can be written as

r(ns)=C(1nS)4/32[3(1nS)1/3+21/3ξ].r(n_{s})=\frac{C(1-n_{S})^{4/3}}{2[3(1-n_{S})^{1/3}+2^{1/3}\xi]}. (37)

Here, we have used Eqs.(28) and (36), respectively.

From Eq.(26) we find that the effective potential V(N)V(N) becomes

V(N)=C1N1/33+ξN1/3,V(N)=\frac{C_{1}\,N^{1/3}}{3+\xi N^{1/3}}, (38)

where C1>0C_{1}>0 is an integration constant (with units of MP4M_{P}^{4}). Note that for the case in which the parameter ξ1/N1/3\xi\gg 1/N^{1/3}, the scalar potential V(N)C1/ξ=V(N)\sim C_{1}/\xi= constant. By using Eq.(27) we obtain that the coupling parameter gg in terms of the number of e-folds NN results

g(N)=g0N4/3,whereg0=β2C12.g(N)=\frac{g_{0}}{N^{4/3}},\;\;\;\;\mbox{where}\,\,\,\,\,g_{0}=\beta^{2}\,C_{1}^{2}. (39)

Here, we note that the coupling g(N)g(N) does not depend of the constant ξ\xi, only of the integration constants β\beta and C1C_{1}, respectively. From the condition in which predominate the Galileon effect, where 𝒜1{\cal{A}}\gg 1, we find a lower bound for the integration constants β\beta and C1C_{1} given by

g0C1κ(27N3+27ξN10/3+9ξ2N11/3+ξ3N4)1/2.g_{0}\,C_{1}\sqrt{\kappa}\gg(27N^{3}+27\xi N^{10/3}+9\xi^{2}N^{11/3}+\xi^{3}N^{4})^{1/2}. (40)

In particular for N=60N=60 and ξ=0.2\xi=-0.2 the lower limit gives g0C1κ1534g_{0}C_{1}\sqrt{\kappa}\gg 1534 and for the value ξ=0.2\xi=0.2 corresponds to g0C1κ3419g_{0}C_{1}\sqrt{\kappa}\gg 3419. Also, in the special case in which the potential at the end of inflation V(N=60)=1011Mp4V(N=60)=10^{-11}M_{p}^{4}, together with ξ=0.2\xi=-0.2 we have C15×1012Mp4C_{1}\simeq 5\times 10^{-12}M_{p}^{4}, then we find that the lower limit for the integration constant β\beta is given by βMp11/2𝒪(1018)\beta M_{p}^{11/2}\gg\,\mathcal{O}(10^{18}). For the case α=0.2\alpha=0.2, we have C19.7×1012Mp4C_{1}\simeq 9.7\times 10^{-12}M_{p}^{4} and then βMp11/22×1018𝒪(1018)\beta M_{p}^{11/2}\gg 2\times 10^{18}\sim\,\mathcal{O}(10^{18}).

On the other hand, from Eq.(23) is easy to find that the relation between NN and ϕ\phi is given by lineal equation i.e., NϕN\propto\phi. Thus, the reconstruction of the coupling parameter gg as a function of the scalar field is given by g(ϕ)ϕ4/3g(\phi)\propto\phi^{-4/3}. This result suggests that again the parameter g(ϕ)g(\phi) has a behavior power law with a negative power. In this point, we mention that in order to obtain the lineal relation NϕN\propto\phi, we establish from Eq.(23) the condition (VN/(κ2gV2))1/3=(V_{N}/(\kappa^{2}gV^{2}))^{1/3}= constant. Thus, combing Eqs.(24) and (25) together with the attractor nSn_{S} given by Eq.(28), we obtained the tensor to scalar ratio rr given by Eq.(36) (or the potential Eq.(38)).

This methodology can be used for any function F(N)F(N), such that (VN/(κ2gV2))1/3=F(N)(V_{N}/(\kappa^{2}gV^{2}))^{1/3}=F(N). In this sense, the Eq.(23) takes of form F(N)𝑑N=𝑑ϕ\int F(N)dN=\int d\phi, being possible to choose any function F(N)F(N) in order to obtain an analytical and invertible solution for N=N(ϕ)N=N(\phi). Subsequently, we solve the differential equation for the variable r(N)r(N) (or V(N)V(N) or also g(N)g(N)), by combining Eqs.(24), (25) and (28) for a specific function F(N)F(N). In particular, for the case in which the function F(N)=F(N)=constant, we have NϕN\propto\phi, and then we found that the tensor to scalar ratio r(N)r(N) is given by Eq.(36), and hence the form of the potential V(N)V(N) and the coupling g(N)g(N) correspond to the equations (38) and (39), respectively.

Refer to caption
Refer to caption
Figure 1: The evolution of the ratio g(N)/g0g(N)/g_{0} versus the number of e-folding NN (left panel) from Eq.(39). In the right panel, the plot shows the tensor to scalar ratio rr versus the spectral index nSn_{S}, from Eq.(37) for three values of the parameter ξ\xi. Here, we show the two-marginalized constraints on the curve r=r(nS)r=r(n_{S}) and as the BICEP2 and Keck Array Collaborations data places stronger limits on the rnSr-n_{S} plane Ob2 .

In Fig.1 we show the evolution of the ratio g(N)/g0g(N)/g_{0} on the number of e-folds NN (left panel) from Eq.(39) and as it decays in terms of the number of e-folds NN (or equivalently ϕ\phi, since NϕN\propto\phi). In the right panel, we show the dependence of the inflationary parameters r=r(nS)r=r(n_{S}) from Eq.(37). Here, we show the two-dimensional marginalized constraints (68%\% and 95%\% C.L.) on the tensor to scalar ratio r=r(ns)r=r(n_{s}), derived from observational dataOb2 . In the plot we have used Eq.(37), together with three different values of the parameter ξ\xi. We note that for values of ξ>0.241\xi>-0.241 the relation given by Eq.(37) is well supported by the BICEP2 and Keck Array Collaborations in which r<0.07r<0.07 and nS0.967n_{S}\simeq 0.967 Ob2 .

V Conclusions

In this paper we have investigated the reconstruction in the model of G-inflation, from the cosmological parameters such as the scalar spectral index nS(N)n_{S}(N) and the tensor to scalar ratio r(N)r(N), in which NN denotes the number of e-folding. By assuming the domination of the Galilean term and under the slow roll approximation, we have developed a general formalism of reconstruction in which the function G(ϕ,X)=g(ϕ)XG(\phi,X)=g(\phi)\,X .

Under this general analysis we have found from the attractor point nS(N)n_{S}(N) and r(N)r(N), integrable expressions for the effective potential V(N)V(N) and the coupling parameter g(N)g(N), respectively. Curiously, in the context of the reconstruction we have found that the tensor to scalar ratio rr is similar to the one obtained in the GR, in which rr depends only of the potential and its derivative, and it does not depend of the coupling gg.

As the simplest example and in order to find the reconstructions of the potential V(ϕ)V(\phi) and the coupling parameter g(ϕ)g(\phi), we have assumed the standard attractors given by nS1=2/Nn_{S}-1=-2/N and r=α/Nr=\alpha/N, for large NN. From these attractors, we have applied our general formalism and also we have found that both, the potential V(N)V(N) and the coupling parameter g(N)g(N) present a power law relation with NN. From the condition 𝒜1{\cal{A}}\gg 1 in which predominate the Galileon term, we have obtained a lower bound for the integration constant β\beta given by βα3/4N5C6α4C(CV02κ1/3)3/4\beta\gg\frac{\alpha^{3/4}N^{\frac{5C-6\alpha}{4C}}}{(CV_{0}^{2}\kappa^{1/3})^{3/4}}. In particular for the specific cases in which N=60N=60, α=4\alpha=4 and V0/Mp4𝒪(1012)V_{0}/M_{p}^{4}\sim\mathcal{O}(10^{-12}), we have obtained that the lower limit for the constant β\beta results βMp11/2𝒪(1018)\beta M_{p}^{11/2}\sim\mathcal{O}(10^{18}). In this context, from the standard attractor point, we have found that the reconstruction of the scalar potential and the coupling parameter are given by Eqs.(34) and (35). Here, we have determined that the parameter g(ϕ)g(\phi) decays as g(ϕ)ϕγ2g(\phi)\propto\phi^{-\gamma_{2}} with γ2\gamma_{2} a positive constant.

Other, important attract in GR corresponds to r=(N(1+ξN))1r=(N(1+\xi N))^{-1} together ns1=2/Nn_{s}-1=-2/N, here we could not find an analytical expression for N=N(ϕ)N=N(\phi), in order to obtain the reconstruction of V(ϕ)V(\phi) and g(ϕ)g(\phi), respectively.

Another, analytical reconstruction in the model of G-inflation corresponds to the tensor to scalar ratio r(N)r(N) given by Eq.(36). Assuming that the function F(N)=F(N)= constant, we have find a simple relation between the number NN and ϕ\phi given by NϕN\propto\phi and then Eq.(36). Here, we have obtained that the parameter gg in terms of the scalar field is given by g(ϕ)ϕ4/3g(\phi)\propto\phi^{-4/3} and the effective potential V(ϕ)ϕ1/3/(3+ξϕ1/3)V(\phi)\propto\phi^{1/3}/(3+\xi^{\prime}\phi^{1/3}). Also, in particular we have find that for the case in which the parameter ξ1/N1/3\xi\gg 1/N^{1/3}, the scalar potential V(ϕ)V(\phi) becomes a constantR1 . In fact, by assuming the condition 𝒜(N)1{\cal{A}}(N)\gg 1, we have obtained a lower bound for the integration constants given by Eq.(40).

Finally in this paper, we have not addressed the reconstruction to another functions G(ϕ,X)G(\phi,X) in the action, such as G(ϕ,X)=g(ϕ)XnG(\phi,X)=g(\phi)\,X^{n} or g1(ϕ)X+g2(ϕ)X2+..g_{1}(\phi)X+g_{2}(\phi)X^{2}+.. or other. In this sense, we hope to return to this point in the near future.

Acknowledgements.
The author thanks Prof. Marco Olivares for useful discussions on the rnSr-n_{S} plane. This work was supported by Proyecto VRIEA-PUCV N0 039.309/2018.

References

  • (1) A.A. Starobinsky, Phys. Lett. B 91, 99 (1980).
  • (2) A. Guth , Phys. Rev. D 23, 347 (1981).
  • (3) A.D. Linde, Phys. Lett. B 108, 389 (1982); A.D. Linde, Phys. Lett. B 129, 177 (1983).
  • (4) K. Sato, Mon. Not. Roy. Astron. Soc.  195, 467 (1981).
  • (5) V.F. Mukhanov and G.V. Chibisov , JETP Letters 33, 532 (1981); S. W. Hawking,Phys. Lett. B 115, 295 (1982).
  • (6) A. Guth and S.-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982); A. A. Starobinsky, Phys. Lett. B 117, 175 (1982).
  • (7) D. Larson et al., Astrophys. J. Suppl. 192, 16 (2011); C. L. Bennett et al., Astrophys. J. Suppl. 192, 17 (2011); N. Jarosik et al., Astrophys. J. Suppl. 192, 14 (2011).
  • (8) G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl.  208, 19 (2013).
  • (9) P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys.  571, A16 (2014); P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys.  571, A22 (2014).
  • (10) P. A. R. Ade et al. [Planck Collaboration], Astron. Astrophys.  594, A20 (2016).
  • (11) P. A. R. Ade et. al. [BICEP2 and Keck Array Collaborations], Phys. Rev. Lett. 116, 031302 (2016).
  • (12) E. Di Valentino et al. [CORE Collaboration], JCAP 1804, 017 (2018); F. Finelli et al. [CORE Collaboration], JCAP 1804, 016 (2018).
  • (13) G. W. Horndeski, Int. J. Theor. Phys. 10, 363 (1974).
  • (14) T. Kobayashi, M. Yamaguchi and J. Yokoyama, Prog. Theor. Phys.  126, 511 (2011).
  • (15) C. Deffayet, O. Pujolas, I. Sawicki and A. Vikman, JCAP 1010, 026 (2010).
  • (16) T. Kobayashi, M. Yamaguchi and J. Yokoyama, Phys. Rev. Lett. 105, 231302 (2010).
  • (17) B. P. Abbott et al. [LIGO Scientific and Virgo Collaborations], Phys. Rev. Lett.  119, no. 16, 161101 (2017).
  • (18) B. P. Abbott et al. [LIGO Scientific and VIRGO Collaborations], Living Rev. Rel.  19, 1 (2016); B. P. Abbott et al. [LIGO Scientific and Virgo and Fermi-GBM and INTEGRAL Collaborations], Astrophys. J.  848, no. 2, L13 (2017).
  • (19) B. P. Abbott et al., Astrophys. J.  848, no. 2, L12 (2017).
  • (20) A. Nicolis, R. Rattazzi, E. Trincherini, Phys. Rev. D 79, 064036 (2009).
  • (21) K. Kamada, T. Kobayashi, M. Yamaguchi and J. Yokoyama, Phys. Rev. D 83, 083515 (2011).
  • (22) K. Kamada, T. Kobayashi, T. Kunimitsu, M. Yamaguchi and J. Yokoyama, Phys. Rev. D 88, no. 12, 123518 (2013).
  • (23) K. Kamada, T. Kobayashi, T. Takahashi, M. Yamaguchi and J. Yokoyama, Phys. Rev. D 86, 023504 (2012).
  • (24) S. Hirano, T. Kobayashi and S. Yokoyama, Phys. Rev. D 94, no. 10, 103515 (2016).
  • (25) S. Unnikrishnan and S. Shankaranarayanan, JCAP 1407, 003 (2014).
  • (26) R. Herrera, JCAP 1705, no. 05, 029 (2017).
  • (27) H. Bazrafshan Moghaddam, R. Brandenberger and J. Yokoyama, Phys. Rev. D 95, no. 6, 063529 (2017).
  • (28) M. He, J. Liu, S. Lu, S. Zhou, Y. F. Cai, Y. Wang and R. Brandenberger, JCAP 1612, no. 12, 040 (2016).
  • (29) Y. F. Cai, J. O. Gong, S. Pi, E. N. Saridakis and S. Y. Wu, Nucl. Phys. B 900, 517 (2015)
  • (30) J. Ohashi and S. Tsujikawa, JCAP 1210, 035 (2012).
  • (31) T. Kobayashi, Phys. Rev. D 81, 103533 (2010).
  • (32) S. Nesseris, A. De Felice, and S. Tsujikawa, Phys. Rev. D 82, 124054 (2010); A. De Felice and S. Tsujikawa, JCAP 1202, 007 (2012); D. Maity and P. Saha, arXiv:1801.08080 [hep-ph]; R. Herrera, N. Videla and M. Olivares, arXiv:1806.04232 [gr-qc].
  • (33) F. Lucchin, S. Matarrese, Phys. Rev. D 32, 1316 (1985); R. Easther, Class. Quantum Grav. 13, 1775 (1996).
  • (34) J. Martin, D. Schwarz, Phys. Lett. B 500, 1-7 (2001).
  • (35) X. z. Li and X. h. Zhai, Phys. Rev. D 67, 067501 (2003).
  • (36) R. Herrera and R. G. Perez, Phys. Rev. D 93, no. 6, 063516 (2016).
  • (37) V. Mukhanov, Eur. Phys. J. C 73, 2486 (2013).
  • (38) T. Chiba, PTEP 2015, no. 7, 073E02 (2015).
  • (39) T. Miranda, J. C. Fabris and O. F. Piattella, JCAP 1709, no. 09, 041 (2017); A. Achúcarro, R. Kallosh, A. Linde, D. G. Wang and Y. Welling, JCAP 1804, no. 04, 028 (2018);S. D. Odintsov and V. K. Oikonomou, Annals Phys.  388, 267 (2018);P. Christodoulidis, D. Roest and E. I. Sfakianakis, arXiv:1803.09841 [hep-th]; P. Carrilho, D. Mulryne, J. Ronayne and T. Tenkanen, arXiv:1804.10489 [astro-ph.CO]; S. D. Odintsov and V. K. Oikonomou, Nucl. Phys. B 929, 79 (2018); S. D. Odintsov and V. K. Oikonomou, Phys. Rev. D 97, no.6, 064005 (2018).
  • (40) R. Kallosh and A. Linde, JCAP 1307, 002 (2013).
  • (41) R. Kallosh and A. Linde, JCAP 1310, 033 (2013).
  • (42) D. Kaiser, Phys. Rev. D 52, 4295–4306 (1995); F Bezrukov and M. Shaposhnikov, Phys Lett B 659, 703, (2008).
  • (43) R. Kallosh, A. Linde, D. Roest, Phys Rev Lett, 112 011303, (2014).
  • (44) R. Herrera, Eur. Phys. J. C 78, no. 3, 245 (2018).
  • (45) Q. G. Huang, Phys. Rev. D 76, 061303 (2007).
  • (46) J. Lin, Q. Gao and Y. Gong, Mon. Not. Roy. Astron. Soc.  459, no. 4, 4029 (2016); Q. Gao, Sci. China Phys. Mech. Astron.  60, no. 9, 090411 (2017).
  • (47) D. Roest, JCAP 1401, 007 (2014).
  • (48) J. Garcia-Bellido and D. Roest, Phys. Rev. D 89, no. 10, 103527 (2014).
  • (49) P. Creminelli et al., Phys. Rev. D 92, no. 12, 123528 (2015).
  • (50) X. Gao and D. A. Steer, JCAP 1112, 019 (2011).
  • (51) A. De Felice and S. Tsujikawa, Phys. Rev. D 84, 083504 (2011).
  • (52) J. Ohashi and S. Tsujikawa, JCAP 1210, 035 (2012).
  • (53) A. De Felice, T. Kobayashi and S. Tsujikawa, Phys. Lett. B 706, 123 (2011); C. Burrage, C. de Rham, D. Seery and A. J. Tolley, JCAP 1101, 014 (2011).