2

## Constraining graviton non-Gaussianity through the CMB bispectra

Valerio De Luca, 1 Gabriele Franciolini, 1 Alex Kehagias, 2 Antonio Riotto, 1, 3 and Maresuke Shiraishi 4

1 D´ epartement de Physique Th´ eorique and Centre for Astroparticle Physics (CAP),

Universit´ de e Gen` eve, 24 quai E. Ansermet, CH-1211 Geneva, Switzerland

Physics Division, National Technical University of Athens, 15780 Zografou Campus, Athens, Greece

3 CERN, Theoretical Physics Department, Geneva, Switzerland

4 Department of General Education, National Institute of Technology,

Kagawa College, 355 Chokushi-cho, Takamatsu, Kagawa 761-8058, Japan

(Dated: September 26, 2019)

Tensor non-Gaussianities are a key ingredient to test the symmetries and the presence of higher spin fields during the inflationary epoch. Indeed, the shape of the three point correlator of the graviton is totally fixed by the symmetries of the de Sitter stage and, in the case of parity conservation, gets contributions only from the ordinary gravity action plus a higher derivative term called the (Weyl) 3 action. We discuss current and future bounds on the three point tensor contribution from the (Weyl) 3 term using cosmic microwave background (CMB) bispectra. Our results indicate that forthcoming experiments, such as LiteBIRD, CMB-S4 and CORE, will detect the presence of the (Weyl) 3 term if M L 4 p 4 ∼ 10 17 r -4 , where L parametrizes the strength of the (Weyl) 3 term and r is the tensor-to-scalar ratio, which corresponds to L glyph[greaterorsimilar] 3 2 . × 10 5 M -1 p , while the current upper limit is M L 4 p 4 = (1 1 . ± 4 0) . × 10 19 r -4 (68%CL).

## I. INTRODUCTION

Current observational evidences suggest that the universe experienced a period of accelerated expansion in its early stages which goes under the name of cosmological inflation [1]. In this phase the spacetime geometry is approximately described by the de Sitter (dS) metric. One of its main unknown features is its particle content, which one could be able to specify looking at signatures in the present day observables with current/future experiments. Additionally, since inflation could have taken place at energies much higher than the ones reachable by collider experiments, it provides a unique opportunity to test high energy physics [2-6].

One very important property of inflation is that, by being described by an approximately de Sitter metric, it enjoys all its background symmetries. As we shall see, those symmetries are so powerful that allow to predict the properties of many observables. During this period, the spacetime isometries form an (approximate) SO(4,1) group which can be identified with the conformal group of a CFT 3 acting on perturbations extending on scales larger than the Hubble radius. One of the signatures of the inflationary epoch that could have reached us is given by the gravitational waves and, in particular, the graviton non-Gaussianities (NGs), which are dictated by symmetry arguments. Thus, any observational confirmation of the graviton NG would pave the way to a deeper understanding of the theory of gravity and particle content during inflation. Indeed, symmetry arguments (better to say the isometries of de Sitter) constrain the three point correlator (dubbed sometime the bispectrum) of the graviton to be conformal invariant. Its shape is fixed (when parity is preserved) by two interacting terms, the ordinary Einstein-Hilbert gravity action and the higherorder derivative action called the (Weyl) 3 term.

Since a sizable contribution of the (Weyl) 3 term has been related to possible causality violation [7], which are supposedly cured by the presence of a tower of higherspin states, an experimental confirmation of a large higher derivative tensor NG would suggest the presence of an infinite tower of higher spin (HS) particles during inflation. Indeed, one particularly interesting setup where the power of symmetries manifests itself is in theories with higher spin fields. The possible signature of the presence of massive (or partially massless) fields during inflation has been extensively studied both theoretically [8-16] and observationally [17-19]. In fact, in the massless case, there are very strong theorems constraining the existence of HS fields in Minkowski spacetime (see no-go theorems [20-22]) which can be avoided in curved spacetime settings [23]. A noticeable example is provided by the Vasiliev theory, which describes the dynamics of interacting massless higher spin particles in curved spacetime [24-26] and admits solutions which are asymptotically de Sitter (or anti de Sitter) [27-29]. In general, the theory contains an infinite tower of interacting massless gauge fields of any spin and, since it accommodates a massless spin-2 particle, the graviton, it is therefore an appealing theory of gravity.

The bispectrum of cosmic microwave background (CMB) is one of the cleanest observables of graviton NGs. So far, signatures of various NG shapes and their detectabilities have been theoretically examined (e.g. Refs. [30-40]). Moreover, some of these shapes have already been tested with the real CMB data [41-45]. No detection of &gt; 2 σ signal yields useful constraints on the background theories (see Ref. [46] for a review).

This paper aims at constraining the graviton NG induced by (Weyl) 3 term by employing the CMB bispectra. CMB signatures of the same NG template have already been analyzed in an early work [31], while how much they are constrained from current and forthcoming CMB ob-

servations remains undiscussed. 1 We here address these issues for the first time.

Through the temperature data analysis, we find no evidence of the graviton NG. This provides an upper bound on M L 4 p 4 as in Eq. (25). On the other hand, our forecasted analysis shows that its minimal detectable size shrinks by some orders of magnitudes by taking B-mode polarization into account. This encourages follow-up investigations in next decadal B-mode experiments such as LiteBIRD [47], CMB-S4 [48], and CORE [49].

The paper is organized as follows. In Sec. II we review the tensor NG shapes as dictated by the dS symmetries and discuss the particular example of the Vasiliev higher spin gravity giving rise to such a scenario. In Sec. III we present the observational test based on the currently available data and the forecasted sensitivity with the next generation of future experiments. The last section is devoted to the conclusions.

## II. GRAVITON NON-GAUSSIANITY IN DE SITTER SPACETIME

During the inflationary epoch, the spacetime is approximately described by a de Sitter metric

$$d s ^ { 2 } = \frac { 1 } { H ^ { 2 } \tau ^ { 2 } } \left ( - d \tau ^ { 2 } + d \vec { x } ^ { 2 } \right ), \text{ \quad \ \ } ( 1 )$$

which is invariant under the symmetry group SO(4,1). In the late-time limit ( τ → 0), the background isometries act as a CFT 3 group containing translations, rotations, dilatations and special conformal transformations. This also goes under the name of dS/CFT 3 correspondence. Under such a symmetry, the structure of the two and three point functions is completely fixed and their latetime limit is naturally captured by the correlator structure in the CFT 3 dual description. In the case of the graviton three point correlator, it can be shown that the symmetries dictate the shape of the bispectrum to be induced by an interacting action [50] composed by the Einstein-Hilbert term

$$S _ { \text{EH} } = \frac { M _ { p } ^ { 2 } } { 2 } \int d \tau d ^ { 3 } x \sqrt { - g } \left [ R - 6 H ^ { 2 } \right ], \quad ( 2 )$$

plus an additional (Weyl) 3 term as 2

$$S _ { W ^ { 3 } } = \frac { M _ { p } ^ { 2 } L ^ { 4 } } { 2 } \int d \tau d ^ { 3 } x \sqrt { - g } W ^ { a b } _ { \ c d } W ^ { c d } _ { \ m n } W ^ { m n } _ { \ a b }. \ ( 3 ) \ \frac { \text{from} } { \text{spin} } \, \text{g} \\ \text{const} \, \text{d} \, \text{d}.$$

1 In Ref. [31], a graviton NG from parity-violating (Weyl) 3 cubic term f ( φ WW ) ˜ 2 has also been studied. This has already been

constrained using the CMB temperature bispectrum in Ref. [42]. One additional possible contribution to the graviton three point function would come from the parity-violating interaction term From now on, we will consider only the parity-conserving

2 ˜ WW 2 . case.

The length scale L determines the strength of the (Weyl) 3 term. The explicit expression for the three point correlator is given in Refs. [50, 51]. The same structure is also found working in the three dimensional CFT 3 dual description [52-54].

One explicit example of the theory of gravity giving rise to such a graviton NG is provided by the Vasiliev theory in de Sitter spacetime, where the interactions between the particles are fully dictated by the higher spin gauge group and the spectrum is composed by an infinite tower of massless higher spin fields. The theory depends only on a free parameter identified with the ratio between the reduced Planck mass and the Hubble rate during inflation, M /H p . The prediction regarding the tree-level latetime cosmological correlation functions of the graviton in such a theory has been worked out in various works [5558], while the full nonlinear result can be computed with the aid of the framework introduced in Ref. [57] making use of the dS/CFT 3 correspondence. In particular, one can analyze the characteristic graviton non-Gaussianity through the three point function predicted by the higher spin gravity and obtained in Ref. [52]. The nice property of the higher spin gravity is that the interactions between the particles are completely fixed, thus also the additional scale L is a clear prediction of the theory. In particular, it is given by [52]

$$H L = \left ( \frac { 8 } { 2 7 0 } \right ) ^ { 1 / 4 } \simeq 0. 4. \text{ \quad \ \ } ( 4 )$$

This means that in the Vasiliev theory the higher derivative corrections to the GR graviton three point function are sizable.

Even though the structure of the graviton three point function is robustly fixed by symmetries, different theories of gravity could lead to the prediction of a different scale L which would impact the possibility of a future experimental confirmation through CMB observations. From now on, we will consider L as a free parameter.

## III. OBSERVATIONAL TESTS OF TENSOR NON-GAUSSIANITIES USING THE CMB BISPECTRA

Here we examine the detectability of a graviton NG from the (Weyl) 3 cubic action (3), as a probe of the higher spin gravity model, using CMB bispectra. Moreover, by constraining the graviton NG with the real CMB data, we test this model.

## A. CMB bispectra

Since the CMB fields are distributed on the 2 sphere, their magnitudes are quantified through the spherical harmonic decomposition. Tensor-mode harmonic coefficients of temperature and E/B-mode polarization fields

( X = T, E, B ) take the form [30, 59]

$$a _ { \ell m } ^ { X } = ( - i ) ^ { \ell } \int \frac { d ^ { 3 } \vec { k } } { 2 \pi ^ { 2 } } \mathcal { T } _ { \ell } ^ { X } ( k ) \sum _ { \lambda = \pm 2 } \left ( \frac { \lambda } { 2 } \right ) ^ { x } \gamma _ { \vec { k } } ^ { ( \lambda ) } - \lambda Y _ { \ell m } ^ { * } ( \hat { k } ), \\ t _ { \wedge \cdots \cdots } \, \nu _ { \cdots \cdots } \, \nu _ { \cdots }$$

$$\text{where } x = \begin{cases} 0 & ( X = T, E ) \\ 1 & ( X = B ) \\ \cdots & \dots & \dots \\ \end{cases} \, \begin{array} { c } \mathbb { \vee } \end{array}, \text{and } \gamma _ { \vec { k } } ^ { ( \pm 2 ) } \text{ is a helicity-state}$$

representation of primordial GW, which is given by

$$\gamma _ { i j } ( \vec { x } ) \equiv \int \frac { d ^ { 3 } \vec { k } } { ( 2 \pi ) ^ { 3 } } e ^ { i \vec { k } \cdot \vec { x } } \sum _ { \lambda = \pm 2 } \gamma _ { \vec { k } } ^ { ( \lambda ) } e _ { i j } ^ { ( \lambda ) } ( \hat { k } ), \quad ( 6 )$$

with the transverse and traceless polarization tensor e ( ± 2) ij ( ˆ ) k obeying e ( λ ) ∗ ij ( ˆ ) k = e ( -λ ) ij ( ˆ ) k = e ( λ ) ij ( -ˆ ) k and e ( λ ) ij ( ˆ ) k e ( λ ′ ) ij ( ˆ ) k = 2 δ λ, -λ ′ . The linear transfer function T X glyph[lscript] ( k ) expresses characteristic tensor-mode features, e.g., the integrated Sachs-Wolfe (ISW) amplification for temperature and the reionization bumps for polarization [60]. The CMB bispectrum then takes the form

$$\text{lariation} \, [ 0 0 ]. \, \text{$line\,\text{$mb$ bispectrum then takes the form} \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \
 \left \langle \prod _ { n = 1 } ^ { 3 } a _ { \ell _ { n } m _ { n } } ^ { X _ { n } } \right \rangle = \left [ \prod _ { n = 1 } ^ { 3 } \left ( - i \right ) ^ { \ell _ { n } } \int _ { 0 } ^ { \infty } \frac { k _ { n } ^ { 2 } d k _ { n } } { 2 \pi ^ { 2 } } \mathcal { T } _ { \ell _ { n } } ^ { X _ { n } } \left ( k _ { n } \right ) \right ] \quad \text{where} \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \	 \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdot \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \

 \text{A helicity-state expression of graviton bispectrum from}$$

Ahelicity-state expression of graviton bispectrum from the (Weyl) 3 cubic action (3) (and also a formula of induced CMB bispectrum) has already been obtained in Ref. [31]. Here, we present a mathematically equivalent but alternative form where contractions between unit vectors and polarization tensors are much simplified, reading

$$\left \langle \prod _ { n = 1 } ^ { 3 } \gamma _ { \vec { k } _ { n } } ^ { ( \lambda _ { n } ) } \right \rangle & = ( 2 \pi ) ^ { 3 } \delta ^ { ( 3 ) } \left ( \sum _ { n = 1 } ^ { 3 } \vec { k } _ { n } \right ) f _ { k _ { 1 } k _ { 2 } k _ { 3 } } \delta _ { \lambda _ { 1 }, \lambda _ { 2 } } \delta _ { \lambda _ { 2 }, \lambda _ { 3 } } \\ & \quad \times e _ { i j } ^ { ( - \lambda _ { 1 } ) } ( \hat { k } _ { 1 } ) e _ { j k } ^ { ( - \lambda _ { 2 } ) } ( \hat { k } _ { 2 } ) e _ { k i } ^ { ( - \lambda _ { 3 } ) } ( \hat { k } _ { 3 } ), \quad ( 8 ) \\ & \quad \quad \text{T h o o }.$$

where

$$f _ { k _ { 1 } k _ { 2 } k _ { 3 } } \equiv 5 7 6 0 \left ( \frac { H } { M _ { p } } \right ) ^ { 8 } M _ { p } ^ { 4 } L ^ { 4 } ( k _ { 1 } + k _ { 2 } + k _ { 3 } ) ^ { - 6 }. \quad ( 9 ) \quad \text{mic} \\ \text{natur}$$

In a similar way to Ref. [31], inserting this into Eq. (7), we compute angular integrals and helicity summations, and simplify the resulting formula by employing the addition rule of angular momenta. The bottom-line formula has a rotational-invariant shape

$$\left \langle \prod _ { n = 1 } ^ { 3 } a _ { \ell _ { n } m _ { n } } ^ { X _ { n } } \right \rangle = \left ( \begin{array} { c c c } \ell _ { 1 } & \ell _ { 2 } & \ell _ { 3 } \\ m _ { 1 } & m _ { 2 } & m _ { 3 } \end{array} \right ) B _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { X _ { 1 } X _ { 2 } X _ { 3 } }, \quad ( 1 0 ) \quad \text{of $ott$} _ { \ell \gtrsim 1$}$$

FIG. 1. Intensity distribution of the CMB temperature bispectrum in 3D glyph[lscript] -space domain. Each axis corresponds to glyph[lscript] 1 , glyph[lscript] 2 , or glyph[lscript] 3 . Red (blue) dense regions correspond to large positive (negative) signal. Here, the bispectrum signal normalized by a constant Sachs-Wolfe bispectrum template [61] is plotted and hence dominant configurations are effectively highlighted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00366v2-figure-b999f892-261b-4e2a-9028-0a687f067b7a.jpg)

where

$$& \text{where} \\ & B _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { X _ { 1 } } X _ { 2 } X _ { 3 } } = \frac { - ( 8 \pi ) ^ { 3 / 2 } } { 5 } \sqrt { \frac { 7 } { 3 } } ( - i ) ^ { \ell _ { 1 } + \ell _ { 2 } + \ell _ { 3 } } \\ & \times \delta _ { x _ { 1 } + x _ { 2 } + x _ { 3 } + \ell _ { 1 } + \ell _ { 2 } + \ell _ { 3 } } \sum _ { L _ { 1 } L _ { 2 } L _ { 3 } } ( - 1 ) ^ { \frac { L _ { 1 } + L _ { 2 } + L _ { 3 } } { 2 } } h _ { L _ { 1 } L _ { 2 } L _ { 3 } } ^ { 0 } \, 0 \, 0 \, 0 \\ \text{from} \quad & \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }new} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\mathbb{ }n_{1}$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$ } \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\ \text{$\nu$} \\
\end{cases}$$

$$\text{i-} \text{ with } \delta _ { l } ^ { ( e ) } \equiv \begin{cases} 0 & ( l = \text{odd} ) \\ 1 & ( l = \text{even} ) \end{cases} \text{ and}$$

$$( 8 ) \quad & h _ { l _ { 1 } l _ { 2 } l _ { 3 } } ^ { s _ { 1 } s _ { 2 } s _ { 3 } } \equiv \sqrt { \frac { ( 2 l _ { 1 } + 1 ) ( 2 l _ { 2 } + 1 ) ( 2 l _ { 3 } + 1 ) } { 4 \pi } } \left ( \begin{array} { c c c } l _ { 1 } & l _ { 2 } & l _ { 3 } \\ s _ { 1 } & s _ { 2 } & s _ { 3 } \end{array} \right ). \\ & \text{The model command line to run } ( 1 1 \times \text{too a much } \text{$\mu$} \end{array}$$

The angle-averaged bispectrum (11) takes a much simpler form than the previous one in Ref. [31] (although both yield identical result); hence, much more speedy numerical analysis is realized. Note that parity-conserving nature of the graviton bispectrum (8) restricts nonvanishing signal to TTT , TTE TEE TBB EEE , , , , and EBB when glyph[lscript] 1 + glyph[lscript] 2 + glyph[lscript] 3 = even, and TTB TEB EEB , , , and BBB when glyph[lscript] 1 + glyph[lscript] 2 + glyph[lscript] 3 = odd.

Intensity distribution of the temperature bispectrum in 3D glyph[lscript] -space domain is plotted in Fig. 1. One can confirm there that the dominant signal comes from equilateral configurations glyph[lscript] 1 ∼ glyph[lscript] 2 ∼ glyph[lscript] 3 (see Ref. [46] for illustrations of other shapes). However, it is highly suppressed for glyph[lscript] glyph[greaterorsimilar] 100 because of the cease of the ISW amplification.

The other bispectra including E-mode and/or B-mode polarization, that are not described here, are also suppressed at small scales (because of the damping nature of the GW induced polarization field) and hence have similar intensity distributions.

## B. Detectability

In what follows, we focus on estimating a bispectrum amplitude parameter M L 4 p 4 using the CMB bispectra (11). Then, assuming

$$\left ( \frac { H } { M _ { p } } \right ) ^ { 2 } = \frac { \pi ^ { 2 } } { 2 } r A _ { S } \, \quad \ \ \ ( 1 3 ) \, \quad \ \ \ \ t r a \,$$

with A S the scalar amplitude and r the tensor-to-scalar ratio (see e.g. Ref. [1]), we rewrite f k k k 1 2 3 into

$$f _ { k _ { 1 } k _ { 2 } k _ { 3 } } \equiv 3 6 0 \pi ^ { 8 } A _ { S } ^ { 4 } r ^ { 4 } M _ { p } ^ { 4 } L ^ { 4 } ( k _ { 1 } + k _ { 2 } + k _ { 3 } ) ^ { - 6 }. \quad ( 1 4 ) \quad \text{the ir}$$

We also fix A S to be an observed value, 2 2 . × 10 -9 [62, 63], and analyze the detectability of M L 4 p 4 by varying r .

The Fisher matrices computed from auto bispectrum only ( T -only, E -only or B -only) and temperature and E-mode polarization jointly ( T + E ) read, respectively,

$$F _ { X } = \sum _ { \substack { \ell _ { i } \\ 1 } } ( - 1 ) ^ { \ell _ { 1 } + \ell _ { 2 } + \ell _ { 3 } } \frac { \left ( \tilde { B } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { X X X } \right ) ^ { 2 } } { 6 C _ { \ell _ { 1 } } ^ { X X } C _ { \ell _ { 2 } } ^ { X X } C _ { \ell _ { 3 } } ^ { X X } } \, \quad \ \ ( 1 5 ) \quad \text{$\text{$with $C_{i}$}$}$$

$$F _ { T + E } = \frac { 1 } { 6 } \sum _ { X _ { i }, X _ { i } ^ { \prime } = T, E } \sum _ { \ell _ { i } } ( - 1 ) ^ { \ell _ { 1 } + \ell _ { 2 } + \ell _ { 3 } } \tilde { B } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { X _ { 1 } X _ { 2 } X _ { 3 } } ( C _ { \ell _ { 1 } } ^ { - 1 } ) ^ { X _ { 1 } X _ { 1 } ^ { \prime } } \quad \text{$On$th$} \quad \text{and$n$} \quad \text{of the} \quad \text{$one$} \quad \text{$one$} \quad \text{$one$}$$

$$\times ( C _ { \ell _ { 2 } } ^ { - 1 } ) ^ { X _ { 2 } X _ { 2 } ^ { \prime } } ( C _ { \ell _ { 3 } } ^ { - 1 } ) ^ { X _ { 3 } X _ { 3 } ^ { \prime } } \tilde { B } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { X _ { 1 } ^ { \prime } X _ { 2 } ^ { \prime } X _ { 3 } ^ { \prime } }, \quad \ \ ( 1 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where C XX ′ glyph[lscript] is the CMB angular power spectrum, and ˜ B X X X 1 2 3 glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 = B X X X 1 2 3 glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 / M L ( 4 p 4 ). Here, we have assumed that any NG contribution to the covariance matrix is negligibly small; thus, the covariance is diagonalized and given by the power spectrum alone. This would be justified as long as the Fisher matrices do not contain the information at very small angular scales. According to the Planck team [45], the diagonal covariance matrix approximation works for glyph[lscript] glyph[lessorsimilar] 2000 by using Planck component separated temperature and polarization maps. For safety, our Fisher matrix computations are performed for glyph[lscript] ≤ 1000. The expected 1 σ errors on M L 4 p 4 from T -only, E -only, T + E , and B -only are given by ∆ M L 4 p 4 | T = 1 / √ f sky F T , ∆ M L 4 p 4 | E = 1 / √ f sky F E , ∆ M L 4 p 4 | T + E = 1 / √ f sky F T + E and ∆ M L 4 p 4 | B = 1 / √ f sky F B , respectively, where f sky is a sky coverage fraction. The Bmode power spectrum in F B is modeled by the sum of primordial contribution proportional to r , lensing Bmode one, and experimental noise spectrum, i.e., C BB glyph[lscript] = C prim glyph[lscript] ( r ) + C lens glyph[lscript] + N glyph[lscript] . Here, we consider four cleanliness levels of the B-mode data: a perfectly-delensed and noiseless full-sky case (i.e., C lens glyph[lscript] = N glyph[lscript] = 0 and glyph[negationslash]

glyph[negationslash]

glyph[negationslash]

f sky = 1), a nondelensed and noiseless full-sky one (i.e., C lens glyph[lscript] = 0, N glyph[lscript] = 0 and f sky = 1), a perfectly-delensed and noisy partial-sky one (i.e., C lens glyph[lscript] = 0, N glyph[lscript] = 0 and f sky = 0 5), and a nondelensed and noisy partial-sky one . (i.e., C lens glyph[lscript] = 0, N glyph[lscript] = 0 and f sky = 0 5), are taken into . account. Since an achievable level of delensing in the bispectrum analysis with future B-mode data has not been fixed, we simply examine two extreme cases: the perfectly-delensed and nondelensed ones. The noise spectrum is computed with experimental uncertainties due to beam, noise, mask and residual foreground expected in the LiteBIRD project [35, 47, 64, 65]. 3 As was done in Ref. [35], we assume that foregrounds arising from galactic dust emission and synchrotron radiation are subtracted by use of 9 channels planned in LiteBIRD (4089 GHz and 280-402 GHz), achieving 2% residual level in CMB maps. On the other hand, we assume a noiseless full-sky measurement of temperature and E-mode polarization and therefore C TT glyph[lscript] , C EE glyph[lscript] and C TE glyph[lscript] do not contain the information of instrumental noise.

glyph[negationslash]

The left panel of Fig. 2 depicts ∆ M L 4 p 4 as a function of glyph[lscript] max at r = 10 -3 . One can see there that ∆ M L 4 p 4 | T , ∆ M L 4 p 4 | E , ∆ M L 4 p 4 | T + E , and ∆ M L 4 p 4 | B of the nondelensed case saturate at glyph[lscript] &gt; 100. This is because the tensor bispectrum (the numerator of F ) rapidly decays for glyph[lscript] &gt; 100, where the covariance (the denominator of F ) is still large due to the scalar-mode signal in C TT glyph[lscript] , C EE glyph[lscript] and C TE glyph[lscript] , and the lensing signal in C BB glyph[lscript] . The similar behavior is seen in ∆ M L 4 p 4 | B of the perfectly-delensed and noisy case because of the dominance of N glyph[lscript] for glyph[lscript] glyph[greaterorsimilar] 100. On the other hand, ∆ M L 4 p 4 | B of the perfectly-delensed and noiseless case can shrink as glyph[lscript] max increases because of the absence of the lensing signal and noise in C BB glyph[lscript] .

The right panel of Fig. 2 presents ∆ M L 4 p 4 as a function of r at glyph[lscript] max = 1000. In the perfectly-delensed and noiseless case, ∆ M L 4 p 4 | B ∝ r -5 / 2 holds for any r since ˜ B glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 ∝ r 4 and C BB glyph[lscript] = C prim glyph[lscript] ( r ) ∝ r . Regarding the other error bars, since C TT glyph[lscript] ∝ r 0 , C EE glyph[lscript] ∝ r 0 and

$$C _ { \ell } ^ { B B } \simeq \begin{cases} C _ { \ell } ^ { \text{prim} } ( r ) \in r & \text{for large $r$} \\ C _ { \ell } ^ { \text{lens} } + N _ { \ell } & \text{for small $r$} \end{cases}, \quad ( 1 7 )$$

the r dependences are expected to be

$$\Delta M _ { p } ^ { 4 } L ^ { 4 } | _ { T } \in r ^ { - 4 },$$

$$\Delta M _ { p } ^ { 4 } L ^ { 4 } | _ { E } \in r ^ { - 4 },$$

$$\Delta M _ { p } ^ { 4 } L ^ { 4 } | _ { T + E } \in r ^ { - 4 },$$

$$\Delta M _ { p } ^ { 4 } L ^ { 4 } | _ { B } \in \begin{cases} r ^ { - 5 / 2 } & \text{for large $r$} \\ r ^ { - 4 } & \text{for small $r$} \end{cases}. \quad ( 2 1 )$$

These behaviors can be confirmed from each line in the right panel. We find there that, at any level of delensing,

3 Comparable cleanliness level will also be achieved by other proposed CMB missions as CMB-S4 [48] and CORE [49].

FIG. 2. Expected 1 σ errors: ∆ M L 4 p 4 from T -only (blue lines), E -only (yellow lines), T + E (red lines) and B -only (cyan lines), as a function of the maximum multipole number glyph[lscript] max (left panel) and the tensor-to-scalar ratio r (right panel). The results in the left and right panels are computed with r = 10 -3 and glyph[lscript] max = 1000, respectively. Each line style of B -only (cyan line) discriminates the cleanliness level of the B-mode data: a perfectly-delensed and noiseless full-sky case (dot-dashed lines), a nondelensed and noiseless full-sky one (solid lines), a perfectly-delensed and LiteBIRD-like noisy one (dotted lines), and a nondelensed and LiteBIRD-like noisy one (dashed lines). For temperature and E-mode polarization, we assume a noiseless full-sky measurement.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00366v2-figure-b5536f5b-086b-43cd-b188-26e2808209ac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00366v2-figure-d49b02a2-e564-4123-8ac0-478213218186.jpg)

the graviton NG with

$$M _ { p } ^ { 4 } L ^ { 4 } \sim 1 0 ^ { 1 7 } r ^ { - 4 } \quad \ \ ( 2 2 ) \quad \ \ s i a$$

is detectable by next decadal B-mode survey aiming at hunting r &lt; 10 -2 such as LiteBIRD, CMB-S4 and CORE, and a smaller signal could be captured further in the future.

## C. Observational constraints

Finally, we test this model with existing CMB data. We then follow a bispectrum estimation pipeline adopted in previous works on constraining other tensor NG shapes [41-45].

is the bispectrum reconstructed from maps, and 〈 a T glyph[lscript] 1 m 1 a T glyph[lscript] 2 m 2 〉 MC denotes the variance of simulated Gaussian maps. This estimator properly works when the covariance matrix is approximately diagonalized. Such a situation is achieved by a recursive inpainting prefiltering technique [69, 70], which has been implemented in map making for the Planck bispectrum estimation in order to avoid a time-consuming full-inverse covariance weighting process [43, 45]. At the price of minimal ( ∼ 5%) loss of optimality, we adopt this. Note that estimates from Eq. (23) rely on r , which however is singled out because of ̂ M L 4 p 4 ∝ r -4 .

Here, M L 4 p 4 is constrained from a coadded temperature data using WMAP 9-years V and W bands [66, 67]. The Planck data is not adopted for this analysis, while we expect from previous analyses of a very similar equilateral NG shape [42, 43, 45] that almost an identical limit is obtained from the temperature data, and it is slightly tightened by adding E-mode polarization data. We then employ an estimator [68]

$$\widehat { M _ { p } ^ { 4 } L ^ { 4 } } = \frac { 1 } { F _ { T } } \sum _ { \ell _ { i } } ( - 1 ) ^ { \ell _ { 1 } + \ell _ { 2 } + \ell _ { 3 } } \frac { \tilde { B } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { T T T } \mathcal { B } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } ^ { T T T } } { 6 C _ { \ell _ { 1 } } ^ { T T } C _ { \ell _ { 2 } } ^ { T T } C _ { \ell _ { 3 } } ^ { T T } }, \quad ( 2 3 ) \quad \text{an u} \quad \text{map}$$

where

$$\mathcal { B } ^ { T T T } _ { \ell _ { 1 } \ell _ { 2 } \ell _ { 3 } } \equiv & \sum _ { m _ { i } } \left ( \begin{array} { c c c } \ell _ { 1 } & \ell _ { 2 } & \ell _ { 3 } \\ m _ { 1 } & m _ { 2 } & m _ { 3 } \end{array} \right ) \left [ a ^ { T } _ { \ell _ { 1 } m _ { 1 } } a ^ { T } _ { \ell _ { 2 } m _ { 2 } } a ^ { T } _ { \ell _ { 3 } m _ { 3 } } \quad \text{ use } \text{ as } \text{$\sum$} \\ & - \left ( \left \langle a ^ { T } _ { \ell _ { 1 } m _ { 1 } } a ^ { T } _ { \ell _ { 2 } m _ { 2 } } \right \rangle _ { \text{MC} } a ^ { T } _ { \ell _ { 3 } m _ { 3 } } + 2 \text{ perm.} \right ) \right ] \left ( 2 4 \right ) \quad \text{ is giving}$$

The best-fit value and the 1 σ errors on M L 4 p 4 are obtained computing Eq. (23) with observed data and 500 Gaussian simulations, respectively. Experimental features, i.e., beam shape and anisotropic noise, are added to the simulations by means of the WMAP-team methodology [71]. Moreover, both of the observed data and the simulations are masked using the KQ75 mask, leading to f sky = 0 688. . All masked regions are inpainted after the removal of monopole and dipole components. The observed maps and the experimental tools are downloaded from the Lambda website [72].

The above estimation process, unfortunately, comes at an unfeasibly high computational cost since a map-bymap computation of Eq. (23) already requires O ( glyph[lscript] 5 max ) numerical operations. This issue is actually solved if the summation in the 3D glyph[lscript] -space domain is factorized as ∑ glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 A glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 = [ ∑ glyph[lscript] 1 a glyph[lscript] ( 1 )][ ∑ glyph[lscript] 2 b glyph[lscript] ( 2 )][ ∑ glyph[lscript] 3 c glyph[lscript] ( 3 )]. In fact, in the usual scalar-mode bispectrum analysis, the estimator sum is reduced straightforwardly since ˜ B TTT glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 is given by a (partially) separable form [68]. In general

situations where ˜ B TTT glyph[lscript] 1 glyph[lscript] 2 glyph[lscript] 3 has a nonseparable structure (see Refs. [30-40] for tensor-mode NG cases), the modal approach [61, 73-76] gives a solution. There, a nonseparable input bispectrum template is replaced by an approximated separable one generated via a modal decomposition based on finite numbers of separable eigenfunctions. This yields a fast estimator form imposing O ( glyph[lscript] 3 max ) numerical operations. In this paper, the modal approach is adopted because the input bispectrum under examination (11) is obviously nonseparable. We perform the modal decomposition using 450 modal basis templates based on polynomials, and a 99 9% reproduction of the . input bispectrum is achieved.

Prior to working with the real data, we test the validity of our bispectrum estimation pipeline using simulations. We then produce O (10 ) 2 NG maps including nonzero M L 4 p 4 by means of the modal approach [61, 73, 74]. Implementing all WMAP experimental features mentioned above in them, we do a map-by-map computation of Eq. (23), and confirm that the average of output values recovers the input M L 4 p 4 . In a similar manner, we also compute Eq. (23) from O (10 ) Gaussian simulations, 2 and confirm that the standard deviation of output values saturates a Cram´ er-Rao bound, 1 / √ f sky F T (see Fig. 2 for values in a noiseless full-sky measurement).

Since all validation tests are successfully passed as above stated, the real data are implemented in our pipeline. An analysis at a maximum angular resolution, glyph[lscript] max = 500, yields

$$M _ { p } ^ { 4 } L ^ { 4 } = ( 1. 1 \pm 4. 0 ) \times 1 0 ^ { 1 9 } r ^ { - 4 } \ \ ( 6 8 \% \, \text{CL} ). \quad ( 2 5 )$$

This comes from the foreground-cleaned data, while the raw data result in a comparable value. This result is unbiased by primordial scalar-mode bispectrum or secondary one (i.e., ISW-lensing bispectrum [77, 78]) because it has a shape completely different from our template (11). These make our result robust.

## D. Physical interpretations

The prediction of the Vasiliev HS theory for the parameter L is L glyph[similarequal] 0 4 . /H , as reported in Eq. (4), from which it is clear that the signal is enhanced for a smaller value of the Hubble rate H . We can also use the relation between the tensor-to-scalar ratio and H in Eq. (13) as

$$r \simeq 9 \times 1 0 ^ { 7 } \left ( \frac { H } { M _ { p } } \right ) ^ { 2 } \quad \quad ( 2 6 ) \quad \stackrel { \infty } { \text{vat} } \\ \text{the} }$$

to rewrite the condition allowing the signal to be above the upper bound in Eq. (25) as

$$H \gtrsim 2 \times 1 0 ^ { - 3 } M _ { p } \simeq 5 \times 1 0 ^ { 1 5 } \, \text{GeV}, \quad \ ( 2 7 ) \quad \ \.$$

which is already ruled out by current constraints on the Hubble rate H &lt; 6 × 10 13 GeV [79].

With future experiments, using Eq. (22), the observability limit will be set to

$$H \gtrsim 5 \times 1 0 ^ { - 4 } M _ { p } \simeq 1 \times 1 0 ^ { 1 5 } \, \text{GeV}. \quad \ ( 2 8 )$$

In light of this, one can conclude that no observational confirmation of the Vasiliev HS gravity can be expected within the next generation of future experiments. This difficulty in observing such a signature in the specific case of the Vasiliev higher spin gravity theory is due to its specific prediction of the combined product HL glyph[similarequal] 0 4. . One could in principle hope that inflation took place with a lower Hubble scale H , enhancing the effect of the (Weyl) 3 interaction onto the graviton NG, but that, in turn, would decrease the amount of tensor modes reaching us (thus the dependence of the expected 1 σ error on r ) making any observation out of reach even for the planned future experiments.

However, a different theory of gravity predicting a different contribution from the (Weyl) 3 term to the three point function would turn out to be testable, provided that the length scale L is predicted to be at least

$$H L \gtrsim 8,$$

where again we utilized the future detectability bound shown in Eq. (22). Given the present bound on the Hubble rate [79], then this condition translates into L glyph[greaterorsimilar] 3 2 . × 10 5 M -1 p . The nonobservation of such a signature would, on the contrary, put an upper bound on the scale L as L glyph[lessorsimilar] 8 /H .

## IV. CONCLUSIONS

We have investigated the capability of measured CMB spectra to provide information on the strength of the (Weyl) 3 term through the tensor NGs generated during inflation. Such term is predicted to be relevant for the three point graviton correlator since it provides a conformal invariant tensor spectrum, together with the standard GR action.

If the strength of the (Weyl) 3 term is fixed by a length scale L , we have concluded that, with the currently available observations, an upper limit is set at M L 4 p 4 = (1 1 . ± 4 0) . × 10 19 r -4 (68%CL), while the new generation of CMB experiments will be able to detect the presence of the (Weyl) 3 term if M L 4 p 4 glyph[greaterorsimilar] 10 17 r -4 , or HL glyph[greaterorsimilar] 8. Unfortunately the well theoretically motivated higher spin theory ` la a Vasiliev does not satisfy such a bound and therefore will not be testable through the tensors NGs in the forthcoming observations.

