pre-inflationary non gaussianities


2025-10-04

On Pre–Inflationary non Gaussianities
M. Meo and A. Sagnotti
\({}^c\) Scuola Normale Superiore and INFN
Piazza dei Cavalieri, 7
56126 Pisa ITALY
e-mail:
mario.meo@sns.it, sagnotti@sns.it**

Abstract

We explore the three–point amplitude of curvature perturbations in scenarios suggested by high–scale supersymmetry breaking in String Theory, where the inflaton is forced to climb a steep exponential potential. We can do it at the price of some simplifications, and more importantly with some assumptions on the softening effects of String Theory. These suggest a scenario proposed long ago by Gasperini and Veneziano, where the initial singularity is replaced by a bounce, and the resulting analysis rests on a scale \(\Delta\) that leaves some signs in the angular power spectrum of the CMB. The amplitude comprises two types of contribution: the first oscillates around the original result of Maldacena and gives no further prospects to detect a non–Gaussian signal, but the second, which is subtly tied to the turning point at the end of the climbing phase, within the window \(62<N<66\) for the inflationary \(e\)-folds could be compatible with Planck constraints and potentially observable. The amplitudes involving the tensor modes contain only the first type of contribution.

1 Introduction↩︎

The three ten–dimensional non–supersymmetric strings free of tachyons [1][3] are a good laboratory to explore the effects of supersymmetry breaking at high scales in String Theory [4]. They all share a common feature, the emergence of an exponential “tadpole potential” \[V(\varphi)\;= \;\frac{T}{2 \kappa^2}\;e^{2 \,\gamma\,\varphi} \label{Vgamma95int}\tag{1}\] that forbids a flat vacuum, where \(\varphi\) denotes a conveniently rescaled dilaton. This occurs at the torus level for the heterotic model of [1], for which \(\gamma=\frac{5}{3}\), and at the (projective-)disk level for the two orientifold models [5] of [2], [3], for which \(\gamma=1\). A striking change of behavior [6] occurs in the resulting spatially flat cosmologies 1 as the parameter \(\gamma\) in eq. 1 is increased and reaches the value corresponding to the orientifold models, which thus concerns all three cases. The scalar field is then forced to emerge from the initial singularity while climbing up the exponential potential, so that the resulting dynamics can be confined to the weak–coupling regime of String Theory and includes a turning point. If the tadpole potential is supplemented by a milder contribution, for example the Starobinsky potential of [9], this setup can induce the onset of inflationary slow–roll [10].

A pre–inflationary period of fast-roll introduces a low–frequency cut in the primordial power spectrum of curvature perturbations. If this region were accessible to us, the cut would depress the first few multipoles in the CMB angular power spectrum [11], which resonates with the lack of power present in them. The cut can be modeled, as in [11], by deforming the Chibisov–Mukhanov primordial power spectrum [12] into 2 \[P(k) \;= \;A \;\frac{\left(\frac{k}{k^\star}\right)^3}{\left[ \left(\frac{k}{k^\star}\right)^2 \;+ \;\left(\frac{\Delta}{k^\star}\right)^2 \right]^{2 - \frac{n_s}{2}}} \;, \label{power95delta95int}\tag{2}\] where \(k^\star\) is a pivot scale, at the price of introducing the new scale \(\Delta\) into the problem. Some evidence for the scale \(\Delta\) was actually detected in CMB data, and reaches the 3\(\sigma\) level if the masked region around the Galactic plane is widened, thus working within a sky fraction of about 39% [14]. The resulting value is of the order of the cosmic Horizon and translates into an energy scale of about \(10^{12} \,e^{N-60}\) GeV at the beginning of an inflationary period of \(N\) \(e\)-folds. However, other features of the power spectrum that are typical of the climbing scenario appear beyond reach. For example, the detection in CMB data of a pre-inflationary peak reflecting the presence of a turning point was attempted in [14] using a simplified model, but did not yield significant results.

The three–point amplitude contains more information on the early stages of the Universe, since it rests on the Green functions before and during the inflationary phase. From a theoretical point of view, the analysis of new cases has some interest of its own, since it can perhaps add to the many lessons that are being drawn from cosmological amplitudes [15], and the setting we are exploring does contain an interesting novelty. This is the presence of a turning point for the scalar where the slow–roll parameter \(\epsilon\) vanishes, making the three–point amplitude based on the Mukhanov–Sasaki variable apparently singular. Or, if you will, making curvature and scalar perturbations somehow inequivalent. If one insists on working with curvature perturbations, which are closer in spirit to the actual observations, the singularity is resolved by the Schwinger–Keldysh contour, and finite contributions emerge from the neighborhood of the turning point, with peculiar and potentially interesting features.

The initial singularity, which is also present in the cosmologies under scrutiny, is another story. All we shall be able to do concretely is appeal to the softening effects of String Theory, arguing for its eventual disappearance in a full treatment. Still, the way to do this will not be completely arbitrary, and we shall be drawn to a scenario that resonates with the “Pre–Big Bang Cosmology" proposed long ago, in String Theory, by Gasperini and Veneziano [16].

With this important proviso, this paper tries to add something on the important issue of non–Gaussianities, but for brevity, we only discuss in detail the equilateral configuration, which yields the strongest signal. As we shall see, the scale \(\Delta\) introduces oscillations in the dimensionless observable \(f_{NL}\) originally computed by Maldacena in [17] (for reviews, see [18]). Moreover, the contribution from the turning point identifies a narrow window in the number of inflationary \(e\)-folds around \(N=63\), where the complete result is sizable and yet lies within limits that emerged from the Planck collaboration [19], [20]. However, as we have anticipated, all this comes at the price of appealing to the role of String Theory to regulate the pre–inflationary dynamics, in a way that we shall try to motivate and appears somewhat natural to us, but cannot be fully justified.

2 Exponential Potentials and Climbing↩︎

Consider the action principle for Einstein gravity minimally coupled to a real scalar field with a potential \(V\), \[{ \cal S} \;= \;\int d^4 x \sqrt{-g}\left[\frac{R}{2\,\kappa_D^2}\;- \;\frac{1}{2}\,g^{\mu \nu }\partial_{\mu}\phi\,\partial_{\nu}\phi\;-\;V(\phi)\right] \;. \label{2der95act}\tag{3}\] While the climbing mechanism emerged in ten–dimensional strings, here we shall work directly with its possible manifestation in four dimensions, referring mostly to the two ten–dimensional orientifold models where the phenomenon sets it. This step affords a partial justification since, if the internal volume is somehow stabilized, a ten–dimensional tadpole potential that is “critical” in the sense that we shall briefly recall translates into a four–dimensional one that is also critical, as explained in [21].

