Heavy spinning particles from signs of primordial non-gaussianities: beyond the positivity bounds

Within the so-called cosmological collider program, imprints of new particles on primordial non-Gaussianities have been studied intensively. In particular, their non- analytic features in the soft limit provide a smoking gun for new particles at the inflation scale. While this approach is very powerful to probe particles of the mass near the Hubble scale, the signal is exponentially suppressed for heavy particles. In this paper, to enlarge the scope of the cosmological collider, we explore a new approach to probing spins of heavy particles from signs of Wilson coefficients of the inflaton effective action and the corresponding primordial non-Gaussianities. As a first step, we focus on the regime where the de Sitter conformal symmetry is weakly broken. It is well known that the leading order effective operator (∂μ𝜙∂μ𝜙)2 is universally positive as a consequence of unitarity. In contrast, we find that the sign of the six derivative operator (∇μ∂ν𝜙)2(∂ρ𝜙)2 is positive for intermediate heavy scalars, whereas it is negative for intermediate heavy spinning states. Therefore, under the assumption of tree-level UV completion, the sign can be used to probe spins of heavy particles generating the effective interaction. We also study phenomenology of primordial non-Gaussianities thereof.


Introduction
The energy scale of inflation could be as high as 10 13 GeV, hence it would be a phenomenon at the highest energy scale we may explore. In particular, primordial non-Gaussianities can be thought of as an extremely high energy particle collider. For example, their nonanalytic behaviors in the soft limit are associated with on-shell particle creations and thus provide a direct evidence of new particles at the inflation scale [1][2][3][4]. Just like resonance signals at particle colliders, it is very powerful to probe particles with the mass m near the Hubble scale H. Such an approach is dubbed the cosmological collider program and has been studied intensively . For heavy particles, on the other hand, the signal is exponentially suppressed by the Boltzmann factor. For scalars, this factor reads e −π m 2 H 2 − 9 4 , which is as small as 10 −4 even for m = 3H! Therefore, it is desirable to develop a new approach to probing heavy particles to enlarge the scope of the cosmological collider. 1 1 In contrast to ordinary colliders, we cannot build a new cosmological collider with a higher energy H!

JHEP12(2019)107
Now let us recall our history of particle physics, where not only resonance signals, but also detailed studies of low-energy effective interactions have been useful to probe new particles. A typical example is the prediction of weak bosons, where angular dependence of the Fermi interactions played an important role. Following the history, we would like to apply a similar idea to inflaton effective interactions and primordial non-Gaussianities. Recall that the inflaton enjoys a shift symmetry under the slow-roll approximation. Its effective Lagrangian then reads [43] where the dots stand for higher dimensional operators and the cutoff Λ is typically associated to the mass of intermediate heavy states. Remarkably, it is well known that α is always positive in a wide class of theories as a consequence of unitarity and analyticity of scattering amplitudes [44] (see, e.g., [45][46][47][48][49][50][51][52][53][54][55][56][57] for recent applications). While the positivity is an elegant consistency condition on IR effective theories, 2 detailed informations of the UV theory, e.g., spins of heavy states, are obscured at the cost of universality.
In this paper we would like to go beyond the positivity and develop an approach to probing spins of heavy states from signs of inflaton effective interactions and the corresponding primordial non-Gaussianities. As a first step, we focus on the regime where the de Sitter conformal symmetry is weakly broken. 3 In this regime, we may focus on operators with four inflaton fields for the study of primordial three-point and four-point functions: which provides the most general effective Lagrangian up to six derivatives after using treelevel equations of motion to remove operator degeneracy. In contrast to α, the sign of β cannot be fixed by analyticity and unitarity essentially because the corresponding fourpoint amplitude vanishes in the forward limit. In other words we may have a chance to probe spins of heavy states from its sign. By generalizing the analysis in [44] to non-forward amplitudes, we show that 1. β > 0 for intermediate heavy scalars, whereas β < 0 for heavy spinning states with spin l = 2, 4, 6, . . ., as long as four-point amplitudes are bounded as < s 2 in the high energy limit (s is the standard Mandelstam variable).
Note that our conclusion is about intermediate "states," which can be multi-particle states in general. See appendix A for more comments. While the argument is applicable to a wide class of theories, it cannot be applied directly to the effective coupling generated 2 In other words, if experiments find violation of the positivity bounds, we have to change our approach to UV completion in a drastic way, which is also an interesting possibility. 3 The nonlinearity parameter fNL for three-point functions is generically smaller than O(1) in this regime.
At this cost, we may utilize developments on Lorentz invariant four-point amplitudes instead. The same remark applies to the elegant works [4,36] based on the de Sitter invariant four-point functions. In order to study the regime fNL > O(1), we need to incorporate more general higher derivative operators such as (∂µφ∂ µ φ) n (n ≥ 3) or work in the effective field theory of primordial perturbations [58].

