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

**footnotetext: Corresponding author.

Nonminimally Assisted Inflation:
A General Analysis

Sang Chul Hyun    Jinsu Kim    Tatsuki Kodama    Seong Chan Park    and Tomo Takahashi
Abstract

The effects of a scalar field, known as the “assistant field,” which nonminimally couples to gravity, on single-field inflationary models are studied. The analysis provides analytical expressions for inflationary observables such as the spectral index (nsn_{s}), the tensor-to-scalar ratio (rr), and the local-type nonlinearity parameter (fNL(local)f_{\rm NL}^{(\rm local)}). The presence of the assistant field leads to a lowering of nsn_{s} and rr in most of the parameter space, compared to the original predictions. In some cases, nsn_{s} may increase due to the assistant field. This revives compatibility between ruled-out single-field models and the latest observations by Planck-BICEP/Keck. The results are demonstrated using three example models: loop inflation, power-law inflation, and hybrid inflation.

1 Introduction

The latest joint analysis of the Planck results [1] and BICEP/Keck (BK) data [2] has put a strong constraint on the spectral index, 0.958ns0.9750.958\leq n_{s}\leq 0.975 (95% C.L.), and a stringent upper bound on the tensor-to-scalar ratio, r0.036r\leq 0.036 (95% C.L.). As a consequence, a plethora of single-field inflationary models have been ruled out. Notably, the chaotic inflation model [3] with a power-law potential V(ϕ)ϕpV(\phi)\sim\phi^{p}, the power-law inflation model [4] with an exponential potential V(ϕ)exp(λϕ)V(\phi)\sim\exp(-\lambda\phi), the hybrid inflation model [5] with an effectively single-field potential V(ϕ)1+(ϕ/μ)2V(\phi)\sim 1+(\phi/\mu)^{2}, and the loop inflation (also known as spontaneously broken SUSY) model [6] with V(ϕ)1+λlog(ϕ/MP)V(\phi)\sim 1+\lambda\log(\phi/M_{\rm P}) are ruled out as either the tensor-to-scalar ratio rr or the spectral index nsn_{s} is too large to be allowed by the latest bounds.

To save single-field inflationary models, many mechanisms have been proposed. For instance, an introduction of the so-called nonminimal coupling of the inflaton field to gravity of the form ξϕ2R\xi\phi^{2}R [7, 8, 9, 10, 11, 12], where RR is the Ricci scalar, is known to reduce the tensor-to-scalar ratio111 See also, e.g., Refs. [13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26] for the analysis of the effects of the nonminimal coupling to gravity in single-field inflation models for observables such as the spectral index and tensor-to-scalar ratio. . The nonminimally-coupled models predict similar tensor-to-scalar ratio and spectral index values as Starobinsky’s R2R^{2} model [27] as they are approximately equivalent to each other under the slow-roll assumption for a particular combination of the inflaton potential and the nonminimal coupling function222 See, e.g., Refs. [28, 29, 30, 31] for the R2R^{2}-type models in the Palatini formulation. . Alternatively, one may modify the inflaton kinetic term. The α\alpha-attractor model [32, 33] is a well-known example where the kinetic term has a pole. All of the aforementioned models, the nonminimally-coupled model, Starobinsky’s R2R^{2} model, and the α\alpha-attractor model, have one thing in common: The potential in the Einstein frame which can be achieved via the Weyl rescaling features a plateau region in the large-field limit where the inflaton slowly rolls during inflation.

In Ref. [34], an alternative approach has been proposed. Instead of directly modifying the inflaton sector in such a way that the Einstein-frame potential becomes flat, an extra scalar field ss, dubbed an assistant field, is introduced. The assistant field does not directly interact with the inflaton field ϕ\phi, while it talks to gravity through the nonminimal coupling smR\sim s^{m}R with m>0m>0. The assistant field is further assumed to be effectively massless. Therefore, the potential in the original Jordan frame is given only by the original ϕ\phi field. Even though the additional assistant field ss has no direct interaction with the original inflaton field ϕ\phi, due to its nonminimal coupling to the Ricci scalar, the potential in the Einstein frame and the dynamics of inflation become nontrivial. It was shown in Ref. [34] that the chaotic inflation model, when assisted by the assistant field, may be revived and become compatible with the latest observational constraints.

In this work, we generalise the analysis of Ref. [34] to a general inflationary model for the original inflaton field ϕ\phi. We remain agnostic as to the form of the ϕ\phi-field potential. We follow the same nomenclature as Ref. [34] and call the extra scalar field ss the assistant field. The characteristics of the assistant field are: (i) it does not directly couple to the original inflaton field ϕ\phi, (ii) it nonminimally couples to gravity through smRs^{m}R, and (iii) it is effectively massless. Including the assistant field ss and its explicit nonminimal coupling term, we perform a general analytical study about inflationary observables such as the spectral index and the tensor-to-scalar ratio. Our analytical formulae for the spectral index and the tensor-to-scalar ratio can readily be applied to a vast range of inflationary models, and one may easily see whether an otherwise ruled-out model can become revived with the help of the assistant field. We also discuss the non-Gaussianity and running of the spectral index that may become sizeable for multifield models, providing the corresponding analytical formula. To demonstrate how the assistant field affects observables in concrete models, we consider as an example three models for the original ϕ\phi field, namely the loop inflation model, the power-law inflation model, and the hybrid inflation model, and show how the presence of the assistant field may bring the models to the observationally-favoured region.

The rest of the paper is organised as follows. In Sec. 2, we start with a generic action for a multifield inflation model with a nonminimal coupling to gravity and perform the Weyl rescaling to the Einstein frame, setting the notations. Defining properties for the assistant field are introduced, and we set our model with a general potential for the original ϕ\phi field. Section 3 presents the general analytical study for the spectral index, the tensor-to-scalar ratio, and the non-Gaussianity. We show how the inflationary observables are modified in comparison with the original predictions due to the presence of the assistant field. As an application of our general analytical formulae, three examples, the loop inflation model, power-law inflation, and hybrid inflation, are discussed in Sec. 4. We show how these models become compatible with the latest observational data. We conclude in Sec. 5. Appendix A summarises the computation and general behaviour of the running of the spectral index.

2 Model

Let us consider the following generic action for multifield inflation:

S=d4xgJ[f(φi)gJμνRJμν(ΓJ)12GJijgJμνμφiνφjVJ(φi)],\displaystyle S=\int d^{4}x\,\sqrt{-g_{\rm J}}\left[f(\varphi^{i})g_{\rm J}^{\mu\nu}R_{{\rm J}\mu\nu}(\Gamma_{\rm J})-\frac{1}{2}G_{{\rm J}ij}g_{\rm J}^{\mu\nu}\partial_{\mu}\varphi^{i}\partial_{\nu}\varphi^{j}-V_{\rm J}(\varphi^{i})\right]\,, (2.1)

where the subscript J indicates that the action is written in the Jordan frame, f(φi)f(\varphi^{i}) is a function of fields that represents the nonminimal coupling, GJijG_{{\rm J}ij} denotes the kinetic mixing between fields, and VJ(φi)V_{\rm J}(\varphi^{i}) is the Jordan-frame potential333 Multifield inflation models with the nonminimal coupling to gravity have been discussed in a certain range of contexts. See, e.g., Refs. [35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52] for such works. . Here, the Ricci tensor RJμνR_{{\rm J}\mu\nu} is given as a function of the connection ΓJ\Gamma_{\rm J} to capture not only the standard metric formulation, but also the Palatini formulation. In the Palatini formulation, the metric and the connection are a priori independent to each other, while in the metric formulation, the connection is related to the metric, becoming the Levi-Civita connection. The connection is assumed to be torsion free.

Through the Weyl rescaling, gJμνgEμν=Ω2gJμνg_{{\rm J}\mu\nu}\rightarrow g_{{\rm E}\mu\nu}=\Omega^{2}g_{{\rm J}\mu\nu}, the Jordan-frame action (2.1) can be brought to the Einstein frame, denoted by the subscript E. The conformal factor Ω2\Omega^{2} is chosen to be 2f/MP22f/M_{\rm P}^{2} with MPM_{\rm P} being the reduced Planck mass. After the Weyl rescaling, we obtain the Einstein-frame action as follows:

S=d4xgE[MP22gEμνREμν(ΓE)12GEijgEμνμφiνφjVE],\displaystyle S=\int d^{4}x\,\sqrt{-g_{\rm E}}\left[\frac{M_{\rm P}^{2}}{2}g_{\rm E}^{\mu\nu}R_{{\rm E}\mu\nu}(\Gamma_{\rm E})-\frac{1}{2}G_{{\rm E}ij}g_{\rm E}^{\mu\nu}\partial_{\mu}\varphi^{i}\partial_{\nu}\varphi^{j}-V_{\rm E}\right]\,, (2.2)

where VEV_{\rm E} is the Einstein-frame potential given by

VE=VJΩ4=MP4VJ4f2,\displaystyle V_{\rm E}=\frac{V_{\rm J}}{\Omega^{4}}=\frac{M_{\rm P}^{4}V_{\rm J}}{4f^{2}}\,, (2.3)

and

GEij=MP22f(GJij+κ3f,if,jf),\displaystyle G_{{\rm E}ij}=\frac{M_{\rm P}^{2}}{2f}\left(G_{{\rm J}ij}+\kappa\frac{3f_{,i}f_{,j}}{f}\right)\,, (2.4)

with f,if/φif_{,i}\equiv\partial f/\partial\varphi^{i}, represents the field-space metric. The parameter κ\kappa parametrises which framework we work in; κ=0\kappa=0 denotes the Palatini formulation while κ=1\kappa=1 corresponds to the metric formulation. Hereinafter, we omit the subscript E for brevity.

Our interests are the cases where the system is comprised of two fields, namely the original inflaton field ϕ\phi and the assistant field ss. The assistant field is effectively massless and does not directly couple to the ϕ\phi field in the Jordan frame. These properties indicate that GJij=δijG_{{\rm J}ij}=\delta_{ij}, i.e., no kinetic mixing between the two scalar fields in the original Jordan frame, and that the Jordan-frame potential VJV_{\rm J} becomes a function of only the ϕ\phi field, i.e., VJ=VJ(ϕ)V_{\rm J}=V_{\rm J}(\phi). Furthermore, the assistant field nonminimally couples to gravity, i.e., f=f(s)f=f(s). Upon imposing a Z2Z_{2} symmetry, a natural choice for the nonminimal coupling from the dimensional analysis viewpoint would be

f(s)=MP22[1+ξ2(sMP)2],\displaystyle f(s)=\frac{M_{\rm P}^{2}}{2}\left[1+\xi_{2}\left(\frac{s}{M_{\rm P}}\right)^{2}\right]\,, (2.5)

where ξ2\xi_{2} is a dimensionless coupling. Moreover, in order for the assistant field not to significantly alter, but to merely assist the inflationary dynamics, the nonminimal coupling term takes not too large values, i.e., ξ2(s/MP)21\xi_{2}(s/M_{\rm P})^{2}\ll 1. We further note that the ξ2(s/MP)20\xi_{2}(s/M_{\rm P})^{2}\rightarrow 0 limit corresponds to the original ϕ\phi-field inflation. In the current work, we focus on the quadratic nonminimal coupling of the form (2.5). However, our results can easily be extended to f(s)=(MP2/2)[1+ξm(s/MP)m]f(s)=(M_{\rm P}^{2}/2)[1+\xi_{m}(s/M_{\rm P})^{m}] with m>2m>2 in a straightforward manner.

The action in the Einstein frame (2.2) is therefore obtained as follows:

S=d4xg[MP22gμνRμν(Γ)12𝒦1gμνμϕνϕ12𝒦2gμνμsνsV],\displaystyle S=\int d^{4}x\,\sqrt{-g}\,\left[\frac{M_{\rm P}^{2}}{2}g^{\mu\nu}R_{\mu\nu}(\Gamma)-\frac{1}{2}\mathcal{K}_{1}g^{\mu\nu}\partial_{\mu}\phi\partial_{\nu}\phi-\frac{1}{2}\mathcal{K}_{2}g^{\mu\nu}\partial_{\mu}s\partial_{\nu}s-V\right]\,, (2.6)

where

𝒦1\displaystyle\mathcal{K}_{1} =MP22f,\displaystyle=\frac{M_{\rm P}^{2}}{2f}\,, (2.7)
𝒦2\displaystyle\mathcal{K}_{2} =MP22f+κ3MP22(f,sf)2,\displaystyle=\frac{M_{\rm P}^{2}}{2f}+\kappa\frac{3M_{\rm P}^{2}}{2}\left(\frac{f_{,s}}{f}\right)^{2}\,, (2.8)

and the Einstein-frame potential is given by

V=MP44f2(s)VJ(ϕ).\displaystyle V=\frac{M_{\rm P}^{4}}{4f^{2}(s)}V_{\rm J}(\phi)\,. (2.9)

We note that the potential is product-separable, V=K(s)VJ(ϕ)V=K(s)V_{\rm J}(\phi), where K(s)=MP4/(4f2(s))K(s)=M_{\rm P}^{4}/(4f^{2}(s)). Upon canonically normalising the ss field through

(σs)2=𝒦2,\displaystyle\left(\frac{\partial\sigma}{\partial s}\right)^{2}=\mathcal{K}_{2}\,, (2.10)

we may rewrite the action as

S=d4xg[MP22gμνRμν(Γ)12e2bgμνμϕμϕ12gμνμσνσV],\displaystyle S=\int d^{4}x\,\sqrt{-g}\,\left[\frac{M_{\rm P}^{2}}{2}g^{\mu\nu}R_{\mu\nu}(\Gamma)-\frac{1}{2}e^{2b}g^{\mu\nu}\partial_{\mu}\phi\partial_{\mu}\phi-\frac{1}{2}g^{\mu\nu}\partial_{\mu}\sigma\partial_{\nu}\sigma-V\right]\,, (2.11)

where b=b(σ(s))b=b(\sigma(s)) is defined via

e2b𝒦1=MP22f.\displaystyle e^{2b}\equiv\mathcal{K}_{1}=\frac{M_{\rm P}^{2}}{2f}\,. (2.12)

We note that the action (2.11) has been explored in detail in Refs. [53, 54]. For brevity, in the following, we set MP=1M_{\rm P}=1.

3 General Analysis

