-
One of the most successful approaches to explain early Universe phenomena is cosmic inflation [1−6], i.e., the accelerated expansion of the early Universe. This important idea has the fundamental implication that the shortcomings of the standard cosmology can be explained in an elegant manner. In addition, the origin of anisotropies observed in the cosmic microwave background (CMB) radiation itself becomes a natural theory [7−16]. In this context, one of the most remarkable advances in modern physics was the establishment of observational constraints that ruled out many inflationary models if they were not supported by observational data [17−20]. Indeed, the observational value of the spectral index
ns and the analysis of the consistent behavior of this spectral index versus the tensor to scalar ratio r help reduce the number of inflationary models. In fact, recent observational data [17] impose constraints on both parameters: an upper limit on the tensor to scalar ratio,r<0.1 (Planck alone), at a 95% confidence level (CL) and a value of the spectral indexns=0.9649±0.0042 at68 % CL.The most famous statement related to the scenario of inflation is that the Higgs boson of the standard model acts as the inflaton [21−26]. There are two approaches to obtain the field equations from the Lagrangian of this theory, namely the metric and Palatini formalisms. In the original scenario [23], general relativity is based on the metric formulation, where all gravitational degrees of freedom are carried by the metric field and the connection is fixed to be the Levi-Civita one. However, in the Palatini formulation of gravity, metric and connection are two independent variables. Interestingly, both formulations lead to the usual Einstein's field equations of motion in minimally coupled scenarios. However, under non-minimal coupling (NMC), different approaches lead to different predictions even when the Lagrangian density of the theory has the same form. In addition, the assumption of considering a non-minimal coupling to gravity is important to sufficiently flatten the Higgs potential at large field values [23] to match observations. A remarkable difference between the metric and Palatini formalisms arises from observational consequences. Indeed, predictions of Palatini Higgs inflation lead to an extremely small tensor to scalar ratio [27, 28] compared to the metric formalism. Another interesting feature of Palatini Higgs inflation is that it has a higher cutoff scale, above which the perturbation theory breaks down, than the metric theory [29]. For reviews on this topic, please see Ref. [30] for the metric and Ref. [31] for the Palatini Higgs inflation. Furthermore, Palatini Higgs inflation lowers the spectral index for the primordial spectrum of density perturbations and reduces the required number of e-folds to answer important cosmology questions [32]. In this study, we developed an alternative approach to connect the metric and the Palatini Higgs inflation called hybrid metric Palatini Higgs inflation
1 . This hybrid metric Palatini scenario was already studied in [34], where anf(R) Palatini correction to the Einstein-Hilbert Lagrangian was added. This type of hybrid theory typically emerges when perturbative quantization techniques are incorporated to Palatini formalisms [35]. It is connected to non-perturbative quantum geometries in interesting ways [36]. Moreover, the scalar-tensor representation of a metric Palatini formalism was found to be useful in cosmology with respect to local experiments, thereby overcoming any matter instabilities that may appear if the scalar field is only weakly connected to matter. In this regard, wormhole geometries and cosmological and astrophysical applications were examined in [37], demonstrating that accelerating solutions are possible. A dynamical system in a hybrid metric Palatini context was also analyzed in [38].In the present paper, we propose a novel approach to modified gravity in which elements from both theories are combined [39]. Thus, one can avoid shortcomings that emerge in pure metric or Palatini approaches, such as the cosmic expansion and structure formation. This recent formalism is called hybrid metric-Palatini gravity, which adds a Palatini scalar curvature to the Einstein-Hilbert action. The benefit of this type of hybrid metric Palatini is to preserve the advantage of the minimal metric approach while improving the non-minimal coupling from the metric by the Palatini one.
The aim of this work was to study the non-minimally coupled Higgs inflation under the hybrid metric-Palatini approach and check the results in light of observational data [17].
The paper is structured as follows. In Sec. II, from the action, we derive the basic field equations of the inflation model with NMC in a hybrid metric Palatini formalism. In Sec. III, we present the Friedmann equation and apply the slow roll conditions on it. In Secs. IV and V, we analyze cosmological perturbations. In Sec. VI, we consider a Higgs inflation model and check its viability. Finally, we present a summary and conclude the manuscript in Sec. VII.
-
We consider a hybrid Palatini model where the scalar field is non-minimally coupled to gravity. Its action is described by
S=∫d4x√−g(M2p2R+12ξϕ2ˆR+Lϕ(gμν,ϕ)),
(1) where g is the determinant of the metric tensor
gμν ;Mp is the Planck mass; R is the Einstein-Hilbert curvature term, determined by the metric tensorgμν ;ˆR is the Palatini curvature, which depends on the metric tensorgμν and connectionΓαβγ and is considered an independent variableˆR=ˆR(gμν,Γαβγ) [40]; ξ is the coupling constant; andLϕ is the lagrangian density of the scalar field ϕ, which takes the formLϕ=−12∇μϕ∇μϕ−V(ϕ),
(2) where
V(ϕ) is the scalar field potential.The variation of this action with respect to the independent connection gives
∇σ(ξϕ2√−ggμν)=0.
(3) The solution of this equation reveals that the independent connection is the Levi- Civita connection of the conformal metric
ˆgμν=ξϕ2gμν ,ˆΓρμσ=12ˆgλρ(∂μˆgλσ+∂σˆgμλ−∂λˆgμσ)=Γρμσ+ωϕ(δρσ∂μ(ϕ)+δρμ∂σ(ϕ)−gμσ∂ρ(ϕ)),
(4) where
ω=1 corresponds to the Palatini approach andω=0 to the metric one. The curvature tensorˆRμν is expressed in terms of the independent connectionˆΓαβγ [40],ˆRμν=ˆΓαμν,α−ˆΓαμα,ν+ˆΓααλˆΓλμν−ˆΓαμλˆΓλαν,
(5) and using Eq. (4), we can rewrite Eq. (5) as
ˆRμν=Rμν+ωϕ2[4∇μϕ∇νϕ−gμν(∇ϕ)2−2ϕ(∇μ∇ν+12gμν◻)ϕ],
where
Rμν is the curvature tensor in the metric formalism. The scalar curvatureˆR can be expressed in terms of the Einstein-Hilbert curvature asˆR=gμνˆRμν=R−6ωϕ◻ϕ.
(6) Varying the action expressed by Eq. (1) with respect to the metric tensor leads to
(M2p+ξϕ2)Gμν=(1+2ξ−4ξω)∇μϕ∇νϕ−(12+2ξ−ξω)gμν(∇ϕ)2−gμνV(ϕ)+2ξ(1+ω)ϕ[∇μ∇ν−gμν◻]ϕ,
(7) which can be rewritten as
F(ϕ)Gμν=κ2Tμν,
(8) where F denotes a function of ϕ given by
F(ϕ)=1+ξκ2ϕ2,
(9) and
Tμν is the matter energy-momentum tensor, which takes the formTμν=A∇μϕ∇νϕ−Bgμν(∇ϕ)2−gμνV(ϕ)+Cϕ[∇μ∇ν−gμν◻]ϕ,
(10) where
A=(1+2ξ−4ξω) ,B=(12+2ξ−ξω) , andC=2ξ(1+ω) are constants.In the case of
ω=0 , Eq. (7) describes NMC in the metric approach [41]. Meanwhile, in the caseξ=0 , we recover the case of general relativity.Finally, let us take the variation of the action expressed by Eq. (1) with respect to ϕ to obtain the modified Klein Gordon equation [40],
◻ϕ+ξˆRϕ−V,ϕ=0,
(11) where
◻ϕ=1√−g∂ν(√−ggμν∂μϕ) is the D'Alembertien andV,ϕ=dV/dϕ . -
In this section, we assume a homogeneous and isotropic Universe described by a spatially flat Robertson-Walker (RW) metric with the signature (–,+,+,+) [42],
ds2=−dt2+a2(t)(dx2+dy2+dz2),
(12) where
a(t) is the scale factor and t is the cosmic time. The Friedmann equation is obtained by taking the 00 component from Eq. (7),H2=κ23F(ϕ)[(12−3ξω)˙ϕ2+V(ϕ)−6Hξ(1+ω)ϕ˙ϕ],
(13) where
H=˙a/a is the Hubble parameter and a dot denotes the differentiation with respect to cosmic time. Under slow roll conditions,˙ϕϕ<<H and˙ϕ2<<V , and Eq. (13) can be approximated byH2≃κ2V(ϕ)3(1+ξκ2ϕ2).
(14) By replacing
◻ϕ ,ˆR , and R by their expressions, the inflaton field equation Eq. (11) becomes−3H˙ϕ(1−6ξω)+12ξϕH2−V,ϕ≃0.
(15) -
In this section, we derive the scalar cosmological perturbations in detail. We choose the Newtonian gauge, in which the scalar metric perturbations of a RW background are given by [43, 44]
ds2=−(1+2Φ)dt2+a(t)2(1−2Ψ)δijdxidxj,
(16) where
Φ(t,x) andΨ(t,x) are the scalar perturbations, also called Bardeen variables.The perturbed Einstein's equations are given by
δF(ϕ)Gμν+F(ϕ)δGμν=κ2δTμν.
(17) For the perturbed metric expressed by Eq. (16), we obtain the individual components of Eq. (17) in the form
−6ξκ2H2ϕδϕ+F(ϕ)[6H(˙Ψ+HΦ)−2∇2a2Ψ]=κ2δT00,
−2F(ϕ)(˙Ψ+HΦ),i=κ2δT0i,
−6ξκ2ϕδϕ(3H2+2˙H)+6F(ϕ)[(3H2+2˙H)Φ+H(˙Φ+3˙Ψ)+¨Ψ+∇23a2(Φ−Ψ)]=κ2δTii,
F(ϕ)a−2(Ψ−Φ),i,j=κ2δTij.
The perturbed energy momentum tensor
δTμν appearing in Eq. (17) is given by [45]δTμν=(−δρaδq,i−a−1δq,iδpδij+δπij),
(19) where
δρ ,δq , andδp represent the perturbed energy density, momentum, and pressure, respectively. The anisotropic stress tensor is given byδπij=(△ij−13δij△)δπ , where△ij is defined by△ij=δik∂k∂j and△=△ii .Now, let us simplify the calculations and study the evolution of perturbations. Therefore, we decompose the function
ψ(x,t) into its Fourier componentsψk(t) asψ(t,x)=1(2π)3/2∫e−ikxψk(t)d3k,
(20) where k is the wave number. The perturbed equations in Eq. (18) can be expressed as
−ξκ2Hϕδϕ+F(ϕ)[H(˙Ψ+HΦ)+k23a2Ψ]=−κ26δρ,
F(ϕ)(˙Ψ+HΦ)=−κ22aδq,
−ξκ2ϕδϕ(3H2+2˙H)+F(ϕ)[(3H2+2˙H)Φ
+H(˙Φ+3˙Ψ)+¨Ψ−k23a2(Φ−Ψ)]=κ22δp,
F(ϕ)(Ψ−Φ),i,j=κ2a2δπij.
By using the perturbed energy momentum tensor, one can write the perturbed energy density, perturbed momentum, perturbed pressure, and anisotropic stress tensor, respectively, as follows:
−δρ=2(A−B)Φ˙ϕ2−2(A−B)˙ϕδ˙ϕ−V,ϕδϕ+3CH[˙ϕδϕ+ϕδ˙ϕ]+6CH(Ψ+Φ)ϕ˙ϕ−Cϕa−2△δϕ,
aδq=−A˙ϕδϕ−Cϕ(δ˙ϕ−Φ˙ϕ−Hδϕ),
δp=2B(˙ϕδ˙ϕ−Φ˙ϕ2)−Vϕδϕ+2CH˙ϕδϕ+Cϕ[2Hδ˙Φ−2Φ¨ϕ+δ¨ϕ−4HΦ˙ϕ−2˙Ψ˙ϕ−a−2△δϕ],
δπij=a−2Cϕδϕ,i,j.
The perturbed equation of motion for ϕ takes the form
2(A−B)˙ϕδ¨ϕ+[2(A−B)¨ϕ+Vϕ−3C˙Hϕ+6(A−C)H˙ϕ−3CH2ϕ]δ˙ϕ+[Vϕϕ˙ϕ+(A−C)˙ϕk2a2−3CH2˙ϕ−2CHϕk2a2]δϕ=2(A−B)[˙Φ˙ϕ2+2Φ˙ϕ¨ϕ]+6C˙Hϕ(Φ+Ψ)˙ϕ+6CH[(˙Φ+˙Ψ)ϕ˙ϕ+(Φ+Ψ)(˙ϕ2+ϕ¨ϕ)]+Cϕ˙ϕk2a2Φ+6AHΦ˙ϕ2+30CH2Φ˙ϕϕ+18CH2Ψ˙ϕϕ+6CϕH(Φ+Ψ)¨ϕ. (23) Therefore, if we adopt the slow roll conditions at large scales, i.e.,
k≪aH , we can neglect˙Φ ,˙Ψ ,¨Φ , and¨Ψ [46, 47]. In fact, throughout the cosmic history of the Universe, significant scales have primarily existed well beyond the Hubble radius, and they have only recently reentered the Universe. Consequently, it is reasonable to consider large scales as a valid assumption. Indeed, to satisfy the longitudinal post-Newtonian limit, we need to consider thatΔΦ≫a2H2×(Φ,˙Φ,¨Φ) ; similar assumptions are taken for the other gradient terms as well. In the case of plane wave perturbation with wavelength λ, when the conditionλ≪1/H is met,H2Φ becomes much smaller thanΔΦ . For˙Φ to be also negligible, the conditiondlogΦdloga≪1λH2 is required, which is satisfied ifλ≪1/H for perturbation growth. The same arguments may be used for¨Φ and for the metric potential Ψ [46, 48]. Hence, we can rewrite Eq. (23) as(1−6ξω)δ¨ϕ+[V,ϕ˙ϕ+6(1−6ξω)H−6ξ(1+ω)H2ϕ˙ϕ]δ˙ϕ+[V,ϕϕ−6ξ(1+ω)H2]δϕ+6H[(1+4ξ−2ξω)˙ϕ+10ξ(1+ω)Hϕ]Φ=0.
(24) Using Eqs. (21b) and (22b), the scalar perturbation Φ can be expressed in terms of the fluctuation of the scalar field
δϕ asΦ=κ2eff(A˙ϕ−CHϕ)2F(ϕ)Hδϕ,
(25) where
κ2eff=κ2/[1+Cκ22F(ϕ)Hϕ˙ϕ] .We define the comoving curvature perturbation as [49]
R=Ψ−Hρ+paδq.
(26) Hence, by considering the slow roll approximations at large scale, and according to Eq. (21b), one can find that
R=Ψ+H˙ϕ[1+Cκ22F(ϕ)Hϕ˙ϕ]δϕ.
(27) Considering the spatially flat gauge where
Ψ=0 , and according to Eq. (27), a new variable can be defined asδϕΨ=δϕ+˙ϕH[1+Cκ22F(ϕ)Hϕ˙ϕ]Ψ.
(28) Using Eq. (21b) in this gauge, Eq. (24) can be expressed as
(1−6ξω)δ¨ϕΨ+3H[(1−6ξω)−2ξHϕ˙ϕ(ω−2)]δ˙ϕΨ+[V,ϕϕ−6ξωH2−6κ2eff((1+2ξ−4ξω)˙ϕ−2ξ(1+ω)Hϕ)×(1+4ξ−2ξω)˙ϕ+10ξ(1+ω)Hϕ2F(ϕ)]δϕΨ=0.
(29) Introducing the Mukhanov-Sasaki variable
v=aδϕΨ allows rewriting the perturbed equation of motion Eq. (29) asv″−1τ2[ν2−14]v=0,
(30) where the derivative with respect to the conformal time τ is denoted by the prime, and the term ν is
ν=32+ϵ−˜η+˜ζ3+2˜χ,
(31) where we have used the slow roll parametres given by
ϵ=1−H′H2=12κ2(VϕV)2C1,
(32) η=a2Vϕϕ3H2,
(33) ζ=6ξω,
(34) χ=κ2eff((1+2ξ−4ξω)ϕ′−2ξ(1+ω)Hϕ)×(1+4ξ−2ξω)ϕ′+10ξ(1+ω)Hϕ2FH2,
(35) and
˜η=1(1−6ξω)η,
(36) ˜ζ=1(1−6ξω)ζ,
(37) ˜χ=1(1−6ξω)χ.
(38) We have also introduced the correction term to the standard expression as
C1=F(ϕ)(1−6ξω)(1−4ξκ2ϕF(ϕ)VVϕ)(1−2ξκ2ϕF(ϕ)VVϕ).
(39) This term characterizes the effect of NMC (through the constant ξ) and the Palatini approach (through ω).
The solution to Eq. (30) is given by [50]
v=aH√2k3(kaH)3/2−ν.
(40) The power spectrum for the scalar field perturbations reads as [49]
Pδϕ=4πk3(2π)3|va|2,
(41) and the spectral index of the power spectrum is given by [49]
ns−1=dLnPδϕdLnk|k=aH=3−2ν,
(42) which can be expressed in terms of slow roll parametres as
ns=1−2ϵ+2˜η−2˜ζ3−4˜χ.
(43) The power spectrum of the curvature perturbations is defined as [49]
A2s=425PR=4254πk3(2π)3|R|2
(44) =(2H5˙ϕ[1+Cκ22F(ϕ)H˙ϕϕ])2Pδϕ,
(45) and assuming the slow-roll conditions, it becomes
A2s=425(2π)2H4˙ϕ2[1+Cκ22F(ϕ)H˙ϕϕ]2=κ6V375π2V2,ϕC2,
(46) where
C2=(1−6ξω)2F(ϕ)[1+Cκ22F(ϕ)H˙ϕϕ]2V2ϕ(2F,ϕV−FV,ϕ)2,
(47) is a correction to the standard expression of the power spectrum. This correction term depends on NMC and the Palatini approach effect.
-
The tensor to scalar ratio is an important observable parameter in cosmology. Observational data [17] provide an upper limit on this ratio,
r<0.1 , at a 95% confidence level. To introduce this parameter, we need to define the tensor perturbations amplitude as [51]A2T=2κ225(H2π)2,
(48) which, in our model, takes the form
A2T=4κ4V600π2C3,
(49) where the correction term
C3 is defined asC3=1F(ϕ).
(50) Furthermore, we can define the tensor to scalar ratio, which is a useful inflationary parameter, as
r=A2TA2S=12κ2V2ϕV2[1+Cκ22F(ϕ)H˙ϕϕ]2(1−6ξω)2.
(51) -
In this section, as an application, we study a Higgs inflationary model in which we consider that the Higgs boson (the inflaton) is NMC to the gravity within the hybrid metric Palatini approach developed in the previous sections. We also check the viability of the model by comparing our results with observational data [17]. In this case, we consider the quartic potential [52]
V(ϕ)=λ4ϕ4,
(52) where λ is the Higgs self-coupling. During inflation, the number of e-folds is given by [53]
N=∫tFtIHdt=∫ϕ(tF)ϕ(tI)H˙ϕdϕ.
(53) From Eq. (15), we have that
˙ϕ=12ξϕH2−Vϕ3H(1−6ξω),
(54) and we obtain
N=(1−6ξω)κ28[ϕ2(tI)−ϕ2(tF)],
(55) where the subscript I and F represent the crossing horizon and end of inflation, respectively. Considering
ϕ2(tI)≫ϕ2(tF) , we obtainϕ2(tI)=8Nκ2(1−6ξω).
(56) Figure 1 depicts the variation of the number of e-folds, N, versus the scalar field for a Higgs self-coupling
λ=0.13 [21] and a coupling constantξ=10−3.5 . This figure shows that for an appropriate range of N, i.e.,50<N<70 , we obtain a large field whereκϕ≫20 .Figure 1. (color online) Plot of the number of e-folds versus the scalar field ϕ for
ξ=10−3.5 andλ=0.13. The slow roll parameter defined in Eq. (32) becomes
ϵ=8κ2ϕ2(1−6ξω)(1−ξκ2ϕ22F),
(57) Figure 2 represents the evolution of the correction term
C1 as a function of the coupling constant ξ. Note that the effect of the Palatini parameter ω onC1 begins from an approximate value of10−4 . Note also that, forξ=0 , the correction term reduces to one, and the standard expression of the slow roll parameter is recovered. Forξ≠0 andω=0 , we recover the slow roll parameter expression in the case of NMC within the metric approach.Figure 2. (color online) Variation of the correction term
C1 as a function of the coupling constant forN=45 .The spectral index of the power spectrum given by Eq. (43) can be expressed as
ns=1−16κ2ϕ2(1−6ξω)(1−ξκ2ϕ22F)+2(1−6ξω)[12Fκ2ϕ2−2ξω−κeff((1−4ξω+2ξ)˙ϕ−2ξ(1+ω)Hϕ)(1−2ξω+4ξ)H˙ϕ+10ξ(1+ω)H2ϕFH3].
(58) Figures 3(a) and 3(b) illustrate the variation of
ns against the number of e-folds N and against the scalar field forN=45 , respectively, forλ=0.13 and for different values of the coupling constant ξ, i.e.,10−3.5,10−4,0, and−10−4 . The gray horizontal bound in both figures represents the limits for the spectral index imposed by Planck data. We conclude that the predictions ofns are consistent with the observational data forξ=10−4 andξ=10−3.5 .Figure 3. (color online) Evolution of
ns against the number of e-folds (a) and against the scalar field (b) for different values of the coupling constant ξ andλ=0.13 .From Eqs. (46) and (49), we can obtain the power spectrum of the amplitudes of the curvature and tensor perturbations as
A2s=λκ6ϕ64800π2C2,
(59) A2T=λκ4ϕ4600π2C3,
(60) respectively.
The behavior of
C2 is shown in Fig. 4. We present this term versus the coupling constant ξ in the cases of the hybrid Palatini metric formalism (blue curve) and metric formalism (green curve). The effect of the Palatini parameter ω onA2s emerges fromξ=5×10−3 .Figure 4. (color online) Variation of the correction term
C2 versus the coupling constant for a number of e-foldsN=45 .The correction term
C3 is plotted as a function of ξ in Fig. 5. Note that the effect of the Palatini parameter emerges from a value ofξ=10−2 .Figure 5. (color online) Variation of the correction term
C3 against the coupling constant for a number of e-foldsN=45 .From Eq. (51), the tensor to scalar ratio can be obtained as
r=8κ2ϕ2[1+Cκ22F(ϕ)H˙ϕϕ]2(1−6ξω)2.
(61) Figure 6 shows the evolution of r versus the number of e-folds N for
λ=0.13 and for selected values of the coupling constant ξ. Note that r lies within the bounds imposed by observational data [17] in the appropriate range of N for the selected values of ξ.Figure 6. (color online) Variation of the tensor to scalar ratio r as a function of the number of e-folds N for different values of the coupling constant.
Figure 7 shows the
(ns,r) plane for different values of the coupling constant ξ in the range of the number of e-folds30≤N≤90 with the constraints from the Planck TT, TE, EE+LowE+lensing (gray contour) as well as Planck TT, TE, EE+lowE+lensing+BK14 data (red contour). Note also thatns−r predictions for the case whereξ≤0 are ruled out at95% confidence level contour according to the current observational data [17]. Furthermore, forξ=10−3.5 , observational parameters lie within68% CL contour for a range of the number of e-folds40.9≤N≤47 (low-N scenario). In addition, we obtain the central value of the index spectralns=0.9649 with a small value of tensor to scalar ratior=0.022 forN=43.43 . Forξ=10−4 , the results are inside the68% CL contour for the range67.8≤N≤86 (high-N scenario). However,N=75.41 givesns=0.9649 andr=0.013 . Thus, we can conclude that NMC in the framework of hybrid metric Palatini can ensure successful Higgs inflation. In the literature, it was reported that NMC in the Pure Palatini formalism requires a large value of ξ and results in an extremely small value of tensor to scalar ratior∼10−12 [54−56]. Therefore, the hybrid model may be an effective approach to solve this issue by increasing the value of r, making it comparable with the corresponding values predicted by the original metric approach. Then, it may be probed by future experiments [57, 58] where the value of the tensor to scalar ratio is on the order ofr∼10−2 . -
In this study, we investigated a cosmological model where the field is non-minimally coupled with gravity in the hybrid metric Palatini approach.
We also analyzed the cosmological perturbations to determine the different parameters during the inflationary period. As previously mentioned, the existence of correction terms to the standard background and perturbative parameters represents the impact of the Palatini approach and the non-minimal coupling between the scalar field and the Ricci scalar.
We applied our model by comprehensively developing a non-minimally coupled inflationary model driven by the Higgs field with a quartic potential within the slow-roll approximation.
We checked our results by plotting the evolution of different inflationary parameters versus the constraints provided by the observational data.
We found that perturbed parameters such as the tensor to scalar ratio and the scalar spectral index are compatible with the observational data for an appropriate range of the number of e-folds for different values of ξ, as shown in Figs. 3 and 6.
We plotted the different correction terms to the standard case versus the coupling constant. We showed that they depend on NMC and the Palatini effect.
Finally, for further checking the consistency of our model, we compared our theoretical predictions with observational data [17] by plotting the Planck confidence contours in the plane of
ns−r (Fig. 6). The results show that the predicted parameters are in good agreement with the Planck data for two values of the NMC constant:ξ=10−3.5 andξ=10−4 .