## ACKNOWLEDGMENTS

V.D.L., G.F. and A.R. are supported by the Swiss National Science Foundation (SNSF), project The Non-Gaussian Universe and Cosmological Symmetries ,

project number: 200020-178787. M.S. is supported by JSPS Grant-in-Aid for Research Activity Start-up Grant Number 17H07319 and JSPS Grant-in-Aid for EarlyCareer Scientists Grant Number 19K14718. M. S. also ac-

- [1] D. H. Lyth and A. Riotto, Phys. Rept. 314 , 1 (1999), arXiv:hep-ph/9807278 [hep-ph].
- [2] X. Chen and Y. Wang, Phys. Rev. D81 , 063511 (2010), arXiv:0909.0496 [astro-ph.CO].
- [3] X. Chen and Y. Wang, JCAP 1004 , 027 (2010), arXiv:0911.3380 [hep-th].
- [4] A. Kehagias and A. Riotto, Fortsch. Phys. 63 , 531 (2015), arXiv:1501.03515 [hep-th].
- [5] N. Arkani-Hamed and J. Maldacena, (2015), arXiv:1503.08043 [hep-th].
- [6] N. Bartolo and G. Orlando, JCAP 1707 , 034 (2017), arXiv:1706.04627 [astro-ph.CO].
- [7] X. O. Camanho, J. D. Edelstein, J. Maldacena, and A. Zhiboedov, JHEP 02 , 020 (2016), arXiv:1407.5597 [hep-th].
- [8] L. Bordin, P. Creminelli, M. Mirbabayi, and J. Norea, JCAP 1609 , 041 (2016), arXiv:1605.08424 [astroph.CO].
- [9] H. Lee, D. Baumann, and G. L. Pimentel, JHEP 12 , 040 (2016), arXiv:1607.03735 [hep-th].
- [10] A. Kehagias and A. Riotto, Fortsch. Phys. 65 , 1700023 (2017), arXiv:1701.05462 [hep-th].
- [11] A. Kehagias and A. Riotto, JCAP 1707 , 046 (2017), arXiv:1705.05834 [hep-th].
- [12] D. Baumann, G. Goon, H. Lee, and G. L. Pimentel, JHEP 04 , 140 (2018), arXiv:1712.06624 [hep-th].
- [13] G. Franciolini, A. Kehagias, and A. Riotto, JCAP 1802 , 023 (2018), arXiv:1712.06626 [hep-th].
- [14] L. Bordin, P. Creminelli, A. Khmelnitsky, and L. Senatore, JCAP 1810 , 013 (2018), arXiv:1806.10587 [hep-th].
- [15] N. Arkani-Hamed, D. Baumann, H. Lee, and G. L. Pimentel, (2018), arXiv:1811.00024 [hep-th].
- [16] G. Goon, K. Hinterbichler, A. Joyce, and M. Trodden, (2018), arXiv:1812.07571 [hep-th].
- [17] N. Bartolo, A. Kehagias, M. Liguori, A. Riotto, M. Shiraishi, and V. Tansella, Phys. Rev. D97 , 023503 (2018), arXiv:1709.05695 [astro-ph.CO].
- [18] G. Franciolini, A. Kehagias, A. Riotto, and M. Shiraishi, Phys. Rev. D98 , 043533 (2018), arXiv:1803.03814 [astroph.CO].
- [19] A. Moradinezhad Dizgah, G. Franciolini, A. Kehagias, and A. Riotto, Phys. Rev. D98 , 063520 (2018), arXiv:1805.10247 [astro-ph.CO].
- [20] S. Weinberg, Phys. Rev. 135 , B1049 (1964).
- [21] S. R. Coleman and J. Mandula, Hadrons and their Interactions: Current and Field Algebra, Soft Pions, Supermultiplets, and Related Topics , Phys. Rev. 159 , 1251 (1967).
- [22] S. Weinberg and E. Witten, Phys. Lett. 96B , 59 (1980).
- [23] X. Bekaert, N. Boulanger, and P. Sundell, Rev. Mod. Phys. 84 , 987 (2012), arXiv:1007.0435 [hep-th].
- [24] M. A. Vasiliev, Phys. Lett. B243 , 378 (1990).
- [25] M. A. Vasiliev, Phys. Lett. B285 , 225 (1992).
- [26] M. A. Vasiliev (2000) arXiv:hep-th/0002183 [hep-th].