In four–dimensional spatially flat cosmologies of the form \[ds^2 \;=\; - \;e^{2B(\xi)}d\xi^2 \;+ \;e^{2A(\xi)}d\vec{x} \cdot d \vec{x}\;,\] which encode the two key properties of homogeneity and isotropy, with potentials that never vanish one can make the gauge choice [7] \[V(\phi)\,e^{2B}\;= \;\frac{M^2}{2\,\kappa^2} \;. \label{gauge95choice}\tag{4}\] Combining it with the redefinitions \[\tau \;= \;\xi \,M \sqrt{\frac{3}{2}} \;, \qquad \varphi \;= \; \kappa \,\phi\, \sqrt{\frac{3}{2}} \;, \qquad A \;= \;\frac{1}{3}\,a\] leads to a neat form for the background equations in an expanding Universe, \[\begin{align} && \frac{da}{d\tau} \;=\;\sqrt{1\;+\;\left(\frac{d\varphi}{d\tau}\right)^2} \;, \nonumber \\ &&\frac{d^2\varphi}{d\tau^2}\;+\;\frac{d\varphi}{d\tau}\sqrt{1\;+\;\left(\frac{d\varphi}{d\tau}\right)^2}\;+\;\frac{1}{2V}\frac{\partial V}{\partial \varphi}\left[1\;+\;\left(\frac{d\varphi}{d\tau}\right)^2\right] \;= \;0 \;, \label{eqs95gen} \end{align}\tag{5}\] where the driving force is due to the logarithmic derivative of the potential \(V\).

A peculiar behavior emerges when these cosmologies are driven by exponential potentials, which can be parametrized as \[V(\varphi)\;= \;\frac{T}{2 \kappa^2}\;e^{2 \,\gamma\,\varphi} \;. \label{Vgamma}\tag{6}\] In the orientifold models of [2], [3], which will be our main focus in this paper, these potentials would reflect an overall positive brane–orientifold tension \(T\), which can be conveniently identified with the scale \(M^2\) in eq. 4 , while \(\gamma \geq 0\) characterizes the driving force on the scalar field due to \(V(\phi)\) 3. Note that in four dimensions, the exponential potential takes the form \[\frac{T}{2 \kappa^2}\;e^{\kappa \, \gamma\, \sqrt{6} \, \phi} \;,\] when expressed in terms of the canonically normalized field \(\phi\).

The “parametric time” \(\tau\) is related to the cosmic time according to \[dt_c \;= \;\sqrt{\frac{2}{3}} \;e^{-\,\gamma\,\varphi} \;\frac{d\tau}{M} \;, \label{tctaugamma}\tag{7}\] and to the conformal time according to \[d\eta \;= \;\sqrt{\frac{2}{3}} \;e^{-\, \frac{a}{3} \,-\, \gamma\,\varphi}\; \;\frac{d\tau}{M} \;. \label{etataugamma}\tag{8}\] With an exponential potential, the second of eqs. 5 reduces to \[\frac{d^2\varphi}{d \tau^2}\;+\;\frac{d\varphi}{d \tau}\sqrt{1\;+ \;\left(\frac{d\varphi}{d \tau}\right)^2}\;+ \;\gamma\left[1\;+ \;\left(\frac{d\varphi}{d \tau}\right)^2\right] \;= \;0\;,\] and for \(\gamma<1\) there are two classes of expanding solutions, with quite different scalar dynamics.

  • In “descending” cosmologies the scalar field \(\varphi\) emerges from the initial singularity from large positive values, and then decreases (one can thus say that \(\varphi\) “descends” the potential) during the subsequent cosmological evolution.

  • In “climbing” cosmologies the scalar field emerges from large negative values of \(\varphi\), and “climbs up” the potential until it reaches a turning point before it starts to descend. In String Theory, this solution has the virtue of involving a bounded string coupling, since \(\varphi\) is naturally related to the dilaton [21].

Both solutions eventually approach the limiting behavior described by \[\varphi(\tau)\;=\;\varphi_0 \;- \;\frac{\gamma \,\tau}{\sqrt{1-\gamma^2}}\;, \qquad a(\tau)\;=\;\frac{\tau}{\sqrt{1-\gamma^2}} \;,\] which is the well–known Lucchin–Matarrese attractor [22] in the present gauge. In conformal coordinates, this limiting behavior is described by \[ds^2 \;= \; \left[\frac{\sqrt{6(1-\gamma^2)}}{M(1-3\gamma^2)} \;\left(-\eta\right) \;e^{-\gamma \varphi_0}\right]^{\frac{2}{1-3\gamma^2}}\Big(-\;d\eta^2\;+ \;d \vec{x} \cdot d\vec{x}\Big)\;,\] while \[\varphi(\eta) \;= \; \frac{\varphi_0}{1-3\gamma^2}\;+ \;\frac{3\gamma}{1-3\gamma^2}\log \left[\frac{M(1-3\gamma^2)}{\sqrt{6(1-\gamma^2)}}(-\eta)\right] \;.\]

Note that for \(\gamma \ll 1\) \[\varphi(\eta) \;\simeq \;\varphi_0 \;+ \;3 \gamma \,\log\left[\frac{M}{\sqrt{6}} \left(- \eta\right) \right]\;, \qquad e^{\frac{1}{3}\,a(\eta)} \;\simeq \;\sqrt{6} \;e^{-\gamma \varphi_0} \frac{1}{M \left(- \eta\right)} \;, \label{aphidS}\tag{9}\] where we have retained contributions up to \({\cal O}(\gamma)\), so that when the system approaches the Lucchin–Matarrese attractor one is describing a de Sitter phase with \[{\cal H} \;=\;\frac{M}{\sqrt{6}} \;e^{\gamma \varphi_0} \label{HM} \;.\tag{10}\] For both classes of solutions the preceding constant value persists until \(\tau={\cal O}\left(\frac{1}{\gamma^2}\right)\), and then \({\cal H}(\tau)\) starts to decrease.

The preceding solutions give rise to power–law inflation for \(\gamma< \frac{1}{\sqrt{3}} \simeq 0.58\), but the situation changes drastically as \(\gamma \to 1\). Both the descending solution and the Lucchin–Matarrese attractor then disappear, while the leftover climbing solution takes the simple form \[\varphi(\tau)\;= \;\widehat{\varphi}_0 \;+ \;\frac{1}{2}\left(\log \tau \;- \;\frac{\tau^2}{2}\right) \;, \qquad a(\tau)\;= \;\frac{1}{2}\left(\log \tau \;+ \;\frac{\tau^2}{2}\right)\;. \label{critical95beh}\tag{11}\]

3 Climbing, \(\Delta\)and the Primordial Power Spectrum↩︎

For \(\gamma \geq 1\) only climbing solutions exist, but here we shall confine our attention to a hard exponential term with \(\gamma=1\), which is directly related to the non–supersymmetric ten–dimensional orientifolds. This case will be central in our discussion, while the behavior for \(\gamma < 1\) will be just indicative of what happens in slow-roll phases.

The climbing behavior is forced by the “hard exponential”, and so it also occurs in the presence of more general potentials of the form \[V \;= \;\frac{T}{2 \kappa^2}\;e^{2 \,\varphi} \;+ \;v(\varphi) \;, \label{vmild}\tag{12}\] with \(v(\varphi)\) a milder non–singular contribution, and can provide a reason for the onset of inflation. This is true, in particular, for the “double–exponential potentials” \[V \;= \;\frac{T}{2 \kappa^2}\;e^{2 \,\varphi} \;+ \;v_0 \,e^{2 \,\gamma\,\varphi} \;, \label{2exp}\tag{13}\] with \(v_0>0\) and \(\gamma< \frac{1}{\sqrt{3}}\). Despite its simplicity, this setting is not exactly solvable, but the results recalled in the previous section capture its main features, since one can work with \[V \;= \;v_0 \,e^{2 \,\gamma\,\varphi} \label{2exp1}\tag{14}\] during most of the evolution, where the “hard” exponential has sub–dominant effects, and with \[V \;= \;\frac{T}{2 \kappa^2}\;e^{2 \,\varphi} \label{2exp2}\tag{15}\] near the turning point, where the “hard exponential” dominates. Integrable variants of the potential 13 exist, and in particular the simple exact solution in the presence of \[V \;= \;\frac{T}{2 \kappa^2}\left[ e^{\frac{2}{\gamma} \,\varphi} \;+ \;e^{{2}{\gamma} \,\varphi}\right] \;. \label{integrable}\tag{16}\] displays very neatly the transition from an early climbing phase to the onset of inflation [21].

