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

Double inflation via non-minimally coupled spectator

Mio Kubota    Kin-ya Oda    Stanislav Rusak    and Tomo Takahashi
Abstract

We argue that double inflation may occur when a spectator field is non-minimally coupled to gravity. As a concrete example, we study a two-field inflationary model where the initial spectator field is non-minimally coupled to gravity while the initial inflaton field is minimally coupled. The non-minimal coupling results in the growth of the spectator field which, in turn, drives the second stage of inflation in a significant region of parameter space. The isocurvature fluctuations originating from the spectator field source adiabatic ones, and hence the spectator non-minimal coupling can modify the inflationary predictions for the spectral index and the tensor-to-scalar ratio even though the initial inflaton field is minimally coupled to gravity. We explicitly show that quadratic chaotic inflation can become viable by the introduction of the spectator non-minimal coupling.

OCHA-PP-370

1 Introduction

The Universe is believed to have experienced a period of nearly exponential expansion, dubbed inflation, in its very early stage. This inflationary period not only helps address the flatness and horizon problems but also generically provides a mechanism to generate primordial density fluctuations necessary for structure formation. A simple possibility is that a single scalar field, called inflaton, is responsible for both driving the inflationary expansion and also generating the primordial fluctuations. Since the nature of the primordial fluctuations depends on models of the inflaton, cosmological observations such as the cosmic microwave background (CMB) can probe these models, and indeed, data from the Planck satellite in combination with other observations, especially CMB B-mode polarization experiments such as BICEP/Keck, provide constraints on these models [1, 2].

Though having only a single field responsible for inflation is attractive, the simplest inflationary potentials are now severely constrained by cosmological data and many have been ruled out. Therefore multi-field models have also been extensively discussed in the literature. These include cases where the inflaton is only responsible for the inflationary expansion while another scalar field, whose energy density is negligible during inflation, and irrelevant to the inflationary dynamics (referred to as a spectator field), is responsible for generating primordial fluctuations. Examples of such models include the curvaton model [3, 4, 5] and modulated reheating scenarios [6, 7]. In high energy theories such as string theory, the existence of multiple scalar fields is generic, which also motivates a multi-field inflation model. In addition, when a spectator field is responsible for primordial fluctuations, the predictions for primordial power spectra are modified and, in some cases, even if the model is excluded by the data as a single-field inflation model, it can be resurrected in the framework of spectator field models; see, e.g., [8, 9, 10, 11, 12, 13, 14] for the curvaton model and [15] for the modulated reheating scenario.

Another extension of a single-field inflation model is to introduce a non-minimal coupling to gravity. A well-known example is the Higgs inflation [16, 17] where the inflaton is identified with the Standard Model Higgs field with non-minimal coupling to gravity.111 There have been early works along this line. See [18, 19, 20, 21] for such early studies. The potential for the Higgs field is quartic during inflation in the pure Standard Model, and therefore the pure Standard Model Higgs inflation is excluded by observational data. However, by introducing a non-minimal coupling, the model becomes consistent with cosmological observations in terms of the spectral index and the tensor-to-scalar ratio (the suppression of the tensor-to-scalar ratio has already been discussed in [21]). Inflationary models with non-minimal coupling have become an extensively studied research area as they generally predict modified inflationary observables; see, e.g., [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34]. In addition, non-minimal couplings generically arise radiatively in quantum field theory in curved space-time [35, 36] even if they do not exist at tree level. Therefore inflation models with non-minimal coupling to gravity are well motivated theoretically.

Given that both multi-field inflationary models and non-minimal couplings are well motivated, it is tempting to consider multi-field models with non-minimal couplings; see, e.g., [37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47] for works along this line. In this paper, we consider a multi-field model in which there exist two scalar fields: a minimally coupled field Φ\Phi and non-minimally coupled one SS, with Φ\Phi being the inflaton and SS a spectator initially. As will be shown in this paper, due to the spectator non-minimal coupling, its amplitude can grow during the inflationary phase driven by the Φ\Phi field until eventually the SS field acquires super-Planckian field values, which, in turn, drives a second inflationary phase after the end of the first inflationary epoch driven by the Φ\Phi field.222 Even in minimally coupled models, when the spectator field is assumed to have a super-Planckian amplitude, the second inflationary phase can be driven by the spectator field (these kinds of models are sometimes called “inflating curvaton models”) [8, 9, 11, 48, 49].

We also discuss primordial fluctuations in such models. Since isocurvature fluctuations generated by the spectator field can source the adiabatic ones, the existence of the spectator non-minimal coupling can also modify primordial power spectra from its minimally coupled counterpart. In particular, the predictions for the spectral index nsn_{s} and the tensor-to-scalar ratio rr are modified. This issue has already been discussed for single-field non-minimally coupled models (see, e.g., [22, 23, 24, 25, 27, 28, 30, 31, 32, 33, 34]). We show in this paper that the spectator non-minimal coupling can also relax the constraints on nsn_{s} and rr so that even inflation models that are excluded in the single-field minimally coupled framework, such as the quadratic chaotic inflation, can become viable due to the existence of the spectator non-minimal coupling.

The structure of this paper is the following. In the next section, we describe the setup for the multi-field inflation model with the non-minimal coupling we consider in this paper. Then in Section 3, we discuss the background evolution, particularly focusing on the growth of the spectator field during the initial inflationary phase driven by Φ\Phi. We also argue that in our setup, a second inflationary phase is somewhat generic due to this growth originating from the existence of the spectator non-minimal coupling. In Section 4, we discuss primordial scalar and tensor power spectra in the model and how they are modified by changing the strength of the non-minimal coupling and the masses of the fields. We also investigate how the inflationary predictions are modified to relax the constraints on inflation models. The final section is devoted to conclusions and discussion.

2 Setup

We consider two scalar fields Φ\Phi and SS that play the role of inflaton and spectator, respectively, at the early stages of inflation. For simplicity, we assume that Φ\Phi is minimally coupled to gravity, while SS non-minimally at the leading order such that the action in the Jordan frame is given by

SJ=d4xg(12Mpl2f(S)R12μΦμΦ12μSμSVJ(Φ,S)),S_{\rm J}=\int\mathrm{d}^{4}x\sqrt{-g}\left(\frac{1}{2}M_{\rm pl}^{2}f(S)R-\frac{1}{2}\partial_{\mu}\Phi\,\partial^{\mu}\Phi-\frac{1}{2}\partial_{\mu}S\,\partial^{\mu}S-V_{\text{J}}(\Phi,S)\right), (2.1)

where Mpl=1/8πG=2.4×1018GeVM_{\rm pl}=1/\sqrt{8\pi G}=2.4\times 10^{18}\,\text{GeV} is the reduced Planck mass and f(S)f(S) characterizes a non-minimal coupling between SS and RR:

f(S)=1+ξS(SMpl)2,f(S)=1+\xi_{S}\left(\frac{S}{M_{\rm pl}}\right)^{2}, (2.2)