Inflationary observables such as the scalar power spectrum 𝒫ζ\mathcal{P}_{\zeta}, the scalar spectral index nsn_{s}, the tensor-to-scalar ratio rr, and the local-type shape-independent nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)} for the action (2.11) have been obtained in, e.g., Refs. [54, 34] by using the δN\delta N formalism [55, 56, 57, 58, 59] under the slow-roll approximation. We do not repeat the derivation here and simply list the resultant expressions. The spectral index and the tensor-to-scalar ratio are given by

ns\displaystyle n_{s} =12ϵσ2ϵϕ4e2Xu2α2/ϵσ+v2/ϵϕ112ηb+2ϵbu2α2/ϵσ+v2/ϵϕ(uαϵϕϵσvϵσϵϕ)2\displaystyle=1-2\epsilon^{\sigma}_{*}-2\epsilon^{\phi}_{*}-\frac{4e^{-2X}}{u^{2}\alpha^{2}/\epsilon_{*}^{\sigma}+v^{2}/\epsilon_{*}^{\phi}}-\frac{1}{12}\frac{\eta_{*}^{b}+2\epsilon_{*}^{b}}{u^{2}\alpha^{2}/\epsilon_{*}^{\sigma}+v^{2}/\epsilon_{*}^{\phi}}\left(u\alpha\sqrt{\frac{\epsilon_{*}^{\phi}}{\epsilon_{*}^{\sigma}}}-v\sqrt{\frac{\epsilon_{*}^{\sigma}}{\epsilon_{*}^{\phi}}}\right)^{2}
+2u2α2/ϵσ+v2/ϵϕ[u2α2ησσϵσ+v2ηϕϕϵϕ+4uvα+12sbsσϵbϵσv(vϵϕ2uαϵσ)],\displaystyle\quad+\frac{2}{u^{2}\alpha^{2}/\epsilon_{*}^{\sigma}+v^{2}/\epsilon_{*}^{\phi}}\left[u^{2}\alpha^{2}\frac{\eta_{*}^{\sigma\sigma}}{\epsilon_{*}^{\sigma}}+v^{2}\frac{\eta_{*}^{\phi\phi}}{\epsilon_{*}^{\phi}}+4uv\alpha+\frac{1}{2}s_{*}^{b}s_{*}^{\sigma}\sqrt{\epsilon_{*}^{b}\epsilon_{*}^{\sigma}}v\left(\frac{v}{\epsilon_{*}^{\phi}}-\frac{2u\alpha}{\epsilon_{*}^{\sigma}}\right)\right]\,, (3.1)
r\displaystyle r =16e2Xu2α2/ϵσ+v2/ϵϕ,\displaystyle=\frac{16e^{-2X}}{u^{2}\alpha^{2}/\epsilon_{*}^{\sigma}+v^{2}/\epsilon_{*}^{\phi}}\,, (3.2)

where

uϵeσϵeσ+ϵeϕ,vϵeϕϵeσ+ϵeϕ,X2be2b,αe2b2be[1+ϵeϕϵeσ(1e2be2b)].\displaystyle u\equiv\frac{\epsilon_{e}^{\sigma}}{\epsilon^{\sigma}_{e}+\epsilon^{\phi}_{e}}\,,\;\;v\equiv\frac{\epsilon_{e}^{\phi}}{\epsilon^{\sigma}_{e}+\epsilon^{\phi}_{e}}\,,\;\;X\equiv 2b_{e}-2b_{*}\,,\;\;\alpha\equiv e^{2b_{*}-2b_{e}}\left[1+\frac{\epsilon_{e}^{\phi}}{\epsilon_{e}^{\sigma}}\left(1-e^{2b_{e}-2b_{*}}\right)\right]\,. (3.3)

Here, sσsgn(V,σ)s^{\sigma}\equiv{\rm sgn}(V_{,\sigma}) and sbsgn(b,σ)s^{b}\equiv{\rm sgn}(b_{,\sigma}) denote the signs of the potential and bb derivatives with respect to the field σ\sigma, respectively. The subscript * (ee) indicates that the quantities are evaluated at the pivot scale (end of inflation). The slow-roll parameters are defined as

ϵσ12(V,σV)2=12(K,σK)2,ϵϕ12(V,ϕVeb)2=12(VJ,ϕVJ)2e2b,\displaystyle\epsilon^{\sigma}\equiv\frac{1}{2}\left(\frac{V_{,\sigma}}{V}\right)^{2}=\frac{1}{2}\left(\frac{K_{,\sigma}}{K}\right)^{2}\,,\qquad\epsilon^{\phi}\equiv\frac{1}{2}\left(\frac{V_{,\phi}}{V}e^{-b}\right)^{2}=\frac{1}{2}\left(\frac{V_{{\rm J},\phi}}{V_{\rm J}}\right)^{2}e^{-2b}\,,
ησσV,σσV=K,σσK,ηϕϕV,ϕϕVe2b=VJ,ϕϕVJe2b,\displaystyle\eta^{\sigma\sigma}\equiv\frac{V_{,\sigma\sigma}}{V}=\frac{K_{,\sigma\sigma}}{K}\,,\qquad\eta^{\phi\phi}\equiv\frac{V_{,\phi\phi}}{V}e^{-2b}=\frac{V_{{\rm J},\phi\phi}}{V_{\rm J}}e^{-2b}\,,
ηϕσV,ϕσVeb,ϵb8b,σ2,\displaystyle\eta^{\phi\sigma}\equiv\frac{V_{,\phi\sigma}}{V}e^{-b}\,,\qquad\epsilon^{b}\equiv 8b_{,\sigma}^{2}\,, (3.4)

and ηb16b,σσ\eta^{b}\equiv 16b_{,\sigma\sigma}. The local-type, shape-independent, nonlinearity parameter is given by

65fNL(local)\displaystyle-\frac{6}{5}f_{\rm NL}^{\rm(local)} =2eX(u2α2/ϵσ+v2/ϵϕ)2[(1ησσ2ϵσ)u3α3ϵσ+(1ηϕϕ2ϵϕ)v3ϵϕ\displaystyle=\frac{2e^{-X}}{(u^{2}\alpha^{2}/\epsilon_{*}^{\sigma}+v^{2}/\epsilon_{*}^{\phi})^{2}}\Bigg{[}\left(1-\frac{\eta_{*}^{\sigma\sigma}}{2\epsilon_{*}^{\sigma}}\right)\frac{u^{3}\alpha^{3}}{\epsilon_{*}^{\sigma}}+\left(1-\frac{\eta_{*}^{\phi\phi}}{2\epsilon_{*}^{\phi}}\right)\frac{v^{3}}{\epsilon_{*}^{\phi}}
+12sbsσu2vα2ϵσϵbϵσ+(uαϵσvϵϕ)2eX𝒞],\displaystyle\qquad+\frac{1}{2}s_{*}^{b}s_{*}^{\sigma}\frac{u^{2}v\alpha^{2}}{\epsilon_{*}^{\sigma}}\sqrt{\frac{\epsilon_{*}^{b}}{\epsilon_{*}^{\sigma}}}+\left(\frac{u\alpha}{\epsilon_{*}^{\sigma}}-\frac{v}{\epsilon_{*}^{\phi}}\right)^{2}e^{X}\mathcal{C}\Bigg{]}\,, (3.5)

where

𝒞\displaystyle\mathcal{C} ϵeσϵeϕϵe2(ϵeσηeϕϕ+ϵeϕηeσσϵe4ϵeϕϵeσϵe12seσsebϵebϵeσ(ϵeϕ)2ϵe).\displaystyle\equiv\frac{\epsilon_{e}^{\sigma}\epsilon_{e}^{\phi}}{\epsilon_{e}^{2}}\left(\frac{\epsilon_{e}^{\sigma}\eta_{e}^{\phi\phi}+\epsilon_{e}^{\phi}\eta_{e}^{\sigma\sigma}}{\epsilon_{e}}-4\frac{\epsilon_{e}^{\phi}\epsilon_{e}^{\sigma}}{\epsilon_{e}}-\frac{1}{2}s_{e}^{\sigma}s_{e}^{b}\sqrt{\frac{\epsilon_{e}^{b}}{\epsilon_{e}^{\sigma}}}\frac{(\epsilon_{e}^{\phi})^{2}}{\epsilon_{e}}\right)\,. (3.6)

We leave detailed discussion on the running of the spectral index αs\alpha_{s} in Appendix A.

The Einstein-frame potential for the assistant field is given by

K(s)=1(1+ξ2s2)2,\displaystyle K(s)=\frac{1}{(1+\xi_{2}s^{2})^{2}}\,, (3.7)

as can be seen from Eqs. (2.5) and (2.9). Thus, it is straightforward to obtain the slow-roll parameters for the ss field or, equivalently, its canonically-normalised version σ\sigma. They are given by

ϵσ=ϵb=8ξ22s21+ξ2s2(1+6κξ2),ησσ=4ξ2[1+3ξ2s2+4ξ22s4(1+6κξ2)][1+ξ2s2(1+6κξ2)]2,\displaystyle\epsilon^{\sigma}=\epsilon^{b}=\frac{8\xi_{2}^{2}s^{2}}{1+\xi_{2}s^{2}(1+6\kappa\xi_{2})}\,,\quad\eta^{\sigma\sigma}=\frac{4\xi_{2}[-1+3\xi_{2}s^{2}+4\xi_{2}^{2}s^{4}(1+6\kappa\xi_{2})]}{[1+\xi_{2}s^{2}(1+6\kappa\xi_{2})]^{2}}\,, (3.8)

where we have used Eq. (2.10). Note also that ηb=16ξ2(1+ξ2s2)/[1+ξ2s2(1+6κξ2)]2\eta^{b}=-16\xi_{2}(1+\xi_{2}s^{2})/[1+\xi_{2}s^{2}(1+6\kappa\xi_{2})]^{2}. One may regard the ss-field value at the pivot scale, ss_{*}, as an input parameter444 Quantum kick for the assistant field in the de Sitter phase can be estimated as (Δσ)2=𝒦2(Δs)2N2H2/(2π)2N2rAs/8(\Delta\sigma)^{2}=\mathcal{K}_{2}(\Delta s)^{2}\simeq N^{2}H_{*}^{2}/(2\pi)^{2}\simeq N^{2}rA_{s}/8, where As2×109A_{s}\approx 2\times 10^{-9} is the magnitude of the scalar power spectrum at the pivot scale and NN is the number of ee-folds. In the ξ2s21\xi_{2}s^{2}\ll 1 limit, we find Δs4.2×104\Delta s\lesssim 4.2\times 10^{-4} (1.8×1041.8\times 10^{-4}) for r<0.2r<0.2 (0.035) with N=60N=60. As far as a relatively larger value of ss is considered, such fluctuations can safely be neglected. . The ss-field value at the end of inflation is then given by requiring the number of ee-folds, NN, to be NeN_{e}; in this work, we choose Ne=60N_{e}=60 from the end of inflation. The number of ee-folds can be solely given by the assistant field ss,

N=eKK,σ𝑑σ=14ξ2[ln(ses)+3κξ2ln(1+ξ2se21+ξ2s2)].\displaystyle N=\int_{e}^{*}\frac{K}{K_{,\sigma}}\,d\sigma=\frac{1}{4\xi_{2}}\left[\ln\left(\frac{s_{e}}{s_{*}}\right)+3\kappa\xi_{2}\ln\left(\frac{1+\xi_{2}s_{e}^{2}}{1+\xi_{2}s_{*}^{2}}\right)\right]\,. (3.9)

In the Palatini formulation, κ=0\kappa=0, and one may obtain an exact analytical solution for ses_{e} as

se=se4ξ2Ne.\displaystyle s_{e}=s_{*}e^{4\xi_{2}N_{e}}\,. (3.10)

In the metric formulation, on the other hand, such a closed analytical form for the solution does not exist. However, an approximated form may be found. In this paper, we assume that the nonminimal coupling of the assistant field ss is small. Thus, we can take the small nonminimal coupling limit where ξ2s21\xi_{2}s^{2}\ll 1. In this limit, we may approximate the number of ee-folds as

N14ξ2ln(ses)+34ξ2(se2s2).\displaystyle N\approx\frac{1}{4\xi_{2}}\ln\left(\frac{s_{e}}{s_{*}}\right)+\frac{3}{4}\xi_{2}\left(s_{e}^{2}-s_{*}^{2}\right)\,. (3.11)

Requiring N=NeN=N_{e} then gives

sese4ξ2Ne[13ξ22s2(e8ξ2Ne1)].\displaystyle s_{e}\approx s_{*}e^{4\xi_{2}N_{e}}\left[1-3\xi_{2}^{2}s_{*}^{2}\left(e^{8\xi_{2}N_{e}}-1\right)\right]\,. (3.12)
Refer to caption
Figure 1: Comparison between the analytical solution for ses_{e} in the Palatini formulation (dotted lines), the approximated analytical solution for ses_{e} in the metric formulation (dashed lines), and the numerical solution for ses_{e} in the metric formulation (solid lines), for three different ss_{*} values, 0.001 (red), 0.01 (blue), and 0.1 (magenta). We observe that the approximated analytical solution in the metric case starts to deviate from the numerical solution as the small nonminimal coupling assumption breaks down. We also see that the difference between the Palatini and the metric formulations is negligible in the small nonminimal coupling limit.

Figure 1 shows the comparison between the analytical solution for ses_{e} in the Palatini formulation (3.10), the approximated analytical solution for ses_{e} in the metric formulation (3.12), and the numerical solution for ses_{e} found by using N=NeN=N_{e} with Eq. (3.9) in the metric formulation. We stress that the solution (3.10) in the Palatini case is exact. Deviations between the approximated solution and the exact numerical solution in the metric case start to appear when the assumption of the small nonminimal coupling breaks down; we shall use ξ2se20.1\xi_{2}s_{e}^{2}\leq 0.1 as the condition for the small nonminimal coupling in this paper. We also observe that the difference between the Palatini case and the metric case is negligible in the small nonminimal coupling limit; deviations become visible when the nonminimal coupling gets larger.

In order to discuss effects of the assistant field ss on the original ϕ\phi-field inflation model in a general way, we remain as agnostic as possible as to the form for the ϕ\phi-field potential, VJ(ϕ)V_{\rm J}(\phi). In the absence of the assistant field, the spectral index and the tensor-to-scalar ratio for the original ϕ\phi-field inflation model, denoted respectively by ns(0)n_{s}^{(0)} and r(0)r^{(0)}, are given by