a

Figure 1: Typically the scale factor (red, dashed) approaches the attractor (black, dotted) while the scalar field is still in the climbing phase with medium–large values of \(\epsilon\) (green, solid), but an initial singularity is readily encountered at earlier times in the region where \(\epsilon \simeq 3\). The dip in the solid curve for \(\epsilon\) corresponds to the turning point..

For the typical values \(\gamma={\cal O}(0.1)\) that are selected by rough comparisons with CMB data, a de Sitter–like expansion takes place during the final descent. The results collected in [11] actually indicate that the same attractor behavior, with an essentially constant \({\cal H}\), also characterizes the final part of the earlier ascent. We shall thus assign to the scale factor of these systems a de Sitter–like behavior after the turning point, and in particular the de Sitter relation between cosmic and conformal time, \[dt_c \;= \;- \;\frac{d\eta}{{\cal H} \,\eta} \;,\] will be used repeatedly in the following. Details of the scalar dynamics play a key role near the turning point, where demanding an essentially constant value \({\cal H}\) for the Hubble parameter links \({\cal H}(\tau)\) to the microscopic parameters of the preceding section for the “hard” exponential. In that region \(H\) evolves according to \[{H}(\tau)\;= \;\frac{{M}_0}{2\,\sqrt{6}}\left(\frac{1}{\tau}+\tau\right)\tau^{1/2}\;e^{-\tau^2/4}\;e^{\varphi_0} \;,\] in terms of the parametric time \(\tau\), where \(M_0\) is the scale associated to the “hard” exponential, and identifying with \({\cal H}\) the value attained by \(H(\tau)\) at \(\tau=1\), which identifies the turning point, leads to \[{\cal H} \;= \;\frac{M_0}{\sqrt{6}}\;e^{-\frac{1}{4}}e^{{\varphi}_0} \;. \label{HM0}\tag{17}\]

a

Figure 2: The non–singular model for \(\epsilon(\eta)\) (blue, dot-dashed), with a dip at \(\eta_t\), where the climbing scalar inverts its motion. As \(\eta\) decreases, the actual \(\epsilon(\eta)\) (orange, solid) experiences a quick growth to its asymptotic value as the Universe collapses, and the dashed vertical line corresponds to the initial singularity..

A slow–roll phase is usually characterized by a set of slow–roll parameters, and primarily by \[\epsilon \;= \;\frac{3\,\left(\frac{d\varphi}{d \tau}\right)^2}{1 \;+\; \left(\frac{d\varphi}{d \tau}\right)^2} \;, \qquad \delta \;= \;3\,\frac{\frac{d^2\varphi}{d t_c^2}}{\frac{d\varphi}{d t_c} \, \frac{d a}{d t_c}} \;. \label{delta95epsilon}\tag{18}\] As can be clearly seen in Fig. 1, tracing the scalar dynamics backwards one encounters a zero of \(\epsilon\) at \(\eta_t\), which corresponds to the turning point, and a sharp increase at earlier times, which is accompanied by the sudden encounter with the initial singularity. To the best of our knowledge, there is no satisfactory way to deal with the divergent contribution that the singularity introduces in three-point amplitudes, but the higher–derivative corrections present in String Theory are expected to resolve it. We shall therefore assume, with good motivations but with an evident degree of arbitrariness, that the singularity would be resolved in a full string treatment, and that the growth of \(\epsilon\) stops as it quickly recovers a slow–roll value, say \(3 \gamma^2\) again, as sketched in fig. 2. As we shall see, if this picture is somehow correct, the three–point amplitude places intriguing constraints on the earlier dynamics.

The climbing behavior typically introduces a low–frequency cut in the primordial power spectrum of curvature perturbations, which resonates with the apparent lack of power present in the first few CMB multipoles. This feature is generally accompanied by a peak that reflects the end of the climbing phase, whose position and size are, however, non–universal features of the scalar dynamics. Leaving aside these finer details, a simple formula, \[P(k) \;= \;A \;\frac{\left(\frac{k}{k^\star}\right)^3}{\left[ \left(\frac{k}{k^\star}\right)^2 \;+ \;\left(\frac{\Delta}{k^\star}\right)^2 \right]^{2 - \frac{n_s}{2}}} \label{power95delta} \;,\tag{19}\] depending on a single additional parameter, the scale \(\Delta\), aside from \(A\) and the reference scale \(k^\star\), can capture the gross features of the low–frequency cut while also approaching the Chibisov–Mukhanov tilt \[P(k) \;\sim \;A \;\left(\frac{k}{k^\star}\right)^{n_s-1} \label{chib-muk}\tag{20}\] for values of \(k\) above \(\Delta\). The modified power spectrum in eq. 19 is actually determined by an exact solution of the Mukhanov–Sasaki equation \[v''_k(\eta)\;+\;\left(k^2\;+ \;\Delta^2 \;- \;\frac{\nu^2\;- \;\frac{1}{4}}{\eta^2}\right)v_k(\eta)\;= \;0 \;,\] with the potential \[W \;= \;\frac{\nu^2\;- \;\frac{1}{4}}{\eta^2} \;- \;\Delta^2 \;, \label{W95Delta}\tag{21}\] which obtains lowering by \(\Delta^2\) the standard attractor \(W\).

Some evidence for the scale \(\Delta\) was found in [14] in the CMB data. The detection level improves up to the 3\(\sigma\) level as the mask around the Galactic plane is enlarged, and the sharpest result, \[\Delta \;= \;\left (0.351 \pm 0.114\right) \times 10^{- 3}\; \mathrm{Mpc}^{- 1} \;, \label{delta95value}\tag{22}\] is obtained within an open sky fraction \(f_{sky}\simeq 39\%\). The preceding value is comparable to the Cosmic Horizon, and retracing the past history of the Universe one can translate it into an energy scale at the onset of inflation [14] \[\Delta_\mathrm{inf} \;\simeq \;3 \times 10^{14} \;e^{N-60} \;\sqrt{\frac{H_\mathrm{inf}}{\mu_{\mathrm{Pl}}}} \;\mathrm{GeV} \;\simeq \;2 \times 10^{12} \;\;e^{N-60} \;\mathrm{GeV} \;, \label{DeltaH}\tag{23}\] where \(H_\mathrm{inf} \simeq 10^{14}\;\mathrm{GeV}\), \(\mu_{\mathrm{Pl}} \simeq 2.4 \times 10^{18}\;\mathrm{GeV}\) is the reduced Planck energy and \(N\) denotes the number of inflationary \(e\)-folds.

4 \(\Delta\) and the Three–Point Amplitude↩︎