JHEP12(2019)107
Four derivative Six derivative Dominant UV States IR a 2,0 > 0, α > 0 a 2,1 < 0, β > 0 Scalar a 2,1 = 0, β = 0 Exact Cancellation a 2,1 > 0, β < 0 Spin l = 2, 4, 6 . . . Table 1. The sign of α in the low energy EFT (or equivalently the coefficient a 2,0 of s 2 in the low-energy four-point scattering) is universally positive. On the other hand, the sign of β (or equivalently the coefficient a 2,1 of s 2 t in the low-energy four-point scattering) reflects which UV states provide dominant contributions to low-energy scattering. The conclusion summarized in the table applies to any weakly coupled, Lorentz-invariant, analytic and unitary UV completion as well as effective interactions mediated by KK gravitons.
by the exchange of KK gravitons, which is one of the main targets in the cosmological collider program, 4 essentially because it involves gravitational dynamics (see section 2.4 for details). Interestingly, however, we find that 2. the six derivative operator generated by the KK graviton exchange has a negative coefficient β < 0.
This shows that a non-positive coefficient β ≤ 0 is a sign of heavy spinning particles with the spin l = 2, 4, 6 . . . in a wide class of theories in our interests. This is one of the main results in our paper. See table 1 for a summary. We also study phenomenology of primordial non-Gaussianities thereof. The organization of the paper is as follows. In section 2 we first study the relation between the sign of the six derivative operator (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 and the spins of intermediate heavy states. We then study shapes of four-point functions (section 3) and three-point functions (section 4) to clarify if those generated by the six derivative operator are distinguishable from the four derivative ones. We conclude in section 5 with a discussion of future directions.

Spin-dependence of six derivative operators
In this section we clarify the spin-dependence of the six derivative operator. In section 2.  Figure 1. Integration contour and analytic structure: The red lines represent branch cuts associated to multi-particle states generated by loops, whereas the red dots are single poles associated to single-particle states. The left panel is the integration contour used to read off the IR coefficients. By deforming it into the one on the right panel, we may clarify how the UV data is encoded into the IR coefficients.
cannot apply the results in sections 2.1-2.2 in general because of gravitational effects, but it is of great interests in the cosmological collider program. A remark on loops is also given in appendix A.

A basis of general UV amplitudes
Let us consider four-point scattering amplitudes M (s, t) of identical massless scalars and elaborate on the relation between IR coefficients and the UV spectrum. For this purpose, we generalize the analysis in [44] to non-forward amplitudes (see also [47]). First, from the IR point of view, it is convenient to expand the amplitude M (s, t) in Mandelstam variables s and t as where we assumed that gravity is subdominant to neglect massless cuts and the massless graviton pole. The coefficient function of s p can then be evaluated by the contour integral, where the integration contour is defined such that it encircles the origin s = 0 and the integrand is analytic inside the contour except for the origin (see the left panel of figure 1). Also we take t infinitesimal to avoid unphysical poles. Next, we deform the contour as depicted in the right panel of figure 1 by assuming that the amplitude is analytic away from the real axis on the physical sheet. If the amplitude JHEP12(2019)107 is bounded as |M (s, t)| < |s| p at UV, the integral at infinity vanishes to obtain ds 2πi where is an infinitesimal positive constant and m 0 is the mass of the lightest intermediate (massive) state. In general the amplitude contains single poles associated with tree-level exchange and branch cuts associated with loops, so that the integral is given by r.h.s. of (2.3) = n g 2 n P ln 1 + 2t where the first line and the second line correspond to single poles and branch cuts, respectively. P n (z) is the Legendre polynomial, m n and l n are the mass and spin of the intermediate particle labeled by n, and g n characterizes the three-point amplitude of two massless scalars and one massive particle n. ρ(M 2 , l) is a spectral density for the transition amplitude from two massless particles to multi-particle states with the energy M and total angular momentum l. s * is the minimum energy squared of intermediate multi-particle states. For notational simplicity, we write eq. (2.4) as r.h.s. of (2.3) = n g 2 n P ln 1 + 2t where one should understand that n stands for both of the discrete sum for tree-level exchange and the integration (continuous sum) for loops. Note that since we are considering identical external massless scalars, the intermediate state can have an even spin only l n = 0, 2, 4, . . . as required by the exchange symmetry (for simplicity, we use "spins" to denote total angular momenta of multi-particle states as well). See, e.g., [59] for details. We now conclude that the IR coefficients are given in terms of the UV data as ∞ q=0 a p,q t q = n g 2 n P ln 1 + 2t (2.6) Note that this relation is applicable only for p satisfying |M (s, t)| < |s| p at UV. On general grounds, we assume that the UV amplitude is bounded as |M (s, t)| < |s| 2 , which is required by the Froissart-Martin bound [60] for example. 5 Hence, the relation is applicable only for p ≥ 2. It is also instructive to reorganize the amplitude as 5 To be precise, the Froissart-Martin bound is applicable only in gapped theories and in the regime 0 ≤ t < 4m 2 (m is the mass of external particles). However, its use can be justified by turning on a tiny mass of the scalar as a regulator, which does not change the conclusion of our analysis. Also note that the inflaton indeed has a tiny mass which is suppressed due to the slow-roll conditions, even though the mass is negligible as long as we are interested in non-Gaussianities generated by higher derivative operators.

JHEP12(2019)107
Here we introduced α p (t) = ∞ q=0 a pq t q for p = 0, 1, which cannot be determined by the contour deformation argument in the above. Eq. (2.7) provides a basis of general UV amplitudes which reproduce the correct factorization and satisfy the Froissart-Martin bound.

Positivity and the beyond
We then explore imprints of the UV data on the IR coefficients a p,q . First, from eq. (2.6), the IR coefficients read is the shifted factorial (also called the Pochhammer symbol) and we used Let us next take a closer look at the first few orders in the t expansion: First, the leading order coefficients in the t expansion are given by for even p , 0 for odd p , which is nothing but the well-known positivity bounds on the O(s 2n ) coefficients of forward amplitudes [44]. This bound is universal and elegant, but detailed information such as spins of the UV states is obscured at the cost of universality.
Similarly, the next-to-leading order coefficients read for odd p . (2.11) The point here is that the sign of a p,1 for even p depends on the spin l n of the intermediate states. More explicitly, we find that intermediate states with the spin higher (lower) than the critical value l * = ( √ 2p + 3 − 1)/2 gives a positive (negative) contribution to a p,1 . In particular, the critical value is l * ∼ 0.8 for p = 2. In terms of the effective Lagrangian (1.2), it means that β is positive only for intermediate scalars as we discuss shortly. This six derivative operator is the first nontrivial operator in the derivative expansion which may be used to probe spins of intermediate states.
On the other hand, a p,1 for odd p is universally positive.

JHEP12(2019)107
As a next example, let us consider the O(t 2 ) term, for which we find (2.12) For odd p, there exists a critical spin l * = ( √ 2p + 5 − 1)/2 below (above) which a p,q is negative (positive). On the other hand, for even p, there exists a window √ p + 2 − 1 < l n < √ p + 2 inside (outside) which a p,q is negative (positive). Note that it has a width 1, so that at most one value of spin can be inside the window. For example, the window is 1 < l n < 2 for p = 2, hence a 2,2 = 0 for spin 2, whereas a 2,2 > 0 for other even spins. Similarly, spin 2 is inside the window for p = 4, 6.

Higher orders
It is straightforward to go higher in a similar way. First, a p,q is a polynomial in l n of order 2q for even p and order 2q − 2 for odd p. It is positive when the spin l n is sufficiently large for fixed p, q. On the other hand, the contribution from a scalar l n = 0 is simply (−1) p+q (p + 1) q /q!(m 2 n ) p+q+1 , 6 hence the sign depends only on (−1) p+q . The sign for general spin depends on details of the polynomial in l n and the phase structure will be richer for higher q. However, we will not go into more details because such higher order terms are not easy to probe phenomenologically, leaving it for future work.
To summarize, we have elaborated on the relation between the sign of IR coefficients a p,q and the spin of intermediate heavy states. The result can also be rephrased in terms of EFT coefficients as follows: for example, four-point scattering amplitudes can be evaluated from the effective Lagrangian (1.2) as where the dots stand for higher derivative terms negligible at low-energy. We then find (2.14) Therefore, our analysis on a 2,1 implies that the coefficient β of the six derivative operator In other words, detection of negative or vanishing β is a smoking gun of spinning heavy states with the even spin l n = 2, 4, . . .. 7 This is one of our main results.

JHEP12(2019)107
It should be noted again that our results are about the spins (or total angular momenta) of heavy states, which can be multi-particle states when loop effects are dominant. However, in the case of inflation, the inflaton has to enjoy the approximate shift symmetry. As a result, the loop effects are generically suppressed by either weak couplings or slow-roll parameters as we discuss in appendix A.3. In such a case, we can use the sign of effective interactions to probe the spins of the intermediate heavy particles exchanged at the treelevel without swamped by multi-particle states associated with loops.

Open superstring amplitudes
In the previous subsection we have shown that the non-positivity of the coefficient β of the six derivative operator (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 is a signature of higher spin states. However, it does not mean that β ≤ 0 holds in all theories with higher spins because intermediate scalars may dominate to make β positive. Therefore, it is useful to demonstrate that there indeed exists a UV completion with higher spins and a non-positive β. In this subsection we study open superstring theory to provide a concrete example for such a UV completion.
Let us consider four-point scattering of identical massless scalars (an extra-dimensional component of the gauge boson) localized on a D3 brane in the open superstring. The corresponding disk amplitude reads 8 (see, e.g., [61]) We also took the unit in which the open string spectrum is given by m 2 = 0, 1, 2, . . .. At IR, we may expand the amplitude in the Mandelstam variables as where the dots stand for fifth and higher orders in Mandelstam variables. We find that the six derivative operator vanishes, which means that there exists an exact cancellation between scalars and higher spins. It is also instructive to explicitly see how the cancellation happens between scalars and higher spins by comparing the open string amplitude (2.15) with the master formula (2.7) to read off the cubic couplings g n . First, the amplitude (2.15) has single poles associated with particles of the mass squared m 2 = 1, 3, 5, . . .. For example, the residue of the s-channel pole reads  It is a polynomial in t of order n + 1, which consists of contributions from intermediate particles with spin up to n + 1. More explicitly, we expand it by Legendre polynomials as where g nl characterizes cubic coupling of two massless scalars and the massive particle of the mass squared m 2 = n and spin l. 9 For example, the n = 1 sector is given by Note that odd spins do not appear as an intermediate particle because of the exchange symmetry of identical scalars mentioned earlier. We then find that the n = 1 sector gives a negative contribution to the coefficient a 2,1 of s 2 t: where we used eq. (2.11). In other words, the scalar contribution dominates over the higher spin one in this sector. It is straightforward to generalize the argument to general n. g 2 nl for the first several orders are given in table 2. We find that for n = 3, 5, . . . higher spins dominate over scalars to obtain a positive contribution to a 2,1 . In figure 2, for illustration, we provide a plot for the coefficient of obtained after summing up over n from 1 to n max , which shows that the coefficient approaches to 0 asymptotically as we increase n max . To summarize, a 2,1 and therefore the Wilson coefficient β of the six derivative operator (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 vanish in the open superstring, which provides a concrete example for UV completion with higher spins and a non-positive β.

KK graviton
So far we have discussed the spin-dependence of the coefficients of the six-derivative operator based on the Froissart-Martin type bound on top of the unitarity and analyticity of scattering amplitudes. An important assumption there was that gravity is subdominant as we noted below eq. (2.1). Even before looking at the six-derivative operator, it is known that the positivity of the s 2 coefficient is subtle in the presence of gravity because the t-channel graviton exchange diagram behaves as ∼ s 2 t in the forward limit t → 0 and dominates over contributions from massive states [44]. Similarly, at least naively, it is not clear if the spin-dependence of the s 2 t coefficient discussed in the previous subsections is applicable in the presence of gravity. 10 Therefore, it is nontrivial if the sign of the six-derivative operator can be used to probe spins of heavy states in models with the Kaluza-Klein (KK) graviton, where the massless graviton inevitably appears. For this reason, in this subsection we study the six-derivative operator generated by the KK graviton by evaluating the four-point scattering amplitude explicitly. What characterizes the KK graviton χ µν is its coupling to the energy-momentum tensor T µν of other particles. Its 4D Lagrangian reads

JHEP12(2019)107
where we introduced χ = χ µ µ and g parametrizes the coupling of the KK graviton and the energy-momentum tensor. The dots stand for interaction terms with two or more χ, which are not relevant for our purpose. The kinetic operator E αβ µν is defined by Note that the kinetic and mass terms are the standard Fierz-Pauli ones (see, e.g., [62][63][64]). This type of Lagrangian appears, e.g., in the Randall-Sundrum (RS) I scenario [65], which has a discrete spectrum of the KK gravitons. Inflaton in this context has been studied in [66][67][68][69][70][71][72][73] and the oscillatory features of non-Gaussianities associated with the on-shell KK graviton creation were studied recently in [37]. We emphasize that the mass m and the coupling g are model-dependent, but the existence of the interaction χ µν T µν is universal due to the KK graviton nature. We then study the effective interactions of massless scalars mediated by the KK graviton. For this, we use the energy-momentum tensor, where we dropped terms arising from interaction terms because they are not relevant for our argument based on four-point amplitudes. Recall that the propagator of the KK graviton is given by (see, e.g., [62][63][64]) with the projector, Here we used a normalized symmetrizer A (ab) = 1 2 (A ab + A ba ). We then obtain the fourpoint amplitude of identical massless scalars, where u = −(s + t). Note that the amplitude behaves as ∼ s 2 for large s, so that it requires an appropriate UV completion at some scale. On the other hand, the IR expansion reads We find that both of s 2 and s 2 t have a positive coefficient. Interestingly, the sign happens to be the same as the massive spin 2 exchange discussed in section 2.2. Therefore, a nonpositive coefficient β implies a spinning heavy state also in the KK graviton case. We again

JHEP12(2019)107
emphasize that the previous argument cannot directly be applied to the KK graviton case at least naively because of the t-channel graviton pole mentioned earlier, so that a separate analysis was required. 11 To conclude this section, a nonpositive coefficient β of the six derivative operator (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 is a signature of intermediate spinning particles with even spin l = 2, 4, . . .. This conclusion is applicable to both (a) non-gravitational theories which respect the Froissart-Martin bound and (b) models with the KK graviton. See also table 1 for a summary of this section.

De Sitter four-point functions
In the rest of the paper we apply our argument on sign of the six derivative operator to cosmological settings. We start from the effective action of a massless scalar φ with a shift symmetry on exact de Sitter space: where the dots stand for operators with more derivatives and/or φ. As we have discussed, the four-derivative operator has to be always positive α > 0, whereas the sign of β tells us which of intermediate heavy scalars and spinning states are dominant. Here one might wonder if we can apply our argument based on the flat space scattering directly to the de Sitter case. Indeed, the counting of derivatives becomes ambiguous due to curvature couplings. For example, we may write a six-derivative operator, which vanishes on flat space. However, such effects are subdominant to the four-derivative operator α by a factor of H 2 /Λ 2 , hence the use of our flat space results is justified as long as H Λ. It is indeed the case because the de Sitter temperature is around the Hubble scale H and the validity of the EFT description requires H Λ. Note that thermal production of heavy particles are exponentially suppressed by the Boltzmann factor, so that their effects are negligible compared to the above mentioned curvature effect. 11 Indeed, the IR coefficients in eq. (2.29) are different from the ones in eqs. (2.10)-(2.11), even though the signs are the same. It does not imply any contradiction in general, but it could provide an interesting direction of the future studies: as we mentioned at footnote 10, we expect that the spin-dependence (2.11) of the six-derivative term will hold even in the presence of gravity, if the contribution of the Regge states is subdominant. If it is true, it would imply that the general results should hold if the KK graviton coupling g in eq. (2.23) is big enough (or if the KK scale is low enough) to dominate over the Regge states effect. Then, the Lagrangian (2.23) with a big g (or a small m), which is phenomenologically interesting, could be incomplete even as a low-energy EFT because of the mismatch. It would require other couplings such as the dilatonic coupling to resolve the mismatch for example. If the mismatch could never be resolved by such couplings, it would imply that models with the KK graviton with a big g (compared to the gravitational coupling) and a small m (compared to the mass scale of Regge states, i.e., the string scale in stringy UV completion) are in the swampland. We leave such a direction for the future work.

JHEP12(2019)107
Our question is now to clarify if we can distinguish the six-derivative operator from the leading-order four-derivative operator. For this purpose, in the rest of this section, we calculate late-time four-point functions generated by these two operators and discuss difference in the shape. Note that the four-point functions calculated in this section are applicable also to inflation at least in the regime where the special conformal symmetry is weakly broken. We also study inflationary three-point functions in the next section.

A brief review of the in-in formalism
We use the in-in formalism to calculate late-time n-point functions of φ in de Sitter space. With an interaction Hamiltonian H I (τ ), the late-time correlator of an operator O(τ ) is evaluated as (see, e.g, [74][75][76]) where T andT denote the time-ordering and anti-time-ordering operators, respectively. We choose the Bunch-Davies vacuum as the state |0 . For our purpose, the first order in H I (τ ) is enough to work with, so that we have where we assumed O is real. In the conformal time, the second order action reads where H is the (constant) Hubble parameter and the prime denote a derivative in conformal time τ . In the interaction picture, the massless scalar φ can then be quantized as with the standard commutation relations, The mode function follows the equation of motion, and it is normalized as For the Bunch-Davis vacuum, we have a mode function,

Comparison of the two shapes
We then investigate the shape of four-point functions with the Lagrangian (3.1). Since there are no three-point interactions, the interaction Hamiltonian relevant for us is simply a minus of the quartic Lagrangian: where the subscript 4 indicates that it contains four φ. It is straightforward to calculate the four-point function of φ by using eq. (3.4). Schematically, we write where A(k i ) and B(k i ) are associated with the four-derivative and six-derivative operators, respectively. Since the full shape is somewhat complicated to present, our shape analysis focuses on two particular limits of momentum configurations, the equilateral limit k 1 = k 2 = k 3 = k 4 and the flat space limit k 1234 = k 1 + k 2 + k 3 + k 4 → 0, which are useful to distinguish the two shapes.
Equilateral limit. As depicted in figure 3, momentum configurations in the equilateral limit are characterized by k 1 = k 2 = k 3 = k 4 = k, the angle γ between the two isosceles JHEP12(2019)107 triangles, and the angle θ between k 1 and −k 2 up to an overall rotation which does not affect scalar four-point functions. In this language we have k 1 · k 2 = −k 2 cos θ , k 1 · k 3 = k 2 · k 4 = k 2 cos θ + 1 2 cos γ + cos θ − 1 2 , Then, the shape functions A(k i ) and B(k i ) are evaluated as A(k i ) = H 8 256Λ 4 k 9 cos 2 θ(103 cos 2 γ + 309) − 206 cos θ sin 2 γ + 103 cos 2 γ + 173 , (3.14) B(k i ) = H 10 4096Λ 6 k 9 − 1311 cos 3 θ sin 2 γ − cos 2 θ(1111 cos 2 γ + 8577) + 6155 cos θ sin 2 γ − (3733 cos 2 γ + 3155) . We find that the θ-dependence is qualitatively different among the two: the four-derivative and six-derivative operators generate up to the second and the third order harmonics in θ, respectively. Therefore, the coefficient of cos 3θ can be used to probe the six-derivative operator without swamped by the four-derivative one. More generally, the n-th order harmonics in θ is sensitive to the 2n-derivative operators. It is essentially because mode functions of all external scalars and their conformal time derivatives do not depend on the angles θ and γ, and thus all the angular dependence is through spatial derivatives in the effective interactions. In this way, the equilateral configurations parameterized by θ and γ, and especially the θ-dependence of four-point functions is useful to separate out contributions from each order in the derivative expansion.
Flat space limit k 1234 → 0. Another interesting limit is the so-called flat space limit k 1234 = k 1 + k 2 + k 3 + k 4 → 0 with k i = 0, which is analogous to the flat space limit of AdS correlators proposed by Raju [77]. Note that this limit cannot be achieved by physical momentum configurations of dS correlators, hence it requires analytic continuation. By virtue of analytic continuation, we may probe deep inside the horizon τ → −∞ or in other words high-energy scattering on de Sitter. As demonstrated in [4,78], correlators in this limit are proportional to the corresponding flat space amplitudes. In our setup the flat space limit k 1234 → 0 reads where we used shorthand notations k 12 ≡ k 1 + k 2 , k 34 ≡ k 3 + k 4 , and k I = k 1 + k 2 . Note that the leading order of A gives the flat space scattering amplitude ∼ s 2 , if we identify k 12 and k 34 as the energy variable and k I as the three momentum. On the other hands, B gives the flat space scattering amplitude ∼ s 3 . More detailed analysis of the flat space limit is given in appendix B, to which we refer the readers with interests.

JHEP12(2019)107
4 Inflationary three-point functions Next we investigate inflationary three-point functions. We use the same inflaton Lagrangian as section 3 (with an appropriate slow-roll potential) and turn on a time-dependent inflaton backgroundφ(t) to write where δφ(t, x) denotes inflaton fluctuations. We assume thatφ is nearly constant under the slow-roll approximation (the dot denotes a derivative in physical time t as usual). Also, as we mentioned earlier, we focus on the regime the de Sitter conformal symmetry is weakly broken, so that we neglect higher order terms inφ in the following analysis. Under this assumption, the cubic interaction Hamiltonian is given by where the subscript 3 indicates that it contains three δφ. Since the cubic interaction is O(φ), corrections to the quadratic Lagrangian are not relevant for our purpose, so that we use the uncorrected linear equation of motion, and the corresponding mode functions for canonical quantization. When calculating threepoint functions, we may use eq. (4.3) to simplify the cubic Hamiltonian (4.2) as where the second line is a total time derivative that is not relevant in the following calculation. 12 Note that the first term of the β term in the first line is the same as the α one up to an overall constant. Now it is straightforward to calculate three-point functions of the scalar curvature perturbation ζ = −(Hδφ)/φ. We schematically write 5) 12 In the in-in formalism, it is sometimes dangerous to neglect total time derivatives because they may provide a non-vanishing boundary term in general. However, it is easy to see that the boundary terms we dropped vanish at the future boundary τ = 0 and thus we may safely focus on the first line of eq. (4.4). where A 3pt and B 3pt are contributions from the four-derivative and six-derivative operators, respectively. P ζ is the scalar power spectrum given by

JHEP12(2019)107
Next we discuss how to distinguish the three-point functions generated by the sixderivative operator from other contributions. First, as depicted in figure 4, both the leading contribution from the four-derivative operator and the next-to-leading one from the sixderivative operator have a peak at the equilateral configuration. Therefore we need to look at details of the shape. For this purpose, we may use, e.g., the ratio between the equilateral configuration and the folded one: It is also worth clarifying if the shape has degeneracy with other operators we have neglected. As we mentioned, the shape B 3pt is from the operator (δφ ) 3 , which arises from the six-point effective interaction (∂ µ φ∂ µ φ) 3 as well. One would then wonder if (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 and (∂ µ φ∂ µ φ) 3 are degenerate as long as we look at inflationary three-point functions. However, recall that (∂ µ φ∂ µ φ) 3 provides the interaction (δφ ) 3 at the cubic order inφ. To work at this order, we need to take into account corrections to the linear equation of motion (4.3) when simplifying the cubic Hamiltonian, which gives a new operator on top of the ones in eq. (4.4). Hence, the potential degeneracy of (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 and (∂ µ φ∂ µ φ) 3 are resolved by carefully studying higher order terms inφ which we neglected. Besides, these two operators generate different shapes of inflationary four-point functions because the four-point interaction of δφ originating from (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 is Lorentz invariant, but the one from (∂ µ φ) 6 is not.
To end this section, we estimate the nonlinearity parameter f N L and note the limitation of our present work for the near future phenomenology. First, f N L sourced by the fourderivative operator is estimated as so that the observable non-Gaussianity f N L 1 can be achieved only whenφ 2 Λ 4 . However, contributions from the six-point interaction γ Λ 8 (∂ µ φ∂ µ φ) 3 and higher are nonnegligible in the regimeφ 2 Λ 4 as long as we employ the ordinary order-estimation of the effective theory, i.e., if the dimensionless parameter γ is O(1) for example. Therefore, our argument based on the truncation at the four-point level is not applicable to the regime f N L 1 generically, which motivates further studies along the line of our present work.

Outlook
To enlarge the scope of the cosmological collider, we initiated a program decoding imprints of heavy particles from effective interactions of primordial perturbations. In this paper, as a first step in this direction, we studied spin-dependence of four-point effective interactions of the inflaton. It is well known that the sign of the four-derivative operator (∂ µ φ∂ µ φ) 2 is universally positive [44] as a consequence of unitarity and analyticity of scattering amplitudes together with the Froissart-Martin bound. In contrast to the universal positivity

JHEP12(2019)107
bound, we demonstrated that the sign of the six-derivative operator (∇ µ ∂ ν φ) 2 (∂ ρ φ) 2 is positive for intermediate scalars, whereas it is negative for intermediate spinning states. In particular, a non-positive coefficient requires heavy spinning particles with the spin s = 2, 4, . . . in tree-level UV completion. This result applies, e.g., to (1) non-gravitational theories enjoying the Froissart-Martin bound and (2) effective interactions generated by the KK graviton. We thus conclude that we may probe spin of heavy intermediate states from the sign of effective interactions by going beyond the positivity bound.
We also studied phenomenology of primordial non-Gaussianities thereof. Since the sixderivative operator is the next-to-leading order correction in the inflaton effective action, we need to explore how to distinguish it from the leading order correction, i.e., the fourderivative operator. First, we found that they are distinguishable from angular dependence of de Sitter four-point functions, especially in the equilateral limit. On the other hand, both operators generate three-point functions with a peak at the equilateral configuration, so that a detailed analysis of the full shape is required to distinguish the two from three-point functions. We again emphasize that the signals we studied are generated by intermediate off-shell heavy particles, so that it is free from the exponential Boltzmann suppression in contrast to the previous studies in the cosmological collider program [1][2][3][4].
There are various directions to explore along the line of our present work. First of all, the present paper focused on a regime where the de Sitter conformal symmetry is weakly broken to neglect six-point interactions such as (∂ µ φ∂ µ φ) 3 and higher. As we mentioned, primordial non-Gaussianities in this regime are generically small f N L 1 if we employ the usual order-estimation of the effective theory. It would be interesting to explore a concrete UV model which has an observable non-Gaussianity f N L 1, but negligible higher-point effective interactions. A more challenging, but important direction is to extend our analysis to the effective field theory of inflation [58], or equivalently to incorporate the higher-point inflaton effective interactions, which covers the regime f N L 1 in a general context. In this regime, we cannot utilize, e.g., the Lorentz invariance and crossing symmetry anymore (see [45] for related discussions), so that we would need another UV input such as Lorentz symmetry restoration at high energy [80]. It would also be interesting to generalize our work to include external spinning particles, e.g., for applications to graviton non-Gaussianities and phenomena other than inflation. We hope to report our progress in these directions elsewhere in the near future.

A Comments on loops
As we mentioned, the intermediate states generating effective interactions can be multiparticle states associated with loops. It will then be useful to demonstrate which sign of s 2 t appears in typical loop diagrams. For illustration, here we consider one-loop amplitudes with an internal massive scalar/fermion.

A.1 Scalar loop
Let us begin by the following two-scalar model: where σ is the massive particle to be integrated out and the interactions are parameterized by the couplings g 1 and g 2 . Four-point scattering amplitudes of φ at one-loop are then Here we defined the bubble, triangle, and box integrals by , (A.4) where the loop integral over l is already Wick rotated. Also the bubble integral (A.3) is defined in d = 4 − and the UV divergence ∼ 1/ has to be subtracted by a counterterm appropriately. Using the standard Feynman integrals, we may rewrite the integrals as where we introduced

JHEP12(2019)107
It is now easy to evaluate IR coefficients analytically by expanding the integrand in K i . The four-and six-derivative terms of the amplitude then read M (s, t) = 3g 4 1 − 8g 2 1 g 2 m 2 + 6g 2 2 m 4 5760π 2 m 6 s 2 + st + t 2 − 10g 4 1 − 27g 2 1 g 2 m 2 + 18g 2 2 m 4 80640π 2 m 10 st(s + t) . (A.11) Here and in the next subsection we suppress a constant piece, which can be eliminated by the φ 4 counterterm, as well as higher derivative terms. We find that the coefficient of s 2 is positive as required by unitarity. On the other hand, the sign of s 2 t depends on the ratio of g 2 1 and g 2 .

A.2 Fermion loop
We next consider a model of a massless (pseudo-)scalar φ and a fermion ψ: where Γ = 1 (Γ = iγ 5 ) when φ is a scalar (pseudo-scalar). The coefficient of s 2 is positive as required by unitarity, whereas the sign of s 2 t is negative, hence scalar intermediate states dominate over spinning ones.

JHEP12(2019)107
where η •• is the flat space metric and u k i (τ ) denotes the mode function (3.10) of φ. Also spatial derivatives should be understood as spatial momenta k i . Notice here that the integral (B.2) is a linear combination of the integrals, with (spatial) momentum-dependent coefficients, so that the highest power in n is dominant in the flat space limit k 1234 → 0. It implies that the contributions from connections are subleading compared with the normal derivative terms in the flat space limit. We therefore focus on the following part of the integral: . . ∂ µp u * k 1 (τ ))(∂ ν 1 . . . ∂ νq u * k 2 (τ )) × (∂ ρ 1 . . . ∂ ρr u * k 3 (τ ))(∂ σ 1 . . . ∂ σs u * k 4 (τ )) . (B.4) We may further simplify the integral in the flat space limit by using the high energy limit of the mode function, Under this approximation, we find where the last line is nothing but what we encounter in the computation of flat space amplitudes. In summary, the four-point correlator in the flat space limit yields comes from the contribution of external legs u k 1 (0)u k 2 (0)u k 3 (0)u k 4 (0).
The last line gives the flat space amplitude.