ns(0)=16ϵ(0)+2η(0),r(0)=16ϵ(0),\displaystyle n_{s}^{(0)}=1-6\epsilon^{(0)}_{*}+2\eta^{(0)}_{*}\,,\quad r^{(0)}=16\epsilon^{(0)}_{*}\,, (3.13)

under the slow-roll approximation, where the slow-roll parameters are defined as

ϵ(0)12(VJ,ϕVJ)2,η(0)VJ,ϕϕVJ.\displaystyle\epsilon^{(0)}\equiv\frac{1}{2}\left(\frac{V_{{\rm J},\phi}}{V_{\rm J}}\right)^{2}\,,\quad\eta^{(0)}\equiv\frac{V_{{\rm J},\phi\phi}}{V_{\rm J}}\,. (3.14)

Comparing with the slow-roll parameters defined in Eq. (3.4), we see that

ϵϕ=ϵ(0)e2b,ηϕϕ=η(0)e2b.\displaystyle\epsilon^{\phi}=\epsilon^{(0)}e^{-2b}\,,\quad\eta^{\phi\phi}=\eta^{(0)}e^{-2b}\,. (3.15)

Writing ϵ(0)\epsilon^{(0)}_{*} and η(0)\eta^{(0)}_{*} in terms of ns(0)n^{(0)}_{s} and r(0)r^{(0)}, we get

ϵϕ=r(0)16e2b,ηϕϕ=12(ns(0)1+38r(0))e2b.\displaystyle\epsilon^{\phi}_{*}=\frac{r^{(0)}}{16}e^{-2b_{*}}\,,\quad\eta^{\phi\phi}_{*}=\frac{1}{2}\left(n^{(0)}_{s}-1+\frac{3}{8}r^{(0)}\right)e^{-2b_{*}}\,. (3.16)

We may thus substitute Eq. (3.16) into the observables, Eqs. (3.1), (3.2), and (3.5), to get rid of the information on the ϕ\phi field at the pivot scale. There still remains, however, dependence on the ϕ\phi field at the end of inflation through the ϕ\phi-field slow-roll parameters ϵeϕ\epsilon^{\phi}_{e} and ηeϕϕ\eta^{\phi\phi}_{e} or, equivalently, ϵe(0)\epsilon^{(0)}_{e} and ηe(0)\eta^{(0)}_{e}. In the current work, we consider two classes.

3.1 Class I: End of inflation via slow-roll violations

The first class we consider is the case where end of inflation is set by ϵe1\epsilon_{e}\approx 1, i.e., violation of slow roll. Examples include the chaotic inflation model with a power-law potential [3] and the loop inflation model with a logarithmic correction [6]. In this case, one may replace ϵeϕ\epsilon^{\phi}_{e} by 1ϵeσ1-\epsilon^{\sigma}_{e}. Utilising Eq. (3.16), one may express nsn_{s} and rr in terms of ξ2\xi_{2}, ss_{*}, ses_{e}, ns(0)n_{s}^{(0)}, and r(0)r^{(0)}. The spectral index is given by

ns(I)\displaystyle n_{s}^{({\rm I})} =1(1+ξ2s2)r(0)816ξ22s21+ξ2s2(1+6κξ2)32(1+ξ2se2)2(1+ξ2se2(1+6κξ2))2AI\displaystyle=1-\frac{(1+\xi_{2}s_{*}^{2})r^{(0)}}{8}-\frac{16\xi_{2}^{2}s_{*}^{2}}{1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})}-\frac{32(1+\xi_{2}s_{e}^{2})^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}{A_{\rm I}}
+4ξ23(1+ξ2s2)(1ξ22s4(1+6κξ2))r(0)3AI(1+ξ2s2(1+6κξ2))3\displaystyle\quad+\frac{4\xi_{2}^{3}(1+\xi_{2}s_{*}^{2})\left(1-\xi_{2}^{2}s_{*}^{4}(1+6\kappa\xi_{2})\right)r^{(0)}}{3A_{\rm I}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)^{3}}
×[(1+ξ2s2(1+6κξ2))BI4ξ2s32sr(0)(1+ξ2se2(1(86κ)ξ2))]2\displaystyle\quad\quad\times\bigg{[}\frac{\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)B_{\rm I}}{4\xi_{2}s_{*}}-\frac{32s_{*}}{r^{(0)}}\left(1+\xi_{2}s_{e}^{2}(1-(8-6\kappa)\xi_{2})\right)\bigg{]}^{2}
8(1+ξ2s2)2AI[2(1+ξ2se2(1(86κ)ξ2))2r(0)(88ns(0)3r(0))\displaystyle\quad-\frac{8(1+\xi_{2}s_{*}^{2})^{2}}{A_{\rm I}}\bigg{[}\frac{2\left(1+\xi_{2}s_{e}^{2}(1-(8-6\kappa)\xi_{2})\right)^{2}}{r^{(0)}}\left(8-8n_{s}^{(0)}-3r^{(0)}\right)
6ξ2(1+ξ2se2(1(86κ)ξ2))BI1+ξ2s2+ξ2(13ξ2s24ξ22s424κξ23s4)BI2s2(1+ξ2s2)2(1+ξ2s2(1+6κξ2))\displaystyle\quad\quad-\frac{6\xi_{2}\left(1+\xi_{2}s_{e}^{2}(1-(8-6\kappa)\xi_{2})\right)B_{\rm I}}{1+\xi_{2}s_{*}^{2}}+\frac{\xi_{2}\left(1-3\xi_{2}s_{*}^{2}-4\xi_{2}^{2}s_{*}^{4}-24\kappa\xi_{2}^{3}s_{*}^{4}\right)B_{\rm I}^{2}}{s_{*}^{2}(1+\xi_{2}s_{*}^{2})^{2}(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2}))}
128ξ22s2(1+ξ2se2(1(86κ)ξ2))2r(0)(1+ξ2s2)(1+ξ2s2(1+6κξ2))],\displaystyle\quad\quad-\frac{128\xi_{2}^{2}s_{*}^{2}\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)^{2}}{r^{(0)}(1+\xi_{2}s_{*}^{2})\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}\bigg{]}\,, (3.17)

where

AI\displaystyle A_{\rm I} =128(1+ξ2s2)(1+ξ2se2(1(86κ)ξ2))2r(0)+(1+ξ2s2(1+6κξ2))BI2s2,\displaystyle=\frac{128(1+\xi_{2}s_{*}^{2})\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)^{2}}{r^{(0)}}+\frac{\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)B_{\rm I}^{2}}{s_{*}^{2}}\,, (3.18)
BI\displaystyle B_{\rm I} =(1+ξ2se2)(se2s2)+8ξ2se2+2ξ22se2(4s2+3κ(se2s2)).\displaystyle=\left(1+\xi_{2}s_{e}^{2}\right)\left(s_{e}^{2}-s_{*}^{2}\right)+8\xi_{2}s_{e}^{2}+2\xi_{2}^{2}s_{e}^{2}\left(4s_{*}^{2}+3\kappa\left(s_{e}^{2}-s_{*}^{2}\right)\right)\,. (3.19)

The same procedure can be repeated for the tensor-to-scalar ratio. We obtain

r(I)\displaystyle r^{({\rm I})} =128[1+ξ2se2(1+6κξ2)]2(1+ξ2se21+ξ2s2)2\displaystyle=128\big{[}1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\big{]}^{2}\left(\frac{1+\xi_{2}s_{e}^{2}}{1+\xi_{2}s_{*}^{2}}\right)^{2}
×[128r(0)(1+ξ2s2)[1+ξ2se2(18ξ2+6κξ2)]2\displaystyle\quad\times\Bigg{[}\frac{128}{r^{(0)}(1+\xi_{2}s_{*}^{2})}\big{[}1+\xi_{2}s_{e}^{2}(1-8\xi_{2}+6\kappa\xi_{2})\big{]}^{2}
+1+ξ2s2(1+6κξ2)s2(1+ξ2s2)2[(1+ξ2se2(1+6κξ2))(s2se2)8ξ2se2(1+ξ2s2)]2]1.\displaystyle\quad\quad+\frac{1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})}{s_{*}^{2}(1+\xi_{2}s_{*}^{2})^{2}}\big{[}(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2}))(s_{*}^{2}-s_{e}^{2})-8\xi_{2}s_{e}^{2}(1+\xi_{2}s_{*}^{2})\big{]}^{2}\Bigg{]}^{-1}\,. (3.20)

The expressions (3.17) and (3.20) are exact so long as the slow-roll approximation holds. Together with Eq. (3.10) or Eq. (3.12), we can determine the spectral index as well as the tensor-to-scalar ratio with a set of input parameters {κ,ξ2,s,ns(0),r(0)}\{\kappa,\xi_{2},s_{*},n_{s}^{(0)},r^{(0)}\}555 To be precise, both ns(0)n_{s}^{(0)} and r(0)r^{(0)} depend on the assistant field, namely ξ2\xi_{2} and ss_{*}, as the ϕ\phi-field value at the end of inflation changes when we take into account the dynamics of the assistant field. The dependence is, however, negligible as the small nonmiminal coupling limit is taken. .

The local-type nonlinearity parameter for Class I is given by

fNL,(I)(local)\displaystyle f_{\rm NL,(I)}^{({\rm local})} =10(1+ξ2s2)3(1+ξ2se2)3AI2(1+ξ2se2(1+6κξ2))2{BI3(1+ξ2se2(1+6κξ2))3s4(1+ξ2s2)2\displaystyle=-\frac{10(1+\xi_{2}s_{*}^{2})^{3}(1+\xi_{2}s_{e}^{2})}{3A_{\rm I}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}\Bigg{\{}\frac{B_{\rm I}^{3}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{3}}{s_{*}^{4}(1+\xi_{2}s_{*}^{2})^{2}}
+256(1+ξ2se2(1+6κξ2))3(1+ξ2se2(1(86κ)ξ2))3(r(0))2(1+ξ2s2)(88ns(0)r(0))\displaystyle\quad+\frac{256\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{3}\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)^{3}}{(r^{(0)})^{2}(1+\xi_{2}s_{*}^{2})}\left(8-8n_{s}^{(0)}-r^{(0)}\right)
+2BI2(1+ξ2s2(1+6κξ2))(1+ξ2se2(1+6κξ2))3(1+ξ2se2(1(86κ)ξ2))s2(1+ξ2s2)2\displaystyle\quad+\frac{2B_{\rm I}^{2}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{3}\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)}{s_{*}^{2}(1+\xi_{2}s_{*}^{2})^{2}}
512(1+ξ2se2(1(86κ)ξ2))(r(0))2s4(1+ξ2s2)(1+ξ2se2)\displaystyle\quad-\frac{512\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)}{(r^{(0)})^{2}s_{*}^{4}(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})}
×[8ξ2s2se(1+ξ2se2(1(86κ)ξ2))116BIr(0)se(1+ξ2s2(1+6κξ2))]2\displaystyle\quad\times\Bigg{[}8\xi_{2}s_{*}^{2}s_{e}\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)-\frac{1}{16}B_{\rm I}r^{(0)}s_{e}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)\Bigg{]}^{2}
×[(1+ξ2se2(1(86κ)ξ2))2\displaystyle\quad\times\Bigg{[}\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)^{2}
+8ξ2[(1+5ξ2se2+4ξ22se4(1+6κξ2))(1+ξ2se2(1(86κ)ξ2))1+ξ2se2(1+6κξ2)\displaystyle\quad\quad+8\xi_{2}\bigg{[}\frac{\left(1+5\xi_{2}s_{e}^{2}+4\xi_{2}^{2}s_{e}^{4}(1+6\kappa\xi_{2})\right)\left(1+\xi_{2}s_{e}^{2}\left(1-(8-6\kappa)\xi_{2}\right)\right)}{1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})}
2ηe(0)ξ2se2(1+ξ2se2)(1+ξ2se2(1+6κξ2))]]}.\displaystyle\quad\quad-2\eta_{e}^{(0)}\xi_{2}s_{e}^{2}(1+\xi_{2}s_{e}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\bigg{]}\Bigg{]}\Bigg{\}}\,. (3.21)

Compared to the tensor-to-scalar ratio and the spectral index, the nonlinearity parameter contains one more parameter, namely ηe(0)\eta_{e}^{(0)}. From Eq. (3.21), we note that the prefactor associated with ηe(0)\eta_{e}^{(0)} is naturally small in the small nonminimal coupling limit. Consequently, the dependence on ηe(0)\eta_{e}^{(0)} is negligible. Thus, when presenting model-independent results, we shall consider ηe(0)=0.01\eta_{e}^{(0)}=0.01. Different choices of ηe(0)\eta_{e}^{(0)} do not lead to any visible difference.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 2: In the upper panel, we show the dependence of the spectral index ns(I)n_{s}^{({\rm I})} (right panel) and the tensor-to-scalar ratio r(I)r^{({\rm I})} (left panel) on the nonminimal coupling parameter ξ2\xi_{2} for three different values of ss_{*}, 0.001 (blue), 0.01 (red), and 0.1 (magenta), with the choice of {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\}. The maximum values of ξ2\xi_{2} are chosen in such a way that ξ2se2=0.1\xi_{2}s_{e}^{2}=0.1 so that the condition of small nonminimal couplings is met. We observe that both the spectral index and the tensor-to-scalar ratio decrease as the nonminimal coupling parameter increases; this behaviour is the general tendency as shown in Fig. 3. In the lower panel, we see the behaviour of the tensor-to-scalar ratio for r(0)=5×104r^{(0)}=5\times 10^{-4} (blue) and 10310^{-3} (red) with the choice of {ns(0),s}={0.98,0.1}\{n_{s}^{(0)},s_{*}\}=\{0.98,0.1\} (left panel) and the behaviour of the spectral index for ns(0)=0.9n_{s}^{(0)}=0.9 (blue) and 0.850.85 (red) with {r(0),s}={0.1,0.1}\{r^{(0)},s_{*}\}=\{0.1,0.1\} (right panel). When ns(0)n_{s}^{(0)} or r(0)r^{(0)} takes a rather small value with a relatively large value of ss_{*}, we observe the increasing behaviour. However, the change of tensor-to-scalar ratio is minuscule, and the increase of the spectral index is not large enough to go inside the latest Planck-BK bounds. In both the upper and lower panels, we see that the difference between the Palatini formulation (dashed) and the metric formulation (solid) is negligible as expected.