We can now link \(\Delta\) to the conformal time of the turning point that concludes the climbing phase. To this end, let us start by recalling that the combination \[z \;= \; e^\frac{a}{3}\;\sqrt{2 \,\epsilon}\] satisfies \[\frac{z''}{z} \;= \;W \;,\] where \(W\) is the Mukhanov–Sasaki potential, whose role was recalled in Section 3. During the descent, the scale factor has settled on the attractor value \[e^{\frac{1}{3}\,a_{LM}}\;\sim \;(-\eta)^{-\frac{1}{1-3\gamma^2}} \;,\] and, as we have stressed, this is typically the case already before the climbing phase ends. Under this assumption, one can deduce the behavior of \[\sigma \;= \;\sqrt{\epsilon}\] from the linear differential equation \[\sigma''\;+ \;\left(1\,-\,2\nu\right)\frac{\sigma'}{\eta} \;+ \;\Delta^2\,\sigma\;= \;0 \;, \label{sigmadelta}\tag{24}\] where \[\nu \;= \;\frac{3}{2} \frac{1 \;- \;\gamma^2}{1 \;- \;3\,\gamma^2} \;.\]

Demanding that \(\epsilon\) approach \(3\gamma^2\) as \(\eta \rightarrow 0^-\), one can conclude that
\[\sigma(\eta)\; \simeq \;c_1 (-\Delta \eta)^{\nu} J_{\nu}(-\Delta \eta) \;- \; \frac{\pi}{2^{\nu}\Gamma(\nu)}\sqrt{\frac{2\nu-3}{2\nu-1}}(-\Delta \eta)^{\nu} Y_{\nu}(-\Delta \eta)\;. \label{sigma95delta}\tag{25}\] For negative values of \(c_1\), the first zero \(\eta_t\) of \(\sigma\) approaches zero, while for large positive values it approaches \(-\,\frac{4.5}{\Delta}\). We shall leave \(\eta_t\) as a free parameter in the range \(\left[-\,\frac{4}{\Delta},-\,\frac{1}{\Delta}\right]\), where the lower values correspond to stronger bounces of the scalar field against the “critical” exponential wall, while the higher values correspond to milder ones.

Figure 3: Left panel: a typical Mukhanov–Sasaki potential for a climbing scalar. Note that W changes sign at \eta^\star and experiences a short growth as \eta decreases further, before tending to -\infty as the initial singularity is approached. Right-panel: an attractor W with \nu=\frac{3}{2} lowered by \Delta^2, as in eq. 21 .

The evolution can be traced backward until the first zero of \(\sigma\) or \(\epsilon\), which can be identified with the turning point at \(\eta_t\) and links it to the scale \(\Delta\). This result will play an important role in our construction of the three–point amplitude, but now the issue is what happened earlier. Fig. 3 compares the typical behavior of \(W\) for a climbing–scalar cosmology with an attractor \(W\) lowered by \(\Delta^2\) as in eq. 21 .

The simplified shape in the right panel of fig. 3 suffices to model the low-\(k\) behavior of the primordial power spectrum and, as we have seen, can also recover the presence of a turning point. However, the three–point amplitude is built from contributions of the type \[\int_{-\infty}^0 d\eta \,f(\eta) \Big[G_1(k_1,\eta)\, G_2(k_2,\eta)\, G_3(k_3,\eta) \;- \;\mathrm{c.c.} \Big] \;, \label{3pointcc}\tag{26}\] where the \(G\)’s denote Green functions or their first time derivatives, and depends crucially on the earlier history of the system, so more information is needed to compute it. One could try to continue with the actual \(W\) determined by the two–derivative Lagrangian. However, the scalar moves faster and faster for earlier times, so that \(\epsilon\) quickly approaches three, its upper bound, and one encounters an initial singularity at a nearby finite value \(\eta_s\) of the conformal time, where the Universe collapses. This fate should be avoided in a complete theory of gravity, and, most importantly for us, the three–point amplitude suffers there from what one could call an ultraviolet divergence, since some integrands are dominated by \[\frac{\log(\eta\,-\,\eta_s)}{\eta \,-\, \eta_s} ,\] as \(\eta \to \eta_s^+\). This is a problem if one tries to enforce the Bunch–Davis condition there, where the evolution determined by eq 3 starts.

Figure 4: Left panel: an extended model for W that follows the curve of eq. 21 up to \eta_0, where it experiences a sudden jump to the corresponding attractor curve without \Delta, has the virtue of leading to simple mode functions. Right panel: the actual transition could occur smoothly, within a small window of conformal time.

Trying to avoid this problem leads to an interesting scenario, as we can now explain. String Theory is supposed to regulate the divergence, albeit in ways that we do not control at present, and therefore some guesswork is needed. One could conceive, for example, that the higher–derivative corrections drive \(\epsilon\) to recover its slow–roll value \(3 \gamma^2\), as in fig. 2, while the Universe continues a de Sitter contraction. If this were the case, the two–derivative equations would suggest that \(W\) returns to the attractor form without \(\Delta\), as in fig. 4, so that \(\epsilon\) could stop moving. Note, in fact, that only if \(\Delta=0\) is eq. 25 solved by a constant \(\sigma\), so that \(\epsilon\) could stay anchored to the attractor value as \(\eta\) decreases beyond the turning point, consistent with the resolution of the singularity illustrated in fig. 2. This picture can also lead to a simple model, where \(\Delta\) disappears suddenly as in the left panel of fig. 4, and the mode functions can be computed exactly. However, the reader will recognize that, as in one–dimensional Quantum Mechanics, a Bunch–Davies mode function \(H_{3/2}^{(1)}\) on the left would call for contributions involving both \(H_{3/2}^{(1)}\) and \(H_{3/2}^{(2)}\) on the right. The presence of both types of mode functions would introduce in the three-point amplitude what could be termed an infrared divergence, since the problem now arises from the upper end of the integration region. This divergence could not be compensated for by the field redefinitions, which yield only finite results, and would also affect higher–point functions, even if further redefinitions, as in [17], were used to remove terms without time derivatives from the three–point amplitude.

Figure 5: Left panel: a sharp transition between an earlier epoch of compression and the expansion that begins at \eta_0 and evolves into the inflationary scenario. Right panel: a smooth model of the transition.

We can now explore an alternative option that can avoid these difficulties. It is suggested by the growth experienced by \(W\) in fig. 3 for decreasing values of \(\eta\), before the initial singularity takes over. What if this behavior persisted in a full string treatment, as in fig. 5? If \(\epsilon\) quickly recovered a slow–roll value, say again \(3\,\gamma^2\), the two–derivative equations would compel the Universe to expand again for decreasing values of \(\eta\). The initial singularity would therefore be replaced by a bounce at a conformal time \(\eta_0\) below \(\eta_t\), both typically \({\cal O}\left(-\, \frac{1}{\Delta}\right)\), and the earlier history of the Universe would play no role in the subsequent evolution.

If a bounce were induced at \(\eta_0\) by string corrections, the Bunch–Davies condition should naturally be imposed there, at the instant of maximal compression. One might object that the same inflationary evolution would emerge if one started in a Bunch–Davies vacuum near the end of the climbing phase, without appealing to string–induced modifications of the action in eq. 3 . However, the bounce scenario grants this peculiar choice of initial condition a special role and is somehow suggested by the typical \(W\) in the left panel of fig. 3 before the singularity is reached. Moreover, the assumption that the singularity is somehow overcome by string corrections conforms to standard expectations, and is not foreign to other treatments of non-Gaussianities that we are aware of, starting from the original work of [17], where the de Sitter phase continues indefinitely in the past.

A bounce scenario was actually considered long ago in String Theory by Gasperini and Veneziano [16], under the spell of \(T\)-dualities. It resolves most of the problems connected with our setting, so we shall now explore its consequences, although it creates another problem. In fact, the reader will not fail to recognize that the modifications impinge on the unique role of the climbing behavior, which will continue to underlie our considerations but, strictly speaking, is only unavoidable in the two–derivative setting.

After the bounce, the Universe undergoes a slow–roll evolution that is similar to the scenario originally considered in [17] and reviewed in [18], but for three important distinctive features:

  • the attractor curve for \(W\) is lowered by \(\Delta^2\), consistently with the modified power spectrum of eq. 19 , so that the mode functions depend on \[\omega \;= \;\sqrt{k^2 \;+ \;\Delta^2} \;; \label{omegakappa}\tag{27}\]

  • the expected contribution to the three–point amplitude acquires an oscillatory behavior depending on \(\omega\) and on the conformal time \(\eta_0\) where the expansion begins, which as we stated is typically \({\cal O}\left(-\, \frac{1}{\Delta}\right)\);

  • curvature perturbations are actually singular at the turning point \(\eta_t\) but, as we shall see, in a way that is regulated by the Schwinger–Keldysh contour. The resulting contributions have an oscillatory behavior determined by \(\omega\) and \(\eta_t\).

This final setting seems relatively simple, and the analysis we are about to present seems well motivated and safe, but we are aware that subtleties of various types were considered in the literature, in a number of different contexts, in particular in [23]. In our case, the solutions of the Mukhanov–Sasaki equation rest essentially on Bessel functions of order \(\frac{3}{2}\), and for \(\eta>\eta_0\) the initial Bunch–Davis condition selects the mode functions \[v_\omega(\eta)\;= \;- \;\frac{\sqrt{\pi}}{2}\sqrt{-\eta}\;H_{\frac{3}{2}}^{(1)}(-\omega \eta) \;= \; \frac{1}{\sqrt{2 \,\omega}} \left(1 \;- \;\frac{i}{\omega\,\eta} \right)\;e^{-i\,\omega\,\eta} \;, \label{vlargenegeta}\tag{28}\] with \(\omega\) as in eq. 27 . Consequently, for \(\eta > \eta_0\) the Green functions that determine the three–point amplitude are \[G_>(\eta,\omega) \;= \;\frac{{\cal H}^2}{4 \,M_P^2 \, \sqrt{3 \, \gamma^2 \, \epsilon (\eta)} \;\omega^{3}} \;e^{i \, \omega\,\eta} \left(1\;- \;i\, \omega\,\eta \right) \;, \label{green}\tag{29}\] with \(\epsilon(\eta) = 3 \,\gamma^2\) away from the turning point.

a

Figure 6: The Schwinger–Keldysh contour for the proposed dynamics starts slightly above the real axis of the \(\eta\) plane at \(\eta_0\), proceeds along the semicircle above \(\eta_t\), turns back at the origin, proceeds along the semicircle below \(\eta_t\) and finally ends at \(\eta_0\) slightly below the real axis..

The three-point amplitude thus rests on three types of ingredients, which depend on the two conformal times \(\eta_0\) and \(\eta_t\), both \({\cal O}\left(-\,\frac{1}{\Delta}\right)\) and with \(\eta_t > \eta_0\) in our picture:

  1. a “principal–part” contribution from the region \(\eta_0 \leq \eta< 0\), \(\eta \neq \eta_t\), on a Schwinger–Keldysh contour whose two sides are slightly displaced, as fig. 6, around the real axis, where the Green function is as in eq. 29 with \(\epsilon(\eta)=3\,\gamma^2\);

  2. a contribution emerging from the neighborhood of \(\eta_t\), the turning point where \(\epsilon(\eta)\) vanishes, which is regulated by the Schwinger–Keldysh contour;

  3. a contribution arising from the field redefinition performed in [17] and reviewed in [18], which is identical to the original result, up to the replacement of \(k_i\) with \(\omega_i\).

5 The Three–Point Amplitude↩︎

We can now compute the three–point amplitude, proceeding as outlined in the preceding section. The cubic action derived in [17] and reviewed in detail in [18] reads \[\begin{align} \mathcal{S}^{(3)} &=& M_{P}^2\int d^4 x \Big\{\epsilon^2 e^{a}\dot{\zeta}^2\zeta\,+\,\epsilon^2 e^{a/3}\zeta \partial_k \zeta \partial^k \zeta \,-\,2 \epsilon^2 e^{a} \dot{\zeta} \partial_k \zeta \partial^k (\partial^{-2} \dot{\zeta})\nonumber \\ &+& \frac{3\,\epsilon}{2}\, e^{a}\,\frac{d}{dt_c}\left(\frac{\dot{\epsilon}}{\epsilon \dot{a}}\right)\dot{\zeta}\zeta^2\,-\,\frac{\epsilon^2}{2}\;e^{a} \left[\dot{\zeta}^2\zeta - \zeta \partial_k \partial_l (\partial^{-2}\dot{\zeta})\partial^k \partial^l (\partial^{-2}\dot{\zeta})\right]+f(\zeta)\frac{\delta \mathcal{L}}{\delta \zeta}\Big\}\,, \end{align}\] where “dots” denote derivatives with respect to cosmic time, and where the last term contains the field redefinition originally performed in [17], which emphasizes the slow–roll limit of the different residual contributions, with \[\begin{align} f(\zeta)&=&(\delta + \epsilon)\zeta^2 \;+ \;\frac{6}{\dot{a}}\,\dot{\zeta}\zeta \;- \;\frac{9}{2\,\dot{a}^2}\;e^{-\frac{2}{3}a}\left[\partial_k\zeta \partial^k \zeta - \partial^{-2}\partial_k \partial_l (\partial^k \zeta \partial^l \zeta)\right]\nonumber \\ &+&\frac{3\,\epsilon}{\dot{a}}\;\left[\partial_k \zeta \partial^k (\partial^{-2}\dot{\zeta})\,-\,\partial^{-2}\partial_k \partial_l (\partial^k \zeta \partial^l (\partial^{-2} \dot{\zeta}))\right] \;. \end{align}\] Following the convenient classification in [18], we can now write the three–point amplitude as \[\begin{align} {\cal A}_3 &=& \int \frac{d^3 \vec{k}_1}{(2 \pi)^3} \frac{d^3 \vec{k}_2}{(2 \pi)^3} \frac{d^3 \vec{k}_3}{(2 \pi)^3}e^{i \vec{k}_1 \cdot \vec{x}} e^{i \vec{k}_2 \cdot \vec{y}} e^{i \vec{k}_3 \cdot \vec{z}}(2 \pi)^3 \delta^{(3)}(\vec{k}_1 + \vec{k}_2 + \vec{k}_3)\nonumber \\ && \Big\{\langle O_1 \rangle + \langle O_2 \rangle + \langle O_3 \rangle + \langle O_4 \rangle + \langle O_5 \rangle + \langle O_f \rangle + 2\, \mathrm{perms} \Big\}\; , \end{align}\] where the last term, \[\begin{align} \langle {\cal O}_f \rangle &=& \frac{\mathcal{H}^4}{16\,M_P^4} \frac{\delta + \epsilon}{\epsilon^2} \, \frac{1}{ \omega_2^3 \, \omega_3^3} \;, \end{align}\] with \(\delta\) and \(\epsilon\) as in eq. 18 , originates from the field redefinition in [17], and, as we saw in Section 3, \[\omega_{2,3} \;= \;\sqrt{k_{2,3}^2 \;+ \;\Delta^2} \;,\] since in our case the attractor curve is lowered by \(\Delta^2\).

As we have anticipated, the three–point amplitude is determined by two contributions to the Schwinger–Keldysh contour displayed in fig. 6 with very different origins. Those of the first type, which we shall refer to as “principal–part” contributions, reflect the quasi–de Sitter expansion of the Universe with a constant \(\epsilon=3 \,\gamma^2\) that begins at \(\eta_0\) and leave out the turning point at \(\eta_t\), while those of the second type emerge precisely from the turning point. Both types of contributions originate from the five integrals \[\begin{align} \langle O_1 \rangle &=& 2i M_{P}^2 \int_{ \eta_0}^{0} \frac{d \eta\;\epsilon(\eta)^2}{(\mathcal{H} \eta)^2}\left[\, G_>'(1)\, G_>'(2) \,G_>(3)\; - \mathrm{c.c.}\right] \;, \nonumber \\ \langle O_2 \rangle &=& -\;2i M_{P}^2 \int_{ \eta_0}^{0} \frac{d \eta\;\epsilon(\eta)^2}{(\mathcal{H} \eta)^2}\left(\vec{k}_1 \cdot \vec{k}_2\right)\left[ G_>(1)\, G_>(2) \,G_>(3)\; - \mathrm{c.c.} \right] \;, \nonumber \\ \langle O_3 \rangle &=& -\;2i M_{P}^2 \int_{ \eta_0}^{0} \frac{d \eta\;\epsilon(\eta)^2}{(\mathcal{H} \eta)^2}\left[\frac{\vec{k}_1 \cdot \vec{k}_3}{k_1^2} + \frac{\vec{k}_2 \cdot \vec{k}_3}{k_2^2}\right]\left[ \, G_>'(1) \, G_>'(2) \,G_>(3)\; - \mathrm{c.c.}\right] \;, \nonumber \\ \langle O_4 \rangle &=& i M_{P}^2 \int_{ \eta_0}^{0} \frac{d \eta}{(\mathcal{H} \eta)^2} \, \frac{\epsilon (\eta) }{\mathcal{H}}\,\frac{d^2}{dt_c^2}\log \epsilon(\eta)\, \left[ \, G_>'(1)\, G_>(2) \,G_>(3)\; - \mathrm{c.c.}\right]\;, \nonumber \\ \langle O_5 \rangle &=& -\;i M_{P}^2 \int_{ \eta_0}^{0} \frac{d \eta\;\epsilon(\eta)^3}{(\mathcal{H} \eta)^2} \, \left[1 \, - \, \frac{\left(\vec{k}_2\cdot \vec{k}_3\right)^2}{k_2^2\,k_3^2}\right]\left[ \, G_>'(1) \, G_>'(2) \,G_>(3)\; - \mathrm{c.c.}\right] \;. \label{O15} \end{align}\tag{30}\] In these expressions, the scale factor \(e^\frac{a}{3}\) has been replaced by its de Sitter form \(\frac{1}{- {\cal H}\,\eta}\), and the measure has been similarly adapted to the conformal time \(\eta\). Note that \(O_2\) involves products of three Green functions, while in \(O_1\), \(O_3\) and \(O_5\) the two “primed” factors are differentiated with respect to the conformal time \(\eta\). In \(O_4\), which has no “principal–part” contribution since it contains derivatives of \(\epsilon\), only one of the Green functions is differentiated.

5.1 Principal–Part Contributions↩︎

In this section, we collect the “principal–part” contributions, showing for brevity only those for the equilateral configuration, which yields the largest result. These follow from standard integrations, but the reader should note that the \(H_\frac{3}{2}^{(1)}\) behavior of the mode functions, without other terms involving \(H_\frac{3}{2}^{(2)}\), is crucial to avoid a divergent contribution to \(\langle O_2 \rangle\) from the upper end of the integration region. The results for the different amplitudes are \[\begin{align} \langle O_1 \rangle &=& \frac{\mathcal{H}^4}{144 \, \gamma^2 \, {M_p}^4 \omega^6}\;\Big[4 \left(1 \;- \;\cos (3 \eta_0 \, \omega )\right) \;- \;3 \,\eta_0 \, \omega \sin (3 \eta_0 \,\omega ) \Big] \;, \nonumber\\ \langle O_2 \rangle &=& \frac{\mathcal{H}^4\;k^2}{288 \,\gamma^2 \, \eta_0 \, {M_p}^4 \omega^9} \Big[3 \left(\eta_0^2 \,\omega^2\,-\,3\right) \sin (3 \eta_0 \, \omega )\;+ \;\eta_0 \, \omega \left(10 \cos (3 \eta_0 \,\omega )\,+\,17\right)\Big] \;,\nonumber\\ \langle O_3 \rangle &=& \frac{\mathcal{H}^4}{144 \, \gamma^2 \, {M_p}^4 \omega^6}\;\Big[4 \left(1 \, - \, \cos (3 \eta_0 \, \omega )\right) \;- \;3 \,\eta_0 \, \omega \sin (3 \eta_0 \,\omega ) \Big] \;, \nonumber\\ \langle O_4 \rangle &=& 0 \;, \nonumber\\ \langle O_5 \rangle &=& - \;\frac{\mathcal{H}^4}{128 \, {M_p}^4 \omega^6} \; \Big[4 \left(1 \, - \, \cos (3 \eta_0 \, \omega )\right) \;- \;3 \,\eta_0 \, \omega \sin (3 \eta_0 \,\omega ) \Big] \;, \label{easy95contributions} \end{align}\tag{31}\] where, as we saw in Section 3, \[\omega \;= \;\sqrt{k^2 \;+ \;\Delta^2} \;.\]

The contribution from \(\langle O_4 \rangle\) vanishes identically, since it involves derivatives of \(\epsilon\), which is constant away from the turning point in our approximation, while the others are proportional in the equilateral configuration. \(\langle O_5 \rangle\) is subdominant in \(\gamma\) or \(\epsilon\), and for this reason it was ignored in [17]. All these contributions oscillate around the result in [17], which can be recovered if one ignores the trigonometric functions.

5.2 Turning–Point Contributions↩︎

We can now describe the contributions that originate from the neighborhood of the turning point, where the dynamics is generally dominated by the hard exponential. Again, for brevity, we only show the equilateral results. In the neighborhood of the turning point, the integration variable \(\eta\) can be traded for the parametric time \(\tau\), and one can rely on the explicit forms of \(\varphi(\tau)\) and \(a(\tau)\) in eq. 11 , so that \[\epsilon(\tau) \;= \;3 \left(\frac{1 \;- \;\tau^2}{1 \;+ \;\tau^2}\right)^2 \;.\] The microscopic data of the “critical” climbing scalar can be then linked to \(\cal {H}\) as explained in Sections 2 and 3, and in particular making use of eq. 17 . In the two branches of the Schwinger–Keldysh contour, these local contributions are computed relying on \[\frac{1}{z \;\pm \;i\,\xi} \;= \;PP\left[\frac{1}{z} \right] \;\mp \;i\, \pi \, \delta(z) \;,\] and on the residues that emerge from the two portions of the integrands, which have individually (anti) analytic properties there.

There is a subtlety here: the Mukhanov–Sasaki variable and the deviation of the scalar field from its cosmological trajectory are connected by coordinate transformations that are singular at the turning point, so that they are somehow inequivalent in the case under scrutiny. With the former choice and thus focusing on curvature perturbations, which are closer to actual experiments, the turning point leaves tangible and calculable imprints in the three–point amplitude, whose interesting features we can now describe.

Only three of the five amplitudes in eqs. 30 play a role, since the others are not singular at the turning point, and the resulting equilateral contributions read \[\begin{align} \langle O_1 \rangle &=& -\;\frac{3\,\pi\, \mathcal{H}^5 \,{e}^\frac{1}{12}}{16\, \gamma^3 \,\eta_t^2\, {M_p}^4 \omega^9} \; \Big[\eta_t \omega \left(\eta_t^2 \omega ^2-3\right) \sin (3 \eta_t \omega )\;+ \;\left(3 \eta_t^2 \omega ^2-1\right) \cos (3 \eta_t \omega )\Big] \;, \nonumber\\ \langle O_2 \rangle &=& 0 \;,\nonumber\\ \langle O_3 \rangle &=& -\;\frac{3\,\pi\, \mathcal{H}^5 \,{e}^\frac{1}{12}}{16\, \gamma^3 \,\eta_t^2\, {M_p}^4 \omega^9} \; \Big[\eta_t \omega \left(\eta_t^2 \omega ^2-3\right) \sin (3 \eta_t \omega )\;+ \;\left(3 \eta_t^2 \omega ^2-1\right) \cos (3 \eta_t \omega )\Big] \;, \nonumber\\ \langle O_4 \rangle &=& \frac{\pi \mathcal{H}}{288 \,{e}^\frac{1}{4} \gamma ^3 \eta_t^6 {M_p}^4 \omega ^9} \Big\{\eta_t \omega \sin (3 \eta_t \omega ) \Big[\eta_t \Big(9 \eta_t^3 \omega ^4-65 \eta_t \omega^2-90 {e}^\frac{1}{4} \eta_t^2 \mathcal{H}^3 \left(\eta_t^2 \omega ^2-3\right)\nonumber \\ &+& 2 {e}^\frac{1}{6} \eta_t \mathcal{H}^2 \left(\eta_t^2 \omega ^2+15\right) \;+ \;9 {e}^\frac{1}{12} \mathcal{H} \left(\eta_t^4 \omega ^4-5 \eta_t^2 \omega ^2+4\right)\Big)\;+ \;60\Big]\nonumber \\ &+&\cos (3 \eta_t \omega ) \Big[\eta_t \Big(33 \eta_t^3 \omega^4-81 \eta_t \omega^2+90 {e}^\frac{1}{4} \eta_t^2 \mathcal{H}^3 \left(1-3 \eta_t^2 \omega ^2\right)\nonumber \\ &-& 2 {e}^\frac{1}{6} \eta_t \mathcal{H}^2 \left(3 \eta_t^4 \omega^4+12 \eta_t^2 \omega^2-5\right)\;+ \;3 {e}^\frac{1}{12} \mathcal{H} \left(9 \eta_t^4 \omega ^4-17 \eta_t^2 \omega ^2+4\right)\Big) \;+\;20\Big]\Big\} \;, \nonumber\\ \langle O_5 \rangle &=& 0 \;. \label{turn95contributions} \end{align}\tag{32}\] Note that these contributions include terms accompanied by different powers of \(\mathcal{H}\). In view of eq. 23 , this has an interesting effect: it introduces a strong dependence of the three–point amplitude, and of the corresponding observables, on the number of inflationary \(e\)-folds.

6 Equilateral \(f_{NL}\)and Comparisons↩︎

Figure 7: Left panel: f_{PP} (black, solid) compared with f_{NL}^{(M)} (red, dot-dashed), for \eta_0= -\,\frac{2}{\Delta}. Right panel: f_{PP} (black, solid) compared with f_{NL}^{(M)} (red, dot-dashed), for \eta_0= -\,\frac{5}{\Delta}. The results displayed correspond to the choice \epsilon=0.03, and the oscillations are determined by \eta_0.

Our starting point is the definition of the convenient dimensionless quantity \[f_{NL}(k)\;= \;40 \, \frac{\gamma^4 \, M_{p}^4\,\Delta^{6}}{ \, \mathcal{H}^4} \left[\left(\frac{k}{\Delta}\right)^2 + 1\right]^{2 \, \nu} \, \langle \zeta(k) \zeta(k) \zeta(k) \rangle \;, \label{eq:fnl}\tag{33}\] where we have identified the reference scale \(k^\star\) with \(\Delta\) and \[\langle \zeta(k) \zeta(k) \zeta(k) \rangle \;= \;\langle O_1 \rangle + \langle O_2 \rangle + \langle O_3 \rangle + \langle O_4 \rangle + \langle O_5 \rangle + \langle O_f \rangle + 2\, \mathrm{perms} \;.\] We also write \[f_{NL}(k) \;= \;f_{PP}(k) \;+ \;f_{t}(k) \;, \label{ppt}\tag{34}\] and in the following we consider separately the two contributions of Sections 5.1 and 5.2, comparing the former with the original result in [17], which in the same conventions reads \[f_{NL}^{(M)}(k)\;= \;\frac{5}{12}\, \left(\frac{k}{\Delta}\right)^{4 \, \nu - 6} \left[-\,3\, \gamma^4\,+\,17\, \gamma^2\,+\,2\, \delta \right] \;. \label{eq:fnlM}\tag{35}\]

Figure 8: Left panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=61. Right panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=62. The results displayed correspond to the choice \epsilon=0.03, and the oscillations are determined by \eta_t.
Figure 9: Left panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=63. Right panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=64. The results displayed correspond to the choice \epsilon=0.03, and the oscillations are determined by \eta_t.
Figure 10: Left panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=65. Right panel: f_{t} for \eta_t= -\,\frac{2}{\Delta} (blue, solid) and for \eta_t= -\,\frac{4}{\Delta} (orange, dashed) for N=66. The results displayed correspond to the choice \epsilon=0.03, and the oscillations are determined by \eta_t.

Here we include for completeness the subdominant term contributed by \(\langle O_5 \rangle\), which was ignored in the original work. This result further simplifies for the Lucchin–Matarrese attractor, to which we shall refer for definiteness, since in this case \[\delta \;= \;- \;3\,\gamma^2 \;,\] so that the contribution from the field redefinition vanishes, but in general slow-roll scenarios this last relation would not hold.

Fig. 7 compares \(f_{PP}\) to \(f_{NL}^{(M)}\), within the range \(\frac{\Delta}{2}<k<2\,\Delta\), for two choices of the conformal time \(\eta_0\) when the expanding phase begins. The two signals are comparable in size, and the values of \(\eta_0\) mostly affect the frequency of the oscillations around the original result in [17].

Figs. 8,  9 and 10 compare the behavior of \(f_t\), within the range \(\frac{\Delta}{2}<k<2\,\Delta\), for different numbers of \(e\)–folds in the interval \(61 \leq N \leq 66\). The solid curves refer to \(\eta_t= -\,\frac{2}{\Delta}\) and the dashed ones to \(\eta_t= -\,\frac{4}{\Delta}\). Note that for \(63\leq N \leq 65\) the values of \(f_t\) are appreciable but not too large and lie within ranges not excluded by the Planck Collaboration [19] (see also [20] for forecasts related to future experiments) 4. Note how in this case the choice of \(\eta_t\) has significant effects on the size of \(f_t\). For \(N<63\) or \(N>65\) the peaks typically lie beyond the range identified by the Planck collaboration, so there would be a clear preference for a small window in the number of \(e\)-folds, if these results were to play a role in the comparison with future data. This can also be seen from fig. 11, which shows the dependence of \(f_t\) for \(k= 1.2\,\Delta\) (solid line) and for \(k=2\,\Delta\) (dashed line) on the number \(N\) of inflationary \(e\)-folds, and where a plateau is clearly present for \(63\leq N \leq 65\).

In summary, \(f_{PP}\) undergoes small oscillations around the original result in [17], but \(f_t\) can offer some prospects for future detection and has the interesting feature of selecting a very small range of \(e\)-folds where our results are sizable and yet not too large. In the present setup, as we have stressed, \(\eta_0\) should be at least slightly below \(\eta_t\), with both \({\cal O}\left(-\,\frac{1}{\Delta}\right)\). Let us also recall that, according to the analysis in [14], \(k=\Delta\) translates into \(\ell \simeq 11\) in the angular power spectrum of the CMB, which justifies the choice to concentrate our attention on the range \(\frac{\Delta}{2}<k<2\,\Delta\). All preceding results become unreliable for values of \(k\) well beyond \(\Delta\), since the \(k \to \infty\) limit is equivalent to letting \(\Delta \to 0\), which would extend the integrals to the standard range \(- \infty< \eta<0\). This would require the usual contour deformation that makes it move farther and farther away from the real axis for large negative values of \(\eta\), which we avoided within the limited range of interest for \(k\). The modification can be simply implemented for \(\eta_0\), while the way to treat the turning point is less clear to us. Note that, within the limited range for \(k\) that we are focusing on, the primordial effects that we are describing can have some direct bearing on the data, without the complications due to the subsequent evolution of the Universe, as is the case for the power spectrum.

a

Figure 11: The typical behavior of \(f_t\) for \(k=1.2\,\Delta\) (blue, solid) and for \(k=2\,\Delta\) (orange, dashed) and \(\epsilon=0.03\), as the number of \(e\)-folds increases from 60 to 67. Note the plateau present in the region \(63<N<65\). This is due to the combined effect of \(\langle O_{1}\rangle\) and \(\langle O_{3}\rangle\), which decrease in absolute value for increasing values of \(N\), and to the opposite behavior of \(\langle O_4 \rangle\)..

7 Conclusions↩︎

This paper elaborates on a setup for the three–point amplitude of curvature perturbations in scenarios inspired by high–scale supersymmetry breaking in String Theory, which involve a scalar field climbing up a steep exponential potential. The analysis rests on three additions to the original work in [17]. The first is a scale \(\Delta\) [11], [14], which introduces a low–frequency cut in the primordial power spectrum of [12]. If that region were accessible to us, the modification would resonate with the lack of power in the first CMB multipoles, and some evidence for this effect was collected in [14]. The scale \(\Delta\) is also related to the conformal time \(\eta_t\) of the turning point that completes the climbing phase, whose values are typically in the range \(-\,\frac{4}{\Delta}< \eta_t< -\,\frac{1}{\Delta}\). The initial singularity, inevitably present at earlier times in the two–derivative formulation based on eq. 3 , introduces what could be termed an ultraviolet divergence in three–point amplitudes. The string corrections are supposed to overcome this problem, albeit in ways that we can only guess at this time. The tendency that manifests itself in fig. 3 before the singular behavior takes over suggests that in a full treatment \(W\) could approach the real axis from below as \(\eta \to - \,\infty\). If \(\epsilon\) quickly resumed a slow-roll value for \(\eta<\eta_t\), the singularity would be replaced by a bounce at a conformal time \(\eta_0<\eta_t\), typically also \({\cal O}\left(-\,\frac{1}{\Delta}\right)\).

This type of scenario resonates with ideas proposed long ago in String Theory by Gasperini and Veneziano [16], and the Bunch–Davies condition is naturally imposed at \(\eta_0\), the instant of maximum contraction for the Universe. The Schwinger–Keldysh contour of fig. 6 then yields a “principal–part” contribution, which undergoes small oscillations around the original result of [17] at a frequency determined by \(\eta_0\). However, curvature perturbations, which appear close in spirit to the actual measurements, receive additional imprints from the turning point. These contributions can be computed in closed form relying on the “critical” microscopic dynamics of the climbing scalar, which dominates there, and the determination of \(\Delta\) obtained in [14] indicates that the resulting signal could be potentially detectable, and yet compatible with restrictions that emerged from the Planck collaboration [19], [20], at least within the window \(62<N<66\) for the number of inflationary \(e\)-folds.

In analogy with the Feynman prescription for Quantum Field Theory, the \(\epsilon\)-prescription for the Schwinger–Keldysh contour is meant to switch off all contributions in the far past. Here we start from a finite value \(\eta_0\) of the conformal time, which is typically a number of order one in units of \(-\,\frac{1}{\Delta}\), and consequently the amplitudes thus obtained are only reliable for momenta that are comparable to \(\Delta\), but not far beyond. In terms of CMB multipoles, \(k \simeq \Delta\) corresponds to \(\ell \simeq 11\) [14], so the results should not be taken at face value much beyond \(\ell=20\).

As we have stressed, the modifications that we were led to introduce at very early times yield well–defined results but lack a compelling dynamical justification. They also have the virtue of leading to a picture that lies in slow–roll regime. Within their limitations, the present results add some novel features to the scenarios explored for the three–point amplitude of curvature perturbations, relying effectively on a single parameter \(\Delta\) that does not appear foreign to the CMB [14].

The analysis can be extended to amplitudes involving tensor perturbations. These receive no contributions from the turning point and therefore undergo small oscillations around the original results in [17], like the “principal part” contribution \(f_{PP}\) to the three–point function of curvature perturbations of fig. 7.

Acknowledgments↩︎

We are grateful to N. Bartolo, P. Natoli, M. Peloso, and especially to E. Dudas, A. Gruppuso, S. Patil and G.L. Pimentel, for stimulating discussions. This work was supported in part by Scuola Normale and in part by INFN (IS GSS-Pi).

References↩︎

[1]
.
[2]
.
[3]
.
[4]
.
[5]
.
[6]
.
[7]
.
[8]
.
[9]
.
[10]
.
[11]
.
[12]
.
[13]
.
[14]
.
[15]
.
[16]
.
[17]
.
[18]
.
[19]
.
[20]
.
[21]
.
[22]
.
[23]
.

  1. These solutions were previously found in [7] for the potentials emerging from the ten–dimensional models, and then in [8], but the peculiar transition for increasing values of \(\gamma\) was first noted in [6].↩︎

  2. Mechanisms to suppress the first few CMB multipoles were previously considered in other contexts, including in [13].↩︎

  3. The relevant values are \(\gamma=\frac{5}{3}\) for the \(SO(16) \times SO(16)\) model of [1], where \(T\) would reflect the one–loop vacuum energy, and \(\gamma=1\) for the tachyon–free non–supersymmetric orientifolds of [2] and [3].↩︎

  4. Those limits are probably too conservative for the effects under scrutiny, since they were obtained within ranges of \(\ell\) that are much wider than those relevant to the present analysis. The interesting message here is the unusual emergence of upper and lower bounds on \(N\) (see Fig. 11), which could be expanded.↩︎