with ξS\xi_{S} being the non-minimal coupling constant, and VJV_{\text{J}} is the Jordan frame potential. As a simplest non-trivial example, we assume that the potential takes the following form:

VJ(Φ,S)=12mΦ2Φ2+12mS2S2.V_{\text{J}}(\Phi,S)=\frac{1}{2}m_{\Phi}^{2}\Phi^{2}+\frac{1}{2}m_{S}^{2}S^{2}. (2.3)

To follow the evolution of the scalar fields and the cosmic expansion, we work in the Einstein frame. For a system with multiple scalar fields, one can write the action in the Einstein frame as333 For the formalism of a general multi-field inflation, see e.g., [50, 37, 39]. We follow the notation of [50].

S=d4xg[12Mpl2R12GIJμφIμφJV(φI)],S=\int\mathrm{d}^{4}x\sqrt{-g}\left[\frac{1}{2}M_{\mathrm{pl}}^{2}R-\frac{1}{2}G_{IJ}\nabla_{\mu}\varphi^{I}\nabla^{\mu}\varphi^{J}-V(\varphi^{I})\right], (2.4)

where the indices run for I=Φ,SI=\Phi,S with the abbreviation (φΦ,φS)=(Φ,S)\left(\varphi^{\Phi},\varphi^{S}\right)=\left(\Phi,S\right) for the field space. Here, GIJG_{IJ} is the field space metric given by

GΦΦ=f1,GSS=f1+32fSfSf2,GΦS=GSΦ=0,G_{\Phi\Phi}=f^{-1},\qquad G_{SS}=f^{-1}+\frac{3}{2}\frac{f_{S}f_{S}}{f^{2}},\qquad G_{\Phi S}=G_{S\Phi}=0, (2.5)

with fS=f/S=2ξSS/Mpl2f_{S}=\partial f/\partial S=2\xi_{S}S/M_{\mathrm{pl}}^{2}, and V(φI)V(\varphi^{I}) is the Einstein-frame potential related to the Jordan-frame counterpart as

V(φI)=VJ(φI)f2.V(\varphi^{I})=\frac{V_{\text{J}}(\varphi^{I})}{f^{2}}. (2.6)

3 Background Evolution

Let us first discuss the background evolution of the fields. We first focus on the case where the non-minimal coupling is small and therefore the field space metric component of the spectator can be approximated as GSSf1G_{SS}\simeq f^{-1}. We can then introduce a canonically normalized spectator field444 The situation for the spectator is analogous to Higgs inflation in Palatini gravity [51, 52, 53, 54].

χMplξSarcsinh(ξSSMpl)\chi\equiv\frac{M_{\rm pl}}{\sqrt{\xi_{S}}}\,\,{\rm arcsinh}\left(\sqrt{\xi_{S}}\frac{S}{M_{\rm pl}}\right) (3.1)

to aid the intuitive understanding of the dynamics. Though this rescaling introduces some kinetic couplings between the spectator and inflaton fields, in the double inflation regime which we are interested in, such couplings will be unimportant. This is due firstly because of slow-roll dynamics and secondly due to relative smallness of the term ξSS2\xi_{S}S^{2} necessary to achieve desirable duration of inflation as we shall discuss.

The equation of motion for the inflaton is then given by

Φ¨+[3H2ξSMpltanh(ξSMpl1χ)χ˙]Φ˙+mΦ2Φcosh2(ξSMpl1χ)=0.\ddot{\Phi}+\left[3H-\frac{2\sqrt{\xi_{S}}}{M_{\mathrm{pl}}}\tanh\left(\sqrt{\xi_{S}}M_{\mathrm{pl}}^{-1}\chi\right)\dot{\chi}\right]\dot{\Phi}+\frac{m_{\Phi}^{2}\Phi}{\cosh^{2}\left(\sqrt{\xi_{S}}M_{\mathrm{pl}}^{-1}\chi\right)}=0. (3.2)

The presence of the spectator acts to scale down the effective mass of the inflaton and to introduce an extra anti-friction term. For slow-roll dynamics and moderate values of the spectator amplitude, however, these changes are negligible and the dynamics is virtually that of the usual quadratic chaotic inflation. The equation of motion for the canonical spectator, on the other hand, is

χ¨+3Hχ˙2ξSMplsinh(ξSMpl1χ)cosh5(ξSMpl1χ)[mΦ2Φ2mS2Mpl22ξS(1sinh2(ξSMpl1χ))]=0,\ddot{\chi}+3H\dot{\chi}-\frac{2\sqrt{\xi_{S}}}{M_{\mathrm{pl}}}\frac{\sinh(\sqrt{\xi_{S}}M_{\mathrm{pl}}^{-1}\chi)}{\cosh^{5}(\sqrt{\xi_{S}}M_{\mathrm{pl}}^{-1}\chi)}\left[m_{\Phi}^{2}\Phi^{2}-\frac{m_{S}^{2}M_{\mathrm{pl}}^{2}}{2\xi_{S}}\left(1-\sinh^{2}(\sqrt{\xi_{S}}M_{\mathrm{pl}}^{-1}\chi)\right)\right]=0, (3.3)

where we have neglected the term proportional to Φ˙2\dot{\Phi}^{2} as it is negligible during slow-roll.

We can then deduce the dynamics as follows. Initially, the inflaton field dominates and induces a force that amplifies the spectator field. Then, as the amplitude of Φ\Phi decreases past the threshold

Φturn=MplmS2ξSmΦ,\Phi_{\mathrm{turn}}=\frac{M_{\mathrm{pl}}m_{S}}{\sqrt{2\xi_{S}}m_{\Phi}}, (3.4)

the force reverses and the spectator is driven back toward zero as long as its amplitude has not exceeded the critical value

χrunaway=MplξSarcsinh(1)orSrunaway=MplξS.\chi_{\mathrm{runaway}}=\frac{M_{\mathrm{pl}}}{\sqrt{\xi_{S}}}\operatorname{arcsinh}(1)\qquad\mathrm{or}\qquad S_{\mathrm{runaway}}=\frac{M_{\mathrm{pl}}}{\sqrt{\xi_{S}}}. (3.5)

If this value is exceeded, the spectator is amplified indefinitely.

Refer to caption
Figure 1: The effective potential for the initial inflaton Φ\Phi and the initial rescaled spectator χ\chi. In the χ\chi direction, the potential switches from being convex at large Φ\Phi to being concave around the origin of χ\chi at small Φ\Phi, causing χ\chi to first be amplified and then roll back (blue trajectory). If χ\chi amplitude grows past the ridge of the potential, however, it continues to be amplified indefinitely resulting in runaway behaviour (red trajectory).

When the spectator is amplified to a large value, it comes to dominate the energy density, driving the second stage of inflation while it rolls back to zero. Essentially the physics in the regime of interest is described by the effective potential

