On the full quantum trispectrum in multi-field DBI inflation
Shuntaro Mizuno∗ ♯, Frederico Arroja† ♭, and Kazuya Koyama‡♮
arXiv:0907.2439v2 [hep-th] 20 Oct 2009
♯School of Physics and Astronomy, University of Nottingham, University Park,
Nottingham NG7 2RD, UK; Research Center for the Early Universe (RESCEU),
Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan.
♭Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan.
♮Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth PO1 3FX, UK.
(Dated: October 20, 2009)
We compute the leading order connected four-point function of the primordial curvature perturbation coming from the four-point function of the fields in multi-field DBI inflation models. We
confirm that the consistency relations in the squeezed limit and in the counter-collinear limit hold as
in single field models thanks to special properties of the DBI action. We also study the momentum
dependence of the trispectra coming from the adiabatic, mixed and purely entropic contributions
separately and we find that they have different momentum dependence. This means that if the
amount of the transfer from the entropy perturbations to the curvature perturbation is significantly
large, the trispectrum can distinguish multi-field DBI inflation models from single field DBI inflation
−2
and the
models. A large amount of transfer TRS ≫ 1 suppresses the tensor to scalar ratio r ∝ TRS
equi
−2
amplitude of the bispectrum fNL ∝ TRS and so it can ease the severe observational constraints on
the DBI inflation model based on string theory. On the other hand, it enhances the amplitude of
equi
equi 2
2
the trispectrum τNL
∝ TRS
fNL
for a given amplitude of the bispectrum.
I.
INTRODUCTION
Precise measurements of the cosmic microwave background (CMB) anisotropies such as those obtained by the
WMAP satellite [1] provide valuable information on the very early universe. Any theoretical model that attempts to
explain the evolution of the universe before the big bang nucleosynthesis will also have to explain the observed CMB
anisotropies. Even though these anisotropies are almost Gaussian, a small amount of non-Gaussianity is still allowed
by the data [2, 3, 4] and future experiments such as PLANCK [5] might detect this small amount of the primordial
non-Gaussianity. Even though the simplest slow-roll single field inflation models predict that the non-Gaussianity of
the fluctuations should be very difficult to be detected [6], recently, theoretical models which can produce sizeable
non-Gaussianity has been extensively studied by many authors [7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21,
22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53,
54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66].
Among them, the Dirac-Born-Infeld (DBI) inflation model can produce large non-Gaussianity by the fact that the
sound speed of the perturbations can be much smaller than one due to the non-trivial form of the kinetic term [67, 68].
DBI inflation models are also well motivated by string theory [13, 27, 69, 70, 71, 72]. In this model, the inflaton fields
are identified with the positions of a moving D3 brane in the six-dimensional internal space. The dynamics of the D3
brane is described by the DBI action. However, recently it has been pointed out that single field DBI inflation driven
by a mobile D3 brane with large non-Gaussianity might contradict with the current WMAP data. For current and
stringent observational constraints and consequences of DBI-inflation see [73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83].
One way to avoid these constraints is to consider multi-field DBI models [36]. Since the position of the brane in
each compact direction is described by a scalar field, DBI inflation is naturally a multi-field inflationary model [84].
As first pointed out by [85], in multi-field inflation models, the curvature perturbation is modified on large scales
due to the entropy perturbation. Even though there are some works considering multi-field inflationary models with
non-canonical kinetic terms [45, 86, 87, 88], the consistent analysis for the entropy modes in the multi-field DBI
inflation model has started only very recently [35, 36, 37].
In [35, 36, 37], the three-point functions in the small sound speed limit and at leading order in the slow-roll expansion
were obtained and it was shown that in addition to the purely adiabatic three-point function, there exists a mixed
component hQσ (k1 )Qs (k2 )Qs (k3 )i where Qσ and Qs are the adiabatic and the entropy perturbations, respectively.
Since the momentum dependence of the three-point function from the adiabatic modes was shown to be identical
∗
†
‡
shuntaro.mizuno@nottingham.ac.uk
arrojaf@yukawa.kyoto-u.ac.jp
Kazuya.Koyama@port.ac.uk
2
with the mixed component, the shape of the bispectrum of the curvature perturbations remains the same as in the
single-field case, while the amplitude is affected by the entropy perturbation.
Even though previous works on the non-Gaussianity in the multi-field DBI inflation model are mainly limited to the
bispectrum, it is expected that the cosmic microwave background (CMB) trispectrum also provides strong constraints
on early universe models. At the moment, the constraints are rather weak given by |τN L | < 108 [89, 90], where τN L
denotes the size of the trispectrum. However, PLANCK will tighten the constraints significantly reaching |τN L | ∼ 560
[91]. Although these estimations are obtained assuming local type non-Gaussianity, we expect similar constraints can
be obtained for non-Gaussianity from DBI inflation. It is also worth noting that the analysis in the single field DBI
inflation model shows that the trispectrum from the contact interaction diagram is enhanced in the small sound speed
limit as τN L ∼ 1/c4s [28]. Even though the contribution from the scalar exchange interaction was overlooked in this
analysis, recently, it has been confirmed that the same scaling holds for the trispectrum from this interaction [62, 63].
As in the bispectrum case [92], the observational constraints depend on the shape of the wave vectors’ configuration.
Therefore, it is important to calculate the full shape dependence of the trispectrum in the multi-field DBI inflation
model. For the details of the observations of the CMB trispectrum, see [93, 94, 95].
For this purpose, recently we have calculated the quantum trispectrum from the contact interaction and obtained
the simple understanding of the origin of the interaction terms of the multi-field DBI inflation [64]. In this paper, as a
natural continuation, we will calculate the scalar exchange trispectrum and obtain the complete theoretical prediction
for the quantum trispectrum of multi-field DBI inflation. For a related work, where only the purely entropic component
of the quantum scalar exchange trispectrum was calculated, see version three of [55].
The structure of this paper is as follows. In the next section, we shall introduce the model and some basic notation.
In section III we summarise the results obtained in [64] for the contact interaction trispectrum. After calculating the
scalar exchange trispectrum in section IV, the momentum dependence of the trispectra coming from the adiabatic,
mixed, and purely entropic contributions are studied separately in section V. The observational constraints on the
multi-field DBI inflation model are discussed in section VI. Section VII is devoted to the conclusions.
II.
THE MODEL
In this section, we will introduce the multi-field Dirac-Born-Infeld (DBI) inflationary model. We will present the
background evolution equations and define some basic notation.
The multi-field DBI inflation model is described by the following action [96]
Z
i
√ h
1
S=
d4 x −g R + 2P̃ (X̃, φI ) ,
2
q
1
I )X̃ − 1 − V (φI ) ,
1
−
2f
(φ
P̃ (X̃, φI ) = −
(1)
f (φI )
where we have set 8πG = 1, R is the Ricci scalar, φI are the scalar fields (I = 1, 2, ..., Nφ ), f (φI ) and V (φI ) are
functions of the scalar fields determined by string theory configurations and X̃ is defined in terms of the determinant
D ≡ det(δνµ + f GIJ ∂ µ φI ∂ν φJ ) as X̃ = (1 − D)/(2f ). Here GIJ is the metric in the field space. We assume that P̃
is a well behaved function of φI and X̃. It is also shown that X̃ is related to the kinetic term of the scalar fields as
[36, 37]
[I
J]
[I
K]
[I
L]
X̃ = GIJ X IJ − 2f XI XJ + 4f 2 XI XJJ XK − 8f 3 XI XJJ XKK XL ,
1
X IJ ≡ − g µν ∂µ φI ∂ν φJ , XIJ = GIK X KJ ,
2
(2)
(3)
where the brackets denote antisymmetrization. It is worth noting that even though X̃ = X ( = GIJ X IJ ) in the
homogeneous background, this does not hold if we take into account the inhomogeneous components.
In the background, we are interested in flat, homogeneous and isotropic Friedman-Robertson-Walker (FRW) universes described by the line element
ds2 = −dt2 + a2 (t)δij dxi dxj ,
(4)
where a(t) is the scale factor. The Friedman equation and the continuity equation read
3H 2 = E0 ,
(5)
Ė0 = −3H E0 + P̃0 ,
(6)
3
where the Hubble rate is H = ȧ/a, a dot denotes derivative with respect to cosmic time t, E0 is the total energy of
the fields which is given by
E0 = 2X0IJ P̃0,X IJ − P̃0 ,
(7)
and the subscript 0 denotes that the quantity is evaluated in the background.
For this model the speed of propagation of scalar perturbations (“speed of sound”), cs , is given by
!
P̃,X̃
2
.
cs ≡
P̃,X̃ + 2X̃ P̃,X̃ X̃ 0
(8)
Since we are interested in the inflationary background, we assume that the form of f (φI ) and V (φI ) are chosen
so that inflation is realized at least for 60 e-foldings. In order to characterize this background, we define the slow
variation parameters, analogues of the slow-roll parameters, as:
ǫ=−
Ḣ
X0
= 2 ,
H2
H cs
ι=
ǫ̇
,
ǫH
χ=
c˙s
.
cs H
(9)
We should note that these slow variation parameters are more general than the usual slow-roll parameters and that
the smallness of these parameters does not imply that the field is rolling slowly. We assume that the rate of change
of the speed of sound is small (as described by χ) but cs is otherwise free to change between zero and one.
We shall consider perturbations on this FRW background. We decompose the scalar field φI into the background
value φI0 and perturbation QI in the flat gauge as,
φI (x, t) = φI0 (t) + QI (x, t) .
(10)
Furthermore, as was done in [97], we decompose the perturbations into instantaneous adiabatic and entropy perturbations, where the adiabatic direction corresponds to the direction of the background fields’ evolution while the
entropy directions are orthogonal to this. We introduce an orthogonal basis eIn , with n = 1, 2, ..., Nφ, in the field space
so that the orthonormality condition are given by [35]
eIn emI =
1 − c2s
1
δmn −
δm1 δn1 ,
cs
cs
(11)
where the adiabatic basis is defined as
eI1 =
r
cs I
φ̇ .
2X0 0
(12)
Notice that the length of the basis vector eI0 is cs and that of the other basis vectors is 1/cs . If we consider the
two-field case (I = 1, 2), the field perturbations are decomposed into the adiabatic field Qσ and the entropy field Qs
as
QI = Qσ eI1 + Qs eI2 .
(13)
Hereafter, for simplicity, we will concentrate on the two-field case although the extension to more fields is straightforward.
For the single field model, since the comoving curvature perturbation R is given by
√
cs H
R= √
Qσ∗ ,
(14)
2X0 ∗
the power spectrum of the primordial quantum fluctuation was given by [68]
PR∗ (k) =
E02
1
1 H2
,
=
2
36π E0 + P̃0
8π 2 cs ǫ
(15)
where it should be evaluated at the time of horizon crossing cs ∗ k = a∗ H∗ . The spectral index is
nR∗ − 1 =
d ln PR∗ (k)
= −2ǫ − ι − χ.
d ln k
However, as explained in the following sections, this will not be the case for multi-field models.
(16)
4
III.
THE TRISPECTRUM FROM THE CONTACT INTERACTION
In this and next sections, we will derive the four-point functions of the field perturbations at leading order in
the slow-roll expansion and in the small sound speed limit. There are two important tree-level diagrams for the
trispectrum. One is a diagram where the interaction occurs at a point, i.e. a contact interaction diagram and the
other is a diagram where a scalar mode is exchanged. In the single field DBI inflation model, it was recently shown
by [62, 63] that both contributions are comparable. Because it is expected that this is also true for the two-field DBI
inflation model, it is necessary to consider both contributions. In this section, we summarize the calculation of the
four-point functions of the field perturbations at horizon crossing coming from the contact interaction diagram [64].
In the next section, we shall calculate the scalar exchange four-point function.
Using the approximations mentioned above and following the ADM formalism [6, 17, 18, 98], the action up to fourth
order can be calculated as [64]
Z
c2
1
a3
(main)
(17)
=
S(2)
dtd3 x 2 Q̇2σ + Q̇2s − s2 ∂i Qσ ∂ i Qσ + ∂i Qs ∂ i Qs ,
2
cs
a
Z
1
a3
c2
(main)
=
S(3)
(18)
∂i Qs ∂ i Qs − ∂i Qσ ∂ i Qσ Q̇σ − 2 ∂i Qσ ∂ i Qs Q̇s ,
dtd3 x p
Q̇3σ + Q̇σ Q̇2s + s2
2
a
2X0 c7s
"
Z
2c2
1
a3
(main)
3
5Q̇4σ + 6Q̇2σ Q̇2s + Q̇4s − 2s 3Q̇2σ ∂i Qσ ∂ i Qσ − Q̇2σ ∂i Qs ∂ i Qs + 4Q̇σ Q̇s ∂i Qσ ∂ i Qs
=
S(4)
dx dt 5
16
cs X 0
a
+Q̇2s ∂i Qσ ∂ i Qσ + Q̇2s ∂i Qs ∂ i Qs
#
2
2
2
c4s
i
i
j
i
i
+ 4
.
∂i Qσ ∂ Qσ − 2 ∂i Qσ ∂ Qσ ∂j Qs ∂ Qs + 4 ∂i Qσ ∂ Qs + ∂i Qs ∂ Qs
a
(19)
It is worth noting that Eqs. (17), (18) and (19) can also be derived using a simpler method based on a Lorentz boost
from the frame where the brane is at rest to the frame where the brane is moving [64].
Based on these actions, we can obtain the third- and fourth-order interaction Hamiltonian densities as [64]
o
−a3
c2s n
int
2
3
i
i
i
H(3) = p
Q̇σ + Q̇σ Q̇s − 2 Q̇σ ∂i Qσ ∂ Qσ − ∂i Qs ∂ Qs + 2Q̇s ∂i Qσ ∂ Qs
,
(20)
a
2 2X0 c7s
2o
2
c2s n
a3
4
2 2
i
i
int
,
(21)
Q̇
Q̇
+
Q̇
Q̇
+
Q̇
+
∂
Q
∂
Q
∂
Q
∂
Q
H(4)
=
i s
s
i s
s
σ
σ s
s
σ
4X0 c5s
a2
where here, Qσ and Qs are the interaction picture fields. One important remark is that while the cubic part of Hint
is the opposite sign of the cubic Lagrangian density Lint , this is generally not true at fourth order.
From the second-order action, we can solve for the perturbations and quantize them according to the standard
procedures of quantum field theory:
Z
1
Qn (η, x) =
d3 k un (η, k)an (k) + u∗n (η, −k)a†n (−k) eik·x ,
(22)
(2π)3
where an (k) and a†n (−k) are the annihilation and creation operator respectively, which satisfy the usual commutation
relations:
an (k1 ), a†m (k2 ) = (2π)3 δ (3) (k1 − k2 )δnm ,
[an (k1 ), am (k2 )] = a†n (k1 ), a†m (k2 ) = 0 .
(23)
At leading order, solutions for the mode functions are given by
un (η, k) = An
1
(1 + ikcs η) e−ikcs η .
k 3/2
(24)
The two point correlation function is then obtained as
h0|Qn (η = 0, k1 )Qm (η = 0, k2 )|0i = (2π)3 δ (3) (k1 + k2 )PQn
2π 2
δnm ,
k13
(25)
5
where the power spectrum is defined as
PQn =
|An |2
,
2π 2
|Aσ |2 = |As |2 =
H2
≡ N2 ,
2cs
(26)
and it should be evaluated at the time of the sound horizon crossing cs ∗ k1 = a∗ H∗ .
In terms of these quantum operators, the connected four-point correlation function coming from the contact interaction diagram is given by [99]
hΩ|Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 )|ΩiCI =
Z 0
i
h
int
(η) |0i ,
dηh0| Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 ), H(4)
−i
(27)
−∞
int
where Qm on the r.h.s. of the equation are the interaction picture fields and H(4)
is the integral of the fourth-order
R 3 int
int
Hamiltonian density (21), that is, H(4) = d xH(4) .
The purely adiabatic, purely entropic and mixed components are respectively given by
4
X
H8
1
ki )
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiCI = (2π)3 δ (3) (
(−36A1 ) ,
4
6
2X0 cs Πi=1 ki3
i=1
A1 =
Π4i=1 ki2
,
K5
K=
4
X
ki ,
(28)
(29)
i=1
4
X
H8
1
1
ki )
−
hΩ|Qs (0, k1 )Qs (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiCI = (2π)3 δ (3) (
A
2 ,
2X0 c6s Π4i=1 ki3
8
i=1
3(k3 + k4 ) 12k3 k4
k 2 k 2 (k3 · k4 )
+
1+
+ perm.
A2 = 1 2 3
K
K
K2
4
X
H8
−1
1
CI
3 (3)
ki )
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qs (0, k3 )Qs (0, k4 )|Ωi = (2π) δ (
6A1 + A3 ,
2X0 c6s Π4i=1 ki3
2
i=1
2 2
3(k3 + k4 ) 12k3 k4
k k (k3 · k4 )
1+
,
+
A3 = 1 2 3
K
K
K2
(30)
(31)
(32)
(33)
where “perm” denotes twenty three permutations of {k1 , k2 , k3 , k4 }.
Following the analysis of [100], if we describe the conversion of the entropy perturbation into the curvature perturbation by a transfer coefficient TRS , the final comoving curvature perturbation R is expressed in terms of the
adiabatic and entropy field perturbations as
√
√
cs H
cs H
R = Aσ Qσ∗ + As Qs∗ , Aσ = √
, As = TRS √
.
(34)
2X0 ∗
2X0 ∗
This implies that the final comoving curvature perturbation power spectrum is
2
PR∗ ,
PR = 1 + TRS
(35)
where the horizon crossing power spectrum PR∗ was given by Eq. (15).
At leading order, the connected four-point function of R coming from the contact interaction diagram is given by 1
hR(k1 )R(k2 )R(k3 )R(k4 )iCI = A4σ hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iCI
+A2σ A2s hQσ (k1 )Qσ (k2 )Qs (k3 )Qs (k4 )iCI + 5 perm.
+A4s hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iCI ,
1
(36)
There are other contributions to the connected four-point function of R, for example those coming from the five- and six-point functions of
the scalar fields that can be expressed in terms of the power spectrum using Wick’s theorem. In this work, we ignore these contributions
because we are only interested in the intrinsically quantum four-point function. These extra terms are generated by the non-linear
evolution outside the horizon and can be important
6
where “5 perm.” denotes five permutations of the four-momenta. Using Eqs. (28), (30) and (32), it is possible to
show a simple relation among the three different terms in the previous equation, as
hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iCI + hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iCI = hQσ (k1 )Qσ (k2 )Qs (k3 )Qs (k4 )iCI + 5 perm.
(37)
This relation is perhaps less surprising if one notes that the different fourth-order interactions present in Eq. (21)
also obey a similar relation, i.e. the purely adiabatic plus the purely entropic interactions (vertices) are equal to the
mixed vertices.
Using the previous results, the contact interaction four-point function of R can be written as
2
CI
4
hR(k1 )R(k2 )R(k3 )R(k4 )i
= Aσ 1 + TRS hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iCI
2
CI
.
(38)
+TRS hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )i
One immediately sees that if the transfer coefficient is non-zero, one can in principle use the trispectrum to distinguish
single field from multi-field DBI inflation.
IV.
THE TRISPECTRUM FROM THE SCALAR EXCHANGE INTERACTION
In this section, we will calculate the four-point functions of the field perturbations at horizon crossing coming
from a diagram where a scalar mode is exchanged. Within the “interaction picture” formalism [99], to calculate the
four-point function resulting from a correlation established via the exchange of a scalar mode, one needs to evaluate
the following time integrals
hΩ|Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 )|ΩiSE
Z 0
Z η
i
i
hh
int
int
(η̃) |0i,
(η) , H(3)
dη
dη̃h0| Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 ), H(3)
=−
−∞
(39)
−∞
int
is the integral of the third-order
where Qm on the r.h.s. of the equation are the interaction picture fields and H(3)
R 3 int
int
Hamiltonian density (20), that is, H(3) = d xH(3) .
The calculation of these time integrals is rather long and tedious. In Appendix A, we present an elegant and simple
way of evaluating these contributions. This method is based on diagrams and simple rules that make the origin and
physical meaning of all the different terms in the following expressions very clear. The purely adiabatic four-point
7
function coming from the scalar exchange diagrams is identical to the single field one [62], and it is given by
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiSE = (2π)3 δ (3) (K)
H4
3
16X0 c9s (k1 k2 k3 k4 ) 2
"
− 9 F1 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F1 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
×
−3c2s (k3 · k4 ) F3 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F3 (−k1 , −k2 , −k12 , k12 , k3 , k4 )
−2(k34 · k4 ) F3 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F3 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+(k1 · k2 ) F4 (−k12 , k1 , k2 , k3 , k4 , k12 ) − F4 (−k12 , −k1 , −k2 , k3 , k4 , k12 )
!
−2(k12 · k2 ) F4 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F4 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
−c4s (k1 · k2 )(k3 · k4 ) F2 (−k12 , k1 , k2 , k12 , k3 , k4 ) − F2 (−k12 , −k1 , −k2 , k12 , k3 , k4 )
−2(k1 · k2 )(k34 · k4 ) F2 (−k12 , k1 , k2 , k3 , k4 , k12 ) − F2 (−k12 , −k1 , −k2 , k3 , k4 , k12 )
−2(k12 · k2 )(k3 · k4 ) F2 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F2 (−k1 , −k2 , −k12 , k12 , k3 , k4 )
!#
+4(k12 · k2 )(k34 · k4 ) F2 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F2 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+23 permutations of{k1 , k2 , k3 , k4 },
(40)
where K = k1 + k2 + k3 + k4 and the information of the shape of the trispectrum is encoded in the Fi functions
defined by Eqs. (B1), (B2), (B3) and (B4) in Appendix B. In the previous expression, “permutations” denote the
other twenty three terms that result from the permutations of {k1 , k2 , k3 , k4 } in the preceding term.
The mixed component of the four-point function coming from the scalar exchange diagrams is
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiSE + 5 terms = (2π)3 δ (3) (K)
H4
3
16X0 c9s (k1 k2 k3 k4 ) 2
"
− 10 F1 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F1 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+2c2s
×
− (k3 · k4 ) F3 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F3 (−k1 , −k2 , −k12 , k12 , k3 , k4 )
+4(k34 · k4 ) F3 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F3 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
−(k1 · k2 ) F4 (−k12 , k1 , k2 , k3 , k4 , k12 ) − F4 (−k12 , −k1 , −k2 , k3 , k4 , k12 )
!
+4(k12 · k2 ) F4 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F4 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
−2c4s (k1 · k2 )(k3 · k4 ) F2 (−k12 , k1 , k2 , k12 , k3 , k4 ) − F2 (−k12 , −k1 , −k2 , k12 , k3 , k4 )
!#
+4(k12 · k2 )(k34 · k4 ) F2 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F2 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+23 permutations of{k1 , k2 , k3 , k4 }.
(41)
8
Finally, the purely entropic four-point function coming from the scalar exchange diagrams is
hΩ|Qs (0, k1 )Qs (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiSE = (2π)3 δ (3) (K)
H4
3
16X0 c9s (k1 k2 k3 k4 ) 2
"
− F1 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F1 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
×
+c2s (k3 · k4 ) F3 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F3 (−k1 , −k2 , −k12 , k12 , k3 , k4 )
+2(k34 · k4 ) F3 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F3 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+(k1 · k2 ) F4 (−k12 , k1 , k2 , k3 , k4 , k12 ) − F4 (−k12 , −k1 , −k2 , k3 , k4 , k12 )
!
+2(k12 · k2 ) F4 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F4 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
−c4s (k1 · k2 )(k3 · k4 ) F2 (−k12 , k1 , k2 , k12 , k3 , k4 ) − F2 (−k12 , −k1 , −k2 , k12 , k3 , k4 )
+2(k1 · k2 )(k34 · k4 ) F2 (−k12 , k1 , k2 , k3 , k4 , k12 ) − F2 (−k12 , −k1 , −k2 , k3 , k4 , k12 )
+2(k12 · k2 )(k3 · k4 ) F2 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F2 (−k1 , −k2 , −k12 , k12 , k3 , k4 )
!#
+4(k12 · k2 )(k34 · k4 ) F2 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F2 (−k1 , −k2 , −k12 , k3 , k4 , k12 )
+23 permutations of{k1 , k2 , k3 , k4 }.
(42)
The remaining possible four-point functions, hΩ|Qσ Qσ Qσ Qs |ΩiSE and hΩ|Qσ Qs Qs Qs |ΩiSE , are zero at leading order
in the slow-roll expansion and in the small sound speed limit.
As in the previous section, the connected four-point function of the curvature perturbation R is related with the
four-point function of the field perturbation as 2
hΩ|R(k1 )R(k2 )R(k3 )R(k4 )|ΩiSE = A4σ hΩ|Qσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )|ΩiSE
+A2σ A2s hΩ|Qσ (k1 )Qσ (k2 )Qs (k3 )Qs (k4 )|ΩiSE + 5 perms
+A4s hΩ|Qs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )|ΩiSE ,
(43)
Again, as in the case of the contact interaction trispectra, the scalar exchange trispectra obey
hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iSE + hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iSE = hQσ (k1 )Qσ (k2 )Qs (k3 )Qs (k4 )iSE + 5 perm.
(44)
Using this fact Eq. (43) can be simplified to
2
hR(k1 )R(k2 )R(k3 )R(k4 )iSE = A4σ 1 + TRS
hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iSE
2
+TRS
hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iSE
.
(45)
Finally, the total four-point function of R, coming from the quantum four-point functions of the fields, is the sum
2
Where again we ignore terms that can be expressed by copies of the power spectrum by using Wick’s theorem.
9
of the contact interaction and scalar exchange contributions as
hR(k1 )R(k2 )R(k3 )R(k4 )itotal = hR(k1 )R(k2 )R(k3 )R(k4 )iCI + hR(k1 )R(k2 )R(k3 )R(k4 )iSE
2
hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iCI + hQσ (k1 )Qσ (k2 )Qσ (k3 )Qσ (k4 )iSE
= A4σ 1 + TRS
2
,
hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iCI + hQs (k1 )Qs (k2 )Qs (k3 )Qs (k4 )iSE
+TRS
(46)
where the different terms can be found in Eqs. (28), (40), (30) and (42). This constitutes one of the main results of
this work.
From Eq. (46) we can read the exact momentum dependence of the trispectrum which is related with the four-point
function as
hR(k1 )R(k2 )R(k3 )R(k4 )itotal = (2π)3 δ (3) (K)TR (k1 , k2 , k3 , k4 ).
(47)
In the next section we shall study the shape of the trispectrum in more detail.
V.
THE SHAPE OF THE TRISPECTRUM AND THE NON-LINEARITY PARAMETER τNL
In this section, we will study the shape dependence of the different trispectra calculated in the previous sections.
We will consider the so-called equilateral configuration. In the next subsection, we shall discuss some consistency
relations that apply for the multi-field DBI trispectrum. After that, we will compute the non-linearity parameter τN L
and we shall plot the different contributions (pure adiabatic and pure entropy) to the total τN L both in the equilateral
case and in a non-equilateral configuration. It is worth noting that from Eqs. (37) and (44), the contribution from
the mixed component can be given from the other two contributions.
A.
Consistency relations
Here, we discuss whether the consistency relations obtained in the single field DBI inflation model can hold also
for the multi-field DBI inflation model. In the following, as in [62], we will define the non-linearity parameter, which
naively parameterizes the size of the trispectrum as
3
4ǫcs
τN L (k1 , k2 , k3 , k4 ) =
TR (k1 , k2 , k3 , k4 )Π4i=1 ki3
2 )
H 2 (1 + TRS
−3
−3
−1
−3
−3
−3
−3
3 3
3 3
3 3
3 3
3 3
3 3
,
× k1 k2 + k3 k4 k13 + k14 + k1 k4 + k2 k3 k12 + k13 + k1 k3 + k2 k4 k12 + k14
(48)
and we discuss behaviours of τN L .
First, we discuss the consistency relation in the squeezed limit. As first pointed out by Maldacena [6] and Seery et
al. [31], all higher order correlators of R in single field inflation obey consistency relations when one of the momentum
vector is very small. Let us consider the limit k1 → 0, which implies that the corresponding mode R(k1 ) leaves the
horizon well before all other modes leave the horizon. By the time the remaining modes exit the horizon, R(k1 ) will
be frozen as a super-horizon mode and its only effect is to deform the background if the mode with k1 which leaves
the horizon earlier than the others is the adiabatic mode (σ). Then in this case, the trispectrum can be expressed as
[31]
hR(k1 )R(k2 )R(k3 )R(k4 )i → −P̃R (k1 )
d
hR(k2 )R(k3 )R(k4 )i,
Hdt
(49)
where the power spectrum P̃R (k) is given by hR(k1 )R(k2 )i = (2π)3 δ (3) (k1 + k2 )P̃R (k1 ). Using the fact that the bis2
−2 3
−1
,
pectrum is Bζ ∝ c−2
s P̃R , it follows that in the squeezed limit the trispectrum is TR ∝ (ǫcs P̃R ). Using P̃R ∝ (cs ǫ)
one finds that the total non-linearity parameter τN L should be at most of order ǫc−2
(i.e.
under
our
approximation
it
s
should vanish) when any of the momentum vectors goes to zero. From Eqs. (28), (30), (32), (40), (41) and (42), we can
show that τN L scales as O(k12 ), which means the consistency relation in the squeezed limit holds for the trispectrum
in multi-field DBI inflation models.
10
However, if the mode with k1 is the entropy mode (s), this explanation can not be applied. Instead, the fact that at
leading order in slow-roll, the third-order action is written as Eq. (18) and there are no three point functions like (sss)
and (sσσ) turns out to be crucial. Actually, if one of the entropy modes leaves the horizon and becomes classical, the
combination of the remaining three variables are (sss) for pure entropy component and (sσσ) for mixed component,
respectively. As mentioned above, since these three-point correlations are absent at leading order in slow roll, the
trispectrum vanishes in this limit.
Similarly, we can discuss the consistency relation in the counter-collinear limit for the trispectra coming from the
scalar exchange diagram. According to [38], in the limit where the momentum of the scalar mode that is exchanged
goes to zero, one can find a simple relation between the scalar exchange trispectrum and the power spectrum, in
a similar way to Maldacena’s consistency relations [6] (see also [28, 101, 102]) since one can treat this mode as a
background again if this mode is an adiabatic mode. Let us suppose that the momentum of the exchanged adiabatic
particle is k12 and that k12 ≪ k1 ≈ k2 , k3 ≈ k4 . Then the mode associated with this scalar particle will cross the
horizon much before the other ki modes, where i = 1, . . . , 4, and it only rescale the spatial background where the
ki modes exist. Then, in the limit k12 → 0, the following relation should hold if the exchange mode is an adiabatic
mode:
hR(k1 )R(k2 )R(k3 )R(k4 )iSE → (2π)3 δ (3) (K)(nR − 1)2 P̃R (k12 )P̃R (k1 )P̃R (k3 ).
(50)
SE
The previous equation implies that in the counter-collinear limit, the scalar exchange non-linearity parameter τN
L is
2
SE
of order ǫ (at leading order in our approximations this is equivalent to say that τN L should vanish). From Eqs. (40),
4
(41) and (42), the non-linearity parameter τN L is shown to scale as O(k12
). Therefore, the consistency relation in
the counter-collinear limit holds for the trispectrum coming from the scalar exchange diagram generated from the
multi-field DBI inflation model. For the trispectrum coming from the other diagram, since there is no k12 dependence,
this consistency relation holds trivially.
For the pure adiabatic, pure entropy and part of the mixed component, the exchanged mode is always an adiabatic
mode. Thus we can apply the above argument to show that τN L vanishes in this limit. On the other hand, if the
mode with k12 is the entropy mode, which happens for the mixed component, we cannot apply the above arguments.
However, in this case, the four-point function becomes a square of the two-point function of (sσ). As was shown by
Eq. (17), there is no two-point correlation between the adiabatic and entropy mode at leading order in slow-roll. This
explains the consistency relation in the counter-collinear limit for the remaining part of the mixed component.
Now that we have checked the two important consistency relations, in the next subsection, we will concentrate on
a concrete and best studied configuration in single field DBI inflation model, the so-called equilateral configuration
[28, 31, 62, 63].
B.
The shape of the trispectrum and τNL in the equilateral configuration
In this subsection, we shall study the shape of the trispectrum in the limit where all the four momentum vectors have
the same magnitude k. If we denote the angle between ki and kj by θij , then in this configuration, we have cos(θ12 ) =
cos(θ34 ) ≡ cos(θ3 ), cos(θ23 ) = cos(θ14 ) ≡ cos(θ1 ), cos(θ13 ) = cos(θ24 ) ≡ cos(θ2 ) and cos(θ1 ) + cos(θ2 ) + cos(θ3 ) = −1
due to momentum conservation. Contrary to what happens in the case of the bispectrum, the equilateral configuration
conditions do not fix all degrees of freedom (dof) required to describe the shape of the trispectrum and we are left
with two angular dof [31].
In the equilateral configuration, the contact interaction trispectrum is constant and only the scalar exchange
trispectrum depends on the remaining angular dof, that we choose to be θ1 and θ2 . The plots of the shape of
the scalar exchange trispectrum coming from the purely adiabatic component (the first term in Eq. (43)) and the
purely entropy component (the last term in Eq. (43)) can be found in Fig. 1. Using Eq. (45), one can easily see that,
for TRS ≪ 1 the purely adiabatic component dominates while for TRS ≫ 1, the purely entropy component dominates.
Because the momentum dependence of the trispectrum changes depending on the value of TRS , the trispectrum can
distinguish multi-field DBI inflation models from single field DBI inflation models if TRS is significantly large. The
exact analytical expressions for the different contributions in the scalar exchange trispectrum can be found in Appendix
C.
Even though the trispectrum includes all the information, for practical purposes it is more convenient to use another
SE
quantity τN L defined in Eq. (48). In Fig. 2, we plot the non-linearity parameter τN
L for the purely adiabatic and
purely entropy components calculated from the respective terms in the scalar exchange trispectrum, Eq. (43).
CI
In Fig. 3, we plot the non-linearity parameter τN
L for the purely entropy component of the contact interaction
trispectrum. The shape of the surfaces for the purely adiabatic and mixed components are the same as the one for
the purely entropy component. However, the amplitudes are different, in particular, for the purely adiabatic plot,
SE
CI
CI
τN
L is negative. On the r.h.s. of Fig. 3, we plot the sum of the non-linearity parameters τN L and τN L that come
11
FIG. 1: The shape of the scalar exchange adiabatic (l.h.s panel) and pure entropy (r.h.s. panel) trispectra as functions of the
variables cos θ1 and cos θ2 . The amplitude of the trispectra at the point cos θ1 = cos θ2 = −1/3 has been normalized to unity.
SE
FIG. 2: Plots of the non-linearity parameter τNL
calculated from the scalar exchange adiabatic (l.h.s. panel) and pure entropy
2
(r.h.s.panel) trispectra as functions of the variables cos θ1 and cos θ2 . The amplitude of the first plot was rescaled by c4s (1+TRS
)3 ,
4
2
3
4
while the amplitude of the r.h.s. plot was rescaled by cs (1 + TRS ) /TRS .
from purely entropic perturbations. If the transfer coefficient is large this will be the dominant shape in the total
non-linearity parameter.
3
Notice that since the factor kij
is multiplying the trispectrum, in terms of τN L , we confirm the consistency relation
in the counter-collinear limit of the trispectrum coming from pure adiabatic, mixed and pure entropy components.
When cos θ1 = −1 or cos θ2 = −1 or cos θ3 = −1 (on the diagonal lines in the plots of Fig. 2, when cos θ2 = − cos θ1 )
we can see the non-linearity parameters are zero.
In order to quote a number for the value of the non-linearity parameter, we need to fix the remaining two degrees
of freedom (θ1 and θ2 ) and in this way fully specify the momentum vectors’ configuration. In the single field DBI
model, there are good reasons to believe that the maximum of τN L occurs for the equilateral configuration with
SE
cos(θ1 ) = cos(θ2 ) = cos(θ3 ) = −1/3 [62] 3 . In the multi-field case, for some contributions (like for example, τN
L
coming from the pure entropy trispectrum, as seen on the r.h.s. panel of Fig. 2) the non-linearity parameter is not
maximized for this particular angular configuration. However, for the sake of comparison with the single field model,
we shall obtain the non-linearity parameter in that configuration.
3
However, the definite proof that this is the case is still missing
12
CI
FIG. 3: L.h.s. panel: Plot of the non-linearity parameter τNL
calculated from the contact interaction pure entropy trispectrum
CI
SE
as a function of cos θ1 and cos θ2 . R.h.s. panel: Plot of the sum of τNL
and τNL
coming from the pure entropy trispectrum as
4
2
4
a function of cos θ1 and cos θ2 . The amplitude of the plots was rescaled by cs (1 + TRS
)3 /TRS
.
So, for the equilateral configuration, with cos(θ1 ) = cos(θ2 ) = cos(θ3 ) = −1/3, we obtain
τN L CI
4σ ∼ −
0.036
2 )3 ,
+ TRS
c4s (1
τN L CI
σσss ∼
2
0.016TRS
,
2
+ TRS )3
c4s (1
τN L CI
4s ∼
4
0.052TRS
,
2
+ TRS )3
c4s (1
(51)
where the different values of τN L correspond to the three terms in the contact interaction trispectrum coming from
the purely adiabatic, mixed and purely entropic components in Eq. (36), respectively.
Similarly, from the scalar exchange trispectrum we obtain
τN L SE
4σ ∼
0.60
2 )3 ,
4
cs (1 + TRS
τN L SE
σσss ∼
2
0.67TRS
2 )3 ,
c4s (1 + TRS
τN L SE
4s ∼
4
0.071TRS
2 )3 ,
c4s (1 + TRS
(52)
where again the different values of τN L correspond to the three terms in the scalar exchange interaction trispectrum
coming from the purely adiabatic, mixed and purely entropic components in Eq. (43), respectively.
It is worth mentioning that from Eqs. (37) and (44), both for the contact interaction (Eq. (51)) and the scalar
exchange interaction (Eq. (52)), we can obtain the relation
2
τN L4σ +
τN Lσσss = TRS
1
2
TRS
τN L4s .
(53)
Therefore, we conclude that the total non-linearity parameter τN L for the multi-field DBI inflation model in the
equilateral configuration with equal angles between the momentum vectors is
τN L tot
equi ∼
2
4
0.68TRS
0.12TRS
0.56
+
+
2 )3
2 )3
2 )3 .
c4s (1 + TRS
c4s (1 + TRS
c4s (1 + TRS
(54)
It is worth noting that when there is no conversion of entropy perturbation into curvature perturbation, i.e. TRS = 0,
Eq. (54) reproduces the result obtained in the single field DBI inflation model [62]. On the other hand, when the
transfer coefficient is large, i.e. TRS ≫ 1, the final curvature perturbation is mainly of entropy origin and the last
term in Eq. (54) dominates over the others.
C.
The non-linearity parameter τNL in a non-equilateral configuration
As we have mentioned before, the equilateral configuration is well investigated in the DBI inflation model. However,
since there are five degrees of freedom to parameterize the shape of the trispectrum and the expression for the
trispectrum is fairly complicated, non-equilateral configurations have attracted less attention. Because it has not
been shown that the equilateral configuration maximizes the value of the non-linearity parameter, in this subsection
we will consider a non-equilateral configuration.
In our previous work [62], we have studied the deviation from the equilateral configuration by fixing a particular
set of degrees of freedom in the single field DBI inflation model case. Here, we apply the same approach and we study
the behaviour of τN L in a non-equilateral configuration.
13
SE
FIG. 4: Plots of the non-linearity parameter τNL
calculated from the scalar exchange adiabatic (l.h.s. panel) and pure entropy
(r.h.s. panel) trispectra for non-equilateral configurations as functions of the momentum amplitude q and cos θ2 . The amplitude
2
2
4
of the first plot was rescaled by c4s (1 + TRS
)3 , while the amplitude of the r.h.s. plot was rescaled by c4s (1 + TRS
)3 /TRS
.
CI
FIG. 5: Plots of the non-linearity parameter τNL
calculated from the contact interaction adiabatic (l.h.s. panel) and pure
entropy (r.h.s. panel) trispectra for non-equilateral configurations as functions of the momentum amplitude q and cos θ2 . The
2
2
4
amplitude of the first plot was rescaled by c4s (1+TRS
)3 , while the amplitude of the r.h.s. plot was rescaled by c4s (1+TRS
)3 /TRS
.
We consider configurations where the two momentum vectors k3 and k4 have a different magnitude q from the
magnitude of k1 and k2 that we normalized to unity. Then in this configuration, we have cos(θ34 ) ≡ cos(θ3 ), cos(θ23 ) =
cos(θ14 ) ≡ cos(θ1 ), cos(θ13 ) = cos(θ24 ) ≡ cos(θ2 ), cos(θ12 ) = q 2 (1 + cos(θ3 )) − 1. Momentum conservation implies
cos θ3 = −(1/q)(cos θ1 + cos θ2 ) − 1. We will allow one of the angular dof to vary freely, the other angular dof θ1 is
fixed by cos θ1 = −1/3 to the value that gives us the maximum for the single field model. In Fig. 4, we plot the
SE
non-linearity parameter τN
L calculated from the scalar exchange adiabatic (l.h.s. panel) and entropic (r.h.s. panel)
trispectra, respectively. The maximum appears at q = 1 and cos(θ2 ) = −1/3, which supports the importance of the
CI
equilateral configuration. In Fig. 5, we plot the non-linearity parameter τN
L calculated from the contact interaction
adiabatic (l.h.s. panel) and entropic (r.h.s. panel) trispectra, respectively. Comparing Figs. 4 and 5, one can see
that around the maximum the purely adiabatic contribution coming from the contact interaction is smaller than the
corresponding one coming from the scalar exchange interaction. Only for the entropy component, the amplitudes of
the scalar exchange and contact interaction trispectra are comparable. We plot the sum of these two terms in Fig. 6.
If the transfer coefficient is large this will be the dominant shape.
VI.
OBSERVATIONAL CONSTRAINTS
Since we have obtained the full quantum trispectra for each component for the equilateral configuration with
cos(θ1 ) = cos(θ2 ) = cos(θ3 ) = −1/3, in terms of TRS , we are in a position to discuss the relation between TRS and
observable quantities. It was shown in [57] that in the small sound speed limit (c2s ≪ 1), the spectral index of the
14
SE
CI
FIG. 6: Plot of the sum of the non-linearity parameters τNL
and τNL
calculated from the pure entropy trispectra, for nonequilateral configurations, as functions of the momentum amplitude q and cos θ2 . The amplitude of the plot was rescaled by
2
4
c4s (1 + TRS
)3 /TRS
.
scalar perturbations nR in two-field DBI inflation model can be written as
q
equil
3.1|fN
f˙
L |r
−
+ α∗ sin(2Θ) + 2β∗ sin2 Θ,
1 − nR ≃
3
4 cos Θ
Hf
ξ2
χ
ι
1
ξ
2
µs + 2 , Θ = arcsin
, β= − −
with α =
aH
2
2 3H 2
a
(55)
T
p RS
2
1 + TRS
!
,
(56)
where ι and χ are the slow-roll parameters defined in Eq. (9). ξ is the coupling between the adiabatic and entropy
perturbation and µs is the mass of the entropy perturbation, which are model dependent quantities. In order to
obtain Eq. (55), they used the following relations:
PT
= 16ǫcs cos2 Θ,
PR
1
equi
cos2 Θ,
fN
L ≃−
3.1c2s
r=
(57)
(58)
equi
where r is the tensor to scalar ratio and fN
L parameterizes the size of the bispectrum in the equilateral configuration.
˙
In the ultra-violet (UV) DBI model (f > 0) and without entropy transfer, i.e. TRS = 0 (or equivalently sin Θ = 0),
Eq. (55) gives the well known stringent constraint as
4
r> q
(1 − nR ),
equi
3.1|fN
|
L
(59)
equi
by which one obtains r ≥ 10−3 if one uses the observational constraints on fN
L and 1 − nR from WMAP5 data
[3]. On the other hand, since for general compact internal spaces, suggested by string theory, there is a geometrical
constraint that limits r to be r < 10−7 [77], the single field UV DBI inflation model is already excluded by the
observational data.
However, if we take into account the transfer from the entropy perturbation, the last two terms in Eq. (55) are no
equi
2
longer negligible. Furthermore, from Eqs. (57) and (58), both r and fN
L are suppressed by cos Θ which scales as
2
1/TRS for TRS ≫ 1. Because of these factors, the constraint given by Eq. (59) cannot be applied and it is possible
to satisfy the geometrical constraint r < 10−7 if TRS ≫ 1 [57]. In such a case, since the last term of Eq. (54) will
equi
equi
equi 2 2
equi
equi
dominate over the other terms, τN
L can be expressed as τN L ∼ 1.14(fN L ) TRS . Thus for a given fN L , τN L is
2
enhanced for a large transfer TRS ≫ 1, while the global magnitude of τN L is reduced, going like 1/TRS
in this case.
VII.
CONCLUSIONS
In this work, we have obtained the full quantum trispectrum of the primordial curvature perturbation at leading
order in the slow-roll expansion and in the small sound speed limit in multi-field DBI inflation models. Previously,
15
we had obtained the contact interaction trispectrum [64]. In this work we completed the calculation of the quantum
trispectrum including the scalar exchange trispectrum.
We studied the momentum dependence of the trispectrum by separating it into the purely adiabatic, mixed and
purely entropic components. For general configurations, as in the single field model, we confirmed that the consistency
relations in the squeezed and in the counter-collinear limit hold. If a mode that crosses the horizon much before the
other modes is an adiabatic one, these can be explained by the usual argument that the mode only changes the
background (just like in the single field models). However, if the mode that crosses the horizon first is an entropic
one, the consistency relations hold because of the absence of the two-point and three-point correlations among the
modes remaining under the horizon scale. Essentially, the absence of these correlations is attributed to the particular
properties of the DBI action discussed in [64].
In the equilateral configuration, we examined the momentum dependence of the scalar exchange trispectra as
functions of the two remaining angular degrees of freedom. We found that the shape of the adiabatic, mixed, and
purely entropic trispectra are different. This means that if the amount of the transfer from the entropy perturbation to
the final curvature perturbation is significantly large, the trispectrum can distinguish multi-field DBI inflation models
from single field DBI inflation models. In the equilateral configuration with cos θ1 = cos θ2 = cos θ3 = −1/3, we showed
that for the contributions for τN L coming from the adiabatic and mixed trispectra, the scalar exchange diagram is
more important than the contact interaction diagram, while for the pure entropy contribution both diagrams are
equally important. This suggests that for quantitative discussions, the analysis done in Ref. [55] taking into account
only the scalar exchange diagram is insufficient, at least, in this configuration. We have also considered configurations
where the equilateral conditions are not satisfied.
After obtaining the complete predictions for the quantum trispectrum in the multi-field DBI inflation model, we
discussed the observational constraints on the non-linearity parameter, for the equilateral configuration mentioned
above. In single field DBI inflation, some models based on string theory set-ups suffer from a severe constraint related
with the geometrical structure of the compact internal spaces. However, in the presence of the transfer from the
entropy mode, this constraint no longer holds. Furthermore, since both the tensor-to-scalar ratio and non-linearity
−2
for a large transfer from the entropy mode to the curvature
parameter for the bispectrum are suppressed like TRS
perturbation, it is easier to satisfy the corresponding constraints [57]. On the other hand, the trispectrum is enhanced
equi 2
equi
equi
2
for a given amplitude of the bispectrum fN
like τN
L ∝ TRS fN L
L , while the global magnitude of τN L is reduced,
2
going like 1/TRS in this case.
In order to obtain a concrete value for the transfer coefficient TRS , we need a specific model. It would be interesting
to study this mixing in specific string theory motivated models [103].
Note added: On the day this work appeared in the arXiv, the paper [104] was also submitted to the arXiv, it
also calculates the trispectrum in multi-field DBI-inflation. While we have concentrated on the trispectrum coming
from the intrinsically quantum four-point functions of the fields which are generated during inflation, Ref. [104], by
assuming a specific model, showed that the contribution to the trispectrum coming from the super-horizon non-linear
dynamics is also important.
Acknowledgments
We would like to thank Takahiro Tanaka for interesting discussions. SM and FA are supported by JSPS Research
Fellowships. KK is supported by ERC, RCUK and STFC. SM is grateful to the ICG, Portsmouth for their hospitality
when this work was completed.
APPENDIX A: THE DIAGRAMMATIC APPROACH TO THE TRISPECTRUM FROM THE SCALAR
EXCHANGE INTERACTION
The calculation of the integrals in Eq. (39) is rather long and it is easy to forget some terms or permutations. In this
appendix, we present a systematic and simple way to keep track of all the different contributions and permutations
to the integrals in (39). This approach is based on diagrams and simple rules that allows one to read the respective
contributions for the scalar exchange four-point function directly from the diagrams.
To calculate the scalar exchange four-point function we need to evaluate
hΩ|Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 )|ΩiSE
Z 0
Z η
i
i
hh
int
int
(η̃) |0i,
(η) , H(3)
=−
dη
dη̃h0| Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 ), H(3)
−∞
−∞
(A1)
16
int
where on the l.h.s. of the equation Qm are the interaction picture fields and H(3)
is the third-order interaction
Hamiltonian that can be found in Eq. (20) or
int
H(3)
(η) = H1int (η) + H2int (η) + H3int (η) + H4int (η) + H5int (η),
where the different type of interactions are
Z
Z
Z
2
3
H1int (η) ≡
d3 xg1 aQ′σ , H2int (η) ≡ d3 xg2 aQ′σ (∂Qσ )2 , H3int (η) ≡ d3 xg3 aQ′σ Q′s ,
Z
Z
H4int (η) ≡
d3 xg4 aQ′σ (∂Qs )2 , H5int (η) ≡ d3 xg5 aQ′s (∂Qσ )(∂Qs ),
(A2)
(A3)
where here Qσ and Qs are the interaction picture fields. The coupling constants are
g1 = M,
1
.
8X0 c7s
where M = − √
g2 = −c2s M,
g3 = M,
g4 = c2s M,
g5 = −2c2s M,
(A4)
The integrand of Eq. (A1) can be written as
i
i
hh
int
int
(η̃)
(η) , H(3)
Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 ), H(3)
5
X
int
int
int
int
int
int
int
int
F Hi (η)Hj (η̃) + Hj (η̃)Hi (η)F − Hj (η̃)F Hi (η) − Hi (η)F Hj (η̃)
=
(A5)
i,j=1
where we have defined F as F ≡ Qm (0, k1 )Qn (0, k2 )Qp (0, k3 )Qq (0, k4 ).
One can immediately see that at the leading order the only non-zero four-point functions
are
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiSE ,
hΩ|Qs (0, k1 )Qs (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiSE
and
SE
hΩ|Qs (0, k1 )Qs (0, k2 )Qs (0, k3 )Qs (0, k4 )|Ωi . The four-point functions hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qs (0, k4 )|ΩiSE
and hΩ|Qσ (0, k1 )Qs (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiSE are zero at this order. This is a direct consequence of the absence
of the interactions Q3s and Q2σ Qs in the third order interaction Hamiltonian (A2). In terms of diagrams this fact is
clear. With the interactions (A3) one cannot draw diagrams with one and three external entropy legs, respectively.
In order to find the non-zero four-point functions we need the following set of rules. We draw the allowed diagrams
with two of the vertices (A3), where one leg of the first vertex is glued with one leg of the other vertex (entropy legs
can only be glued with entropy legs, similarly for adiabatic legs). The resulting four-point function can be constructed
as:
• Multiply by −2(2π)3 δ (3) (K)N 4 /(k1 k2 k3 k4 )3/2 , where N is the normalization of the mode functions. Multiply
by the coupling constants of the vertices that enter the diagram. The overall minus sign is already present in
(A1). The Dirac delta function ensures momentum conservation. The factor N 4 /(k1 k2 k3 k4 )3/2 comes from the
normalization of the four external legs, described by F . Finally, the factor of two can be understood if one notes
that the contribution to the four-point function coming from the first term in the r.h.s. of Eq. (A5) is real and
the operator in the second term is just the Hermitian conjugate of the first term. Similarly with the third and
fourth terms in the r.h.s. of Eq. (A5).
• To each leg of the vertex with a spatial derivative associate the usual factor −ik.
• The different ways of drawing the same-looking diagram (ignoring the k labels but with the order in the leftto-right of the vertices the same) give a multiplicative symmetry factor.
• The contribution coming from the third (and fourth) term in the r.h.s. of Eq. (A5) is equal to the contribution
coming from the first (and second) term but with the sign of k1 and k2 (but not k12 ), that enter the arguments
of Fi functions, changed to minus sign.
• In the diagrams with vertices H1int or H3int the form factor is made of F1 functions. The definitions of the Fi
functions can be found in Appendix B. In the diagrams with vertices H2int , H4int or H5int the form factor is
made of F2 functions. In the diagrams where the l.h.s. vertex is either one of H1int or H3int and the r.h.s. vertex
is either one of H2int , H4int or H5int , the form factor is made of F3 functions. In the diagrams where the l.h.s.
vertex is either one of H2int , H4int or H5int and the r.h.s. vertex is either one of H1int or H3int , the form factor
is made of F4 functions. The order of the six arguments in Fi is read from the diagrams too. The first three
arguments are read from the l.h.s. vertex and the last three arguments from the r.h.s vertex. For example, for
the diagram made of two H3int vertices, the first and fourth arguments of F2 are the legs with time derivatives.
k1 in the first vertex and k4 in the second vertex.
17
.
σ
σ
k1
_
k 12
k3
k 12
.
σ
k4
k2
.
.
σ
σ
σ
FIG. 7: The two leading order vertices that contribute for the scalar exchange trispectrum in the single field case. Dot over
σ denotes that in that leg the interaction is a time derivative. ∇σ means that the interaction is a spatial derivative. There is
momentum conservation at the vertices, hence at each vertex the momentum vectors have to sum to zero. Overall momentum
conservation implies that k12 = −k34 , hence the labels in the r.h.s. diagram.
.
σ
k1
_
k 12
.
σ
.
.
σ
k3
k2
σ
.
σ
_
k 12
k 12
.
σ
σ
k4
σ
.
σ
k2
σ
k3
k1
.
k 12
σ
k4
σ
FIG. 8: The two diagrams that contribute to the adiabatic four-point function, where the l.h.s. vertex is H1int and the r.h.s.
vertex is H2int . The first, on the l.h.s. panel, we glue two legs with a time and a spatial derivative. The second, on the r.h.s.
panel, we glue two legs with time derivatives.
• At the end, add twenty three permutations of {k1 , k2 , k3 , k4 } to restore the symmetry of the four-point functions
under the exchange of the momentum vectors.
Let us consider an explicit example and calculate the contribution for the adiabatic four-point function of a diagram
where the l.h.s. vertex is H1int and the r.h.s. vertex is H2int . These vertices are schematically depicted in Fig. 7 and
they are the same ones as in the single field case.
There are two distinct diagrams that one can draw, as can be seen in Fig. 8, hence there will be two terms. For
the diagram on the l.h.s. of Fig. 8 we get a symmetry factor six. This is because there are three ways we can glue the
∇σ leg to the σ̇ leg and there are two ∇σ legs. For the diagram on the r.h.s of Fig. 8 we get a symmetry factor three.
This is because there are three ways of gluing the σ̇ leg of the H2int vertex with one of the σ̇ legs of the H1int vertex.
In both diagrams, we have two spatial derivatives so we get a factor (−i)2 . These diagrams will be proportional to
F3 .
Finally, the diagram on the l.h.s. gives
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiSE
−2N 4 g1 g2
⊃ (−i)2 6(2π)3 δ(K)
3 (k12 · k4 ) (F3 (k1 , k2 , −k12 , k3 , k4 , k12 ) − F3 (−k1 , −k2 , −k12 , k3 , k4 , k12 ))
(k1 k2 k3 k4 ) 2
+23 permutations of{k1 , k2 , k3 , k4 },
(A6)
and the diagram in the r.h.s. gives
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qσ (0, k3 )Qσ (0, k4 )|ΩiSE
−2N 4 g1 g2
⊃ (−i)2 3(2π)3 δ(K)
3 (k3 · k4 ) (F3 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F3 (−k1 , −k2 , −k12 , k12 , k3 , k4 ))
(k1 k2 k3 k4 ) 2
+23 permutations of{k1 , k2 , k3 , k4 }
(A7)
Now, let us consider another example and calculate the contribution for the mixed four-point function of a diagram
where the l.h.s. vertex is H4int and the r.h.s. vertex is H5int . These diagrams involves vertices that contain entropy
legs, they are schematically depicted in Fig. 9.
Because the entropy and the adiabatic mode functions are the same, if one substitutes s with σ in the vertices of
Fig. 9 we get the vertices of Fig. 7. Because the Fi functions only depend on the mode functions, that are the same
18
.
.
σ
σ
k1
k1
_
g3
k 12
.s
_
g4
k2
k 12
s
k2
.s
s
.s
k1
_
g5
k 12
s
k2
σ
FIG. 9: The leading order vertices containing entropy modes that contribute for the scalar exchange trispectrum in the multifield case. gi denotes the coupling constants. As before, the dot over s denotes that in that leg the interaction is a time
derivative. ∇s means that the interaction is a spatial derivative. A label s in the leg means that it is a entropy mode. σ
denotes an adiabatic leg.
.
σ
g
4
s
k3
k1
_
k 12
s
k2
s
.
s
k 12
g
.
σ
g
5
_
4
k2
k4
σ
s
σ
k3
k1
k 12
k 12
s
s
g
5
k4
.
s
FIG. 10: The two diagrams that contribute to the mixed four-point function, where the l.h.s. vertex is H4int and the r.h.s.
vertex is H5int .
for the entropy and the adiabatic modes we don’t need to define other Fi functions, and the single field ones are
enough.
There are two diagrams that contribute in this example, as can be seen in Fig. 10, hence there will be two terms.
In both diagrams, we have four spatial derivatives, this gives a factor of (−i)4 . The diagram on the l.h.s. of the figure
has a symmetry factor two, because there are two ways of gluing the ṡ leg of the r.h.s. vertex with one of the two ∇s
legs of the l.h.s. vertex. Similarly, the symmetry factor of the r.h.s. diagram is also two.
The l.h.s. diagram contributes with
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiSE + 5 perms.
−2N 4 g4 g5
⊃ (−i)4 2(2π)3 δ(K)
3 (k2 · (−k12 )) (k3 · k4 ) (F2 (k1 , k2 , −k12 , k12 , k3 , k4 ) − F2 (−k1 , −k2 , −k12 , k12 , k3 , k4 ))
(k1 k2 k3 k4 ) 2
+23 permutations of{k1 , k2 , k3 , k4 },
(A8)
and the diagram in the r.h.s. gives
hΩ|Qσ (0, k1 )Qσ (0, k2 )Qs (0, k3 )Qs (0, k4 )|ΩiSE + 5 perms.
−2N 4 g4 g5
⊃ (−i)4 2(2π)3 δ(K)
3 (k2 · (−k12 )) (k12 · k4 ) (F2 (k1 , k2 , −k12 , k4 , k3 , k12 ) − F2 (−k1 , −k2 , −k12 , k4 , k3 , k12 ))
(k1 k2 k3 k4 ) 2
+23 permutations of{k1 , k2 , k3 , k4 },
(A9)
where 5 perms. denotes the following permutations of the displayed term [1 ↔ 3], [1 ↔ 4], [2 ↔ 3], [2 ↔ 4] and
[1 ↔ 3, 2 ↔ 4].
19
APPENDIX B: DEFINITIONS OF THE Fi FUNCTIONS
In this Appendix, we present the definitions and the analytical expressions of the four Fi functions (with i = 1, . . . , 4)
that appear in the main text.
Because there are only four different types of double time integrals over the mode functions, we define four Fi
functions (with i = 1, . . . , 4). Their explicit expressions are
0
Z
η
Z
′
′
′
′
′
′
dη̃a(η̃)U ∗ (η, k1 )U ∗ (η, k2 )U ∗ (η, k3 )U ∗ (η̃, k4 )U ∗ (η̃, k5 )U ∗ (η̃, k6 )
−∞
−∞
1
1
A2
N 6 c6s
A
= −4 2 |k1 · · · k6 | 2 3 3 1 + 3 + 6 2 ,
(B1)
H
A C
C
C
√
where N is the normalisation of the mode functions given by N = H/ 2cs , A is defined by the sum of the last
three arguments of the Fi functions as A = k4 + k5 + k6 and C is defined by the sum of all the arguments as
C = k1 + k2 + k3 + k4 + k5 + k6 . The remaining functions are
Z 0
Z η
′
′
F2 (k1 , k2 , k3 , k4 , k5 , k6 ) =
dηa(η)
dη̃a(η̃)U ∗ (η, k1 )U ∗ (η, k2 )U ∗ (η, k3 )U ∗ (η̃, k4 )U ∗ (η̃, k5 )U ∗ (η̃, k6 )
F1 (k1 , k2 , k3 , k4 , k5 , k6 ) =
dηa(η)
−∞
−∞
1
N 6 c2s |k1 k4 | 2
k5 k6
k5 + k6
1
= − 2
1+
+2 2
3
H |k2 k3 k5 k6 | 2 AC
A
A
1
1
k5 k6 (k2 + k3 )
+
k2 + k3 + k5 + k6 + ((k2 + k3 ) (k5 + k6 ) + 2k5 k6 ) + 2
C
A
A2
2
+ 2 k5 k6 + (k2 + k3 ) (k5 + k6 ) + k2 k3
C
1
k2 k3 k5 k6
+ (k2 k3 (k5 + k6 ) + 2k5 k6 (k2 + k3 )) + 2
A
A2
k2 k3 k5 k6
k2 k3 k5 k6
6
+ 24
,
+ 3 k2 k3 (k5 + k6 ) + k5 k6 (k2 + k3 ) + 2
C
A
C4
F3 (k1 , k2 , k3 , k4 , k5 , k6 ) =
=
F4 (k1 , k2 , k3 , k4 , k5 , k6 ) =
Z
0
dηa(η)
−∞
N 6 c4
2 2s
H
Z
Z
′
′
′
−∞
1
3
|k5 k6 | 2
dηa(η)
′
dη̃a(η̃)U ∗ (η, k1 )U ∗ (η, k2 )U ∗ (η, k3 )U ∗ (η̃, k4 )U ∗ (η̃, k5 )U ∗ (η̃, k6 )
|k1 k2 k3 k4 | 2
0
−∞
η
(B2)
Z
k5 k6
k5 k6
k5 + k6
3
k5 k6
1
k
+
k
+
2
+
12
1
+
,
+
+
2
5
6
AC 3
A
A2
C
A
C2
(B3)
η
′
′
′
′
dη̃a(η̃)U ∗ (η, k1 )U ∗ (η, k2 )U ∗ (η, k3 )U ∗ (η̃, k4 )U ∗ (η̃, k5 )U ∗ (η̃, k6 )
−∞
1
N 6 c4 |k1 k4 k5 k6 | 2 1
A A2
k2 + k3
A (k2 + k3 ) + k2 k3
= 2 2s
1+ + 2 +
+2
3
3
H
C
C
C
C2
|k2 k3 | 2 A C
A
A2
+3 3 (A (k2 + k3 ) + 2k2 k3 ) + 12k2 k3 4 ,
C
C
(B4)
where the function U (η, k) is a modified mode function given by
U (η, k) = N
1
3
|k| 2
(1 + ikcs η) e−ikcs η .
(B5)
If the sign of the argument k of U is positive then U is equal to the mode function, if the sign is negative then U
equals the complex conjugate of the mode function.
20
APPENDIX C: THE SCALAR EXCHANGE TRISPECTRUM IN THE EQUILATERAL LIMIT
In this Appendix, we present the scalar exchange trispectrum of the primordial curvature perturbation R in the
equilateral configuration. These analytical expressions can be derived from the general expression for the scalar
exchange trispectrum, Eq. (43). We will use these simpler equations to study the shape dependence of the trispectrum
and to calculate the non-linearity parameter τN L in subsection V B.
At leading order in the small sound speed and in the slow-roll expansion, the scalar exchange trispectrum is
hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE = hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE 4σ
+hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE σσss
+hΩ|R(0, Rk1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE 4s ,
(C1)
where the three different contributions come from the terms in Eq. (43) that are coming from the purely adiabatic,
mixed and purely entropy components.
Explicitly, they are given by
hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE 4σ = (2π)3 δ (3) (K)
1
1
2
2 )3 D1 (k12 /k)P̃R∗ (k)P̃R∗ (k12 ) + 2 perms,
c4s (1 + TRS
(C2)
hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE σσss = (2π)3 δ (3) (K)
hΩ|R(0, k1 )R(0, k2 )R(0, k3 )R(0, k4 )|ΩiSE 4s = (2π)3 δ (3) (K)
2
1
TRS
2
D2 (k12 /k)P̃R∗
(k)P̃R∗ (k12 ) + 2 perms,
2
3
(1 + TRS ) c4s
(C3)
4
TRS
1
2
2 )3 c4 D3 (k12 /k)P̃R∗ (k)P̃R∗ (k12 ) + 2 perms,
(1 + TRS
s
(C4)
where we have defined three Di (x) functions, with i = 1, 2, 3, as
x4 (61952 + 48384x − 10240x2 − 20832x3 − 5184x4 + 1412x5 + 824x6 + 103x7 )
,
29 (x + 2)4
x4 (31232 + 23744x − 6272x2 − 10496x3 − 1248x4 + 1904x5 + 824x6 + 103x7 )
,
D2 (x) =
28 (x + 2)4
x4 (512 − 896x − 2304x2 − 160x3 + 2688x4 + 2396x5 + 824x6 + 103x7 )
D3 (x) =
,
(C5)
29 (x + 2)4
p
p
p
where k12 = k 2(1 + cos θ3 ), k13 = k 2(1 + cos θ2 ) and k14 = k 2(1 + cos θ1 ). k denotes the common amplitude
of all four momentum vectors kj , with j = 1, . . . , 4 and the rescaled power spectrum is P̃R∗ (k) = H 2 /(4ǫcs k 3 ) that
should be evaluated at sound horizon crossing.
It is worth mentioning that from Eqs. (37) and (44), the following relation holds,
D1 (x) =
D2 = D1 + D3 .
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
http://map.gsfc.nasa.gov/.
A. P. S. Yadav and B. D. Wandelt, Phys. Rev. Lett. 100, 181301 (2008), 0712.1148.
WMAP, E. Komatsu et al., Astrophys. J. Suppl. 180, 330 (2009), 0803.0547.
K. M. Smith, L. Senatore and M. Zaldarriaga, JCAP 0909 (2009) 006 [arXiv:0901.2572 [astro-ph]].
http://www.rssd.esa.int/index.php?project=Planck
J. M. Maldacena, JHEP 0305 (2003) 013 [arXiv:astro-ph/0210603].
A. D. Linde and V. F. Mukhanov, Phys. Rev. D56, R535 (1997), astro-ph/9610219.
N. Bartolo, S. Matarrese, and A. Riotto, Phys. Rev. D65, 103505 (2002), hep-ph/0112261.
(C6)
21
[9]
[10]
[11]
[12]
[13]
[14]
[15]
[16]
[17]
[18]
[19]
[20]
[21]
[22]
[23]
[24]
[25]
[26]
[27]
[28]
[29]
[30]
[31]
[32]
[33]
[34]
[35]
[36]
[37]
[38]
[39]
[40]
[41]
[42]
[43]
[44]
[45]
[46]
[47]
[48]
[49]
[50]
[51]
[52]
[53]
[54]
[55]
[56]
[57]
[58]
[59]
[60]
[61]
[62]
[63]
[64]
[65]
[66]
[67]
[68]
[69]
[70]
[71]
F. Bernardeau and J.-P. Uzan, Phys. Rev. D66, 103506 (2002), hep-ph/0207295.
F. Bernardeau and J.-P. Uzan, Phys. Rev. D67, 121301 (2003), astro-ph/0209330.
G. Dvali, A. Gruzinov, and M. Zaldarriaga, Phys. Rev. D69, 023505 (2004), astro-ph/0303591.
P. Creminelli, JCAP 0310, 003 (2003), astro-ph/0306122.
M. Alishahiha, E. Silverstein and D. Tong, Phys. Rev. D 70 (2004) 123505 [arXiv:hep-th/0404084].
A. Gruzinov, Phys. Rev. D71, 027301 (2005), astro-ph/0406129.
N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) [arXiv:astro-ph/0406398].
K. Enqvist, A. Jokinen, A. Mazumdar, T. Multamaki, and A. Vaihkonen, Phys. Rev. Lett. 94, 161301 (2005),
astro-ph/0411394.
D. Seery and J. E. Lidsey, JCAP 0506 (2005) 003 [arXiv:astro-ph/0503692].
D. Seery and J. E. Lidsey, JCAP 0509 (2005) 011 [arXiv:astro-ph/0506056].
A. Jokinen and A. Mazumdar, JCAP 0604, 003 (2006), astro-ph/0512368.
D. H. Lyth, JCAP 0511, 006 (2005), astro-ph/0510443.
M. P. Salem, Phys. Rev. D72, 123516 (2005), astro-ph/0511146.
D. Seery and J. E. Lidsey, JCAP 0701, 008 (2007), astro-ph/0611034.
M. Sasaki, J. Valiviita, and D. Wands, Phys. Rev. D74, 103003 (2006), astro-ph/0607627.
K. A. Malik and D. H. Lyth, JCAP 0609, 008 (2006), astro-ph/0604387.
N. Barnaby and J. M. Cline, Phys. Rev. D73, 106012 (2006), astro-ph/0601481.
L. Alabidi and D. Lyth, JCAP 0608, 006 (2006), astro-ph/0604569.
X. Chen, M. x. Huang, S. Kachru and G. Shiu, JCAP 0701 (2007) 002 [arXiv:hep-th/0605045].
X. Chen, M. x. Huang and G. Shiu, Phys. Rev. D 74 (2006) 121301 [arXiv:hep-th/0610235].
X. Chen, R. Easther, and E. A. Lim, JCAP 0706, 023 (2007), astro-ph/0611645.
L. Alabidi, JCAP 0610, 015 (2006), astro-ph/0604611.
D. Seery, J. E. Lidsey and M. S. Sloth, JCAP 0701, 027 (2007) [arXiv:astro-ph/0610210].
C. T. Byrnes, M. Sasaki, and D. Wands, Phys. Rev. D74, 123519 (2006), astro-ph/0611075.
T. Suyama and M. Yamaguchi, Phys. Rev. D77, 023505 (2008), 0709.2545.
F. Arroja and K. Koyama, Phys. Rev. D 77, 083517 (2008) [arXiv:0802.1167 [hep-th]].
F. Arroja, S. Mizuno and K. Koyama, JCAP 0808, 015 (2008) [arXiv:0806.0619 [astro-ph]].
D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, Phys. Rev. Lett. 101, 061301 (2008) [arXiv:0804.3139 [hep-th]].
D. Langlois, S. Renaux-Petel, D. A. Steer and T. Tanaka, Phys. Rev. D 78, 063523 (2008) [arXiv:0806.0336 [hep-th]].
D. Seery, M. S. Sloth and F. Vernizzi, JCAP 0903, 018 (2009) [arXiv:0811.3934 [astro-ph]].
M. Sasaki, Prog. Theor. Phys. 120, 159 (2008), 0805.0974.
C. T. Byrnes, K.-Y. Choi, and L. M. H. Hall, JCAP 0810, 008 (2008), 0807.1101.
C. T. Byrnes, K.-Y. Choi, and L. M. H. Hall, JCAP 0902, 017 (2009), 0812.0807.
B. Dutta, L. Leblond, and J. Kumar, Phys. Rev. D78, 083522 (2008), 0805.1229.
A. Naruko and M. Sasaki, Prog. Theor. Phys. 121, 193 (2009), 0807.0180.
T. Suyama and F. Takahashi, JCAP 0809, 007 (2008), 0804.0425.
X. Gao, JCAP 0806, 029 (2008) [arXiv:0804.1055 [astro-ph]].
H. R. S. Cogollo, Y. Rodriguez, and C. A. Valenzuela-Toledo, JCAP 0808, 029 (2008), 0806.1546.
Y. Rodriguez and C. A. Valenzuela-Toledo, (2008), 0811.4092.
K. Ichikawa, T. Suyama, T. Takahashi, and M. Yamaguchi, Phys. Rev. D78, 023513 (2008), 0802.4138.
C. T. Byrnes, JCAP 0901, 011 (2009), 0810.3913.
S. Li, Y.-F. Cai, and Y.-S. Piao, Phys. Lett. B671, 423 (2009), 0806.2363.
D. Langlois, F. Vernizzi, and D. Wands, JCAP 0812, 004 (2008), 0809.4646.
C. Hikage, K. Koyama, T. Matsubara, T. Takahashi, and M. Yamaguchi, (2008), 0812.3500.
M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama, and F. Takahashi, JCAP 0811, 019 (2008), 0808.0009.
Q. G. Huang, JCAP 0811 (2008) 005 [arXiv:0808.1793 [hep-th]].
X. Gao and B. Hu, JCAP 0908, 012 (2009) [arXiv:0903.1920 [astro-ph.CO]].
Y. F. Cai and H. Y. Xia, Phys. Lett. B 677 (2009) 226 [arXiv:0904.0062 [hep-th]].
D. Langlois, S. Renaux-Petel and D. A. Steer, JCAP 0904, 021 (2009) [arXiv:0902.2941 [hep-th]].
X. Gao, (2009), 0904.4187.
Q. G. Huang, JCAP 0905, 005 (2009) [arXiv:0903.1542 [hep-th]].
Q. G. Huang, JCAP 0906, 035 (2009) [arXiv:0904.2649 [hep-th]].
J. Khoury and F. Piazza, JCAP 0907, 026 (2009) [arXiv:0811.3633 [hep-th]].
F. Arroja, S. Mizuno, K. Koyama and T. Tanaka, Phys. Rev. D 80, 043527 (2009) [arXiv:0905.3641 [hep-th]].
X. Chen, B. Hu, M. x. Huang, G. Shiu and Y. Wang, JCAP 0908, 008 (2009) [arXiv:0905.3494 [astro-ph.CO]].
S. Mizuno, F. Arroja, K. Koyama and T. Tanaka, Phys. Rev. D 80, 023530 (2009) [arXiv:0905.4557 [hep-th]].
C. T. Byrnes and G. Tasinato, JCAP 0908, 016 (2009) [arXiv:0906.0767 [astro-ph.CO]].
X. Gao, M. Li and C. Lin, arXiv:0906.1345 [astro-ph.CO].
C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. B 458 (1999) 209 [arXiv:hep-th/9904075].
J. Garriga and V. F. Mukhanov, Phys. Lett. B 458 (1999) 219 [arXiv:hep-th/9904176].
E. Silverstein and D. Tong, Phys. Rev. D 70 (2004) 103505 [arXiv:hep-th/0310221].
X. Chen, Phys. Rev. D 71 (2005) 063506 [arXiv:hep-th/0408084].
X. Chen, JHEP 0508 (2005) 045 [arXiv:hep-th/0501184].
22
[72]
[73]
[74]
[75]
[76]
[77]
[78]
[79]
[80]
[81]
[82]
[83]
[84]
[85]
[86]
[87]
[88]
[89]
[90]
[91]
[92]
[93]
[94]
[95]
[96]
[97]
[98]
[99]
[100]
[101]
[102]
[103]
[104]
X. Chen, Phys. Rev. D 72 (2005) 123518 [arXiv:astro-ph/0507053].
S. Kecskemeti, J. Maiden, G. Shiu and B. Underwood, JHEP 0609, 076 (2006) [arXiv:hep-th/0605189].
J. E. Lidsey and D. Seery, Phys. Rev. D75, 043505 (2007), astro-ph/0610398.
D. Baumann and L. McAllister, Phys. Rev. D 75, 123508 (2007) [arXiv:hep-th/0610285].
R. Bean, S. E. Shandera, S. H. Henry Tye and J. Xu, JCAP 0705, 004 (2007) [arXiv:hep-th/0702107].
J. E. Lidsey and I. Huston, JCAP 0707, 002 (2007) [arXiv:0705.0240 [hep-th]].
H. V. Peiris, D. Baumann, B. Friedman, and A. Cooray, Phys. Rev. D76, 103517 (2007), 0706.1240.
T. Kobayashi, S. Mukohyama and S. Kinoshita, JCAP 0801, 028 (2008) [arXiv:0708.4285 [hep-th]].
F. Gmeiner and C. D. White, JCAP 0802 (2008) 012 [arXiv:0710.2009 [hep-th]].
L. Lorenz, J. Martin, and C. Ringeval, JCAP 0804, 001 (2008), 0709.3758.
R. Bean, X. Chen, H. Peiris and J. Xu, Phys. Rev. D 77, 023527 (2008) [arXiv:0710.1812 [hep-th]].
S. Bird, H. V. Peiris and D. Baumann, Phys. Rev. D 80, 023534 (2009) [arXiv:0905.2412 [hep-th]].
D. A. Easson, R. Gregory, D. F. Mota, G. Tasinato and I. Zavala, JCAP 0802 (2008) 010 [arXiv:0709.2666 [hep-th]].
A. A. Starobinsky and J. Yokoyama, arXiv:gr-qc/9502002.
M. x. Huang, G. Shiu and B. Underwood, Phys. Rev. D 77 (2008) 023511 [arXiv:0709.3299 [hep-th]].
D. Langlois and S. Renaux-Petel, JCAP 0804 (2008) 017 [arXiv:0801.1085 [hep-th]].
S. Renaux-Petel and G. Tasinato, JCAP 0901 (2009) 012 [arXiv:0810.2405 [hep-th]].
L. Boubekeur and D. H. Lyth, Phys. Rev. D 73, 021301 (2006) [arXiv:astro-ph/0504046].
L. Alabidi and D. H. Lyth, JCAP 0605, 016 (2006) [arXiv:astro-ph/0510441].
N. Kogo and E. Komatsu, Phys. Rev. D 73, 083007 (2006) [arXiv:astro-ph/0602099].
D. Babich, P. Creminelli and M. Zaldarriaga, JCAP 0408, 009 (2004) [arXiv:astro-ph/0405356].
W. Hu, Phys. Rev. D 64, 083005 (2001) [arXiv:astro-ph/0105117].
T. Okamoto and W. Hu, Phys. Rev. D 66, 063008 (2002) [arXiv:astro-ph/0206155].
P. Creminelli, L. Senatore and M. Zaldarriaga, JCAP 0703, 019 (2007) [arXiv:astro-ph/0606001].
R. G. Leigh, Mod. Phys. Lett. A 4, 2767 (1989).
C. Gordon, D. Wands, B. A. Bassett and R. Maartens, Phys. Rev. D 63 (2000) 023506 [arXiv:astro-ph/0009131].
R. Arnowitt, S. Deser and C. W. Misner, Phys. Rev. 117 (1960) 1595.
S. Weinberg, Phys. Rev. D 72 (2005) 043514 [arXiv:hep-th/0506236].
D. Wands, N. Bartolo, S. Matarrese and A. Riotto, Phys. Rev. D 66 (2002) 043520 [arXiv:astro-ph/0205253].
P. Creminelli and M. Zaldarriaga, JCAP 0410, 006 (2004), astro-ph/0407059.
C. Cheung, A. L. Fitzpatrick, J. Kaplan, and L. Senatore, JCAP 0802, 021 (2008), 0709.0295.
H-Y. Chen, J-O. Gong, K. Koyama and G. Tasinato, in preparation.
S. Renaux-Petel, arXiv:0907.2476 [hep-th].