The Effective Field Theory of nonsingular cosmology: II

TL;DR

This paper clarifies and resolves instabilities in nonsingular cosmological models using Effective Field Theory, introduces a Genesis-inflation model, and explores its implications for primordial perturbations and CMB anomalies.

Contribution

It identifies and cures instabilities in nonsingular cubic Galileon models within EFT, proposing a new Genesis-inflation model with potential observational signatures.

Findings

Resolved instabilities in nonsingular models

Constructed a Genesis-inflation scenario

Predicted a large-scale cutoff in primordial spectrum

Abstract

Based on the Effective Field Theory (EFT) of cosmological perturbations, we explicitly clarify the pathology in nonsingular cubic Galileon models and show how to cure it in EFT with new insights into this issue. With the least set of EFT operators that are capable to avoid instabilities in nonsingular cosmologies, we construct a nonsingular model dubbed the Genesis-inflation model, in which a slowly expanding phase (namely, Genesis) with increasing energy density is followed by slow-roll inflation. The spectrum of the primordial perturbation may be simulated numerically, which shows itself a large-scale cutoff, as the large-scale anomalies in CMB might be a hint for.

Tables

Table 1. Table 1 : Pathology in nonsingular cubic Galileon cosmological models and its cure in EFT by either Q T subscript 𝑄 𝑇 Q_{T} or Q m ~ 4 subscript 𝑄 subscript ~ 𝑚 4 Q_{\tilde{m}_{4}} .

Nonsingular cubic Galileon models
	Initially	$\gamma >0$ throughout
Crossing 0 for $\gamma$?	$\surd$	$\times$
is inevitable (“no-go”)?	$\surd$	$\surd$
Phase with (Pathological phase)	$\dot{\gamma }>0$ around $\gamma \simeq 0$
Curing pathology in EFT
Conditions of $c_s^2>0$	(11)	(17)
Applying $Q_T$	(12)	(18)
Applying $Q_{{\tilde{m}}_4}$	(13)	(19)

Equations

$$S_{\zeta }^{(2)}=\int d^{4}x{a}^3{c}_1\left[{\dot{\zeta }}^2-c_s^2\frac{(\partial \zeta {)}^2}{a^2}\right],$$

$$c_s^2=\left(\frac{{\dot{c}}_3}{a}-c_2\right)/c_1,$$

$$c_{3}\Big{|}_{t_{f}}-c_{3}\Big{|}_{t_{i}}>\int_{t_{i}}^{t_{f}}ac_{2}dt\,.$$

$$S_{\gamma }^{(2)}=\frac{M_p^2}{8}\int d^{4}x{a}^3{Q}_T\left[{\dot{\gamma }}_{ij}^2-c_T^2\frac{({\partial }_k{\gamma }_{ij}{)}^2}{a^2}\right],$$

$$\displaystyle c_{1}=\frac{Q_{T}}{4\gamma^{2}M_{p}^{2}}\Big{[}2M_{p}^{4}Q_{T}\dot{f}H-2M_{p}^{2}Q_{T}\left(2fM_{p}^{2}\dot{H}+\ddot{f}M_{p}^{2}-4M_{2}^{4}\right)$$

$$c_2=f{M}_p^2,$$

$$c_3=\frac{a{M}_p^2}{\gamma }{Q}_T{Q}_{{\tilde{m}}_4},$$

$$\gamma =H{Q}_T-\frac{m_3^3}{2M_p^2}+\frac{1}{2}\dot{f},Q_{{\tilde{m}}_4}=f+\frac{2{\tilde{m}}_4^2}{M_p^2},$$

$$\left(H\gamma +\frac{{\dot{Q}}_T}{Q_T}\gamma +\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}\gamma -c_T^2\frac{{\gamma }^2}{Q_{{\tilde{m}}_4}}-\dot{\gamma }\right)\frac{Q_T{Q}_{{\tilde{m}}_4}}{{\gamma }^2}>0.$$

$$-\dot{\gamma }>0.$$

$$\left(\frac{{\dot{Q}}_T}{Q_T}\gamma +\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}\gamma -\frac{c_T^2{\gamma }^2}{Q_{{\tilde{m}}_4}}-\dot{\gamma }\right)Q_{{\tilde{m}}_4}>0.$$

$$\frac{{\dot{Q}}_T}{Q_T}\gamma >\dot{\gamma },$$

$$\mathrm{or}{Q}_{{\tilde{m}}_4}<0,\mathrm{or}\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}\gamma >\dot{\gamma }+\frac{c_T^2{\gamma }^2}{Q_{{\tilde{m}}_4}}\left(\mathrm{for}{Q}_{{\tilde{m}}_4}⩾0\right).$$

$$Q_{{\tilde{m}}_4}=\frac{\gamma }{a{M}_p^2}\int a\left(c_1{c}_s^2+c_2\right)dt,$$

$$H\gamma -{\gamma }^2-\dot{\gamma }>0.$$