Veff=ξSmΦ2Φ2+mS2Mpl2sinh2(ξSχ/Mpl)2ξScosh4(ξSχ/Mpl),V_{\rm eff}=\frac{\xi_{S}m_{\Phi}^{2}\Phi^{2}+m_{S}^{2}M_{\rm pl}^{2}~{}{\rm sinh}^{2}(\sqrt{\xi_{S}}\chi/M_{\rm pl})}{2\xi_{S}~{}{\rm cosh}^{4}(\sqrt{\xi_{S}}\chi/M_{\rm pl})}, (3.6)

which is depicted in Fig. 1: When Φ\Phi is large (upper left direction), VeffV_{\text{eff}} is convex in the direction of χ\chi, causing it to be amplified. At smaller values of Φ\Phi (lower right direction), VeffV_{\text{eff}} becomes concave in the χ\chi direction, causing it to roll back, and drives the second inflationary stage.

3.1 Numerical analysis

Having gained an analytical understanding of the dynamics we proceed to numerically solve the full equations of motion

𝒟tφ˙I+3Hφ˙+GIJV,J=0\mathcal{D}_{t}\dot{\varphi}^{I}+3H\dot{\varphi}+G^{IJ}V_{,J}=0 (3.7)

which is derived from the action (2.4). Here 𝒟t\mathcal{D}_{t} is the field-space-covariant time derivative: 𝒟tϕ˙I=ϕ¨I+ΓIϕ˙JJKϕ˙K\mathcal{D}_{t}\dot{\phi}^{I}\ =\ddot{\phi}^{I}+\Gamma^{I}{}_{JK}\dot{\phi}^{J}\dot{\phi}^{K}.

Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Refer to caption
Figure 2: Evolution of the initial spectator SS for the cases ξS=0.005\xi_{S}=0.005 (left), 0 (middle), and 0.005-0.005 (right) as a function of the number of ee-folds NN, which is taken to be 0 at the start of our numerical calculation. The values of S0S_{0} in the legend are given in units of MplM_{\rm pl}.

In the lower panels of Fig. 2, we show the evolution of the initial spectator SS as a function of the number of ee-folds NN for ξs=0.005\xi_{s}=0.005 (left). For reference, we also plot the cases ξS=0\xi_{S}=0 (middle) and ξS=0.005\xi_{S}=-0.005 (right). The initial value of Φ\Phi is set to be Φ0=20Mpl\Phi_{0}=20M_{\rm pl} and in this setup, the energy density of Φ\Phi dominates over that of SS at the beginning of the simulation. The initial inflaton amplitude is not important for the second inflationary stage, as long as it ensures initial domination of the Φ\Phi field because at first the dynamics is effectively single-field and quickly approaches an attractor. The masses of Φ\Phi and SS are taken to be mΦ=5×106Mplm_{\Phi}=5\times 10^{-6}M_{\rm pl} and mS=5×107Mplm_{S}=5\times 10^{-7}M_{\rm pl}. The ratio of the masses determines the dynamics, and changing the masses while keeping the ratio fixed merely changes the time scale. In the minimally coupled case ξS=0\xi_{S}=0 (middle), SS slowly rolls down the potential, and when HmSH\sim m_{S} is reached, it starts to oscillate. However, when the non-minimal coupling ξS\xi_{S} is negatively non-zero (right), SS decreases faster than the minimally coupled counterpart and becomes more suppressed. On the other hand, when the non-minimal coupling ξS\xi_{S} is positively non-zero (left), SS grows during the initial inflationary period.

Refer to caption
Figure 3: Evolution of HH as a function of the number of ee-folds NN for the case of ξS=0.005\xi_{S}=0.005 with several values of S0S_{0}. The values of S0S_{0} given in the legend are in units of MplM_{\rm pl}.

Even if the initial value of SS is smaller than MplM_{\rm pl}, SS gets increased to super-Planckian field values due to the existence of the spectator non-minimal coupling. Indeed, due to this behavior, the SS field can dominate the Universe after the end of the first inflationary phase while it is still rolling slowly, and can drive the second inflationary phase. This can be confirmed by the evolution of the Hubble parameter in Fig. 3. As mentioned above, Φ\Phi is initially the dominant component and drives the first period of inflationary expansion. This corresponds to the early phase where the Hubble parameter HH in Fig. 3 almost stays constant. Once inflation ends, HH sharply drops, which happens at N90100N\sim 90\text{--}100 in the figure. However, for the cases of S0=MplS_{0}=M_{\rm pl} and S0=0.5MplS_{0}=0.5M_{\mathrm{pl}}, HH flattens out again, indicating that the second inflationary stage occurs. This second inflation is indeed caused by the initial growth of SS seen in the lower-left panel of Fig. 2.

Refer to caption
Figure 4: Evolution of ϵH\epsilon_{H} as a function of the number of ee-folds NN for the case of ξS=0.005\xi_{S}=0.005 with several values of S0S_{0}. The values of S0S_{0} given in the legend are in units of MplM_{\rm pl}.

Fig. 4 shows the evolution of the slow-roll parameter ϵHH˙/H2\epsilon_{H}\equiv-\dot{H}/H^{2}. As can be seen from the figure, ϵH\epsilon_{H} grows steadily during the Φ\Phi-driven stage until it exceeds unity and then decreases again as the second inflation begins to be driven by the amplified SS field. However, if the initial amplitude is not sufficiently large, as in the case of S0=0.1MplS_{0}=0.1M_{\mathrm{pl}}, the field does not have enough time in the first inflationary stage to become dominant enough to drive the second inflationary stage. On the other hand, if the initial value is too large, it will be amplified beyond the point of no return, as we saw from the analytic discussion. We have confirmed this behavior numerically.

Refer to caption
Figure 5: The number of ee-folds generated during the second inflation N2N_{2} for mass ratio mΦ/mS=10m_{\Phi}/m_{S}=10. In the lower left region (white), there is no second inflation generated by the SS field as it never becomes dominant. In the upper right region (salmon), SS is amplified past the point of no return.

As a result of this dynamics, there is a range of initial spectator amplitude and non-minimal coupling values where the second inflationary stage can be realized. We have performed a numerical scan of the parameter space to determine the viable region which is shown in Fig. 5.555 The mass ratio in the numerical simulation was taken to be mΦ/mS=10m_{\Phi}/m_{S}=10. As mentioned earlier, the dynamics of the second inflationary stage is insensitive to the initial inflaton amplitude and absolute mass values. As seen from the figure, even when the initial amplitude of SS is sub-Planckian, the second inflation can be generated, and its ee-folding number N2N_{2} can also be very large, N260N_{2}\gtrsim 60.666 In fact, it can be tuned to be arbitrarily large on the edge from the runaway as we can choose initial values such that the field ends up balancing on the ridge of the Einstein frame potential. If the field overshoots and ends up in the runaway phase, the second phase of inflation never terminates. Such cases correspond to the upper right region in the figure. The lower white corner corresponds to cases where the SS field never becomes dominant.