The upper panel of Fig. 2 shows the dependence of the spectral index ns(I)n_{s}^{({\rm I})} (right panel) and the tensor-to-scalar ratio r(I)r^{({\rm I})} (left panel) on ξ2\xi_{2} for three different ss_{*} values with the choice of {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\}. The ranges of ξ2\xi_{2} are chosen such that ξ2se2\xi_{2}s_{e}^{2} satisfies the condition ξ2se20.1\xi_{2}s_{e}^{2}\leq 0.1 in order that the analytic formulae work well. One may observe that both the spectral index and the tensor-to-scalar ratio decrease as the nonminimal coupling parameter ξ2\xi_{2} increases. The lower panel of Fig. 2 presents the behaviour of the tensor-to-scalar ratio (left panel) and the spectral index (right panel) for rather small values of r(0)r^{(0)} and ns(0)n_{s}^{(0)}, respectively, with a relatively large value of s=0.1s_{*}=0.1. In this case, we observe the opposite behaviour; both the spectral index and the tensor-to-scalar ratio increase as the nonminimal coupling parameter ξ2\xi_{2} increases. However, the change of the tensor-to-scalar ratio is minuscule, and the increase of the spectral index is not large enough to be allowed by the latest Planck-BK bounds. For all the cases presented, the difference between the Palatini and metric formulations is negligible as the nonminimal coupling is small. For most of the parameter space, the general tendency is that both the spectral index as well as the tensor-to-scalar ratio become lowered as we increase ξ2\xi_{2}. Therefore, we may bring models that originally predict large {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\} to the observationally-favoured region as demonstrated in Fig. 3 on the nsn_{s}rr plane.

Refer to caption
Figure 3: Behaviour of the spectral index ns(I)n_{s}^{({\rm I})} and the tensor-to-scalar ratio r(I)r^{({\rm I})} for {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} (red), {0.95,0.05}\{0.95,0.05\} (blue), {1.05,0.06}\{1.05,0.06\} (cyan), and {1,0.03}\{1,0.03\} (magenta) with s=103s_{*}=10^{-3} (solid) and 10110^{-1} (dashed). Only the metric cases are presented here as there is little difference between the Palatini and the metric formulations. The latest Planck-BK 1-sigma (2-sigma) bound is depicted as the black solid (dashed) line, while the Planck-only 1-sigma (2-sigma) bound is presented with the grey solid (dashed) line. Both the spectral index and the tensor-to-scalar ratio decrease as the nonminimal coupling parameter ξ2\xi_{2} takes a larger value as indicated by arrows; the maximum value of ξ2\xi_{2} is chosen such that ξ2se2=0.1\xi_{2}s_{e}^{2}=0.1, i.e., until the small nonminimal coupling condition is valid. Thus, one may save some otherwise ruled-out models; models that originally predict large spectral index and/or tensor-to-scalar ratio can be brought to the observationally-favoured region.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 4: In the upper panel, the evolution of the nonlinearity parameter in terms of the nonminimal coupling parameter is shown for Class I. Two values of ss_{*} are considered, 0.001 (left panel) and 0.1 (right panel) with various choices of {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} (red), {0.95,0.05}\{0.95,0.05\} (blue), {1.05,0.06}\{1.05,0.06\} (cyan), and {1.0,0.03}\{1.0,0.03\} (magenta). We observe that the nonlinearity parameter initially shows an increasing behaviour. In the lower panel, we present the predictions in the fNL,(I)(local)f_{\rm NL,(I)}^{({\rm local})}ns(I)n_{s}^{({\rm I})} plane. The shaded region corresponds to the latest bounds on the spectral index. In the case of s=0.001s_{*}=0.001, the nonlinearity parameter for the {ns(0),r(0)}={1.05,0.06}\{n_{s}^{(0)},r^{(0)}\}=\{1.05,0.06\} case becomes slightly larger than the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. In the case of s=0.1s_{*}=0.1, however, all the constraints are safely satisfied. For all cases, ηe(0)=0.01\eta_{e}^{(0)}=0.01 is chosen. However, the dependence on ηe(0)\eta_{e}^{(0)} is weak as the prefactor is small in the small nonminimal coupling limit. Only the metric formulation is presented as there is little difference between the metric and the Palatini formulations.

The behaviour of the nonlinearity parameter as we increase the nonminimal coupling parameter ξ2\xi_{2} is shown in the upper panel of Fig. 4. There is little difference between the metric and the Palatini formulations, and thus, only the metric formulation is presented. Four different choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\} are considered, {0.98,0.1}\{0.98,0.1\}, {0.95,0.05}\{0.95,0.05\}, {1.05,0.06}\{1.05,0.06\}, and {1.0,0.03}\{1.0,0.03\}, and two values of ss_{*} are chosen, 0.001 and 0.1. We observe that fNL(local)f_{\rm NL}^{({\rm local})} initially increases as we increase ξ2\xi_{2}, and then it decreases. Moreover, we note that the fNL(local)f_{\rm NL}^{({\rm local})} value becomes larger when the ss_{*} value is small which is an expected behaviour from the analytical expression (3.21). In the lower panel of Fig. 4, the predictions are shown in the fNL(local)f_{\rm NL}^{({\rm local})}nsn_{s} plane together with the observationally-allowed bound for the spectral index. For the case of s=0.001s_{*}=0.001, the nonlinearity parameter for the {ns(0),r(0)}={1.05,0.06}\{n_{s}^{(0)},r^{(0)}\}=\{1.05,0.06\} case becomes slightly larger than the Planck 2-sigma bound [60], 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. In the case of s=0.1s_{*}=0.1, however, all the constraints are safely satisfied. It is interesting to note that constraints on fNL(local)f_{\rm NL}^{\rm(local)} can be improved as |fNL(local)|<𝒪(0.1)𝒪(1)|f_{\rm NL}^{\rm(local)}|<{\cal O}(0.1)-{\cal O}(1) in future galaxy surveys and 21 cm line of neutral hydrogen observations [61, 62, 63, 64, 65, 66, 67, 68, 69]. Hence, the currently allowed parameter range can even further be probed by such future observations of non-Gaussianity.

3.2 Class II: End of inflation via a separate sector

The second class we consider includes scenarios where end of inflation is provided by a separate sector other than the ϕ\phi and ss fields. Examples include the power-law inflation with an exponential potential [4] and the hybrid inflation model [5]. In the case of power-law inflation, the slow-roll parameters in the original ϕ\phi-field model become constant. Thus, inflation does not end via the standard slow-roll violation, and a separate end-of-inflation mechanism is needed. In the hybrid inflation model, the end of inflation is achieved by the so-called waterfall field666 There exists a parameter space where end of inflation via the breaking of the slow-roll condition is possible. This case falls into Class I discussed in Sec. 3.1. . In the presence of the assistant field ss, one may use the ss field to end inflation via the slow-roll violation, ϵe1\epsilon_{e}\approx 1, in which case Class I applies777 When the assistant field is responsible for end of inflation, ss_{*} (or ξ2\xi_{2}) becomes a given quantity rather than a free, input parameter. . The violation of slow-roll via a nonminimal coupling can also be realised in some single-field nonminimal inflation models [23]. In this subsection, we do not assume that the assistant field is responsible for end of inflation. Instead, we assume that there is a separate mechanism that ends inflation without substantially modifying the observables.

As we do not have the condition ϵe1\epsilon_{e}\approx 1, we are no longer able to replace ϵeϕ\epsilon_{e}^{\phi}. In other words, information on the ϕ\phi field at the end of inflation or, equivalently, ϵe(0)\epsilon_{e}^{(0)}, is required. One exception is when the ϕ\phi-field slow-roll parameters are approximately constant, as in the power-law inflation case; in this special case, we have ϵe(0)ϵ(0)=r(0)/16\epsilon_{e}^{(0)}\approx\epsilon_{*}^{(0)}=r^{(0)}/16, and thus, we do not need any information on the ϕ\phi-field value at the end of inflation. We aim to present a general analysis without relying on a specific model. Therefore, we remain agnostic with respect to the original ϕ\phi-field model and treat ϵe(0)\epsilon_{e}^{(0)} as an independent parameter. Choosing a specific model for the ϕ\phi field would correspond to choosing a specific value for ϵe(0)\epsilon_{e}^{(0)}.

From Eq. (3.1), we obtain the spectral index as follows:

ns(II)\displaystyle n_{s}^{({\rm II})} =1(1+ξ2s2)r(0)816ξ22s21+ξ2s2(1+6κξ2)\displaystyle=1-\frac{(1+\xi_{2}s_{*}^{2})r^{(0)}}{8}-\frac{16\xi_{2}^{2}s_{*}^{2}}{1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})}
32(1+ξ2se2)2AII(1+ξ2se2(1+6κξ2))2[8ξ22se2+ϵe(0)(1+ξ2se2)(1+ξ2se2(1+6κξ2))]2\displaystyle\quad-\frac{32(1+\xi_{2}s_{e}^{2})^{2}}{A_{\rm II}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}\Big{[}8\xi_{2}^{2}s_{e}^{2}+\epsilon^{(0)}_{e}(1+\xi_{2}s_{e}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\Big{]}^{2}
+r(0)ξ2(1+ξ2s2)(1+ξ2se2)2(1ξ22s4(1+6κξ2))12AIIs2(1+ξ2s2(1+6κξ2))\displaystyle\quad+\frac{r^{(0)}\xi_{2}(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})^{2}\left(1-\xi_{2}^{2}s_{*}^{4}(1+6\kappa\xi_{2})\right)}{12A_{\rm II}s_{*}^{2}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}
×[BII1+ξ2se2(1+6κξ2)128ϵe(0)ξ2s2r(0)(1+ξ2s2(1+6κξ2))]2\displaystyle\quad\quad\times\left[\frac{B_{\rm II}}{1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})}-\frac{128\epsilon^{(0)}_{e}\xi_{2}s_{*}^{2}}{r^{(0)}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}\right]^{2}
16(1+ξ2s2)2AII[(1+ξ2se2)2r(0)(ϵe(0))2(88ns(0)3r(0))\displaystyle\quad-\frac{16(1+\xi_{2}s_{*}^{2})^{2}}{A_{\rm II}}\Bigg{[}\frac{(1+\xi_{2}s_{e}^{2})^{2}}{r^{(0)}}\left(\epsilon^{(0)}_{e}\right)^{2}\left(8-8n_{s}^{(0)}-3r^{(0)}\right)
3BIIϵe(0)ξ2(1+ξ2se2)2(1+ξ2s2)(1+ξ2se2(1+6κξ2))64ξ22s2(1+ξ2se2)2r(0)(1+ξ2s2)(1+ξ2s2(1+6κξ2))(ϵe(0))2\displaystyle\quad\quad-\frac{3B_{\rm II}\epsilon^{(0)}_{e}\xi_{2}(1+\xi_{2}s_{e}^{2})^{2}}{(1+\xi_{2}s_{*}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)}-\frac{64\xi_{2}^{2}s_{*}^{2}(1+\xi_{2}s_{e}^{2})^{2}}{r^{(0)}(1+\xi_{2}s_{*}^{2})\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}\left(\epsilon^{(0)}_{e}\right)^{2}
+BII2ξ2(1+ξ2se2)2(1ξ2s2(3+4ξ2s2(1+6κξ2)))2s2(1+ξ2s2)2(1+ξ2s2(1+6κξ2))(1+ξ2se2(1+6κξ2))2],\displaystyle\quad\quad+\frac{B_{\rm II}^{2}\xi_{2}(1+\xi_{2}s_{e}^{2})^{2}\left(1-\xi_{2}s_{*}^{2}(3+4\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2}))\right)}{2s_{*}^{2}(1+\xi_{2}s_{*}^{2})^{2}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}\Bigg{]}\,, (3.22)

where

AII\displaystyle A_{\rm II} =128(1+ξ2s2)(1+ξ2se2)2r(0)(ϵe(0))2+(1+ξ2se2)2(1+ξ2s2(1+6κξ2))BII2s2(1+ξ2se2(1+6κξ2))2,\displaystyle=\frac{128(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})^{2}}{r^{(0)}}\left(\epsilon^{(0)}_{e}\right)^{2}+\frac{(1+\xi_{2}s_{e}^{2})^{2}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)B_{\rm II}^{2}}{s_{*}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}\,, (3.23)
BII\displaystyle B_{\rm II} =8ξ2se2+ϵe(0)(se2s2)(1+ξ2se2(1+6κξ2)).\displaystyle=8\xi_{2}s_{e}^{2}+\epsilon^{(0)}_{e}\left(s_{e}^{2}-s_{*}^{2}\right)\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\,. (3.24)

Similarly, we find the tensor-to-scalar ratio as