$$\frac{\dot{\gamma }}{H\gamma }\simeq \frac{\dot{H}}{H^2}\sim (-t{)}^{n-1}\gg 1,$$

$$\left(\frac{{\dot{Q}}_T}{Q_T}\gamma +\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}\gamma -\frac{c_T^2{\gamma }^2}{Q_{{\tilde{m}}_4}}+H\gamma -\dot{\gamma }\right)Q_{{\tilde{m}}_4}>0.$$

$$\frac{{\dot{Q}}_T}{Q_T}>c_T^2\gamma -H+\frac{\dot{\gamma }}{\gamma },$$

$$\mathrm{or}\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}<\frac{c_T^2\gamma }{Q_{{\tilde{m}}_4}}-H+\frac{\dot{\gamma }}{\gamma }\left(\mathrm{initially}{Q}_{{\tilde{m}}_4}<0\right).$$

$${\left({\dot{Q}}_T/Q_T\right)}^{-1}\frac{\dot{\gamma }}{\gamma }=n/p<1.$$

$$\frac{{\dot{Q}}_{{\tilde{m}}_4}}{Q_{{\tilde{m}}_4}}>\frac{c_T^2\gamma }{Q_{{\tilde{m}}_4}}-H+\frac{\dot{\gamma }}{\gamma }$$