4 Evolution of fluctuations and power spectra

In this section, we discuss the evolution of fluctuations in our setup. First, we briefly summarize how we calculate the perturbations in the model. We use the formalism of [50, 37, 39] for decomposing the fluctuations into adiabatic and entropic components as follows. A scalar field φI\varphi^{I} can be divided into the background part φ¯I\bar{\varphi}^{I} and its perturbation QIQ^{I} as

φI=φ¯I+QI.\varphi^{I}=\bar{\varphi}^{I}+Q^{I}\,. (4.1)

We can define a new basis {enI}(n=σ,s)\{e_{n}^{I}\}~{}(n=\sigma,{s}) such that the first vector points along the background trajectory, corresponding to the adiabatic mode whose basis vector can be written as

eσIφ˙IGJKφ˙Jφ˙K=φ˙Iσ˙,e_{\sigma}^{I}\equiv\frac{\dot{\varphi}^{I}}{\sqrt{G_{JK}\dot{\varphi}^{J}\dot{\varphi}^{K}}}=\frac{\dot{\varphi}^{I}}{\dot{\sigma}}\,, (4.2)

where we have defined σ˙GIJφ˙Iφ˙J\dot{\sigma}\equiv\sqrt{G_{IJ}\dot{\varphi}^{I}\dot{\varphi}^{J}}. The remaining isocurvature direction is then perpendicular to the background trajectory. The perturbations QIQ^{I} in the new basis are written as QI=QnenIQ^{I}=Q^{n}e_{n}^{I} with n=σ,sn=\sigma,s.

The equations of motion for the adiabatic and isocurvature modes QσQ_{\sigma} and QsQ_{s} are

Q¨σ+3HQ˙σ+(k2a2+μσ2)Qσ\displaystyle\ddot{Q}_{\sigma}+3H\dot{Q}_{\sigma}+\left(\frac{k^{2}}{a^{2}}+\mu_{\sigma}^{2}\right)Q_{\sigma} =ddt(ΞQs)(H˙H+V,σσ˙)ΞQs\displaystyle=\frac{\mathrm{d}}{\mathrm{d}t}\Big{(}\Xi Q_{s}\Big{)}-\left(\frac{\dot{H}}{H}+\frac{V_{,\sigma}}{\dot{\sigma}}\right)\Xi Q_{s} (4.3)
and
Q¨s+3HQ˙s+(k2a2+μs2)Qs\displaystyle\ddot{Q}_{s}+3H\dot{Q}_{s}+\left(\frac{k^{2}}{a^{2}}+\mu_{s}^{2}\right)Q_{s} =Ξ[Q˙σ+(H˙Hσ¨σ˙)Qσ],\displaystyle=-\Xi\left[\dot{Q}_{\sigma}+\left(\frac{\dot{H}}{H}-\frac{\ddot{\sigma}}{\dot{\sigma}}\right)Q_{\sigma}\right], (4.4)

respectively, where μσ\mu_{\sigma}, μs\mu_{s}, and Ξ\Xi are given by

μσ2\displaystyle\mu_{\sigma}^{2} =\displaystyle= V;σσ1a3ddt(a3σ˙2H)V,s2σ˙2,\displaystyle V_{;\sigma\sigma}-\frac{1}{a^{3}}\frac{\mathrm{d}}{\mathrm{d}t}\left(\frac{a^{3}\dot{\sigma}^{2}}{H}\right)-\frac{V_{,{s}}^{2}}{\dot{\sigma}^{2}}, (4.5)
μs2\displaystyle\mu_{s}^{2} =\displaystyle= V;ss+12σ˙2~V,s2σ˙2,\displaystyle V_{;{ss}}+\frac{1}{2}\dot{\sigma}^{2}\tilde{\mathcal{R}}-\frac{V_{,{s}}^{2}}{\dot{\sigma}^{2}}, (4.6)
Ξ\displaystyle\Xi =\displaystyle= 2V,sσ˙,\displaystyle-\frac{2V_{,{s}}}{\dot{\sigma}}\,, (4.7)

in which V;σσV_{;\sigma\sigma} and V;ssV_{;{ss}} are defined as projections of 𝒟I𝒟JV\mathcal{D}_{I}\mathcal{D}_{J}V onto the appropriate directions and ~\tilde{\mathcal{R}} is the Ricci curvature scalar of the field manifold defined by GIJG_{IJ}. The comoving curvature perturbation in the flat gauge is given as

=Hσ˙Qσ.\mathcal{R}=\frac{H}{\dot{\sigma}}Q_{\sigma}\,. (4.8)

The evolution of {\cal R} can be calculated by solving the equations of motion (4.3) and (4.4) along with the background ones. The primordial (scalar) power spectrum 𝒫{\cal P_{R}} is computed by evaluating {\cal R} at the end of inflation and given by

𝒫ζ(k)\displaystyle{\cal P}_{\zeta}(k) =\displaystyle= k32π2|𝐤|2,\displaystyle\frac{k^{3}}{2\pi^{2}}|{\cal R_{\mathbf{k}}}|^{2}\,, (4.9)

while its scale dependence is characterized by the spectral index

ns1dlog𝒫ζdlogk.n_{s}-1\equiv\frac{d\log{\cal P_{\zeta}}}{d\log k}\,. (4.10)

We also discuss primordial gravitational waves generated during inflation. Gravitational waves correspond to the tensor perturbations whose power spectrum can be calculated with the standard procedure (see, e.g., [55]) as

𝒫T(k)=k32π2|vk|2,{\cal P}_{T}(k)=\frac{k^{3}}{2\pi^{2}}|v_{k}|^{2}\,, (4.11)

where vkv_{k} is the Fourier component of vij=Mplhij/(2a)v_{ij}=M_{\rm pl}h_{ij}/(2a) in a plane wave expansion, with hijh_{ij} being the tensor perturbation in the FRW metric. vkv_{k} satisfies the following equation of motion:

vk¨+Hvk˙+(k2a2H2a¨a)vk=0.\ddot{v_{k}}+H\dot{v_{k}}+\left(\frac{k^{2}}{a^{2}}-H^{2}-\frac{\ddot{a}}{a}\right)v_{k}=0. (4.12)

To quantify the size of the tensor power spectrum, the tensor-to-scalar ratio rr is commonly adopted, which can be compared with observational constraints. rr is defined at some reference scale k=kk=k_{\ast} as

r𝒫T𝒫ζ|k=k.r\equiv\left.\frac{{\cal P}_{T}}{{\cal P}_{\zeta}}\right|_{k=k_{\ast}}\,. (4.13)

4.1 Numerical analysis