r(II)\displaystyle r^{({\rm II})} =128(1+ξ2se2)2(8ξ22se2+ϵe(0)(1+ξ2se2)(1+ξ2se2(1+6κξ2)))2(1+ξ2se2(1+6κξ2))2\displaystyle=\frac{128(1+\xi_{2}s_{e}^{2})^{2}\left(8\xi_{2}^{2}s_{e}^{2}+\epsilon^{(0)}_{e}(1+\xi_{2}s_{e}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\right)^{2}}{\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}
×[(1+ξ2se2)2(1+ξ2s2(1+6κξ2))s2(1+ξ2se2(1+6κξ2))2[8ξ2se2+ϵe(0)(se2s2)(1+ξ2se2(1+6κξ2))]2\displaystyle\quad\times\Bigg{[}\frac{(1+\xi_{2}s_{e}^{2})^{2}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}{s_{*}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}\Big{[}8\xi_{2}s_{e}^{2}+\epsilon^{(0)}_{e}(s_{e}^{2}-s_{*}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\Big{]}^{2}
+128(1+ξ2s2)(1+ξ2se2)2r(0)(ϵe(0))2]1.\displaystyle\quad\quad+\frac{128(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})^{2}}{r^{(0)}}\left(\epsilon^{(0)}_{e}\right)^{2}\Bigg{]}^{-1}\,. (3.25)

Note that these expressions are exact so long as the slow-roll approximation holds. Together with the solution for ses_{e}, Eq. (3.10) or Eq. (3.12), the spectral index (3.22) as well as the tensor-to-scalar ratio (3.25) can be computed with a set of input parameters {κ,ξ2,s,ϵe(0),ns(0),r(0)}\{\kappa,\xi_{2},s_{*},\epsilon_{e}^{(0)},n_{s}^{(0)},r^{(0)}\}. Class II is more general than Class I in the sense that the expressions for Class I, Eq. (3.17) and Eq. (3.20), can be obtained by taking ϵe(0)=(1ϵeσ)e2be\epsilon_{e}^{(0)}=(1-\epsilon_{e}^{\sigma})e^{2b_{e}}. The special case of the constant slow-roll parameter can be captured by setting ϵe(0)=r(0)/16\epsilon_{e}^{(0)}=r^{(0)}/16.

The local-type nonlinearity parameter for Class II is given by

fNL,(II)(local)\displaystyle f_{\rm NL,(II)}^{({\rm local})} =10(1+ξ2s2)(1+ξ2se2)3AII2(ϵe(0)(1+ξ2se2)+8ξ22se21+ξ2se2(1+6κξ2))\displaystyle=-\frac{10(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})}{3A_{\rm II}^{2}}\left(\epsilon_{e}^{(0)}(1+\xi_{2}s_{e}^{2})+\frac{8\xi_{2}^{2}s_{e}^{2}}{1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})}\right)
×[256(ϵe(0))3(1+ξ2s2)(1+ξ2se2)3(r(0))2(88ns(0)r(0))\displaystyle\quad\times\Bigg{[}\frac{256(\epsilon_{e}^{(0)})^{3}(1+\xi_{2}s_{*}^{2})(1+\xi_{2}s_{e}^{2})^{3}}{(r^{(0)})^{2}}\left(8-8n_{s}^{(0)}-r^{(0)}\right)
+2BII2ϵe(0)(1+ξ2se2)3(1+ξ2s2(1+6κξ2))s2(1+ξ2se2(1+6κξ2))2+BII3(1+ξ2se2)3s4(1+ξ2se2(1+6κξ2))3]\displaystyle\quad\quad+\frac{2B_{\rm II}^{2}\epsilon_{e}^{(0)}(1+\xi_{2}s_{e}^{2})^{3}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}{s_{*}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}}+\frac{B_{\rm II}^{3}(1+\xi_{2}s_{e}^{2})^{3}}{s_{*}^{4}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{3}}\Bigg{]}
+20ϵe(0)se2(1+ξ2s2)2(1+ξ2se2)43AII2(1+ξ2se2(1+6κξ2))3(ϵe(0)(1+ξ2se2)+8ξ22se21+ξ2se2(1+6κξ2))1\displaystyle\quad+\frac{20\epsilon_{e}^{(0)}s_{e}^{2}(1+\xi_{2}s_{*}^{2})^{2}(1+\xi_{2}s_{e}^{2})^{4}}{3A_{\rm II}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{3}}\left(\epsilon_{e}^{(0)}(1+\xi_{2}s_{e}^{2})+\frac{8\xi_{2}^{2}s_{e}^{2}}{1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})}\right)^{-1}
×(BII(1+ξ2s2(1+6κξ2))s2(1+ξ2se2(1+6κξ2))128ξ2ϵe(0)r(0))2[(ϵe(0))2(1+ξ2se2)(1+ξ2se2(1+6κξ2))2\displaystyle\quad\times\left(\frac{B_{\rm II}\left(1+\xi_{2}s_{*}^{2}(1+6\kappa\xi_{2})\right)}{s_{*}^{2}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)}-\frac{128\xi_{2}\epsilon_{e}^{(0)}}{r^{(0)}}\right)^{2}\Big{[}(\epsilon_{e}^{(0)})^{2}(1+\xi_{2}s_{e}^{2})\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)^{2}
+8ξ2ϵe(0)(1+5ξ2se2+4ξ22se4(1+6κξ2))16ξ22se2ηe(0)(1+ξ2se2(1+6κξ2))].\displaystyle\quad\quad+8\xi_{2}\epsilon_{e}^{(0)}\left(1+5\xi_{2}s_{e}^{2}+4\xi_{2}^{2}s_{e}^{4}(1+6\kappa\xi_{2})\right)-16\xi_{2}^{2}s_{e}^{2}\eta_{e}^{(0)}\left(1+\xi_{2}s_{e}^{2}(1+6\kappa\xi_{2})\right)\Big{]}\,. (3.26)

Similar to the case of Class I, the nonlinearity parameter contains one more parameter compared to the tensor-to-scalar ratio and the spectral index, namely ηe(0)\eta_{e}^{(0)}. As the prefactor associated with ηe(0)\eta_{e}^{(0)} is naturally small in the small nonminimal coupling limit, the dependence on ηe(0)\eta_{e}^{(0)} is negligible. Thus, when presenting model-independent results, we shall consider ηe(0)=0.01\eta_{e}^{(0)}=0.01. Different choices of ηe(0)\eta_{e}^{(0)} do not lead to any visible difference.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 5: In the upper panel, we show the dependence of the spectral index ns(II)n_{s}^{({\rm II})} and the tensor-to-scalar ratio r(II)r^{({\rm II})} on the nonminimal coupling parameter ξ2\xi_{2} for different choices of {s,ee(0)}\{s_{*},e_{e}^{(0)}\} with {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\}. The maximum values of ξ2\xi_{2} are given via the small nonminimal coupling condition, ξ2se2=0.1\xi_{2}s_{e}^{2}=0.1. We observe that, similar to Class I, both the spectral index and the tensor-to-scalar ratio decrease as ξ2\xi_{2} increases. Also similar to Class I, the difference between the Palatini formulation (dashed) and the metric formulation (solid) is negligible. In the lower panel, we see the behaviour of the tensor-to-scalar ratio for two different values of r(0)r^{(0)} and three different values of ϵe(0)\epsilon_{e}^{(0)} with {ns(0),s}={0.98,0.1}\{n_{s}^{(0)},s_{*}\}=\{0.98,0.1\} (left panel) and the behaviour of the spectral index for different choices of {ns(0),ϵe(0)}\{n_{s}^{(0)},\epsilon_{e}^{(0)}\} with {r(0),s}={0.1,0.1}\{r^{(0)},s_{*}\}=\{0.1,0.1\} (right panel). As the difference between the Palatini and the metric formulations is negligible, we only present the metric case in the lower panel. When ϵe(0)\epsilon_{e}^{(0)} is chosen to be relatively large, the behaviour is similar to Class I. On the other hand, when ϵe(0)\epsilon_{e}^{(0)} takes a very small value, the tensor-to-scalar ratio decreases initially and then increases again. Furthermore, unlike Class I, the increase of the spectral index is strong enough to go inside the allowed bounds. The same behaviour is observed for different sets of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\}; see also Fig. 6.

In the upper panel of Fig. 5, we present how the spectral index ns(II)n_{s}^{({\rm II})} and the tensor-to-scalar ratio r(II)r^{({\rm II})} behave as we change ξ2\xi_{2} for different sets of {s,ϵe(0)}\{s_{*},\epsilon_{e}^{(0)}\} with a specific choice of {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\}. The maximum values of ξ2\xi_{2} are chosen such that the small nonminimal coupling condition holds, i.e., ξ2se2=0.1\xi_{2}s_{e}^{2}=0.1. The result is similar to Class I; both the spectral index and the tensor-to-scalar ratio decrease as the nonminimal coupling parameter ξ2\xi_{2} increases. Furthermore, we see that the difference between the Palatini and metric formulations is negligible. Thus, models that originally predict large {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\} and that belong to Class II can as well be brought to the observationally-favoured region. In the lower panel of Fig. 5, we show the behaviour of the tensor-to-scalar ratio (left panel) and the spectral index (right panel) for rather small r(0)r^{(0)} and ns(0)n_{s}^{(0)} values, respectively, when ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2}, 5×1045\times 10^{-4}, and 10510^{-5}. When ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2}, the results are similar to Class I; the tensor-to-scalar ratio as well as the spectral index show an increasing behaviour, but the change is not significant. On the other hand, when ϵe(0)\epsilon_{e}^{(0)} takes a very small value, the tensor-to-scalar ratio initially decreases and then increases again. The spectral index initially increases and then decreases again, which remains to be similar to Class I, but, unlike Class I, the increase of the spectral index is strong enough to bring the prediction to the allowed bounds. This is a stark difference between Class I and Class II: Not only can models that originally predict large spectral index and/or tensor-to-scalar ratio be brought to the observationally-favoured region, but models that predict small spectral index values may also be revived and become compatible with the latest Planck-BK bounds in the presence of the assistant field. These behaviours are observed for different choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\} as demonstrated in Fig. 6 on the nsn_{s}rr plane.

Refer to caption
Refer to caption
Figure 6: Left panel: Behaviour of the spectral index ns(II)n_{s}^{({\rm II})} and the tensor-to-scalar ratio r(II)r^{({\rm II})} for {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} (red), {0.95,0.05}\{0.95,0.05\} (blue), {1.05,0.06}\{1.05,0.06\} (cyan), and {1,0.03}\{1,0.03\} (magenta) for s=0.001s_{*}=0.001 (solid) and 0.10.1 (dashed) with ϵe(0)=0.01\epsilon_{e}^{(0)}=0.01. Right panel: Behaviour of ns(II)n_{s}^{({\rm II})} and r(II)r^{({\rm II})} for {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} (red), {0.95,0.05}\{0.95,0.05\} (blue), {1.05,0.06}\{1.05,0.06\} (cyan), {1.02,0.001}\{1.02,0.001\} (magenta), and {0.92,0.002}\{0.92,0.002\} (brown) for ϵe(0)=105\epsilon_{e}^{(0)}=10^{-5} (solid) and 10310^{-3} (dashed) with s=0.1s_{*}=0.1. As the difference between the Palatini and the metric formulations is negligible, only the metric cases are presented. The latest Planck-BK 1-sigma (2-sigma) bound is depicted as the black solid (dashed) line, while the Planck-only 1-sigma (2-sigma) bound is presented with the grey solid (dashed) line. When ϵe(0)\epsilon_{e}^{(0)} is relatively large, both the spectral index and the tensor-to-scalar ratio decrease as the nonminimal coupling parameter ξ2\xi_{2} takes a larger value as indicated by arrows. Thus, similar to Class I, models that originally predict large spectral index and/or tensor-to-scalar ratio can be brought to the observationally-favoured region. When ϵe(0)\epsilon_{e}^{(0)} takes a very small value, the spectral index initially shows an increasing behaviour for small ns(0)n_{s}^{(0)} cases, thereby entering the observationally-allowed bounds. Thus, unlike Class I, models that originally predict small spectral index may as well become revived for small ϵe(0)\epsilon_{e}^{(0)} values.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 7: In the upper panel, the evolution of the nonlinearity parameter in terms of the nonminimal coupling parameter is shown for Class II. Two values of ss_{*} are considered, 0.001 (left panel) and 0.1 (right panel) with various choices of {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} (red), {0.95,0.05}\{0.95,0.05\} (blue), {1.05,0.06}\{1.05,0.06\} (cyan), and {1.0,0.03}\{1.0,0.03\} (magenta). Furthermore, two values of ϵe(0)\epsilon_{e}^{(0)} are considered, 10210^{-2} (solid) and 10510^{-5} (dashed). Similar to Class I, we observe that the nonlinearity parameter initially shows an increasing behaviour, and then it decreases. In the lower panel, we present the predictions in the fNL,(II)(local)f_{\rm NL,(II)}^{({\rm local})}ns(II)n_{s}^{({\rm II})} plane. The shaded region corresponds to the latest bounds on the spectral index. In the case of s=0.001s_{*}=0.001 with ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2}, the nonlinearity parameter for the {ns(0),r(0)}={1.05,0.06}\{n_{s}^{(0)},r^{(0)}\}=\{1.05,0.06\} case becomes slightly larger than the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. When ϵe(0)=105\epsilon_{e}^{(0)}=10^{-5}, only a narrow parameter range is allowed for the {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} case, while the rest cases are ruled out. In the case of s=0.1s_{*}=0.1, however, all the constraints are safely satisfied for both ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2} and 10510^{-5}. For all cases, ηe(0)=0.01\eta_{e}^{(0)}=0.01 is chosen. However, the dependence on ηe(0)\eta_{e}^{(0)} is weak as the prefactor is small in the small nonminimal coupling limit. Only the metric formulation is presented as there is little difference between the metric and the Palatini formulations.

The behaviour of the nonlinearity parameter as we increase the nonminimal coupling parameter ξ2\xi_{2} is shown in the upper panel of Fig. 7. Only the metric formulation is presented as the difference between the metric and the Palatini formulations is negligible. Four different choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\} are considered, {0.98,0.1}\{0.98,0.1\}, {0.95,0.05}\{0.95,0.05\}, {1.05,0.06}\{1.05,0.06\}, and {1.0,0.03}\{1.0,0.03\}, and two values of ss_{*} are chosen, 0.001 and 0.1. We additionally considered two scenarios of ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2} and ϵe(0)=105\epsilon_{e}^{(0)}=10^{-5}. Similar to Class I, we observe that fNL(local)f_{\rm NL}^{({\rm local})} initially increases as we increase ξ2\xi_{2}, and then it decreases. Moreover, the fNL(local)f_{\rm NL}^{({\rm local})} value tends to be larger when the ss_{*} value is small. In the lower panel of Fig. 7, the predictions are shown in the fNL(local)f_{\rm NL}^{({\rm local})}nsn_{s} plane. In the case of s=0.001s_{*}=0.001 with ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2}, the nonlinearity parameter for the {ns(0),r(0)}={1.05,0.06}\{n_{s}^{(0)},r^{(0)}\}=\{1.05,0.06\} case becomes slightly larger than the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. When ϵe(0)=105\epsilon_{e}^{(0)}=10^{-5}, only a narrow parameter range is allowed for the {ns(0),r(0)}={0.98,0.1}\{n_{s}^{(0)},r^{(0)}\}=\{0.98,0.1\} case, while the rest cases are ruled out. In the case of s=0.1s_{*}=0.1, however, all the constraints are safely satisfied for both ϵe(0)=102\epsilon_{e}^{(0)}=10^{-2} and 10510^{-5}. The currently allowed parameter range can even further be probed by future observations of non-Gaussianity such as galaxy surveys and 21 cm line of neutral hydrogen observations [61, 62, 63, 64, 65, 66, 67, 68, 69] as the constraints on fNL(local)f_{\rm NL}^{\rm(local)} can be improved as |fNL(local)|<𝒪(0.1)𝒪(1)|f_{\rm NL}^{\rm(local)}|<{\cal O}(0.1)-{\cal O}(1).

4 Examples