$$\displaystyle S=\int d^{4}x\sqrt{-g}\Big{[}{M_{p}^{2}\over 2}R+M_{p}^{2}g_{1}(\phi)X+g_{2}(\phi)X\Box\phi+g_{3}(\phi)X^{2}-M_{p}^{4}V(\phi)$$

$$f=1,$$

$$\Lambda (t)=M_p^{4}V-\frac{1}{2}{g}_2{\dot{\varphi }}^2\left(3H\dot{\varphi }+\ddot{\varphi }\right)+\frac{1}{4}{g}_3{\dot{\varphi }}^4,$$

$$c(t)=\frac{M_p^2}{2}{g}_1{\dot{\varphi }}^2-\frac{1}{2}{g}_2{\dot{\varphi }}^2\left(3H\dot{\varphi }-\ddot{\varphi }\right)+\frac{1}{2}{g}_{2,\varphi }{\dot{\varphi }}^4+\frac{1}{2}{g}_3{\dot{\varphi }}^4,$$

$$M_2^4(t)=-\frac{1}{4}{g}_2{\dot{\varphi }}^2\left(3H\dot{\varphi }+\ddot{\varphi }\right)+\frac{1}{4}{g}_{2,\varphi }{\dot{\varphi }}^4+\frac{1}{2}{g}_3{\dot{\varphi }}^4,$$

$$m_3^3(t)=-g_2{\dot{\varphi }}^3,$$

$$m_4^2=0,$$

$${\tilde{m}}_4^2\ne 0.$$

$$3H^2{M}_p^2=\frac{M_p^2}{2}{g}_1{\dot{\varphi }}^2-3g_{2}H{\dot{\varphi }}^3+\frac{1}{2}{g}_{2,\varphi }{\dot{\varphi }}^4+\frac{3}{4}{g}_3{\dot{\varphi }}^4+M_p^{4}V,$$

$$\dot{H}{M}_p^2=-\frac{M_p^2}{2}{g}_1{\dot{\varphi }}^2+\frac{3}{2}{g}_{2}H{\dot{\varphi }}^3-\frac{1}{2}{g}_2{\dot{\varphi }}^2\ddot{\varphi }-\frac{1}{2}{g}_{2,\varphi }{\dot{\varphi }}^4-\frac{1}{2}{g}_3{\dot{\varphi }}^4,$$

$$+M_p^2{V}_{\varphi },$$

$$e^{2\varphi }=\frac{3f_3}{2M_p^2{f}_1}{\dot{\varphi }}^2.$$

$$\dot{\varphi }=\frac{1}{(-t)},t<0.$$

$$H=\frac{f_3-2f_2}{12M_p^2}\frac{1}{(-t{)}^3}$$

$$a(t)=e^{\int Hdt}=\exp \left(\frac{f_3-2f_2}{24M_p^2{t}^2}\right)\simeq 1+\left(\frac{f_3-2f_2}{24M_p^2{t}^2}\right),$$

Full text

I Introduction

Inflation is still being eulogized for its simplicity and also criticized for its past-incompleteness [1][2]. A complete description of the early universe requires physics other than only implementing inflation.

To the best of current knowledge, the inflation scenario will be past-complete, only if it happens after a nonsingular bounce which is preceded by a contraction [3][4][5], or a slow expansion phase (namely, Genesis phase) with increasing energy density [6][7][8]. These two possibilities will be called bounce-inflation and Genesis-inflation, respectively. Besides being past-complete, a bounce-inflation or Genesis-inflation scenario may explain the probable large-scale anomalies in cosmological microwave background (CMB) [4][9]. The nonsingular Quintom bounce [10][11][12] (see also [13]), the ekpyrotic universe [14][15], the Genesis scenario [16][17][18][19], and the slow expansion scenario [20][21][22][23] have acquired intensive attention. In classical nonsingular (past-complete) cosmologies, the Null Energy Condition (NEC) must be violated for a period.

The ghost-free bounce models [24][25][26][27][28] have been obtained in cubic Galileon [29] and full Horndeski theory [30][31][32]. However, recently, it has been proved in Ref. [33] that there exists a “no-go” theorem for cubic Galileon, i.e., the gradient instability ($c_{s}^{2}<0$ with $c_{s}$ being sound speed of the scalar perturbations) is inevitable in the corresponding models. See also [34] for the extension to the full Horndeski theory, and [35] for the attempts to avoid “no-go” in Horndeski theory. Relevant studies can also be found in [36][37][38].

Recently, in Ref. [39] (see also [40]), we dealt with this issue in the framework of the Effective Field Theory (EFT) [41][42][43][44], which has proved to be a powerful tool. In EFT, the quadratic action of scalar perturbation could always be written in the form (see [39] for detailed derivations)

[TABLE]

where we have neglected higher-order spatial derivatives of the scalar perturbation $\zeta$, the sound speed squared of scalar perturbation

[TABLE]

with the coefficients $c_{1}$, $c_{2}$ and $c_{3}$ being time dependent parameters in general, and $c_{1}>0$ is needed to avoid the ghost instability. The condition for avoiding gradient instability is $c_{s}^{2}\sim{\dot{c}_{3}/a}-c_{2}>0$, which is usually integrated as

[TABLE]

The condition of satisfying the inequality is to have $c_{3}$ cross 0, which is hardly possible in models based on the cubic Galileon [33][34]. However, we found that it can easily be satisfied by applying the EFT operator $R^{(3)}\delta g^{00}$ (with $R^{(3)}$ being the 3-dimensional Ricci scalar on spacelike hypersurface and $\delta g^{00}=g^{00}+1$), so that the gradient instability can be cured.

Though the integral approach (3) is simple and efficient, some details of curing the pathology might actually be missed. In this paper, based on the EFT, using a “non-integral approach”, we revisit the nonsingular cosmologies. We begin straightly with (2), and clarify the origin of pathology and show how to cure it in EFT with new insights into what is happening (Sec. II). To have practice in this clarification, we build a stable model of the Genesis-inflation scenario by using the $R^{(3)}\delta g^{00}$ operator (Sec. III). As a supplementary remark, we discuss a dilemma in the Genesis scenario (Sec. IV).

II Re-proof of the “No-Go” and its avoidance in EFT

The EFT is briefly introduced in Appendix A. In the unitary gauge, the quadratic action of tensor perturbation $\gamma_{ij}$ is (see [39] for the derivation of Eqs. (4) to (8))

[TABLE]

where $Q_{T}=f+{2m_{4}^{2}\over M_{p}^{2}}>0$, $c_{T}^{2}=f/Q_{T}>0$, $f$ and $m_{4}^{2}$ are coefficients defined in the EFT action (63).

The quadratic action of the scalar perturbation $\zeta$ is given by Eq. (1) with

[TABLE]

where $M_{2}^{4}$, $m_{3}^{3}$ and $\tilde{m}_{4}^{2}$ are coefficients defined in the EFT action (63) and they could be time dependent in general.

Only if $c_{1}>0$ and $c_{s}^{2}>0$, the model is free from ghost and gradient instabilities, respectively. In nonsingular cosmological models based on the cubic Galileon [26][27][28], $c_{1}>0$ is not hard to obtain, as can be seen from Eq. (5), since the cubic Galileon contributes the ${m_{3}^{3}(t)\over 2}\delta K\delta g^{00}$ operator in EFT. However, since $c_{3}$ is also affected by ${m_{3}^{3}(t)\over 2}\delta K\delta g^{00}$ through $\gamma$, $c_{s}^{2}<0$ is actually inevitable, as will be demonstrated in the following.

Since $c_{1}>0$, the requirement of $c_{s}^{2}>0$ equals

[TABLE]

Here, $Q_{T}\neq Q_{\tilde{m}_{4}}$ is required, which cannot be embodied by the Horndeski theory [30][31][32]. Thus whether $c_{s}^{2}>0$ or not is controlled by the parameter set ($H$, $\gamma$, $Q_{T}$, $c_{T}^{2}$, $Q_{\tilde{m}_{4}}$).

In the following, with condition (9), we will re-prove the “no-go” theorem for the cubic Galileon, and clarify how to cure it in EFT. Different from the proof in [33][39][40], the re-proof is directly based on the derivative inequality instead of integrating it, which we called “non-integral approach”. We assume that after the beginning of the hot “big bang” or inflation, $\gamma=H>0$, $\dot{\gamma}<0$ and $Q_{\tilde{m}_{4}}=1$.

II.1 Case I: initially $\gamma<0$

Since initially $\gamma<0$, $\gamma$ has to cross 0 from $\gamma<0$ to $\gamma>0$ at $t_{\gamma}$. The analysis below is also applicable for all cases with $\gamma$ crossing 0 from $\gamma<0$ to $\gamma>0$.

In the ekpyrotic and bounce models, initially $\gamma=H<0$. In the Genesis model [16] and slow expansion model [21], $H>0$ during the Genesis, but actually $\gamma=H-{m^{3}_{3}\over 2M_{p}^{2}}<0$, as discussed in Sec. IV. Both belong to the Case I.

In the cubic Galileon case*, $f=Q_{T}=Q_{\tilde{m}_{4}}=1$. Around $t_{\gamma}$, condition (9) is*

[TABLE]

We see that $c_{s}^{2}<0$ is inevitable around $t_{\gamma}$, since ${\dot{\gamma}}>0$. Thus the nonsingular models based on cubic Galileon is pathological, as first proved by LMR in [33].

In the EFT case*, around $t_{\gamma}$, condition (9) requires*

[TABLE]

We might have $c_{s}^{2}>0$, only if (considering only the case where only one of $Q_{T}$ and $Q_{\tilde{m}_{4}}$ is modified while the unmodified one is unity) around $\gamma=0$

[TABLE]

[TABLE]

In solution (12), at $t_{\gamma}$, $\gamma=0$ suggests $Q_{T}=0$. Here, since $\gamma=0$ at $t_{\gamma}$, $c_{1}\sim{Q_{T}/\gamma^{2}}$ diverges. One possibility of removing this divergence is that $\gamma\sim(t-t_{\gamma})^{p}$ and $Q_{T}\sim(t-t_{\gamma})^{n}$ around $t_{\gamma}$, with $n\geqslant 2p$ and $p$, $n$ being constants. In Ijjas and Steinhardt’s model [35], $\gamma\sim t-t_{\gamma}$ while $Q_{T}\sim(t-t_{\gamma})^{2}$, which belongs to this case.

In the bounce model based on the cubic Galileon, Eq. (8) gives $\gamma=H-{m_{3}^{3}\over 2M_{p}^{2}}\neq H$. Generally, the NEC is violated when ${\dot{H}}>0$, while the period of $c_{s}^{2}<0$ corresponds to the phase with $\gamma\simeq 0$ and ${\dot{\gamma}}>0$, these two phases do not necessarily coincide, see Eq. (8). As pointed out by Ijjas and Steinhardt [35], it is the sign’s change of $\gamma$ that causes the pathology. Here, we reconfirmed this point.

In solution (13), if $Q_{\tilde{m}_{4}}>0$, at $t_{\gamma}$, $\gamma=0$ suggests $Q_{\tilde{m}_{4}}=0$; while if $Q_{\tilde{m}_{4}}<0$, since $Q_{\tilde{m}_{4}}=1$ eventually, $Q_{\tilde{m}_{4}}$ must cross 0 at $t_{\tilde{m}_{4}}$ (generally $t_{\tilde{m}_{4}}\neq t_{\gamma}$), at which $\dot{Q}_{\tilde{m}_{4}}\gamma>c_{T}^{2}\gamma^{2}$ must be satisfied. In both cases, ${Q}_{\tilde{m}_{4}}=0$ is required, as proposed by Cai et.al [39] and Creminelli et.al [40].

We see again the details of $Q_{\tilde{m}_{4}}$ crossing 0. In both the Genesis model and the bounce model, initially ${Q}_{\tilde{m}_{4}}=1$, so if $Q_{\tilde{m}_{4}}<0$ around $t_{\gamma}$, ${Q}_{\tilde{m}_{4}}$ must cross 0 twice. Thus it seems that $\dot{Q}_{\tilde{m}_{4}}\gamma>c_{T}^{2}\gamma^{2}$ is hard to implement. However, with (2) and (7), one always could solve ${Q}_{\tilde{m}_{4}}$ for any given $c_{s}^{2}$,

[TABLE]

*where $Q_{T}=1$.

II.2 Case II: $\gamma>0$ throughout

Since $\gamma>0$ throughout, we must have ${\dot{\gamma}}\geqslant 0$ during some period initially555Of course, in Case II, we could also have $\dot{\gamma}<0$ during some period, but what we focus on is the period (i.e., $\dot{\gamma}\geqslant 0$) where pathologies appear., otherwise $\gamma$ will diverge in the infinite past.

In the cubic Galileon case*, condition (9) is*

[TABLE]

In the bounce model, $H<0$ in the contracting phase, and in the Genesis model, $H\sim 0$ in the Genesis phase, both suggest $H\gamma-\gamma^{2}-{\dot{\gamma}}<0$ 666In the case where $\gamma$ grows from 0 initially, (15) is also obeyed no more.. Thus $c_{s}^{2}<0$ is inevitable in the corresponding phases, so the nonsingular models based on the cubic Galileon is pathological.

We see the Genesis model in the cubic Galileon version again in detail. During the slow expansion (Genesis phase), $H\sim{1/(-t)^{n}}$ with the constant $n>1$. Thus

[TABLE]

which implies $H\gamma\ll{\dot{\gamma}}$. Thus with (15), we see that $c_{s}^{2}<0$ is inevitable in the slow expansion phase. It seems that if $n=1$, $H\gamma\ll{\dot{\gamma}}$ might be avoided. However, when $n=1$, we have $H=p/(-t)$ and $a\sim{1/(-t)^{p}}$ with constant $p$, thus $a\rightarrow 0$ in the infinite past. From (16), we see that $c_{s}^{2}>0$ requires $p={H^{2}/{\dot{H}}}>1$. Therefore, the universe is singular, or from another point of view, it is geodesically incomplete since the affine parameter of the graviton geodesics $\int_{t_{i}}^{t_{f}}adt$ is finite for $p>1$ when $t_{i}\rightarrow-\infty$.

In the EFT case*, condition (9) requires*

[TABLE]

*We might have $c_{s}^{2}>0$, only if (considering only the case where either $Q_{T}$ or $Q_{\tilde{m}_{4}}$ is modified) *

[TABLE]

[TABLE]

Generally, $-H{\gamma}+{\dot{\gamma}}>0$, as in Genesis model and bounce model. Thus the solution (18) suggests ${\dot{Q}}_{T}>0$, so that we will have $Q_{T}=0$ in infinite past. Thus based on (12) and (18), it seems that though the pathology can be cured by applying $Q_{T}$, $Q_{T}=0$ is inevitable. A model with (18) has been proposed by Kobayashi **[34]** ($Q_{T}\sim{1\over(-t)^{p}}$, $p>n>1$). During the Genesis $\gamma\sim H\sim{1/(-t)^{n}}$, $n>1$, (17) is

[TABLE]

Initially, $Q_{T}\sim{1\over(-t)^{p}}=0$.

In solution (19), $Q_{\tilde{m}_{4}}$ must cross 0 at $t_{\tilde{m}_{4}}$ to $Q_{\tilde{m}_{4}}>0$, as pointed out by Cai et.al **[39]** and Creminelli et.al **[40]**. Around $t_{\tilde{m}_{4}}$, $\dot{Q}_{\tilde{m}_{4}}>c_{T}^{2}\gamma$ must be satisfied.

In (17), if $Q_{\tilde{m}_{4}}>0$ throughout,

[TABLE]

is obtained. Thus, similar to (18), we have $Q_{\tilde{m}_{4}}=0$ (which definitely requires $\gamma=0$) in the infinite past. In the Genesis model, $Q_{\tilde{m}_{4}}\sim 1/(-t)^{p}$ and $\gamma\sim 1/(-t)^{n}$ with $p>n$, since $\dot{Q}_{\tilde{m}_{4}}/{Q}_{\tilde{m}_{4}}>\dot{\gamma}/\gamma$. However, $p>n$ indicates $\dot{Q}_{\tilde{m}_{4}}<\gamma$ in the infinite past ($Q_{\tilde{m}_{4}}=0$), which violates the inequality (21). Thus $Q_{\tilde{m}_{4}}>0$ throughout seems unworkable.

III Application to Genesis-inflation

In this section, we will build a nonsingular model with the solution (19), in which the slow-roll inflation is preceded by a Genesis phase. A Genesis phase is a slowly expanding phase originating form the Minkowski vacuum with a drastic violation of NEC, i.e., $\epsilon\ll-1$, thus the energy density is increasing with the expansion of the universe and hence is free from the initial singularity **[16]**[17]** (see also **[20]**). As will be shown below, our model cannot only get rid of the pathology of instability, but also give rise to a flat spectrum with interesting features at large scales.

III.1 The setup of the model

The action of the model is

[TABLE]

where $X=-\nabla_{\mu}\phi\nabla^{\mu}\phi/2$, $\Box\phi=\nabla_{\mu}\nabla^{\mu}\phi$, and $\phi$ is a dimensionless scalar field, so dimensionless are $g_{1}(\phi)$, $g_{2}(\phi)$, $g_{3}(\phi)$ and $V(\phi)$.

Mapped into the EFT action (63), (22) corresponds to

[TABLE]

We can get the background equations

[TABLE]

where “${}_{,\phi}=d/d\phi$” and “${}_{,\phi\phi}=d^{2}/d\phi^{2}$”.

Initially, the universe is slowly expanding (in the Genesis phase), $H\simeq 0$. We set $V=0$, $g_{1}=-f_{1}e^{2\phi}$, $g_{2}=f_{2}$ and $g_{3}=f_{3}$, see e.g. Ref. **[7]**, with $f_{1,2,3}$ being dimensionless constants. Thus with Eq. (30), we have $\frac{M_{p}^{2}}{2}g_{1}\dot{\phi}^{2}+\frac{3}{4}g_{3}\dot{\phi}^{4}=0$, which suggests

[TABLE]

*The solution is *

[TABLE]

Eq. (III.1) reads $\dot{H}={f_{3}-2f_{2}\over 4M_{p}^{2}}\dot{\phi}^{4}$. Thus we get

[TABLE]

after the integration. In principle, there could be a constant, i.e., $H={f_{3}-2f_{2}\over 12M_{p}^{2}}{1\over(-t)^{3}}+const$, however, in that case we will have $H\approx const$ initially, which is geodesically incomplete, see also **[45]**.

Additionally, from Eq. (35), we have

[TABLE]

while we set $a(-\infty)=1$.

During inflation, we set $g_{1}=1$ and $g_{2}=g_{3}=0$, since we require that the inflationary phase is controlled by a simple slow-roll field777The behaviors of these $g_{i}$’s in the two phases can easily be matched together by making use of some shape functions [5][46]..

III.2 The primordial perturbation and its spectrum

In the unitary gauge, the quadratic action of the scalar perturbation is presented in the form of Eq. (1). The coefficients $c_{i}$ are (substituting Eqs. (23) to (29) into (5)(6)(7))

[TABLE]

*where *

[TABLE]

The sound speed squared $c_{s}^{2}$ of scalar perturbation is defined in Eq. (2). Here, when $\tilde{m}_{4}^{2}\equiv 0$ or $Q_{\tilde{m}_{4}}=1$, the sound speed squared of the scalar perturbation is reduced to

[TABLE]

It is easy to see that $c_{s0}^{2}=1$ for inflation, since $g_{2}=g_{3}=0$, but not for Genesis. However, using the operator ${\tilde{m}_{4}^{2}(t)\over 2}R^{(3)}\delta g^{00}$, we could always set $c_{s}^{2}=1$ in the Genesis phase, which requires $\tilde{m}_{4}^{2}=-{2M_{p}^{2}(f_{2}+f_{3})\over 4f_{2}+f_{3}}$. This suggests $Q_{\tilde{m}_{4}}=-{3f_{3}\over 4f_{2}+f_{3}}$ is a constant at $|-t|\gg 1$, which is consistent with the solution (19).

The equation of motion of $\zeta$ is

[TABLE]

with $u=z\zeta$, $z=\sqrt{2a^{2}c_{1}}$, the prime denotes the derivative with respect to the conformal time $\tau=\int dt/a$. The initial state is the Minkowski vacuum, thus $u={1\over\sqrt{2c_{s}k}}e^{-ic_{s}k\tau}$ for $\zeta$ modes deep inside the horizon. The power spectrum of $\zeta$ is

[TABLE]

In the following, we will analytically estimate the spectrum of the scalar perturbation. We set $c_{s}^{2}=1$ throughout for simplicity, which could be implemented by using $Q_{\tilde{m}_{4}}(t)$, as will be illustrated by the numerical simulation.

In the Genesis phase, substituting Eqs. (34), (35) into (5), we have

[TABLE]

Thus

[TABLE]

Then it is straightforward to obtain ${z^{\prime\prime}\over z}\approx\frac{\left(f_{3}-2f_{2}\right){}^{2}}{72M_{p}^{4}\tau^{6}}\approx{0\over\tau^{2}}$, where $\tau=\int{1\over a}dt\approx t$. Thus the solution of Eq. (42) is

[TABLE]

*where $C_{11}$ and $C_{12}$ are functions of $k$, $H^{(1)}_{\nu}$ and $H^{(2)}_{\nu}$ are the first and the second kind Hankel functions of *$\nu-$th order, respectively. The initial condition $u={1\over\sqrt{2k}}e^{-ik\tau}$ indicates

[TABLE]

In the inflation phase, $c_{1}=\epsilon M_{p}^{2}$, thus $z=\sqrt{2\epsilon a^{2}M_{p}^{2}}$. We set $\epsilon\ll 1$ as a constant during inflation. Then ${z^{\prime\prime}/z}\approx(2+3\epsilon)/\tau^{2}$. The solution of Eq. (42) is

[TABLE]

with $\nu_{2}\approx 3/2+\epsilon$.

We require that $u_{1}(\tau_{m})=u_{2}(\tau_{m})$ and $u_{1}^{\prime}(\tau_{m})=u_{2}^{\prime}(\tau_{m})$, with $\tau_{m}$ approximately corresponding to the beginning time of inflation phase, and we obtain

[TABLE]

The power spectrum of $\zeta$ is given by

[TABLE]

where $P_{\cal R}^{inf}=\frac{H^{2}_{inf}}{8\pi^{2}M_{p}^{2}\epsilon}\cdot\left({k\over aH}\right)^{3-2\nu_{2}}$ is the power spectrum of scalar perturbation modes that exit horizon during inflation. We see that for the perturbation modes exiting horizon in the Genesis phase, $-k\tau_{m}\ll 1$, $|C_{21}-C_{22}|^{2}\simeq(-k\tau_{m})^{2}$, thus $P_{\cal R}\sim k^{2}$ is strong blue-tilted, while for the perturbation modes exiting horizon in the inflation phase, $-k\tau_{m}\gg 1$, $|C_{21}-C_{22}|^{2}\simeq 1$, thus $P_{\cal R}\sim k^{3-2\nu_{2}}=k^{-2\epsilon}$ is flat with a slightly red tilt.

Tensor perturbation is unaffected by the $R^{(3)}\delta g^{00}$ operator. Its quadratic action is given in Eq. (4) with $Q_{T}=1$ and $c_{T}^{2}=1$. The spectrum of primordial GWs can be calculated similarly, see also Ref. **[6]**. Since $z_{T}^{\prime\prime}/z_{T}=a^{\prime\prime}/a$, we have

[TABLE]

where $P_{T}^{inf}=\frac{2H^{2}_{inf}}{\pi^{2}M_{p}^{2}}\cdot\left({k\over aH}\right)^{3-2\nu_{2}}$ is the power spectrum of tensor perturbation modes that exit horizon during inflation. Thus the spectrum of primordial GWs has a shape similar to that of the scalar perturbation.

III.3 Numerical simulation

In the numerical calculation, we set

[TABLE]

with $f_{1,2,3}$, $q_{1,2,3,4}$, $\phi_{0,1,2}$ and $\lambda$ being dimensionless constants. When $\phi\ll\phi_{0}$, we have $g_{1}=-f_{1}e^{2\phi}$, $g_{2}=f_{2}$ and $g_{3}=f_{3}$, which brings a Genesis phase (36), while $\phi\gg\phi_{0}$, we have $g_{1}=1$ and $g_{2}=g_{3}=0$, the slow-roll inflation will occur with $V(\phi)\sim\phi^{2}$. When $\phi\ll\phi_{2}$, $V(\phi)\approx 0$, while $\phi\gg\phi_{2}$, $V(\phi)\approx{\lambda\over 2}(\phi-\phi_{1})^{2}$. We do not require $\phi_{0}=\phi_{2}$ but $\phi_{0}>\phi_{2}$.

We start the simulation at $t_{i}\ll-1$, and we set

[TABLE]

and

[TABLE]

We show the evolution of $\phi$ and $\dot{\phi}$ in Fig. 1, and the evolution of $a$, $H$ and $\epsilon$ in Fig. 2. In Fig. 3(a), $c_{1}$ is plotted, and $c_{1}>0$ is satisfied. In Fig. 3(b), we see that $\gamma$ does not cross 0, which implies that, in the Genesis phase, $c_{s0}^{2}<0$ (see Fig. 4(a)), as proved in Sec. IV. By including the operator $R^{(3)}\delta g^{00}$, we could have $c_{s}^{2}>0$, and so cure the gradient instability. The spectrum of the scalar perturbation can be simulated numerically, which is plotted in Fig. 5. The spectrum obtained has a cutoff at large scale $k<k_{*}$ and is nearly scale-invariant for $k>k_{*}$, as displayed in Eq. (51).

IV The dilemma of $\gamma$ in the Genesis scenario

In the Genesis scenario based on the cubic Galileon **[16]**, (also **[21]**), we have

[TABLE]

during the Genesis, where $f_{2}<0$.

In Ref. **[16]**, $f_{3}=-f_{2}$, which suggests $\gamma={f_{2}\over 4M_{p}^{2}}\dot{\phi}^{3}<0$. Thus if a hot “big bang” or inflation ($\gamma=H>0$) starts after the Genesis phase, $\gamma$ must cross 0 at $t_{\gamma}$ ($c_{s0}^{2}<0$* around $t_{\gamma}$, which may be cured by applying $Q_{\tilde{m}_{4}}$). It is obvious that when $\gamma=0$, $c_{1}$ in (1) will be divergent. Though this divergence might not be a problem, it will affect the numerical simulation for perturbations [47][48]**, unless $Q_{T}/\gamma^{2}$ is finite at $t_{\gamma}$, as in Ijjas and Steinhardt’s model [35].*

In the model of **[7]**, the Genesis is followed by Galileon inflation **[49]**. Though $f_{3}=-f_{2}$ and $\gamma={f_{2}\over 4M_{p}^{2}}\dot{\phi}^{3}<0$ in the Genesis phase, one might also have $\gamma<0$ for Galileon inflation, since $g_{2}\neq 0$ in (22) during inflation. Thus it seems that $\gamma$ might not necessarily cross 0. However, after inflation, $\gamma$ crossing 0 is still inevitable.

In our model, the Genesis is followed by the slow-roll inflation, $\gamma=H>0$ for inflation. To not cross 0, initially $\gamma$ must satisfy $\gamma>0$. In the Genesis phase, this suggests $f_{3}>-4f_{2}$. Thus we will have $\gamma>0$ throughout. However, for the cubic Galileon, the expense is

[TABLE]

during the Genesis. Here, this pathology is cured in EFT by applying (19).

V Conclusion

Based on the EFT of cosmological perturbations, we revisit the nonsingular cosmologies, using the “non-integral approach”. By doing this, we could have a clearer understanding of the pathology in nonsingular Galileon models and its cure in EFT.

We clarify the application of the operator ${\tilde{m}_{4}^{2}}R^{(3)}\delta g^{00}/2$ in EFT, which is significant for curing the gradient instability. We show that if $Q_{\tilde{m}_{4}}<0$ around $\gamma=0$ is adopted to cure the gradient instability, in solution (13) (with $\gamma<0$ and ${Q}_{\tilde{m}_{4}}=1$ initially), ${Q}_{\tilde{m}_{4}}$ must cross 0 twice; while in solution (19) (with $\gamma>0$ throughout), initially $Q_{\tilde{m}_{4}}<0$ must be satisfied, $Q_{\tilde{m}_{4}}$ will cross 0 to $Q_{\tilde{m}_{4}}>0$ at $t_{\tilde{m}_{4}}$, and crosses 0 only once. Thus at a certain time, ${Q}_{\tilde{m}_{4}}$ meeting 0 is required, as pointed out first by Cai et.al **[39]**, and also by Creminelli et.al **[40]**.

We also clarify that in the bounce model with $\gamma<0$ initially, $c_{s}^{2}<0$ will occur in the phase with $\gamma\simeq 0$ and ${\dot{\gamma}}>0$, while the NEC is violated when ${\dot{H}}>0$ (bounce phase), these two phases do not necessarily coincide. As pointed out by Ijjas and Steinhardt **[35]**, it is the sign’s change of $\gamma$ that causes $c_{s}^{2}<0$. Here, we verify this point. In Genesis model **[16]**[7]**, and also **[21]**, the case is similar, as discussed in Sec. IV.

The nonsingular model with the solution (19) ($\gamma>0$* throughout) has not been studied before. In Sec. III, we design such a model, in which a slow expansion phase (namely, the Genesis phase) is followed by slow-roll inflation. Under the unitary gauge, since ${\dot{\gamma}}>0$ and $\gamma>0$ (not crossing 0), the evolution of primordial perturbations can be simulated numerically. The simulation displays that the spectrum acquires a large-scale cutoff, as expected in Ref. [6].*

We conclude that, based on EFT, not only a stable nonsingular cosmological scenario may be built without getting involved in unknown physics, but also the phenomenological possibilities of its implementation are far richer than expected (see also recent **[50]**[51]** for the higher spatial derivative operators).

Acknowledgments**

We thank Youping Wan for helpful comments. The work of YSP is supported by NSFC, Nos. 11575188, 11690021, and also supported by the Strategic Priority Research Program of CAS, Nos. XDA04075000, XDB23010100. The work of T. Q. is supported by NSFC under Grant Nos. 11405069 and 11653002.

Appendix A EFT of cosmological perturbations

With the ADM line element, we have

[TABLE]

and $\sqrt{-g}=N\sqrt{h}$, where $N_{i}=h_{ij}N^{j}$. The unit one-form tangent vector is defined as $n_{\nu}=n_{0}(dt/dx^{\mu})=(-N,0,0,0)$ and $n^{\nu}=g^{\mu\nu}n_{\mu}=({1/N},-{N^{i}/N})$, which satisfies $n_{\mu}n^{\mu}=-1$. On the hypersurface, the induced 3-dimensional metric is $H_{\mu\nu}=g_{\mu\nu}+n_{\mu}n_{\nu}$, thus

[TABLE]

The extrinsic curvature is $K_{\mu\nu}\equiv{1\over 2}{\cal L}_{n}H_{\mu\nu}$, where ${\cal L}_{n}$ is the Lie derivative with respective to $n^{\mu}$. The induced 3-dimensional Ricci scalar $R^{(3)}$ associated with $H_{\mu\nu}$ is

[TABLE]

Without higher-order spatial derivatives, the EFT reads **[39]**

[TABLE]

where $\delta g^{00}=g^{00}+1$, $\delta K_{\mu\nu}=K_{\mu\nu}-H_{\mu\nu}H$ with $H$ being the Hubble parameter. The coefficient set $(f,c,\Lambda,M_{2},m_{3}$, $m_{4},{\tilde{m}}_{4})$ specifies different theories and could be time-dependent in general 888Different conventions of the nomenclatures of these coefficients were adopted during the development of the EFT of cosmological perturbations (see e.g., [41][42][43][44]). Here, we follow the convention used in Refs. [43][44].. A particular subset $(m_{4}={\tilde{m}}_{4})$ of EFT (63) is the Horndeski theory. $S_{m}[g_{\mu\nu},\psi_{m}]$ is the matter part, which is minimally coupled to the metric $g_{\mu\nu}$.

To obtain the quadratic actions for scalar and tensor perturbations, we will work in the unitary gauge, thus we set

[TABLE]

Then we follow the standard method first used by Maldacena **[52]**, it is straightforward (though tedious) to obtain the quadratic actions of scalar perturbation $\zeta$ and tensor perturbation $\gamma_{ij}$, as exhibited in Eqs. (1) and (4), respectively (see **[39]** for detailed derivations).