In our numerical analysis, we solve the fluctuation equations of motion (4.3), (4.4), and (4.12) altogether. We evolve the fluctuations from an initial time inside the horizon, through the first inflationary stage, to the end of the second inflation. We identify the reference scale kk_{\ast} that has exited the horizon 60 ee-folds before the end of the second stage of inflation. We normalize the masses to match the observed spectrum amplitude 𝒫2.2×109\mathcal{P}_{\mathcal{R}}\simeq 2.2\times 10^{-9} at the reference scale, with the mass ratio being fixed.

Following [56, 57, 58], we solve the evolution of fluctuations twice: once with Qs=0Q_{s}=0 and the initial Bunch-Davies vacuum amplitude Qσ=1/2kQ_{\sigma}=1/\sqrt{2k} and once more with the opposite choice. This ensures the absence of spurious cross-correlations in power spectra. We have also solved the evolution of fluctuations in the original basis (Φ,S)(\Phi,S) and verified that the results agree as long as this technique is employed.

4.2 Scalar and tensor power spectra

In this section, we discuss the predictions for scalar and tensor power spectra in the model. To obtain the spectra we scan through the momentum space, solving the evolution of fluctuations as outlined above for each kk-mode. We focus specifically on the case where the second inflation is driven by the SS field which has grown to reach a super-Planckian value due to the existence of the non-minimal coupling.777 For numerical studies of power spectra in double inflation with two quadratic potentials in the minimally coupled case, see e.g. [59, 58].

Refer to caption
Refer to caption
Figure 6: Scalar and tensor power spectra (right and left panels, respectively) for several values of S0S_{0}. ξS\xi_{S} is fixed by requiring that the second inflation lasts for N2=10N_{2}=10, which gives ξS=0.005\xi_{S}=0.005, 0.011, 0.016, and 0.022 for S0=MplS_{0}=M_{\rm pl}, 101Mpl10^{-1}M_{\rm pl}, 102Mpl10^{-2}M_{\rm pl}, and 103Mpl10^{-3}M_{\rm pl}, respectively. The ratio of the masses is taken to be mΦ/mS=10m_{\Phi}/m_{S}=10. The value of mΦm_{\Phi} (or mSm_{S}) is fixed by requiring that the scalar power spectrum at log(k/kt)=50\log(k_{\ast}/k_{t})=-50 is normalized to 𝒫=2.2×109{\cal P}_{\cal R}=2.2\times 10^{-9}.

In Fig. 6, we show scalar and tensor power spectra for several values of S0S_{0}. The value of ξS\xi_{S} is chosen such that the second inflationary stage lasts N2=10N_{2}=10 ee-folds for each S0S_{0}. The resultant values are respectively ξS=0.005\xi_{S}=0.005, 0.011, 0.016, and 0.022 for S0=MplS_{0}=M_{\rm pl}, 101Mpl10^{-1}M_{\rm pl}, 102Mpl10^{-2}M_{\rm pl}, and 103Mpl10^{-3}M_{\rm pl}. We have fixed the mass ratio to be mΦ/mS=10m_{\Phi}/m_{S}=10. The power spectra are plotted as a function of log(k/kt)\log(k/k_{t}), where ktk_{t} is the wave number of the mode that exited the horizon at the beginning of the second inflation. The size of mΦm_{\Phi} is fixed by requiring that the scalar power spectrum at the reference scale kk_{\ast}, given as log(k/kt)=50\log(k_{\ast}/k_{t})=-50 in the current choice, is normalized to 𝒫=2.2×109{\cal P}_{\cal R}=2.2\times 10^{-9}.

Both scalar and tensor power spectra are suppressed for the modes k>ktk>k_{t} having exited the horizon during the second inflation. This is mainly because the first inflation, during which the modes k<ktk<k_{t} exited the horizon, is driven by Φ\Phi whose mass is assumed to be larger than that of SS. Although the shape of the power spectra generally varies depending on the value of S0S_{0} (and ξS\xi_{S}), the main overall feature remains the same regardless of the values of the parameters.

Regarding the behavior of the tensor power spectrum 𝒫T{\cal P}_{T}, effects of the spectator non-minimal coupling ξS\xi_{S} can be understood by noticing that, when ξS\xi_{S} gets larger (or S0S_{0} gets smaller to keep N2=10N_{2}=10), the adiabatic fluctuations are more sourced by the isocurvature ones, in which the intrinsic adiabatic mode should give less contribution to the total one and the mass of Φ\Phi is lowered to satisfy the normalization condition for 𝒫{\cal P_{R}}. This results in the decrease of the Hubble parameter during inflation. Since the tensor power spectrum is given as 𝒫T(Hinf/Mpl)2{\cal P}_{T}\propto(H_{\rm inf}/M_{\rm pl})^{2} with HinfH_{\rm inf} being the Hubble parameter during inflation, the amplitude of the tensor power spectrum is lowered as S0S_{0} decreases. Furthermore, sourcing of isocurvature fluctuations is more prominent on large scales since the sourcing continues longer for the modes which exited the horizon earlier. Since we fix the number of the second inflation to be N2=10N_{2}=10, the case with smaller S0S_{0} experiences more sourcing of isocurvature fluctuation, and thus 𝒫{\cal P_{R}} for smaller S0S_{0} gets more red-tilted.

Refer to caption
Refer to caption
Figure 7: Scalar and tensor power spectra (right and left panels, respectively) for several values of the ratio mΦ/mS=1m_{\Phi}/m_{S}=1, 5, 10, and 15. We take S0=0.1MplS_{0}=0.1M_{\rm pl} and ξS0.01\xi_{S}\simeq 0.01, such that the number of ee-folds during the second inflation is N2=10N_{2}=10. The value of mΦm_{\Phi} is fixed by requiring that the scalar power spectrum is normalized at log(k/kt)=50\log(k_{\ast}/k_{t})=-50 as Pζ=2.2×109P_{\zeta}=2.2\times 10^{-9}.

In Fig. 7, scalar and tensor power spectra for several values of the ratio mΦ/mSm_{\Phi}/m_{S} are shown. The values of S0S_{0} is fixed at S0=0.1MplS_{0}=0.1M_{\rm pl} and ξS\xi_{S} is once again fixed by requiring that the number of ee-folds during the second inflation is N2=10N_{2}=10. However, in practice, ξS0.01\xi_{S}\simeq 0.01 for all cases shown. As can be seen from the figure, the ratio mΦ/mSm_{\Phi}/m_{S} mainly affects the power spectra for modes with k>ktk>k_{t} which exited the horizon in the second inflation driven by SS. Since mΦm_{\Phi} is effectively fixed by the normalization to the observed spectrum, increasing the ratio amounts to decreasing mSm_{S}. Therefore, for the modes k>ktk>k_{t}, the situation is the same as the one for a single-field quadratic chaotic inflation without a non-minimal coupling, and lowering mSm_{S} means that the scalar power spectrum is suppressed. The tensor mode also gets decreased since the Hubble rate during the second inflation is lowered. Although both power spectra 𝒫{\cal P_{R}} and 𝒫T{\cal P}_{T} seem to be almost unchanged for k<ktk<k_{t}, when the mass ratio is decreased to be close to unity, adiabatic perturbations are more sourced by isocurvature ones, which indicates that 𝒫{\cal P_{R}} gets more red-tilted. Furthermore, this means that the tensor-to-scalar ratio is suppressed since only the scalar mode is enhanced. This trend is more clearly seen in terms of the spectral index and the tensor-to-scalar ratio as we discuss below.