As an application of our general analysis presented in Sec. 3, we consider three examples: the loop inflation model which belongs to Class I888 Chaotic inflation model is also categorised as Class I, and it has already been studied in Ref. [34]. and the power-law inflation and hybrid inflation models which belong to Class II. These models are ruled out by the latest Planck-BK results. We show that, with the help of the assistant field, these three models may become compatible with the latest observations. When presenting results, we only choose the metric formulation as there exists little difference between the metric and the Palatini formulations.

4.1 Loop inflation – an example for Class I

Refer to caption
Figure 8: Predictions for the loop inflation model. The original model, depicted by the magenta line, is ruled out by the latest Planck [1] (grey) and Planck-BK [2] (black) bounds, due to the largeness of the spectral index. By varying the nonminimal coupling parameter ξ2\xi_{2}, the effect of the assistant field is presented for α=0.1\alpha=0.1 and α=1\alpha=1 with two values of ss_{*}, 0.0010.001 (solid) and 0.10.1 (dashed). The cut in the case of {α,s}={0.1,0.1}\{\alpha,s_{*}\}=\{0.1,0.1\} is due to the fact that the maximum value of ξ2\xi_{2} is chosen such that ξ2se2=0.1\xi_{2}s_{e}^{2}=0.1. In the presence of the assistant field, the spectral index gets suppressed. Thus, the model becomes compatible with the latest observational constraints.

Loop inflation is described by the following potential:

VJ(ϕ)=Λ4(1+αlogϕ),\displaystyle V_{\rm J}(\phi)=\Lambda^{4}(1+\alpha\log\phi)\,, (4.1)

where α\alpha is the coefficient of one-loop correction term [6]. As we are interested in phenomenological aspects of the model, we take α\alpha as a free parameter rather than focusing on its origin. For the potential (4.1), we obtain the slow-roll parameters (3.14) as

ϵ(0)=α22ϕ(1+αlogϕ)2,η(0)=αϕ2(1+αlogϕ).\displaystyle\epsilon^{(0)}=\frac{\alpha^{2}}{2\phi(1+\alpha\log\phi)^{2}}\,,\quad\eta^{(0)}=-\frac{\alpha}{\phi^{2}(1+\alpha\log\phi)}\,. (4.2)

The field value of the inflaton ϕ\phi at the end of inflation, ϕe\phi_{e}, characterised by ϵ(0)(ϕ=ϕe)=1\epsilon^{(0)}(\phi=\phi_{e})=1, is then given by

ϕe=12[W0(e1/α2)]1,\displaystyle\phi_{e}=\frac{1}{\sqrt{2}}\left[W_{0}\left(\frac{e^{1/\alpha}}{\sqrt{2}}\right)\right]^{-1}\,, (4.3)

where W0W_{0} is the zero-branch of Lambert’s WW function, and we have assumed α>0\alpha>0. Since end of inflation may be achieved via the slow-roll violation, the loop inflation model belongs to Class I. Imposing 60 ee-folds, i.e., 60=e(VJ/VJ,ϕ)𝑑ϕ60=\int_{e}^{*}(V_{\rm J}/V_{{\rm J},\phi})d\phi, we obtain the field value at the pivot scale. We depict the original prediction for the spectral index ns(0)n_{s}^{(0)} and the tensor-to-scalar ratio r(0)r^{(0)} in Fig. 8 by varying α\alpha (magenta line). The original loop inflation model is clearly ruled out by the latest observational constraints.

Using the analytical expressions (3.17) and (3.20), the effect of the assistant field ss is presented in Fig. 8 for α=0.1\alpha=0.1 and α=1\alpha=1 with two values of ss_{*}, 0.0010.001 (solid line) and 0.10.1 (dashed line). We see that the spectral index nsn_{s} decreases as the nonminimal coupling parameter ξ2\xi_{2} increases. Hence, we may rescue the original loop inflation model. For example, in the case of α=0.1\alpha=0.1, the model becomes compatible with the Planck-BK 2-sigma bounds for the range of ξ2=(2.12.3)×102\xi_{2}=(2.1-2.3)\times 10^{-2} for s=0.001s_{*}=0.001 and ξ2=(1.21.4)×102\xi_{2}=(1.2-1.4)\times 10^{-2} for s=0.1s_{*}=0.1. In the case of α=1\alpha=1, the corresponding ranges for the nonminimal coupling parameter are ξ2=(1.92.1)×102\xi_{2}=(1.9-2.1)\times 10^{-2} for s=0.001s_{*}=0.001 and ξ2=(1.11.3)×102\xi_{2}=(1.1-1.3)\times 10^{-2} for s=0.1s_{*}=0.1. For each case, we present in Fig. 9 the field trajectory during inflation that yields ns=0.965n_{s}=0.965.

Refer to caption
Figure 9: Field trajectories for the loop inflation model. Two choices of ss_{*}, 0.001 (solid) and 0.1 (dashed), are considered with α=0.1\alpha=0.1 (red) and 11 (blue). For each case, we have chosen the value of ξ2\xi_{2} in such a way that the spectral index becomes 0.965.
Refer to caption
Figure 10: Nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)} in terms of the spectral index nsn_{s} for the loop inflation model. Two choices of ss_{*} are considered, 0.001 (left) and 0.1 (right), and the cases with α=0.1\alpha=0.1 and 11 are shown in each panel. The shaded region represents the Planck 2-sigma bound on the spectral index, ns=[0.955,0.976]n_{s}=[0.955,0.976]. We note that, for both the s=0.001s_{*}=0.001 and s=0.1s_{*}=0.1 cases, the nonlinearity parameter is well within the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. We further observe that, for a relatively large value of ss_{*} such as the s=0.1s_{*}=0.1 case, the nonlinearity parameter tends to be tiny.

Figure 10 shows the nonlinearity parameter fNL(local)f_{\rm NL}^{(\rm local)}, obtained by using Eq. (3.21), as a function of the spectral index nsn_{s} for the loop inflation model. Similar to Fig. 8, two choices of ss_{*}, 0.001 and 0.1, are considered. The shaded grey region indicates the Planck 2-sigma bound on the spectral index. We note that, for both the s=0.001s_{*}=0.001 and s=0.1s_{*}=0.1 cases, the nonlinearity parameter is well within the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. We further observe that, for a relatively large value of ss_{*} such as the s=0.1s_{*}=0.1 case, the nonlinearity parameter tends to be tiny.

Refer to caption
Refer to caption
Figure 11: Running of the spectral index αs\alpha_{s} in terms of the spectral index nsn_{s} for the loop inflation model. Similar to Fig. 10, two choices of ss_{*} are considered with α=0.1\alpha=0.1 and 11 for both cases. We observe that αs\alpha_{s} shows an increasing behaviour as nsn_{s} decreases, and that a larger value of ss_{*} results in smaller αs\alpha_{s}. See also Appendix A.

Finally, we present the prediction for the running of the spectral index αs\alpha_{s}, discussed in Appendix A, in Fig. 11. We have considered two choices of ss_{*} and two values for the parameter α\alpha as before. We observe that the running of the spectral index shows an increasing behaviour as the spectral index decreases. We further see that the running of the spectral index tends to be smaller when a larger value is considered for ss_{*}.

4.2 Power-law inflation – an example for Class II

Refer to caption
Figure 12: Predictions for the power-law inflation model. The original model, depicted by the magenta line, is ruled out by the latest Planck (grey) and Planck-BK (black) bounds. By varying the nonminimal coupling parameter ξ2\xi_{2}, the effect of the assistant field is presented for λ=0.1\lambda=0.1, 0.05, and 0.01 with two choices of ss_{*}, 0.001 (solid) and 0.1 (dashed). In the presence of the assistant field, the spectral index as well as the tensor-to-scalar ratio get suppressed. Thus, the power-law inflation model becomes compatible with the latest observational constraints.

Power-law inflation [4] is described by the potential

VJ(ϕ)=Λ4exp(λϕ).\displaystyle V_{\rm J}(\phi)=\Lambda^{4}\exp(-\lambda\phi)\,. (4.4)

As the potential is given by an exponential function, the slow-roll parameters (3.14) become constants,

ϵ(0)=12λ2,η(0)=λ2.\displaystyle\epsilon^{(0)}=\frac{1}{2}\lambda^{2}\,,\quad\eta^{(0)}=\lambda^{2}\,. (4.5)

Therefore, inflation does not end via slow-roll violations, calling for a separate sector for end of inflation. The power-law inflation model thus belongs to Class II with a special property of ϵe(0)=ϵ(0)\epsilon_{e}^{(0)}=\epsilon_{*}^{(0)}. The spectral index and the tensor-to-scalar ratio are given by

ns(0)=1λ2,r(0)=8λ2.\displaystyle n_{s}^{(0)}=1-\lambda^{2}\,,\quad r^{(0)}=8\lambda^{2}\,. (4.6)

The original predictions, ns(0)n_{s}^{(0)} and r(0)r^{(0)}, are depicted in Fig. 12 as a magenta line for various values of λ\lambda. We see that the original model is ruled out by the latest Planck-BK bounds as either the spectral index or the tensor-to-scalar ratio is too large.

From the general analysis performed in Sec. 3.2, one expects that the inclusion of the assistant field may bring the original predictions to the observationally-favoured region; see, e.g., Fig. 6. Using the analytical expressions (3.22) and (3.25), the spectral index and the tensor-to-scalar ratio in the presence of the assistant field can be computed. In Fig. 12, we present the predictions in the nsn_{s}rr plane for λ=0.1\lambda=0.1, 0.05, and 0.01 with two choices of ss_{*}, 0.1 and 0.001, by varying the nonminimal coupling parameter ξ2\xi_{2}. With the help of the assistant field, the power-law inflation model may thus become compatible again with the latest Planck-BK constraints. For example, in the case of λ=0.01\lambda=0.01, the predictions for the spectral index and the tensor-to-scalar ratio enter the Planck-BK 22-sigma bound for the range of ξ2=(1.051.11)×102\xi_{2}=(1.05-1.11)\times 10^{-2} for s=0.001s_{*}=0.001 and ξ2=(4.55.6)×103\xi_{2}=(4.5-5.6)\times 10^{-3} for s=0.1s_{*}=0.1. In the case of λ=0.05\lambda=0.05, the corresponding range for the nonminimal coupling parameter is ξ2=(1.291.36)×102\xi_{2}=(1.29-1.36)\times 10^{-2} for s=0.001s_{*}=0.001 and ξ2=(6.06.9)×103\xi_{2}=(6.0-6.9)\times 10^{-3} for s=0.1s_{*}=0.1. When λ=0.1\lambda=0.1, although nsn_{s} and rr get smaller by increasing ξ2\xi_{2}, and in particular, nsn_{s} can be brought to the observationally allowed range, the suppression of rr is not enough to be allowed by the Planck-BK constraints. Field trajectories for aforementioned parameter choices are shown in Fig. 13 with ϕ=0\phi_{*}=0 as an example.

Refer to caption
Figure 13: Field trajectories for the power-law inflation model. Two choices of ss_{*}, 0.001 (solid) and 0.1 (dashed), are considered with λ=0.01\lambda=0.01 (red), 0.050.05 (blue), and 11 (green). For each case, we have chosen the value of ξ2\xi_{2} in such a way that the spectral index becomes 0.965. For the demonstration of field trajectories, ϕ=0\phi_{*}=0 is chosen.
Refer to caption
Figure 14: Nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)} in terms of the spectral index nsn_{s} for the power-law inflation model. Two choices of ss_{*} are considered, 0.001 (left) and 0.1 (right). The shaded region represents the Planck 2-sigma bound on the spectral index, ns=[0.955,0.976]n_{s}=[0.955,0.976]. In the case of λ=0.1\lambda=0.1 with s=0.001s_{*}=0.001, the nonlinearity parameter may become slightly larger than the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. On the other hand, for a relatively large value of ss_{*} such as the s=0.1s_{*}=0.1 case, the nonlinearity parameter tends to be tiny.

Figure 14 shows the nonlinearity parameter fNL(local)f_{\rm NL}^{(\rm local)}, obtained by using Eq. (3.26), as a function of the spectral index nsn_{s} for the power-law inflation model. Similar to Figs. 12, two choices of ss_{*}, 0.001 and 0.1, are considered. The shaded grey region indicates the Planck 2-sigma bound on the spectral index. In the case of λ=0.1\lambda=0.1 with s=0.001s_{*}=0.001, the nonlinearity parameter may become slightly larger than the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. On the other hand, for a relatively large value of ss_{*} such as the s=0.1s_{*}=0.1 case, the nonlinearity parameter tends to be tiny.

Refer to caption
Refer to caption
Figure 15: Running of the spectral index αs\alpha_{s} in terms of the spectral index nsn_{s} for the power-law inflation model. Similar to Fig. 14, two choices of ss_{*} are considered with three different values for λ\lambda. We observe that αs\alpha_{s} initially increases and then decreases as nsn_{s} decreases, and that a larger value of ss_{*} results in smaller αs\alpha_{s}. See also Appendix A.

Finally, the running of the spectral index αs\alpha_{s}, discussed in Appendix A, is presented in Fig. 15. Again, two choices of ss_{*} and three values for the parameter λ\lambda are considered. We observe that, as the spectral index decreases, the running of the spectral index first shows an increasing behaviour, and then it decreases. We further see that the running of the spectral index tends to be smaller when a larger value is considered for ss_{*}. We note that this behaviour aligns with the general analysis presented in Appendix A.

4.3 Hybrid inflation – an example for Class II

The scalar potential for the hybrid inflation model [5] is given by

V(ϕ,ψ)=12m2ϕ2+λ4(ψ2Δ2)2+λ2ϕ2ψ2,\displaystyle V(\phi,\psi)=\frac{1}{2}m^{2}\phi^{2}+\frac{\lambda^{\prime}}{4}(\psi^{2}-\Delta^{2})^{2}+\frac{\lambda}{2}\phi^{2}\psi^{2}\,, (4.7)

where ψ\psi is the waterfall field, mm is a mass of ϕ\phi, λ\lambda and λ\lambda^{\prime} are dimensionless coupling constants, and Δ\Delta is a constant parameter that has the dimension of mass. The inflationary trajectory mainly goes through the valley of the potential along ψ=0\psi=0. Thus, during inflation, the hybrid inflation model can effectively be described by a single-field model with the inflaton potential given by

