Solid Consistency Lorenzo Bordin,a,b Paolo Creminelli,c Mehrdad Mirbabayi,c,d Jorge Noren˜ae a SISSA, via Bonomea 265, 34136, Trieste, Italy b INFN, National Institute for Nuclear Physics, Via Valerio 2, 34127 Trieste, Italy 7 c Abdus Salam International Centre for Theoretical Physics 1 Strada Costiera 11, 34151, Trieste, Italy 0 2 d Stanford Institute for Theoretical Physics, Stanford University, Stanford, CA 94305, USA r a M e Instituto de F´ısica, Pontificia Universidad Cato´lica de Valpara´ıso, Casilla 4059, Valpara´ıso, Chile 9 ] O Abstract C We argue that isotropic scalar fluctuations in solid inflation are adiabatic in the super-horizon limit. . h Duringthesolidphasethisadiabaticmodehaspeculiarfeatures: constantenergy-densityslicesandco- p - movingslicesdonotcoincide, andtheircurvatures,parameterizedrespectively byζ and ,bothevolve o R r in time. The existence of this adiabatic mode implies that Maldacena’s squeezed limit consistency re- t s lation holds after angular average over the long mode. The correlation functions of a long-wavelength a [ spherical scalar mode with several short scalar or tensor modes is fixed by the scaling behavior of the 2 correlators of short modes, independently of the solid inflation action or dynamics of reheating. v 2 8 3 1 Introduction and Results 4 0 . 1 0 Consistency relations (CRs) in single-field inflation are a consequence of adiabaticity: a long mode is 7 locally unobservable and its effect can be removed by a coordinate redefinition [1, 2]. In the presence 1 : of additional fields, long-wavelength relative fluctuations (entropy modes) can be locally observed and v i CRs are violated. This common lore is challenged when one considers models of inflation with a X different symmetry structure that cannot be described in the framework of the EFT of inflation [3]. r a In this paper we focus on the case of solid inflation [4, 5] where the “stuff” that drives inflation has the same symmetry as an ordinary solid. Here the situation is different from the usual case. In solids there is a single scalar excitation: the longitudinal phonon. However this mode is not adiabatic: the perturbation is anisotropic and this anisotropy is locally observable even at very long wavelengths. The absence of adiabaticity suggests at first sight that one cannot derive any CR. This conclusion is too quick. The existence of CRs also with this different symmetry structure can be seen in this way. If one considers an isotropic perturbation of the solid, i.e. a dilation or compression, this will be adiabatic since in solids there is a unique relation between the pressure and 1 the energy density, p(ρ). Since the solid experiences all states of compression as the universe expands, this perturbation cannot be locally distinguished from the unperturbed evolution. We are going to verify this statement in Section 2 showing that an isotropic superposition of linear scalar modes is indeed adiabatic. This adiabatic mode is not standard: the two variables ζ and do not coincide and R they are both time dependent. This stems from the fact that the solids do not admit curved FRW solution, but only flat ones. The existence of adiabatic modes imply CRs for the variable ζ. This does not happen for R since the diffeomorphism which removes the long mode cannot be written in terms of in a model- R independent way. In Section 3 we are going to verify the CRs in various cases, both when the short modes are inside the Hubble radius and outside. The conclusion is that, in models with the symmetry pattern of solid inflation, after reheating correlation functions satisfy the usual CRs once an average over therelative orientation between longand shortmodeshas beendone. We discussthe implications of this and open questions in Section 4. Before proceeding, let us recall, following [5], that the dynamics of the solid is described in terms of three scalar fields φI which parametrise the position of the elements of the solid: φI = xI+πI. The action can be written in terms of SO(3)-invariant objects built out of the matrix BIJ ∂ φI∂µφJ. µ ≡ One can choose these invariants to be ([...] indicates a trace) [B2] [B3] X [B] , Y , Z . (1) ≡ ≡ [B]2 ≡ [B]3 So the action, at lowest order in derivatives and including gravity, is M2 S = d4x√ g PlR+F(X,Y,Z) . (2) − 2 Z (cid:20) (cid:21) 2 Long Isotropic Perturbation 2.1 Proof of adiabaticity We want to show that, in solid inflation, a scalar perturbation becomes adiabatic once we average over thesolidangle,i.e.wetakeasuperpositionofscalarFouriermodeswhichisisotropic. Wehavetoprove that this kind of perturbation, in the long-wavelength limit, can be brought back to the unperturbed solution via a suitable diffeomorphism. We start from the so-called Spatially Flat Slicing Gauge (SFSG) which is defined as the gauge where the spatial part of the metric is only perturbed by tensor modes. For the rest of this Section we will only be interested in scalar modes: g = a2δ , φI = xI +πI. (3) ij ij The triplet πI consists of a scalar, π , plus a transverse vector, πI, which we neglect in the following. L T 2 The constraint equations give the lapse N = 1+δN and the longitudinal part of the shift N [5]:1 L a2H˙ π˙ H˙π /H 3a2H˙π˙ /k2 +H˙π /H L L L L δN = − , N = − . (4) − kH 1 3H˙a2/k2 L 1 3H˙a2/k2 − − In SFSG, the gauge-invariant variable ζ is given by ∂π ζ = + (π2), (5) 3 O where ∂π ∂ πi.2 Therefore, by performing the following time diff, i ≡ 1 x0 x0+ξ0(t,x), with ξ0(t,x) = ∂π, (6) → 3H we go from SFSG to the ζ gauge, defined by the condition δρ = 0, where ρ is the energy density. The − spatial part of the metric now reads g = a(t)2 (1+2ζ(t,x)) δ , while one can write π as a function ij ij L of ζ. Using eq.s (4) and (6), one can verify that the expression of δN in ζ gauge is − k d π 1 L δN = . (7) −3dt H 1 3H˙ a2 (cid:16) (cid:17) − k2 In the limit k/aH 0 the time-dependence of π is slow-roll suppressed, see eq. (15), therefore L → δN 0 on super-horizon scales. → To reproduce the unperturbed FRW solution one has to eliminate the perturbation in the scalar fields φi. This can be done by a redefinition of the spatial coordinates xi xi+ξi(t,x), with ξi(t,x) = πi(t,x). (8) → − Since now the scalars are unperturbed, it is natural to call this Unitary Gauge (UG). In UG the shift vanishes on super-horizon scales d π H L N = . (9) L −dt H 1 3H˙ a2 (cid:16) (cid:17) − k2 In this gauge, for long wavelength, δN = N = πi = 0. However the spatial part of the metric is still L perturbed g = a(t)2 (1+2ζ(t,x)δ +∂ ∂ χ(t,x)), with χ(t,x) = 6 ∂−2ζ(t,x). (10) ij ij i j − Theperturbation is purely anisotropic, i.e. the volume is not perturbedbecauseof thegauge condition δρ = 0. Therefore if one considers a spherically symmetric superposition of scalar modes, the metric perturbations in eq. (10) average to zero d2kˆ d2kˆ (2ζ δ k k χ ) = (2ζ δ 6kˆ kˆ ζ )= 0 . (11) k ij i j k k ij i j k 4π − 4π − Z Z 1We use the notation: πi = √∂i 2πL+πTi and analogously for Ni. 2Notice that when we expand−∇around the background φI = xI one has ∂iφI = δiI, so that there is no distinction between capital and lower-case spatial indeces. 3 This shows that a spherically symmetric superposition of scalar modes is adiabatic. Notice that, if one works in ζ gauge, the transformation that eliminates a long-wavelength mode − and goes back to FRW is a rescaling of the spatial coordinates: this is quite similar to the standard case of single-field inflation. However here the rescaling is time-dependent and adiabaticity requires an average over directions. In the next Section we will see that the adiabaticity gives rise to CRs: the only difference with the standard case is that they hold only after the spherical average. (The time dependence of the rescaling is immaterial because one is usually interested in correlation functions at equal time.) Instead of using time-slices with δρ = 0, one could use slices that are orthogonal to the 4-velocity of the solid. The perturbation of the spatial part of the metric is called ζ in the first case and R in the second. Contrary to the usual case, in Solid Inflation the variables ζ and differ even on R super-horizon scales. At linear level 1 ζ˙+ǫHζ = . (12) R ǫH 1+k2/3a2H2ǫ Since the two slicings do not coincide, one needs a time-diff to go from one to the other. This is the difference of the two time diff.s to go from SFSG to ζ gauge and to gauge respectively: − R− ζ ζ˙ δtR→ζ = δtζ −δtR = H − HR ≃−ǫH2 , (13) where the last equation holds on super-horizon scales. The property that δN vanishes in ζ gauge on − largescales, eq.(7),willnotholdin gauge. Thismeansthattogofrom gaugetotheunperturbed R− R− FRW one has to supplement the rescaling of spatial coordinates with the time diff eq. (13). As we will discuss in the next Section this implies that the CR for will contain an extra piece: the time diff R induces a piece involving the time-derivative of the short modes. 2.2 Squeezed vs Super-Squeezed regimes Intakingthelong-wavelength limitk 0oneendsupintheregimek aHǫ1/2. However oneexpects → ≪ that the adiabaticity arguments above hold whenever k is comfortably outside the Hubble radius and in particular also in the intermediate regime aH k aHǫ1/2. This is indeed the case. For example ≫ ≫ in this regime the expression for the lapse in ζ gauge is − d ζ 1 d ζ k2 1 k2 k k δN = < = . (14) | k| (cid:12)(cid:12)dt (cid:18)H(cid:19) 1−3H˙ ka22(cid:12)(cid:12) (cid:12)dt (cid:18)H(cid:19) a2H23H˙ (cid:12) O(cid:18)a2H2(cid:19) (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) (cid:12) (Notice that in the ineq(cid:12)(cid:12)uality above we are(cid:12)(cid:12)not(cid:12) assuming that the(cid:12)term proportional to H˙ in the denominator dominates.) The same argument works for the shift: also in the intermediate regime the physical difference with the unperturbed solution are suppressed when the mode is superhorizon. Therefore we expect the CR to hold both in the intermediate and in the super-squeezed, k ≪ aHǫ1/2,regime. However,togetanalyticalresultsoneisforcedtoexpandthesolutionoftheconstraints 4 in different ways in the two regimes and therefore one has to assume one of the two regimes. This also applies to the expression of the wavefunction, which can be expanded in the two limits to give [5] 1+ic kτ + 1c2k2τ2 eicLkτ , if c kτ ǫ, π (τ,k) Bk L 3 L | L | ≥ (15) L ≃ (Bk(cid:0)1+ǫ(1+c2L)log(−cL(cid:1)kτ) (−cLkτc)−5s/2−η/2−ǫ, if |cLkτ| ≤ ǫ, with (cid:2) (cid:3) 3 H 1 = . (16) Bk −2M c5/2ǫ1/2k5/2 Pl L 2.3 Adiabatic modes and the Weinberg theorem Solid Inflation is an interesting exception to many general theorems on cosmological perturbations. Weinberg [10, 11] showed that, under quite general assumptions, one can always find an adiabatic mode which features identical and time-independent ζ and on super-horizon scales. In the Solid R case, ζ and are neither equal (see eq. (12)) nor time-independent (both ζ and have a slow-roll R R suppressed time-dependence on super-horizon scales). Of course by linearity these properties are not changed by the spherical average. In the original paper on Solid Inflation [5] (see also [9]) the authors addressed the issue of why a scalar Fourier mode does not comply with Weinberg analysis. The point is that the solid supportsa large anisotropic stress, so that even in the long-wavelength limit the stress energy tensor remains anisotropic and thus locally distinguishable from the unperturbedsolution: the mode is not adiabatic. The theorem [10, 11] assumes the decay of the anisotropic stress for k 0. → Here we are considering a spherical average of Fourier modes and the problem takes a somewhat different form. Indeed, as we discussed, the perturbation is now adiabatic, since the anisotropic stress averages to zero. However, this adiabatic mode is still different from the one of Weinberg. One can still check that the assumptions of the theorem do not hold: the 0i component of Einstein equations is not regular for k 0, since δu diverges in that limit. This does not allow to continue a homogeneous → perturbation to a physical one at finite momentum. This however looks rather technical. What is the physical reason why the adiabatic mode we are considering is different from the standard case? Why doesn’t adiabaticity ensure that ζ is constant? In the standard case, the conservation of ζ and therelation ζ = can beunderstoodfrom a linearized R version of a spatially curved FRW. Neglecting short-wavelength perturbations and working at linear order, the curvature of a constant ρ slice is given by 4 (3)R = ∂2ζ. (17) −a2 Since for a curved FRW the spatial curvature κ = a2 (3)R/6 is constant, super-horizon ζ fluctuations better be time-independent. Moreover for a curved FRW the surfaces of constant density are perpen- dicular to the 4-velocity so we need ζ = . The adiabatic mode we are discussing in Solid Inflation R does not have these properties and this is related to the fact that one does not have curved FRW solution in this model. 3 This can be understood in terms of symmetries: the internal symmetries 3 Curved FRW solutions are allowed if we change the internal metric in the Lagrangian, see [12]. However, 5 of the φI is isomorphic to the symmetries of flat Euclidean space. This allows to write flat FRW solutions, but not spatially curved solutions, which have a different group of isometries.4 3 Angle-averaged Consistency Relations In ζ gauge a long mode averaged over the direction can be removed by a rescaling of the spatial − coordinates. The derivation of the CR is very similar to the standard case, apart from the required angular average. In the case of the scalar 3-point function we get d2qˆ dlogk3P (k) ′ ζ 4π hζqζkζ−k−qiq≪k = − dlogk Pζ(q)Pζ(k), (18) Z where here and in the following the prime indicates that a momentum-conserving delta function, (2π)3δ( k ), is dropped. We stress that this result, in the limit of exact scale-invariance when the i i RHS of eq. (18) vanishes, was already discussed in [8]. P 3.1 Check of the consistency relation for ζζζ h i Let us check the CR eq. (18). In Solid Inflation the quadratic action is (ǫ) while the cubic action O is (ǫ0). Thus the 3-point function is slow-roll enhanced, f ζ3 / ζ2 2 = (ǫ−1) [5]. Since NL O ∼ O the tilt of the 2-point function outside the horizon is (ǫ), a non-trivial (in the sense of non-zero) O (cid:10) (cid:11) (cid:10) (cid:11) verification of eq. (18) in this regime would require taking into account corrections to the leading bispectrum at second-order in slow-roll. This is quite challenging. We content ourselves with the first order correction: at this order the two sides of eq. (18) should vanish when all the modes are outside the horizon. When the short modes are inside the horizon the scale-dependence of the spectrum is not slow-roll suppressed and the LHS of eq. (18) should thus be non-zero.5 The check will be done in the regime k aHǫ1/2 for all the modes. ≫ To do this we compute the cubic Lagrangian in SFSG up to (ǫ), calculate the bispectrum and O then transform to ζ gauge. The (ǫ) corrections to the bispectrum of three super-horizon modes − O was studied in detail in [7]. Thus, we skip most of the technical steps. However, we identify a missing term in [7] that is important for the CR to work. The in-in calculation up to this order consists of the sum of three pieces. Schematically, ζζζ (3) π(τ,k) + (3) π(τ,k) + (3) π(τ,k) , (19) h i∼ LO(1)× O(1) LO(1)× O(ǫ) LO(ǫ)× O(1) where L(O3()ǫn) and π(τ,k)O(ǫn) are respectively the cubic Lagrangian and the wavefunctions evaluated at nth order in slow-roll. The leading cubic Lagrangian, (3) , was calculated in [5] (eq. (D.2)). In LO(1) our argument for constancy of ζ in the standard case is based on the fact that curvature is a free parameter of the backgroundsolution. In the models discussed in [12] curvature is uniquely fixed in terms of energy density, so we don’t expect ζ to be conserved. 4This also implies the usual curvature problem takes a somewhat different flavour in this class of models. 5For a discussion of CRs in standard inflation when the short modes are inside the horizon see [6]. 6 Fourier space and in the squeezed limit it reduces to 8 L(O3()1) squeezed = −81FY 1−3cos2(θ) πL,qπL,kπL,−q−k , (20) (cid:12) (cid:12) (cid:0) (cid:1) where θ is the relative ang(cid:12)le between the long and the short modes. The above expression gives zero when one takes the angular average. This means there is no contribution (ǫ−1) to the LHS of O eq. (18). Notice that the cancellation after angular average holds independently of the explicit form of the wavefunctions. Therefore, the second term of eq. (19) vanishes and only the last term is relevant for checking the CR. Cubic scalar Lagrangian at (ǫ). At first look, expanding eq. (D.1) of [5] up to first order in O slow-roll seems like a formidable task. Since Y and Z in (1) are defined in such a way that they start from second order in perturbations, one has to expand 1 (3) = F δX(3)+F δX(1)δX(2)+ F (δX(1))3+F δX(1)δY(2)+(F δY(3)) +(Y Z), LO(ǫ) X XX 6 XXX XY Y O(ǫ) ↔ (21) where subscripts on F denote partial derivatives. However, there are several simplifications [7]. As we will see, one only needs the SFSG deformation matrix BIJ at zeroth order in slow-roll parameters. This allows neglecting Ni and δN (which are slow-roll suppressed in the regime we are considering): 1 BIJ = δIJ +∂IπJ +∂JπI +∂ πI∂ πJ π˙Iπ˙J + (ǫ) (22) a2 k k − O (cid:0) (cid:1) Therefore, δX terminates at quadratic order, apart from slow-roll suppressed corrections 3 2 1 a2 δX = δ[B] = ( ∂π+ ∂ π ∂ π π˙2)+ (ǫ). (23) a2 3 3 i j i j − 3 i O Given that derivatives of F with respect to X are slow-roll suppressed, there is no contribution from F δX to (ǫ) cubic Lagrangian. Another simplification is that if at some order in perturbations the X O corrections to [Bn] involve at most m n of the BIJ factors, then ≤ [Bn] [Bm] δ = δ , for all n m . (24) [B]n [B]m ≥ (cid:18) (cid:19) (cid:18) (cid:19) This implies that δY(2) = δZ(2), (25) and since slow-roll corrections to BIJ start at (π2) O (3) (3) δY = δZ . (26) O(ǫ) O(ǫ) Therefore, the following combinations appear in (21) (F +F )δY(3) , (F +F )δX(1)δY(2). (27) Y Z O(ǫ) XY XZ However F +F = (ǫ) and F +F = (ǫ2), so these terms are negligible. Using Y Z XY XZ O O a4 2a6 F = ǫF, F = ǫF, F = 3M2H2, (28) XX − 9 XXX 27 − pl 7 one gets 2 8 2 (3) = ǫM2 H2a3 (∂π)∂ πk∂ πk (∂π)3 ǫM2 H2a2(∂π)π˙2 LO(ǫ) Pl 3 j j − 27 − 3 Pl i (cid:20) (cid:21) 4 + (F +F )a2 (∂π)π˙2 3∂ πjπ˙iπ˙j . (29) 27 Y Z i − i (cid:2) (cid:3) The terms with time derivatives in this equation are absent from eq. (35) of [7]. Note that the appearance of the combination F + F on the second line is a consequence of (24) since time- Y Z derivatives appear in BIJ starting from quadratic order. The angular average of this term is zero, hence it does not contribute to CR while the last term on the first line does contribute. Field redefinition. Since we are interested in the bispectrum of ζ at f = (1), we need to find NL O the relation between ζ and π at quadratic order and to zeroth order in ǫ. We start from BIJ in SFSG given in (22). The last term can be neglected because in the squeezed limit at least one of the two π’s will be out of the horizon with a slow-roll suppressed time evolution. The assumption of spherical symmetry simplifies the expression, δIJ 2 1 BIJ = δIJX(t) = 1+ (∂π)+ (∂π)2 . (30) a2 3 9 (cid:18) (cid:19) Now we perform the time diffeomorphism that leads to the ζ gauge (the analogue of eq. (6) but now − at second order), where X(t) takes its unperturbed value X(t+ξ0(t,x);x) = X¯(t)= a−2. (31) At non-linear order the spatial part of the ζ gauge metric is defined as g = a2e2ζδ . Therefore up ij ij − to quadratic order in π, we obtain 1 1 1 ζ = Hξ0 = ∂π (∂π)2+ (∂π˙)(∂π)+ ((∂π)3). (32) 3 − 18 9H O One can now put together the in-in computation, using the Lagrangian in the first line of eq. (29), with the definition of ζ, eq. (32), to get d2qˆ 1 d2qˆ 1 4π hζqζkζ−q−ki′q≪k = 27 4π h(∂π)q(∂π)k(∂π)−k−qi′q≪k −2Pζ(q)Pζ(k)+ HPζ(q)P˙ζ(k). (33) Z Z When allthemodesareoutsidethehorizon thelast termon theRHSof eq.(33)is slow-roll suppressed and can be neglected. A straightforward calculation shows that the other two terms cancel each other confirming that d2qˆ ′ 4π hζqζkζ−k−qiq≪k = O(ǫ), (34) Z as implied by eq. (18).6 When the short modes are inside the horizon the cancellation between the first two terms on the RHS of eq. (33) still holds, but now the last term is non-negligible since the 6 It is challengingto performthe calculationatnext orderin slow-roll,to test the CR. The cancellationafter angular average in eq. (33) can be seen only after the explicit in-in integral (in contrast to what happens at leading order, eq. (20)). This means that, at higher order, we should compute time integrals involving the wavefunctions at (ǫ). O 8 time dependence is not slow-roll suppressed in this regime. One has d2qˆ 1 dlogk3P (k,τ) 4π hζqζkζ−k−qi′q≪k = HPζ(q)P˙ζ(k) = − dlogζk Pζ(q)Pζ(k,τ), (35) Z where in the last passage we used that the power spectrum is of the form P = k−3f(kτ) as dictated ζ by scale invariance, up to corrections of order slow-roll. The CR eq. (18) is verified. 3.2 Check of the consistency relation for ζγγ h i One novel feature of our analysis is that the (spherically averaged) adiabatic modes of solid inflation feature a time-dependent ζ outside the horizon. Since this time dependence arises at (ǫ) it would be O nice to check the CR at this order. As discussed this is quite challenging for ζζζ , but it is doable for h i ζγγ , where the scalar mode is taken to be long. The leading term in ζγγ is (1) [8], so one just h i h i O needs to do the calculation including the first-order slow-roll corrections. The leading Lagrangian has the form (see eq. (A.7) of [8]) 1 (3) (∂π)γ γ +γ γ ∂kπi. (36) LO(ǫ0) ∝ −3 ij ij ij jk It averages to zero in the squeezed limit πLq→0 independently of the wavefunctions. Thus we do not need to consider the slow-roll corrections to the wavefunctions. We go directly to the computation of (3) . LO(ǫ) Cubic Lagrangian at (ǫ). Thereare two terms which contribute to (3) . Thefirstarises from O LO(ǫ) the expansion of the function F(X,Y,Z), while the second is from the Einstein-Hilbert action. The expansion of F gives F (3) F δX +F δZ + XXδX2+F δXδY +F δXδZ (37) LO(ǫ) ⊃ X Z 2 XY XZ where, 1 δX = a−2 2(∂π) 2γ ∂ πj + γ 2+γ γ ∂ πi , (38) ij i ij ij jk k − 2 (cid:20) (cid:21) 4 1 2 2 δY = γ ∂ πj + γ 2+ γ γ ∂ πi γ 2(∂π) , (39) ij i ij ij jk k ij −9 9 3 − 9 (cid:20) (cid:21) 4 1 8 8 δZ = γ ∂ πj + γ 2+ γ γ ∂ πi γ 2(∂π) . (40) ij i ij ij jk k ij −9 9 9 − 27 (cid:20) (cid:21) One gets 1 (3) = 3ǫa2M2 H2 γ γ ∂ πj γ2(∂π) . (41) LO(ǫ) − Pl ij jk k − 3 ij (cid:20) (cid:21) This gives zero after the angular average. The contribution which arises from the Einstein-Hilbert action (plus the appropriate boundary terms) is (3) a3 N (3)R+N−1(E Eij E2) LO(ǫ) ⊃ ij − O(ǫ)πγγ h i a2 a3 = δN γ′ 2+(∂ γ )2 γ′ ∂ γ Nk. (42) − 8 ij l ij − 4 ij k ij h i 9 The in-in calculation can be done separately in the intermediate regime aH q aHǫ1/2 and in the ≫ ≫ super-squeezed regime q aHǫ1/2. 7 In both regimes the in-in computation of the 3-point function ≪ gives 1 d2qˆ (∂π) γsγs = 2ǫP (q)P (k). (43) 3 4π q k −k q≪k ζ γ Z (cid:10) (cid:11) Tensor modes at quadratic order in perturbations. Thefinalexpressionof thebispectrum is given once one considers the contribution coming from the time diff that has to be performed to go from SFSG to ζ gauge. This changes the tensor perturbations at quadratic level (see eq. (A.8) of − [1]). The interesting part (for us) is 1 (∂π) γ = γ + γ˙ , (44) ζ π π H 3 where γ and γ denote tensor perturbations respectively in ζ gauge and SFSG. This adds a contri- ζ π − bution to the bispectrum: ζγ γ = ζγ γ + 2 ζ∂πγ γ˙ , which is given to leading order by h ζ ζi h π πi 3H h π πi 1 d P (q) P (k) = 2ǫ(1+c2)P (q)P (k). (45) H ζ dt γ − L ζ γ Eqs. (43) and (45) give d2qˆ dlogk3P (k) ζ γsγs = 2c2ǫP (q)P (k) = γ P (q)P (k). (46) 4π q k −k q≪k − L ζ γ − dlogk ζ γ Z (cid:10) (cid:11) This is exactly what is predicted by the CR. In conclusion, this computation confirms that CRs after spherical average hold even at slow roll order, i.e. when the time dependence of the long mode cannot be neglected. 3.3 Consistency relation for ? R There are two differences if one wants to use the variable instead of ζ. First, as discussed above, R starting from gauge one needs an extra time diff to map a long mode into an unperturbed FRW. R− This changes the CR and introduces a time derivative of the short mode 2-point function. Moreover on super-horizon scale, in both squeezed and super-squeezed regimes, the relation between and ζ R depends on c and therefore is non-universal, L c2ζ. (47) R ≃ − L Hence, the spatial rescaling (8), which is determined in terms of ζ, will be model-dependent when written in terms of . This means that in this gauge we are not going to be able to write an explicit R CR, i.e. a model independent relation among correlation functions. 7Notice thateq.(42), evaluatedinthe super-squeezedlimit is (1). But δN, NL π˙L whichhasa slow-roll O ∝ suppressed time dependence and so the final expression of the bispectrum is (ǫ). Notice also that, naively, O the term proportional to the shift seems to cancel when one performs the angular average. However in the super-squeezed limit this term would give a bispectrum 1/q4 for q 0 (see the expression for Ni in eq. (4)). ∝ → This leading behaviour cancels, even before taking the angular average. The subleading contribution has the correct 1/q3 dependence and does contribute to the angular average. 10