4.3 Spectral index and tensor to scalar ratio

We calculate nsn_{s} and rr at the reference scale corresponding to the mode which exited the horizon at log(k/kt)=50\log(k/k_{t})=-50. Since we calculate the power spectrum 𝒫{\cal P}_{\cal R} numerically, the spectral index, defined in Eq. (4.10), is obtained by numerically evaluating the derivative of 𝒫{\cal P}_{\cal R} with respect to kk. For the tensor-to-scalar ratio rr, defined in Eq. (4.13), we take the ratio of 𝒫{\cal P_{R}} and 𝒫T{\cal P}_{T} at the above mentioned reference scale.

Refer to caption
Figure 8: Spectral index nsn_{s} vs. the tensor-to-scalar ratio rr predicted in the model varying ξSS0\sqrt{\xi_{S}}S_{0}. The values of ξS\xi_{S} and S0S_{0} are chosen to keep N2=10N_{2}=10 unchanged. The mass ratio is fixed as mΦ/mS=10m_{\Phi}/m_{S}=10. Observational constraints on nsn_{s} and rr from Planck [1] and Planck+BICEP/Keck 2018 [2] are also depicted.

In Fig. 8, the predictions for nsn_{s} and rr varying ξSS0\sqrt{\xi_{S}}S_{0} are shown. Since the combination of ξSS2\xi_{S}S^{2} appears in the function f(S)f(S) which represents a non-minimal coupling in the action, this is the main measure of the effect, though the additional difference can be introduced by the kinetic mixing. In practice, we vary the value of S0S_{0} and fix ξS\xi_{S} to keep the duration of the second inflation at N2=10N_{2}=10. As already discussed earlier, as S0S_{0} (or ξSS0\sqrt{\xi_{S}}S_{0}) decreases, 𝒫{\cal P_{R}} gets red-tilted and 𝒫T{\cal P}_{T} is suppressed, which explains the behavior of the predictions in the nsn_{s}rr plane. Notice that the limit of ξSS00\sqrt{\xi_{S}}S_{0}\rightarrow 0 does not correspond to a minimally coupled case because here we fix the number of ee-folds during the second inflation to be N2=10N_{2}=10. In this case, even when ξSS0\sqrt{\xi_{S}}S_{0} is assumed to be small, we need to take a large value of ξS\xi_{S} even for a small value of S0S_{0}, which is the reason why nsn_{s} gets more red-tilted for a smaller ξSS0\sqrt{\xi_{S}}S_{0}.

Refer to caption
Figure 9: Dependence of nSn_{S} and rr on the mass ratio mΦ/mSm_{\Phi}/m_{S} with S0=1.5MplS_{0}=1.5M_{\rm pl}. ξS\xi_{S} is chosen such that the duration of the second inflation is N2=10N_{2}=10.

The dependence of nsn_{s} and rr on the mass ratio mΦ/mSm_{\Phi}/m_{S} is shown in Fig. 9. The value of S0S_{0} is fixed as S0=1.5MplS_{0}=1.5M_{\rm pl} and the number of ee-folds during the second inflation is again taken to be 10, which determines the value of ξS\xi_{S}. When the ratio is mΦ/mS7m_{\Phi}/m_{S}\gtrsim 7, ξS\xi_{S} is almost unchanged and takes the value ξS0.0036\xi_{S}\simeq 0.0036. As already mentioned, when mΦ/mSm_{\Phi}/m_{S} is decreased to be close to unity, adiabatic perturbations are more sourced by the isocurvature ones, which induces more scale-dependence and thus the spectral index becomes smaller. When the sourcing from isocurvature fluctuations is stronger, 𝒫{\cal P_{R}} is enhanced, and the tensor-to-scalar ratio decreases. Note that larger values of mΦ/mSm_{\Phi}/m_{S} are favorable since as the ratio increases, the predictions in the nsn_{s}rr plane move to the center of the current cosmological constraints. Although the tensor-to-scalar ratio is generically suppressed in our setup, the scalar power spectrum becomes more red-tilted when more adiabatic fluctuations are sourced by the isocurvature mode, and thus the mass ratio should be taken not to be inconsistent with the measurement of nsn_{s}.

Refer to caption
Figure 10: Dependence of nsn_{s} and rr on duration of second inflation N2N_{2}. The initial value of SS is fixed as S0=4MplS_{0}=4M_{\rm pl} and ξS\xi_{S} is chosen such that the value of N2N_{2} depicted in the color bar is realized. The mass ratio is taken to be mΦ/mS=10m_{\Phi}/m_{S}=10.

Finally, we show the predictions of our model in the nsn_{s}rr plane, varying the number of ee-folds during the second inflation. In the figure, the mass ratio is fixed to be mΦ/mS=10m_{\Phi}/m_{S}=10 and the initial value of SS to be S0=4MplS_{0}=4M_{\rm pl}. The value of ξS\xi_{S} is then fixed by requiring a specific value of N2N_{2}. To have a large N2N_{2}, we need to have a bigger ξS\xi_{S} since more growth of SS is necessary to have a longer period of the second inflation. Thus, when N2N_{2} is larger, the effects of isocurvature field SS becomes bigger (more adiabatic perturbations are sourced by SS). Therefore large N2N_{2} indicates that 𝒫{\cal P_{R}} is more red-tilted and rr is more suppressed.

In Figs. 8, 9, and 10, observational constraints from Planck [1] as well as from Planck+BICEP/Keck 2018 [2] are superimposed. As can be seen, the spectral index and the tensor-to-scalar ratio are predicted to be well within the current cosmological constraints for a range of model parameters, indicating that a quadratic chaotic inflation model can be rescued by introducing a spectator field non-minimally coupled to gravity.

5 Conclusion and Discussion

In this paper, we have investigated the two-field inflation model where the field SS that initially plays the role of spectator comes to dominate and causes the second inflation due to the non-minimal coupling to gravity. We have first studied the evolution of the background quantities of the fields and have shown that the second inflationary phase is generically realized due to the existence of the spectator non-minimal coupling ξS\xi_{S}. When ξS\xi_{S} is positive, SS field can grow during the first inflationary phase driven by Φ\Phi and once SS acquires a super-Planckian field value, it can generate the second inflation. However, when ξS\xi_{S} or the initial amplitude is too large a runaway behavior can result where the field keeps growing and the second inflation does not terminate. However, in a sizable region of parameter space, the SS field drives a finite amount of inflation and then returns to zero.