knowledges the Center for Computational Astrophysics, National Astronomical Observatory of Japan, for providing the computing resources of Cray XC50.

- [27] M. A. Vasiliev, Fortsch. Phys. 35 , 741 (1987), [Yad. Fiz.45,1784(1987)].
- [28] C. Iazeolla, E. Sezgin, and P. Sundell, Nucl. Phys. B791 , 231 (2008), arXiv:0706.2983 [hep-th].
- [29] R. Aros, C. Iazeolla, J. Norea, E. Sezgin, P. Sundell, and Y. Yin, JHEP 03 , 153 (2018), arXiv:1712.02401 [hep-th].
- [30] M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Prog. Theor. Phys. 125 , 795 (2011), arXiv:1012.1079 [astro-ph.CO].
- [31] M. Shiraishi, D. Nitta, and S. Yokoyama, Prog. Theor. Phys. 126 , 937 (2011), arXiv:1108.0175 [astro-ph.CO].
- [32] M. Shiraishi, JCAP 1206 , 015 (2012), arXiv:1202.2847 [astro-ph.CO].
- [33] M. Shiraishi, JCAP 1311 , 006 (2013), arXiv:1308.2531 [astro-ph.CO].
- [34] M. Shiraishi, A. Ricciardone, and S. Saga, JCAP 1311 , 051 (2013), arXiv:1308.6769 [astro-ph.CO].
- [35] M. Shiraishi, C. Hikage, R. Namba, T. Namikawa, and M. Hazumi, Phys. Rev. D94 , 043506 (2016), arXiv:1606.06082 [astro-ph.CO].
- [36] P. D. Meerburg, J. Meyers, A. van Engelen, and Y. Ali-Hamoud, Phys. Rev. D93 , 123511 (2016), arXiv:1603.02243 [astro-ph.CO].
- [37] G. Domnech, T. Hiramatsu, C. Lin, M. Sasaki, M. Shiraishi, and Y. Wang, JCAP 1705 , 034 (2017), arXiv:1701.05554 [astro-ph.CO].
- [38] H. W. H. Tahara and J. Yokoyama, PTEP 2018 , 013E03 (2018), arXiv:1704.08904 [astro-ph.CO].
- [39] N. Bartolo, G. Orlando, and M. Shiraishi, JCAP 1901 , 050 (2019), arXiv:1809.11170 [astro-ph.CO].
- [40] R. Kothari and D. Nandi, (2019), arXiv:1901.06538 [astro-ph.CO].
- [41] M. Shiraishi and T. Sekiguchi, Phys. Rev. D90 , 103002 (2014), arXiv:1304.7277 [astro-ph.CO].
- [42] M. Shiraishi, M. Liguori, and J. R. Fergusson, JCAP 1501 , 007 (2015), arXiv:1409.0265 [astro-ph.CO].
- [43] P. A. R. Ade et al. (Planck), Astron. Astrophys. 594 , A17 (2016), arXiv:1502.01592 [astro-ph.CO].
- [44] M. Shiraishi, M. Liguori, and J. R. Fergusson, JCAP 1801 , 016 (2018), arXiv:1710.06778 [astro-ph.CO].
- [45] Y. Akrami et al. (Planck), (2019), arXiv:1905.05697 [astro-ph.CO].
- [46] M. Shiraishi, Front. Astron. Space Sci. 6 , 49 (2019), arXiv:1905.12485 [astro-ph.CO].
- [47] M. Hazumi et al. (LiteBIRD), Proc. SPIE Int. Soc. Opt. Eng. 8442 , 844219 (2012).
- [48] K. N. Abazajian et al. (CMB-S4), (2016), arXiv:1610.02743 [astro-ph.CO].
- [49] J. Delabrouille et al. (CORE), JCAP 1804 , 014 (2018), arXiv:1706.04516 [astro-ph.IM].
- [50] J. M. Maldacena and G. L. Pimentel, JHEP 09 , 045 (2011), arXiv:1104.2846 [hep-th].
- [51] J. M. Maldacena, JHEP 05 , 013 (2003), arXiv:astroph/0210603 [astro-ph].