VJ(ϕ)=Λ4[1+(ϕμ)2],\displaystyle V_{\rm J}(\phi)=\Lambda^{4}\left[1+\left(\frac{\phi}{\mu}\right)^{2}\right]\,, (4.8)

where Λ=(λ/4)1/4Δ\Lambda=(\lambda^{\prime}/4)^{1/4}\Delta and μ=λ/2Δ2/m\mu=\sqrt{\lambda^{\prime}/2}\Delta^{2}/m. The slow-roll parameters (3.14) are given by

ϵ(0)=2(ϕ/μ)2μ2[1+(ϕ/μ)2]2,η(0)=2μ2[1+(ϕ/μ)2].\displaystyle\epsilon^{(0)}=\frac{2(\phi/\mu)^{2}}{\mu^{2}[1+(\phi/\mu)^{2}]^{2}}\,,\quad\eta^{(0)}=\frac{2}{\mu^{2}[1+(\phi/\mu)^{2}]}\,. (4.9)

We are interested in the parameter space where end of inflation is governed by the waterfall phase, not by slow-roll violations. We thus focus on the μ>1/2\mu>1/\sqrt{2} region so that the slow-roll parameters remain to be smaller than unity for any value of ϕ\phi. In this case, the model falls into Class II. We shall view ϕe\phi_{e}, the field value of the inflaton at the end of inflation, as a free parameter, together with μ\mu. Once ϕe\phi_{e} and μ\mu are fixed, we use 60 ee-folds, i.e., 60=e(VJ/VJ,ϕ)𝑑ϕ60=\int_{e}^{*}(V_{\rm J}/V_{{\rm J},\phi})d\phi, to find the field value at the pivot scale and thus the values of the spectral index and the tensor-to-scalar ratio. The original predictions are shown in Fig. 16 for ϕe=0.1μ\phi_{e}=0.1\mu and ϕe=0.01μ\phi_{e}=0.01\mu (magenta lines) by varying μ\mu. We see that the original model is ruled out by the latest observational constraints.

Refer to caption
Figure 16: Predictions for the hybrid inflation model. Two values of ϕe\phi_{e} are considered: ϕe=0.1μ\phi_{e}=0.1\mu (left panel) and ϕe=0.01μ\phi_{e}=0.01\mu (right panel). The original model, depicted by the magenta lines, is ruled out by the latest Planck (grey) and Planck-BK (black) bounds. By varying the nonminimal coupling parameter ξ2\xi_{2}, the effect of the assistant field is presented for μ={7.5,10,15,30}\mu=\{7.5,10,15,30\} when ϕe=0.1μ\phi_{e}=0.1\mu (left panel) and μ={2,6,8,14}\mu=\{2,6,8,14\} when ϕe=0.01μ\phi_{e}=0.01\mu (right panel) with two choices of ss_{*}, 0.001 (solid) and 0.1 (dashed). In the presence of the assistant field, the model may become compatible with the latest observational constraints.
Refer to caption
Refer to caption
Figure 17: Field trajectories for the hybrid inflation model. Two choices of ss_{*}, 0.001 (solid) and 0.1 (dashed), are considered. The left panel shows the ϕe=0.1μ\phi_{e}=0.1\mu case with μ=7.5\mu=7.5 (red), 1010 (magenta), 1515 (green), and 3030 (blue), whereas the right panel shows the ϕe=0.01μ\phi_{e}=0.01\mu case with μ=2\mu=2 (red), 66 (blue), 88 (green), and 1414 (magenta). For each case, we have chosen the value of ξ2\xi_{2} in such a way that the spectral index becomes 0.965, except for {s,ϕe,μ}={0.001,0.01μ,2}\{s_{*},\phi_{e},\mu\}=\{0.001,0.01\mu,2\} in which ns=0.965n_{s}=0.965 cannot be reached for any values of ξ2\xi_{2}, and hence, ξ2=0.0001\xi_{2}=0.0001 is chosen for illustration.

In Fig. 16, using the analytical expressions (3.22) and (3.25), the spectral index nsn_{s} and the tensor-to-scalar ratio rr in the presence of the assistant field are presented for μ={7.5,10,15,30}\mu=\{7.5,10,15,30\} with s={0.001,0.1}s_{*}=\{0.001,0.1\} for the case of ϕe=0.1μ\phi_{e}=0.1\mu and μ={2,6,8,14}\mu=\{2,6,8,14\} with s={0.001,0.1}s_{*}=\{0.001,0.1\} for the case of ϕe=0.01μ\phi_{e}=0.01\mu. As expected, the inclusion of the assistant field brings the original predictions to the observationally-favoured region by suppressing the spectral index as well as the tensor-to-scalar ratio. For example, in the case of ϕe=0.1μ\phi_{e}=0.1\mu and μ=10\mu=10, the prediction becomes compatible with the Planck-BK 2-sigma bounds for the range of ξ2=(1.01.1)×102\xi_{2}=(1.0-1.1)\times 10^{-2} for s=0.001s_{*}=0.001 and ξ2=(4.35.3)×103\xi_{2}=(4.3-5.3)\times 10^{-3} for s=0.1s_{*}=0.1. In the case of ϕe=0.01μ\phi_{e}=0.01\mu and μ=6\mu=6, the corresponding ranges for the nonminimal coupling parameter are ξ2=(5.35.9)×103\xi_{2}=(5.3-5.9)\times 10^{-3} and ξ2=(2.24.3)×103\xi_{2}=(2.2-4.3)\times 10^{-3} for s=0.001s_{\ast}=0.001 and s=0.1s_{\ast}=0.1, respectively. Field trajectories for all the parameter choices are shown in Fig. 17.

Refer to caption
Figure 18: Nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)} in terms of the spectral index nsn_{s} for the hybrid inflation model. Two choices of ss_{*} are considered, 0.001 (left) and 0.1 (right). In addition, ϕe=0.1μ\phi_{e}=0.1\mu (solid lines) and ϕe=0.01μ\phi_{e}=0.01\mu (dashed lines) are considered. The shaded region represents the Planck 2-sigma bound on the spectral index, ns=[0.955,0.976]n_{s}=[0.955,0.976]. In the case of s=0.001s_{*}=0.001, the nonlinearity parameter may become large and incompatible with the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. On the other hand, for a relatively large value of ss_{*} such as the s=0.1s_{*}=0.1 case, the nonlinearity parameter tends to be tiny.

Figure 18 shows the nonlinearity parameter fNL(local)f_{\rm NL}^{(\rm local)}, obtained by using Eq. (3.26), as a function of the spectral index nsn_{s} for the hybrid inflation model. Similar to Fig. 16, two choices of ss_{*}, 0.001 and 0.1, are considered with ϕe=0.1μ\phi_{e}=0.1\mu and ϕe=0.01μ\phi_{e}=0.01\mu. The shaded grey region indicates the Planck 2-sigma bound on the spectral index. In the case of s=0.001s_{*}=0.001, the nonlinearity parameter may become large. For example, in the cases of {ϕe,μ}={0.1μ,7.5}\{\phi_{e},\mu\}=\{0.1\mu,7.5\}, {0.1μ,10}\{0.1\mu,10\}, and {0.01μ,6}\{0.01\mu,6\}, the nonlinearity parameter is incompatible with the Planck 2-sigma bound, 11.1<fNL(local)<9.3-11.1<f_{\rm NL}^{({\rm local})}<9.3. On the other hand, the nonlinearity parameter tends to be tiny in the s=0.1s_{*}=0.1 case.

Refer to caption
Refer to caption
Figure 19: Running of the spectral index αs\alpha_{s} in terms of the spectral index nsn_{s} for the hybrid inflation model. Similar to Fig. 18, two choices of ss_{*} are considered with two different values for ϕe\phi_{e}. We observe that αs\alpha_{s} initially increases and then decreases as nsn_{s} decreases, and that a larger value of ss_{*} results in smaller αs\alpha_{s}. See also Appendix A.

The running of the spectral index αs\alpha_{s}, discussed in Appendix A, is presented in Fig. 19. Similar to Fig. 18, two choices of ss_{*} and two values for ϕe\phi_{e} are considered. The general tendency of the behaviour of αs\alpha_{s} remains the same as those for the previous two example models. In other words, the running of the spectral index shows an increasing behaviour before it decreases again, and it tends to take a smaller value when ss_{*} takes a larger value. For the detailed discussion of the running of the spectral index, see Appendix A.

5 Conclusion

We have investigated the impact of the assistant field on single-field inflationary models. The assistant field is characterised by being nonminimally coupled to gravity, effectively massless, and having no direct coupling to the original inflaton field. Without specifying a potential for the inflaton field, we performed a general analysis and presented analytical expressions for inflationary observables such as the spectral index nsn_{s}, the tensor-to-scalar ratio rr, and the local-type nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)}, for both the metric and Palatini formulations.

For Class I, where the end of inflation is achieved through slow-roll violations, the assistant field reduces both nsn_{s} and rr, potentially reviving many single-field models that were previously ruled out by observations. For Class II, where the end of inflation is determined by a separate sector, a small ϵe(0)\epsilon_{e}^{(0)} value may bring a small nsn_{s} into the observationally-favoured region. The compatibility of the nonlinearity parameter fNL(local)f_{\rm NL}^{\rm(local)} with observations is also shown for both Class I and Class II.

Our results are demonstrated using three example models: loop inflation (Class I), power-law inflation, and hybrid inflation (Class II). These models were previously ruled out due to large values of nsn_{s} and/or rr. Our findings show that the presence of the assistant field can bring the predictions of nsn_{s} and rr into the observationally-favoured region, making the models compatible with the current observational bounds.

It is worth noting that not all single-field models can be revived with the help of the assistant field. For instance, models that originally predict a small nsn_{s} and belong to either Class I or Class II with a large ϵe(0)\epsilon_{e}^{(0)} cannot be revived. The impact of a massive assistant field on such models remains an open question and will be the subject of future study.

Acknowledgments

This work was supported by National Research Foundation grants funded by the Korean government (NRF-2021R1A4A2001897) and (NRF-2019R1A2C1089334) (S.C.P.), and JSPS KAKENHI Grants No. 19K03874 (T.T.).

Appendix A Running of the Spectral Index

The running of the scalar spectral index in a general multi-field inflation model is given by [70]

αs\displaystyle\alpha_{s} dnsdlnk4ϵ22ϵ˙H+2N,iN,jGmnN,mN,n(4ϵwij+2wiwjkkwijφ˙0kH;k)(ns1)2,\displaystyle\equiv\frac{dn_{s}}{d\ln k}\simeq 4\epsilon^{2}-2\frac{\dot{\epsilon}}{H}+2\frac{N_{,i}N_{,j}}{G^{mn}N_{,m}N_{,n}}\left(4\epsilon w^{ij}+2w^{i}{}_{k}w^{jk}-w^{ij}{}_{;k}\frac{\dot{\varphi}^{k}_{0}}{H}\right)-(n_{s}-1)^{2}\,, (A.1)

evaluated at the horizon crossing, which we denote by the subscript *, where

wij=u(i;j)+13Rm(ij)nφ˙0mφ˙0nH2,ui=V,i3H2,\displaystyle w_{ij}=u_{(i;j)}+\frac{1}{3}R_{m(ij)n}\frac{\dot{\varphi}^{m}_{0}\dot{\varphi}^{n}_{0}}{H^{2}}\,,\qquad u_{i}=-\frac{V_{,i}}{3H^{2}}\,, (A.2)

with φ0i\varphi^{i}_{0} being the background fields. For the two-field model under consideration, one finds that

ϵ˙H|\displaystyle\frac{\dot{\epsilon}}{H}\bigg{|}_{*} 2ϵϕ(ηϕϕ2ϵϕ)2ϵσ(ησσ2ϵσ)+ϵϕϵσϵbsbsσ,\displaystyle\approx-2\epsilon^{\phi}_{*}\left(\eta^{\phi\phi}_{*}-2\epsilon^{\phi}_{*}\right)-2\epsilon^{\sigma}_{*}\left(\eta^{\sigma\sigma}_{*}-2\epsilon^{\sigma}_{*}\right)+\epsilon^{\phi}_{*}\sqrt{\epsilon^{\sigma}_{*}\epsilon^{b}_{*}}s^{b}_{*}s^{\sigma}_{*}\,,

and

2N,iN,jGmnN,mN,n(4ϵwij+2wiwjkkwijφ˙kH;k)\displaystyle 2\frac{N_{,i}N_{,j}}{G^{mn}N_{,m}N_{,n}}\left(4\epsilon w^{ij}+2w^{i}{}_{k}w^{jk}-w^{ij}{}_{;k}\frac{\dot{\varphi}^{k}}{H}\right) 2(N,ϕ)2Aϕϕ+2(N,σ)2Aσσ+4N,ϕN,σAϕσe2b(N,ϕ)2+(N,σ)2,\displaystyle\approx\frac{2(N_{,\phi})^{2}A^{\phi\phi}+2(N_{,\sigma})^{2}A^{\sigma\sigma}+4N_{,\phi}N_{,\sigma}A^{\phi\sigma}}{e^{-2b}(N_{,\phi})^{2}+(N_{,\sigma})^{2}}\,,

where

Aϕϕ\displaystyle A^{\phi\phi} =4ϵwϕϕ+2e2b(wϕϕ)2+2(wϕσ)2+ebsϕ2ϵϕwϕϕ+;ϕsσ2ϵσwϕϕ,;σ\displaystyle=4\epsilon w^{\phi\phi}+2e^{2b}(w^{\phi\phi})^{2}+2(w^{\phi\sigma})^{2}+e^{-b}s^{\phi}\sqrt{2\epsilon^{\phi}}w^{\phi\phi}{}_{;\phi}+s^{\sigma}\sqrt{2\epsilon^{\sigma}}w^{\phi\phi}{}_{;\sigma}\,,
Aσσ\displaystyle A^{\sigma\sigma} =4ϵwσσ+2e2b(wϕσ)2+2(wσσ)2+ebsϕ2ϵϕwσσ+;ϕsσ2ϵσwσσ,;σ\displaystyle=4\epsilon w^{\sigma\sigma}+2e^{2b}(w^{\phi\sigma})^{2}+2(w^{\sigma\sigma})^{2}+e^{-b}s^{\phi}\sqrt{2\epsilon^{\phi}}w^{\sigma\sigma}{}_{;\phi}+s^{\sigma}\sqrt{2\epsilon^{\sigma}}w^{\sigma\sigma}{}_{;\sigma}\,,
Aϕσ\displaystyle A^{\phi\sigma} =4ϵwϕσ+2e2bwϕϕwϕσ+2wϕσwσσ+ebsϕ2ϵϕwϕσ+;ϕsσ2ϵσwϕσ,;σ\displaystyle=4\epsilon w^{\phi\sigma}+2e^{2b}w^{\phi\phi}w^{\phi\sigma}+2w^{\phi\sigma}w^{\sigma\sigma}+e^{-b}s^{\phi}\sqrt{2\epsilon^{\phi}}w^{\phi\sigma}{}_{;\phi}+s^{\sigma}\sqrt{2\epsilon^{\sigma}}w^{\phi\sigma}{}_{;\sigma}\,,