We also investigated the power spectra for both the scalar and tensor modes in the model. As discussed in Sec. 4, adiabatic fluctuations are more sourced by isocurvature ones as ξSS2\xi_{S}S^{2} (or ξS\xi_{S}) increases. Due to this isocurvature sourcing augmented by the spectator non-minimal coupling, the scale-dependence of the scalar power spectrum is modified. Although the tensor power spectrum itself is not directly affected by this isocurvature sourcing, the overall amplitude can change through the requirement of the normalization condition imposed on the scalar power spectrum. But due to the change of the scalar mode, the tensor-to-scalar ratio rr can also be affected. We explicitly showed nsn_{s} and rr are modified and can fit well within the current observational constraints. With this setup, the quadratic chaotic inflation model can become viable by assuming a spectator with non-minimal coupling to gravity.

Multi-field dynamics is common to high energy models and non-minimal coupling to gravity is a generic feature of quantum field theory in curved space-time. Therefore, it is necessary to consider the interplay of such features in the early Universe. As we have shown, even simple dynamics can have a significant impact on cosmological predictions and change the picture considerably from the single-field case.

Acknowledgements

This work is supported in part by the Communidad de Madrid “Atracción de Talento investigador” Grant No. 2017-T1/TIC-5305 and the MICINN (Spain) project PID2019-107394GB-I00 (SR) as well as by JSPS KAKENHI Grant Numbers 19H01899 (KO), 21H01107 (KO), 17H01131 (TT), and 19K03874 (TT) and MEXT KAKENHI Grant Number 19H05110 (TT).