- [52] D. Anninos, V. De Luca, G. Franciolini, A. Kehagias, and A. Riotto, JCAP 1904 , 045 (2019), arXiv:1902.01251 [hep-th].
- [53] H. Osborn and A. C. Petkou, Annals Phys. 231 , 311 (1994), arXiv:hep-th/9307010 [hep-th].
- [54] S. Giombi, in Proceedings, Theoretical Advanced Study Institute in Elementary Particle Physics: New Frontiers in Fields and Strings (TASI 2015): Boulder, CO, USA, June 1-26, 2015 (2017) pp. 137-214, arXiv:1607.02967 [hep-th].
- [55] D. Anninos, T. Hartman, and A. Strominger, Class. Quant. Grav. 34 , 015009 (2017), arXiv:1108.5735 [hepth].
- [56] D. Anninos, F. Denef, and D. Harlow, Phys. Rev. D88 , 084049 (2013), arXiv:1207.5517 [hep-th].
- [57] D. Anninos, F. Denef, R. Monten, and Z. Sun, (2017), arXiv:1711.10037 [hep-th].
- [58] N. Bobev, T. Hertog, and Y. Vreys, JHEP 11 , 140 (2016), arXiv:1610.01497 [hep-th].
- [59] M. Shiraishi, S. Yokoyama, D. Nitta, K. Ichiki, and K. Takahashi, Phys. Rev. D82 , 103505 (2010), arXiv:1003.2096 [astro-ph.CO].
- [60] J. R. Pritchard and M. Kamionkowski, Annals Phys. 318 , 2 (2005), arXiv:astro-ph/0412581 [astro-ph].
- [61] J. R. Fergusson, M. Liguori, and E. P. S. Shellard, Phys. Rev. D82 , 023502 (2010), arXiv:0912.5516 [astroph.CO].
- [62] N. Aghanim et al. (Planck), (2018), arXiv:1807.06209 [astro-ph.CO].
- [63] Y. Akrami et al. (Planck), (2018), arXiv:1807.06211 [astro-ph.CO].
- [64] T. Matsumura et al. , Journal of Low Temperature Physics (2016), 10.1007/s10909-016-1542-8.
- [65] T. Matsumura et al. , (2013), 10.1007/s10909-013-09961, [J. Low. Temp. Phys.176,733(2014)], arXiv:1311.2847
- [astro-ph.IM].
- [66] C. L. Bennett et al. (WMAP), Astrophys. J. Suppl. 208 , 20 (2013), arXiv:1212.5225 [astro-ph.CO].
- [67] G. Hinshaw et al. (WMAP), Astrophys. J. Suppl. 208 , 19 (2013), arXiv:1212.5226 [astro-ph.CO].
- [68] E. Komatsu, D. N. Spergel, and B. D. Wandelt, Astrophys.J. 634 , 14 (2005), arXiv:astro-ph/0305189 [astroph].
- [69] H. F. Gruetjen, J. R. Fergusson, M. Liguori, and E. P. S. Shellard, Phys. Rev. D95 , 043532 (2017), arXiv:1510.03103 [astro-ph.CO].
- [70] M. Bucher, B. Racine, and B. van Tent, JCAP 1605 , 055 (2016), arXiv:1509.08107 [astro-ph.CO].
- [71] E. Komatsu et al. (WMAP), Astrophys. J. Suppl. 180 , 330 (2009), arXiv:0803.0547 [astro-ph].
- [72] http://lambda.gsfc.nasa.gov .
- [73] J. R. Fergusson, M. Liguori, and E. P. S. Shellard, JCAP 1212 , 032 (2012), arXiv:1006.1642 [astro-ph.CO].
- [74] M. Shiraishi, M. Liguori, and J. R. Fergusson, JCAP 1405 , 008 (2014), arXiv:1403.4222 [astro-ph.CO].
- [75] J. R. Fergusson, Phys. Rev. D90 , 043533 (2014), arXiv:1403.7949 [astro-ph.CO].
- [76] M. Shiraishi, M. Liguori, J. R. Fergusson, and E. P. S. Shellard, JCAP 1906 , 046 (2019), arXiv:1904.02599 [astro-ph.CO].
- [77] W. Hu, Phys. Rev. D62 , 043007 (2000), arXiv:astroph/0001303 [astro-ph].
- [78] A. Lewis, A. Challinor, and D. Hanson, JCAP 1103 , 018 (2011), arXiv:1101.2234 [astro-ph.CO].
- [79] Y. Akrami et al. (Planck), (2018), arXiv:1807.06211 [astro-ph.CO].
Valerio De Luca Gabriele Franciolini Alex Kehagias Antonio Riotto Maresuke Shiraishi
Constraining graviton non-Gaussianity through the CMB bispectra
Tensor non-Gaussianities are a key ingredient to test the symmetries and the
presence of higher spin fields during the inflationary epoch. Indeed, the shape
of the three point correlator of the graviton is totally fixed by the
symmetries of the de Sitter stage and, in the case of parity conservation, gets
contributions only from the ordinary gravity action plus a higher derivative
term called the (Weyl)$^3$ action. We discuss current and future bounds on the
three point tensor contribution from the (Weyl)$^3$ term using cosmic microwave
background (CMB) bispectra. Our results indicate that forthcoming experiments,
such as LiteBIRD, CMB-S4 and CORE, will detect the presence of the (Weyl)$^3$
term if $M_p^4 L^4 \sim 10^{17} r^{-4}$, where $L$ parametrizes the strength of
the (Weyl)$^3$ term and $r$ is the tensor-to-scalar ratio, which corresponds to
$L\gtrsim 3.2 \times 10^5 M_p^{-1}$, while the current upper limit is $M_p^4
L^4 = (1.1 \pm 4.0) \times 10^{19} r^{-4}$ (68%CL).