and

N,ϕ|\displaystyle N_{,\phi}\Big{|}_{*} =12sϕϵϕϵeϕϵee2beb,\displaystyle=\frac{1}{\sqrt{2}}\frac{s^{\phi}_{*}}{\sqrt{\epsilon^{\phi}_{*}}}\frac{\epsilon^{\phi}_{e}}{\epsilon_{e}}e^{2b_{e}-b_{*}}\,,
N,σ|\displaystyle N_{,\sigma}\Big{|}_{*} =12sσϵσ(1ϵeϕϵee2be2b).\displaystyle=\frac{1}{\sqrt{2}}\frac{s^{\sigma}_{*}}{\sqrt{\epsilon^{\sigma}_{*}}}\left(1-\frac{\epsilon^{\phi}_{e}}{\epsilon_{e}}e^{2b_{e}-2b_{*}}\right)\,.

From the definition of wijw_{ij} in Eq. (A.2), we find that

wϕϕ|\displaystyle w^{\phi\phi}\Big{|}_{*} =e2b[2ϵϕηϕϕ12sbsσϵbϵσ+124ϵσ(ηb+2ϵb)],\displaystyle=e^{-2b_{*}}\left[2\epsilon^{\phi}_{*}-\eta^{\phi\phi}_{*}-\frac{1}{2}s^{b}_{*}s^{\sigma}_{*}\sqrt{\epsilon^{b}_{*}\epsilon^{\sigma}_{*}}+\frac{1}{24}\epsilon^{\sigma}_{*}\left(\eta^{b}_{*}+2\epsilon^{b}_{*}\right)\right]\,,
wϕσ|\displaystyle w^{\phi\sigma}\Big{|}_{*} =12ebsϕϵϕ[sbϵb112sσϵσ(ηb+2ϵb)],\displaystyle=\frac{1}{2}e^{-b_{*}}s^{\phi}_{*}\sqrt{\epsilon^{\phi}_{*}}\left[s^{b}_{*}\sqrt{\epsilon^{b}_{*}}-\frac{1}{12}s^{\sigma}_{*}\sqrt{\epsilon^{\sigma}_{*}}\left(\eta^{b}_{*}+2\epsilon^{b}_{*}\right)\right]\,,
wσσ|\displaystyle w^{\sigma\sigma}\Big{|}_{*} =2ϵσησσ+124ϵϕ(ηb+2ϵb),\displaystyle=2\epsilon^{\sigma}_{*}-\eta^{\sigma\sigma}_{*}+\frac{1}{24}\epsilon^{\phi}_{*}\left(\eta^{b}_{*}+2\epsilon^{b}_{*}\right)\,,

and

wϕϕ|;ϕ\displaystyle w^{\phi\phi}{}_{;\phi}\Big{|}_{*} =12ϵϕsϕeb[ξϕϕϕ+6ηϕϕϵϕ8(ϵϕ)2+12ϵbϵϕ124sbsσϵbϵσϵϕ(ηb+2ϵb)],\displaystyle=\frac{1}{\sqrt{2\epsilon^{\phi}_{*}}}s^{\phi}_{*}e^{-b_{*}}\left[-\xi^{\phi\phi\phi}_{*}+6\eta^{\phi\phi}_{*}\epsilon^{\phi}_{*}-8(\epsilon^{\phi}_{*})^{2}+\frac{1}{2}\epsilon^{b}_{*}\epsilon^{\phi}_{*}-\frac{1}{24}s^{b}_{*}s^{\sigma}_{*}\sqrt{\epsilon^{b}_{*}\epsilon^{\sigma}_{*}}\epsilon^{\phi}_{*}\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\right]\,,
wϕϕ|;σ\displaystyle w^{\phi\phi}{}_{;\sigma}\Big{|}_{*} =12ϵbe2b{13sbϵσξb+sbϵb(ηϕϕ2ϵϕ)+12sbϵb[112ϵσηb(ηϕϕ2ϵσ)]\displaystyle=\frac{1}{\sqrt{2\epsilon^{b}_{*}}}e^{-2b_{*}}\bigg{\{}\frac{1}{3}s^{b}_{*}\epsilon^{\sigma}_{*}\xi^{b}_{*}+s^{b}_{*}\epsilon^{b}_{*}\Big{(}\eta^{\phi\phi}_{*}-2\epsilon^{\phi}_{*}\Big{)}+\frac{1}{2}s^{b}_{*}\epsilon^{b}_{*}\left[\frac{1}{12}\epsilon^{\sigma}_{*}\eta^{b}_{*}-\Big{(}\eta^{\phi\phi}_{*}-2\epsilon^{\sigma}_{*}\Big{)}\right]
+124sσϵbϵσ[2(ηb+2ϵb)(ησσ2ϵσ)3ηb]},\displaystyle\qquad\qquad\qquad+\frac{1}{24}s^{\sigma}_{*}\sqrt{\epsilon^{b}_{*}\epsilon^{\sigma}_{*}}\left[2\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\Big{(}\eta^{\sigma\sigma}_{*}-2\epsilon^{\sigma}_{*}\Big{)}-3\eta^{b}_{*}\right]\bigg{\}}\,,
wσσ|;ϕ\displaystyle w^{\sigma\sigma}{}_{;\phi}\Big{|}_{*} =1242sϕϵϕeb[(ηb+2ϵb)(2ηϕϕ4ϵϕ+sbsσϵbϵσ)12ϵb],\displaystyle=\frac{1}{24\sqrt{2}}s^{\phi}_{*}\sqrt{\epsilon^{\phi}_{*}}e^{b_{*}}\left[\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\Big{(}2\eta^{\phi\phi}_{*}-4\epsilon^{\phi}_{*}+s^{b}_{*}s^{\sigma}_{*}\sqrt{\epsilon^{b}_{*}\epsilon^{\sigma}_{*}}\Big{)}-12\epsilon^{b}_{*}\right]\,,
wσσ|;σ\displaystyle w^{\sigma\sigma}{}_{;\sigma}\Big{|}_{*} =12ϵσsσ{ξσσσ+2ϵσ(3ησσ4ϵσ)112sbsσϵσϵbϵϕ[(ϵb)24ξb]},\displaystyle=\frac{1}{\sqrt{2\epsilon^{\sigma}_{*}}}s^{\sigma}_{*}\left\{-\xi^{\sigma\sigma\sigma}_{*}+2\epsilon^{\sigma}_{*}\Big{(}3\eta^{\sigma\sigma}_{*}-4\epsilon^{\sigma}_{*}\Big{)}-\frac{1}{12}s^{b}_{*}s^{\sigma}_{*}\sqrt{\frac{\epsilon^{\sigma}_{*}}{\epsilon^{b}_{*}}}\epsilon^{\phi}_{*}\left[\Big{(}\epsilon^{b}_{*}\Big{)}^{2}-4\xi^{b}_{*}\right]\right\}\,,
wϕσ|;ϕ\displaystyle w^{\phi\sigma}{}_{;\phi}\Big{|}_{*} =1482{sbϵb[48(ηϕϕ2ϵϕ)24(ησσ2ϵσ)+(ϵϕϵσ)(ηb+2ϵb)]\displaystyle=\frac{1}{48\sqrt{2}}\bigg{\{}s^{b}_{*}\sqrt{\epsilon^{b}_{*}}\left[48\Big{(}\eta^{\phi\phi}_{*}-2\epsilon^{\phi}_{*}\Big{)}-24\Big{(}\eta^{\sigma\sigma}_{*}-2\epsilon^{\sigma}_{*}\Big{)}+\Big{(}\epsilon^{\phi}_{*}-\epsilon^{\sigma}_{*}\Big{)}\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\right]
+2sσϵσ[6ϵσ(ηϕϕ2ϵϕ)(ηb+2ϵb)]},\displaystyle\qquad\qquad\qquad+2s^{\sigma}_{*}\sqrt{\epsilon^{\sigma}_{*}}\left[6\sqrt{\epsilon^{\sigma}_{*}}-\Big{(}\eta^{\phi\phi}_{*}-2\epsilon^{\phi}_{*}\Big{)}\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\right]\bigg{\}}\,,
wϕσ|;σ\displaystyle w^{\phi\sigma}{}_{;\sigma}\Big{|}_{*} =1242sϕϵϕeb{3(ηb2ϵb)(ηb+2ϵb)(ησσ2ϵσ)\displaystyle=\frac{1}{24\sqrt{2}}s^{\phi}_{*}\sqrt{\epsilon^{\phi}_{*}}e^{-b_{*}}\bigg{\{}3\Big{(}\eta^{b}_{*}-2\epsilon^{b}_{*}\Big{)}-\Big{(}\eta^{b}_{*}+2\epsilon^{b}_{*}\Big{)}\Big{(}\eta^{\sigma\sigma}_{*}-2\epsilon^{\sigma}_{*}\Big{)}
+12sbsσϵσϵb[2(ϵb)2ϵbηb16ξb]}.\displaystyle\qquad\qquad\qquad\qquad\qquad+\frac{1}{2}s^{b}_{*}s^{\sigma}_{*}\sqrt{\frac{\epsilon^{\sigma}_{*}}{\epsilon^{b}_{*}}}\left[2\Big{(}\epsilon^{b}_{*}\Big{)}^{2}-\epsilon^{b}_{*}\eta^{b}_{*}-16\xi^{b}_{*}\right]\bigg{\}}\,.

Here, we have defined

ξϕϕϕV,ϕϕϕV,ϕV2e4b,ξσσσV,σσσV,σV2,ξb8b,σσσb,σ.\displaystyle\xi^{\phi\phi\phi}\equiv\frac{V_{,\phi\phi\phi}V_{,\phi}}{V^{2}}e^{-4b}\,,\quad\xi^{\sigma\sigma\sigma}\equiv\frac{V_{,\sigma\sigma\sigma}V_{,\sigma}}{V^{2}}\,,\quad\xi^{b}\equiv 8b_{,\sigma\sigma\sigma}b_{,\sigma}\,.

Combining all the individual terms, together with Eq. (3.1), into Eq. (A.1) gives the slow-roll-approximated analytical expression for the running of the scalar spectral index. As the calculation is straightforward and the resultant expression is rather long, we do not present the final expression here.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 20: In the upper panel, the evolution of the running of the spectral index is shown for Class I in terms of the nonminimal coupling parameter ξ2\xi_{2}. Two values of ss_{*} are considered, 0.001 (left) and 0.1 (right) with various choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\}. Furthermore, we have fixed αs(0)=104\alpha_{s}^{(0)}=10^{-4}. In the lower panel, we present the predictions in the αs(I)\alpha_{s}^{({\rm I})}ns(I)n_{s}^{({\rm I})} plane. The shaded region corresponds to the latest bounds on the spectral index. We observe that the running first grows and then eventually decreases as ξ2\xi_{2} increases. Only the metric formulation is considered.
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 21: In the upper panel, the evolution of the running of the spectral index is shown for Class II in terms of the nonminimal coupling parameter ξ2\xi_{2}. Two values of ss_{*} are considered, 0.001 (left) and 0.1 (right) with various choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\}. Furthermore, two different values, 10210^{-2} and 10510^{-5}, are chosen for ϵe(0)\epsilon_{e}^{(0)}. For all cases, we have fixed αs(0)=104\alpha_{s}^{(0)}=10^{-4}. In the lower panel, we present the predictions in the αs(II)\alpha_{s}^{({\rm II})}ns(II)n_{s}^{({\rm II})} plane. The shaded region corresponds to the latest bounds on the spectral index. We again observe that αs\alpha_{s} grows and then decreases as ξ2\xi_{2} increases.

Figure 20 shows the behaviour of the running of the spectral index as a function of the nonminimal coupling parameter (upper panel) and the spectral index (lower panel) for Class I. Two values of ss_{*}, 0.001 (left) and 0.1 (right), are considered, together with various choices of {ns(0),r(0)}\{n_{s}^{(0)},r^{(0)}\}. For all cases, we have fixed αs(0)\alpha_{s}^{(0)} to be 10410^{-4}. We observe that the running of the spectral index first grows as the nonminimal coupling parameter increases. Eventually, the running of the spectral index decreases. This behaviour is similar to that of the nonlinearity parameter. For Class II, the behaviour of the running of the spectral index remains the same as presented in Fig. 21. On top of the various choices for ss_{*}, ns(0)n_{s}^{(0)}, and r(0)r^{(0)}, we have considered two values, 10210^{-2} and 10510^{-5}, for ϵe(0)\epsilon_{e}^{(0)}.

The latest Planck experiment constrains the running of the spectral index as, e.g., 0.0158αs0.0012-0.0158\leq\alpha_{s}\leq-0.0012 (Planck TT,TE,EE+lowEB+lensing) at 68% C.L. On the other hand, once the running of the running of the spectral index is taken into account, the constraint becomes

0.001αs0.025,(Planck TT+lowE+lensing)\displaystyle 0.001\leq\alpha_{s}\leq 0.025\,,\quad\text{(Planck TT+lowE+lensing)} (A.3)

or 0.008αs0.012-0.008\leq\alpha_{s}\leq 0.012 (Planck TT,TE,EE++lowE++lensing), both at 68% C.L. [1]. For all cases we have considered, αs0.015\alpha_{s}\lesssim 0.015. We also note that the running of the spectral index tends to be smaller when ss_{*} takes a larger value, giving αs0.007\alpha_{s}\lesssim 0.007, which is well within the allowed bounds obtained with the running of the running, while they are larger than the bounds where the running of the running is not taken into account.

Actually, the current constraint on the running is rather weak. It can be well improved in future observations such as the 21 cm line of neutral hydrogen [71, 72, 73] and galaxy surveys with CMB [74, 75, 76], which may be useful to differentiate some cases that give the same predictions on nsn_{s} and rr.

References