References

  • [1] Planck Collaboration, Y. Akrami et al., Planck 2018 results. X. Constraints on inflation, Astron. Astrophys. 641 (2020) A10, [arXiv:1807.06211].
  • [2] BICEP, Keck Collaboration, P. A. R. Ade et al., Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season, Phys. Rev. Lett. 127 (2021), no. 15 151301, [arXiv:2110.00483].
  • [3] K. Enqvist and M. S. Sloth, Adiabatic CMB perturbations in pre - big bang string cosmology, Nucl. Phys. B626 (2002) 395–409, [hep-ph/0109214].
  • [4] D. H. Lyth and D. Wands, Generating the curvature perturbation without an inflaton, Phys. Lett. B524 (2002) 5–14, [hep-ph/0110002].
  • [5] T. Moroi and T. Takahashi, Effects of cosmological moduli fields on cosmic microwave background, Phys. Lett. B522 (2001) 215–221, [hep-ph/0110096]. [Erratum: Phys. Lett.B539,303(2002)].
  • [6] G. Dvali, A. Gruzinov, and M. Zaldarriaga, A new mechanism for generating density perturbations from inflation, Phys. Rev. D69 (2004) 023505, [astro-ph/0303591].
  • [7] L. Kofman, Probing string theory with modulated cosmological fluctuations, astro-ph/0303614.
  • [8] D. Langlois and F. Vernizzi, Mixed inflaton and curvaton perturbations, Phys. Rev. D 70 (2004) 063522, [astro-ph/0403258].
  • [9] T. Moroi, T. Takahashi, and Y. Toyoda, Relaxing constraints on inflation models with curvaton, Phys. Rev. D 72 (2005) 023502, [hep-ph/0501007].
  • [10] T. Moroi and T. Takahashi, Implications of the curvaton on inflationary cosmology, Phys. Rev. D 72 (2005) 023505, [astro-ph/0505339].
  • [11] K. Ichikawa, T. Suyama, T. Takahashi, and M. Yamaguchi, Non-Gaussianity, Spectral Index and Tensor Modes in Mixed Inflaton and Curvaton Models, Phys. Rev. D 78 (2008) 023513, [arXiv:0802.4138].
  • [12] K. Enqvist and T. Takahashi, Mixed Inflaton and Spectator Field Models after Planck, JCAP 10 (2013) 034, [arXiv:1306.5958].
  • [13] V. Vennin, K. Koyama, and D. Wands, Encyclopædia curvatonis, JCAP 11 (2015) 008, [arXiv:1507.07575].
  • [14] T. Fujita, M. Kawasaki, and S. Yokoyama, Curvaton in large field inflation, JCAP 09 (2014) 015, [arXiv:1404.0951].
  • [15] K. Ichikawa, T. Suyama, T. Takahashi, and M. Yamaguchi, Primordial Curvature Fluctuation and Its Non-Gaussianity in Models with Modulated Reheating, Phys. Rev. D 78 (2008) 063545, [arXiv:0807.3988].
  • [16] D. S. Salopek, J. R. Bond, and J. M. Bardeen, Designing Density Fluctuation Spectra in Inflation, Phys. Rev. D 40 (1989) 1753.
  • [17] F. L. Bezrukov and M. Shaposhnikov, The Standard Model Higgs boson as the inflaton, Phys. Lett. B659 (2008) 703–706, [arXiv:0710.3755].
  • [18] B. L. Spokoiny, Inflation and generation of perturbations in broken symmetric theory of gravity, Phys. Lett. 147B (1984) 39–43.
  • [19] T. Futamase and K.-i. Maeda, Chaotic Inflationary Scenario in Models Having Nonminimal Coupling With Curvature, Phys. Rev. D39 (1989) 399–404.
  • [20] R. Fakir and W. G. Unruh, Improvement on cosmological chaotic inflation through nonminimal coupling, Phys. Rev. D41 (1990) 1783–1791.
  • [21] E. Komatsu and T. Futamase, Complete constraints on a nonminimally coupled chaotic inflationary scenario from the cosmic microwave background, Phys. Rev. D59 (1999) 064029, [astro-ph/9901127].
  • [22] A. Linde, M. Noorbala, and A. Westphal, Observational consequences of chaotic inflation with nonminimal coupling to gravity, JCAP 1103 (2011) 013, [arXiv:1101.2652].
  • [23] L. Boubekeur, E. Giusarma, O. Mena, and H. Ramírez, Does Current Data Prefer a Non-minimally Coupled Inflaton?, Phys. Rev. D 91 (2015) 103004, [arXiv:1502.05193].
  • [24] T. Tenkanen, Resurrecting Quadratic Inflation with a non-minimal coupling to gravity, JCAP 12 (2017) 001, [arXiv:1710.02758].
  • [25] R. Z. Ferreira, A. Notari, and G. Simeon, Natural Inflation with a periodic non-minimal coupling, JCAP 11 (2018) 021, [arXiv:1806.05511].
  • [26] M. Artymowski, Z. Lalak, and K.-Y. Oda, Hill-climbing dark inflation, arXiv:1807.06830.
  • [27] I. Antoniadis, A. Karam, A. Lykkas, T. Pappas, and K. Tamvakis, Rescuing Quartic and Natural Inflation in the Palatini Formalism, JCAP 03 (2019) 005, [arXiv:1812.00847].
  • [28] T. Takahashi and T. Tenkanen, Towards distinguishing variants of non-minimal inflation, JCAP 04 (2019) 035, [arXiv:1812.08492].
  • [29] R. Jinno, K. Kaneta, K.-y. Oda, and S. C. Park, Hillclimbing inflation in metric and Palatini formulations, Phys. Lett. B 791 (2019) 396–402, [arXiv:1812.11077].
  • [30] M. Shokri, F. Renzi, and A. Melchiorri, Cosmic Microwave Background constraints on non-minimal couplings in inflationary models with power law potentials, Phys. Dark Univ. 24 (2019) 100297, [arXiv:1905.00649].
  • [31] T. Takahashi, T. Tenkanen, and S. Yokoyama, Violation of slow-roll in nonminimal inflation, Phys. Rev. D 102 (2020), no. 4 043524, [arXiv:2003.10203].
  • [32] Y. Reyimuaji and X. Zhang, Natural inflation with a nonminimal coupling to gravity, JCAP 03 (2021) 059, [arXiv:2012.14248].
  • [33] D. Y. Cheong, S. M. Lee, and S. C. Park, Reheating in Models with Non-minimal Coupling in metric and Palatini formalisms, arXiv:2111.00825.
  • [34] T. Kodama and T. Takahashi, Relaxing inflation models with non-minimal coupling: A general study, arXiv:2112.05283.
  • [35] N. Birrell and P. Davies, Quantum Fields in Curved Space. Cambridge Monographs on Mathematical Physics. Cambridge Univ. Press, Cambridge, UK, 2, 1984.
  • [36] L. E. Parker and D. Toms, Quantum Field Theory in Curved Spacetime: Quantized Field and Gravity. Cambridge Monographs on Mathematical Physics. Cambridge University Press, 8, 2009.
  • [37] D. I. Kaiser and A. T. Todhunter, Primordial Perturbations from Multifield Inflation with Nonminimal Couplings, Phys. Rev. D 81 (2010) 124037, [arXiv:1004.3805].
  • [38] J. White, M. Minamitsuji, and M. Sasaki, Curvature perturbation in multi-field inflation with non-minimal coupling, JCAP 07 (2012) 039, [arXiv:1205.0656].
  • [39] D. I. Kaiser, E. A. Mazenc, and E. I. Sfakianakis, Primordial Bispectrum from Multifield Inflation with Nonminimal Couplings, Phys. Rev. D 87 (2013) 064004, [arXiv:1210.7487].
  • [40] R. N. Greenwood, D. I. Kaiser, and E. I. Sfakianakis, Multifield Dynamics of Higgs Inflation, Phys. Rev. D 87 (2013) 064021, [arXiv:1210.8190].
  • [41] D. I. Kaiser and E. I. Sfakianakis, Multifield Inflation after Planck: The Case for Nonminimal Couplings, Phys. Rev. Lett. 112 (2014), no. 1 011302, [arXiv:1304.0363].
  • [42] J. White, M. Minamitsuji, and M. Sasaki, Non-linear curvature perturbation in multi-field inflation models with non-minimal coupling, JCAP 09 (2013) 015, [arXiv:1306.6186].
  • [43] K. Schutz, E. I. Sfakianakis, and D. I. Kaiser, Multifield Inflation after Planck: Isocurvature Modes from Nonminimal Couplings, Phys. Rev. D 89 (2014), no. 6 064044, [arXiv:1310.8285].
  • [44] S. Karamitsos and A. Pilaftsis, Frame Covariant Nonminimal Multifield Inflation, Nucl. Phys. B 927 (2018) 219–254, [arXiv:1706.07011].
  • [45] L.-H. Liu and T. Prokopec, Non-minimally coupled curvaton, JCAP 06 (2021) 033, [arXiv:2005.11069].
  • [46] T. Modak and K.-y. Oda, Echoes of 2HDM inflation at the collider experiments, Eur. Phys. J. C 80 (2020), no. 9 863, [arXiv:2007.08141]. [Erratum: Eur.Phys.J.C 81, 518 (2021)].
  • [47] S. M. Lee, T. Modak, K.-y. Oda, and T. Takahashi, The R2R^{2}-Higgs inflation with two Higgs doublets, Eur. Phys. J. C 82 (2022), no. 1 18, [arXiv:2108.02383].
  • [48] K. Dimopoulos, K. Kohri, D. H. Lyth, and T. Matsuda, The inflating curvaton, JCAP 03 (2012) 022, [arXiv:1110.2951].
  • [49] K. Enqvist, T. Sawala, and T. Takahashi, Structure formation with two periods of inflation: beyond PLaIn Λ\LambdaCDM, JCAP 10 (2020) 053, [arXiv:1905.13580].
  • [50] D. Langlois and S. Renaux-Petel, Perturbations in generalized multi-field inflation, JCAP 04 (2008) 017, [arXiv:0801.1085].
  • [51] F. Bauer and D. A. Demir, Inflation with Non-Minimal Coupling: Metric versus Palatini Formulations, Phys. Lett. B 665 (2008) 222–226, [arXiv:0803.2664].
  • [52] S. Rasanen and P. Wahlman, Higgs inflation with loop corrections in the Palatini formulation, JCAP 11 (2017) 047, [arXiv:1709.07853].
  • [53] L. Jarv, A. Racioppi, and T. Tenkanen, Palatini side of inflationary attractors, Phys. Rev. D97 (2018), no. 8 083513, [arXiv:1712.08471].
  • [54] R. Jinno, M. Kubota, K.-y. Oda, and S. C. Park, Higgs inflation in metric and Palatini formalisms: Required suppression of higher dimensional operators, JCAP 03 (2020) 063, [arXiv:1904.05699].
  • [55] J. E. Lidsey, A. R. Liddle, E. W. Kolb, E. J. Copeland, T. Barreiro, and M. Abney, Reconstructing the inflation potential : An overview, Rev. Mod. Phys. 69 (1997) 373–410, [astro-ph/9508078].
  • [56] Z. Lalak, D. Langlois, S. Pokorski, and K. Turzynski, Curvature and isocurvature perturbations in two-field inflation, JCAP 07 (2007) 014, [arXiv:0704.0212].
  • [57] C. van de Bruck and M. Robinson, Power Spectra beyond the Slow Roll Approximation in Theories with Non-Canonical Kinetic Terms, JCAP 08 (2014) 024, [arXiv:1404.7806].
  • [58] M. Braglia, D. K. Hazra, L. Sriramkumar, and F. Finelli, Generating primordial features at large scales in two field models of inflation, JCAP 08 (2020) 025, [arXiv:2004.00672].
  • [59] B. Feng and X. Zhang, Double inflation and the low cmb quadrupole, Phys. Lett. B 570 (2003) 145–150, [astro-ph/0305020].