Generic 3-point statistics with tensor modes in light of Swampland and Trans-Planckian Censorship Conjecture

Recently proposed Swampland Criteria (SC) and Trans-Planckian Censorship Conjecture (TCC) question the theoretical consistency of inflationary cosmology. While SC can be in tension with the recent observational bound on tensor-to-scalar ratio which is a consequence of single field consistency relation, the requirement of TCC translates to the fact that inflation can not produce a detectable signal for stochastic gravitational wave, raising a question on future prospects of detection of Primordial Gravitational Waves (PGW). But the constraint can be relaxed significantly by considering Non Bunch Davies (NBD) initial states, that in turn brings back the observational relevance of PGW via its 2-point function. Also with NBD states the tension between SC and observational bound on tensor-to-scalar ratio vanishes. Focusing on these theoretical importance of NBD states, in this article we develop consistent 3-point statistics with tensor modes for all possible correlators (auto and mixed) for NBD initial states in a generic, model independent framework of Effective Field Theory of inflation and analyze the properties of the correlators in the light of TCC. We also construct the templates of the corresponding nonlinearity parameters fNL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{NL}$$\end{document} for different shapes of relevance and investigate if any of the 3-point correlators could be of interest for future CMB missions. Our analysis reveals that the prospects of detecting the tensor auto correlator are almost nil whereas the mixed correlators might be relevant for future CMB missions.


Introduction
Nowadays theoretical consistency of any viable model of inflation can be tested against recently proposed Swampland Criteria (SC) [1], and Trans-Planckian Censorship Conjecture (TCC) [2]. Even though the success of Inflationary Cosmology lies in its indigenous ability to solve the puzzles with Standard Big Bang Cosmology and to provide the seeds for Large Scale Structures (LSS) we see today, as well as in its profound consistency with the highly precise Cosmic Microwave Background (CMB) data, the latest being the Planck 2018 data [3], recently it was revealed that the single field, slow roll inflationary paradigm itself is plagued with serious theoretical inconsistency [1]. In brief, in order to embed any model of inflation consistently in a Quantum Gravity theory, it has to satisfy a set of criteria, collectively called the Swampland Criteria, that states that (i) the field excursion by scalar fields in field space is bounded by φ < O(1)M pl ; and (ii) the potential of a scalar field which rolls and dominates the energy density of the universe has to satisfy the condition (1) and being the first slow roll parameter and M pl is the Planck mass. Now if we consider the single field consistency relation r = 16 = 8c 2 with r being the tensor-to-scalar ratio, we can see that the requirement of SC of c ∼ O (1) can ruin the recent observational bound r < 0.064 1 [3]. The second constraint on the theory of inflation comes from recently conjectured Trans-Planckian Censorship Conjecture [2]. In order to produce the seeds of perturbations in LSS, the quantum vacuum fluctuations exit the horizon and get frozen during inflation only to re-enter the horizon at a later stage of cosmological evolution. However, one can see that if inflation lasts more than the minimum period required, all the observable modes has to have a length scale smaller than Planck scale, which once again leads to inconsistent theory of Quantum Gravity. This problem is known as Trans-Planckian Problem [6]. The only wayout seems to be coming through a conjecture, called the Trans-Planckian Censorship Conjecture [2], that states that the length scales smaller than Planck scale can never exit the Hubble Horizon. As a consequence this conjecture places a tight theoretical constraint on the tensor-to-scalar ratio r < 10 −30 and so a constraint on the first slow roll parameter can be written as, < 10 −31 from the single field consistency relation r = 16 [7].
Looking at the constraint on put by the SC and TCC, it may appear that they are in tension with each other, and together they may serve as a no go theorem for a consistent single field slow roll model of inflation. However, it was revealed of late that TCC can be derived from the first SC [8]. Also the second SC can be refined using the first SC and Bousso's entropy bound [9] and can be restated as, a field theory will not be in Swampland if (1) [10,11]. Consequently, one can have a valid inflationary scenario consistent with the refined SC if one satisfies the following condition: the second slow roll parameter η ∼ O(1) with the first slow roll parameter still satisfying the observational bound 1. The constraint on η from Planck observation [3] is given by η ∼ −0.02 which can be considered as close to O(1) as argued in [12] . As we have already stated, the only requirement from TCC is 1. Hence, to summarize the theoretical constraints, any single field slow roll inflationary theory would be consistent with both the refined SC and TCC if one satisfies 1 and η ∼ O(1) simultaneously.
Even though one takes care of the theoretical perspective and the possible wayout as discussed above while inflationary model building, on the backdrop this leads to an altogether novel problem from observational point of view, namely, the future prospects of detection of primordial gravitational waves (PGW). As stated earlier TCC requires a very tiny 1 In the case of SC the first slow roll parameter has to satisfy the condition √ 2 > c ∼ O(1). If we consider ∼ 0.006 respecting the bounds of Planck 2018 then c ∼ 0.11 which can also be considered to be O (1). But in [4,5] it was shown that the lowest value of c ∼ 0.8 and this value contradicts the observational bound on tensor-to-scalar ratio as r = 8 c 2 /2. value for r < 10 −30 . This tiny value is beyond the scope of any future observation, let alone the next generation CMB missions like CMB-S4 [13], LiteBIRD [14,15], COrE [16] etc.
This bound is further tightened to r < 10 −47 if one assumes that the pre-inflationary era was radiation dominated [17]. Although the bound can slightly be relaxed to r < 10 −10 if one assumes that the equation of state of inflation changes from −1 to −(1/3) after a couple of e-folding [18], this is still much below the scope of future CMB missions. Thus, those theoretical constraints put a serious question on the prospects of PGW to act as the smoking gun of inflationary paradigm [19], namely, whether or not it can act as the 'Holy Grail' of inflationary cosmology in the light of these theoretical constraints.
However, with further progress of theoretical analysis, it was recently found that one can bypass the SC within the standard single field inflationary framework if one considers the initial states as follows: the tensor modes in Non Bunch Davies state (NBD) and scalar modes in Bunch Davies (BD) state [4,20]. Subsequently, it was also revealed that if the initial states of both the scalar and tensor modes are NBD, one can relax the TCC bound on r as well, thereby revitalising the possibility of producing a detectable PGW by vacuum fluctuation of inflaton [21] orders close to the range of future CMB missions. Also pre-inflation scenarios suggest that one should set NBD vacuum as the initial states [22]. Thus, the above prescription of considering NBD states for both scalar and tensor modes to obtain a detectable PGW signal from a consistent theoretical framework brings back the relevance of PGW as a probe of inflationary cosmology. The effect of NBD on the scalar perturbations has been studied to a considerable extent [23][24][25][26][27][28][29][30], mostly with theoretical motivation. However, in the light of PGW, the motivation to explore NBD states gets an observational relevance as well. It is important to note that other than initial states modification the upper bound on r can be relaxed by other mechanisms as well e.g in [31] it is shown that by considering a oscillatory stage after inflation the upper bound on r can be significantly relaxed as r ≤ 10 −8 .
In this article our primary intention is to develop consistent 3-point statistics with tensor modes with NBD initial states keeping in mind the above discussion on its relevance in the context of SC and TCC. We would primarily focus on TCC to study the properties of 3-point functions and for any possible enhancement in signal strength we would focus on the parameter space where TCC bound is relaxed due to NBD states. To this end, we would bring forth different correlators relevant for the non-Gaussian behavior of PGW followed by the templates for the corresponding nonlinearity parameters f N L for different relevant shapes considering NBD initial states for both scalar and tensor modes; and study their prospects of detection in future CMB missions. We believe the present analysis is important for couple of reasons. First, the higher order correlation functions can be very sensitive to the choice of initial state (BD/NBD) and hence can be an important probe for the initial vacuum as well. Secondly, the shape of the auto correlation of tensor modes can shed some light on the source of PGW. Thirdly, the mixed correlator of scalar and tensor modes can probe additional properties of PGW which will be discussed in the next sections. To name one, correlators of tensor-scalar-scalar type can serve as a probe of spatial diffeomorphism breaking during inflation [32]. From observational point of view, this mixed correlator also reflects on the quadrupole moment and might be of interest for future CMB missions. Last but not the least, these non-Gaussian properties serve as additional probe for PGW along with the 2-point function r . Thus, detection of any of the correlators would give us a hint on PGW and initial state. Some preliminary studies on the prospects of tensor non-Gaussianities have been reported by the present authors in [33,34] and also by some other authors [19,29,[35][36][37]. However, almost all of them consider BD vacuum. Here we would like to extend the analysis for NBD states and explore the prospects of both auto and mixed correlators.
As of now, the observational bounds on tensor non-Gaussianities are not too tight. The current constraints on the amplitude of tensor bispectrum in equilateral limit is given by Planck 2018 as: f T N L = 800 ± 1100 [38]. However, for the reasons mentioned above, the hunt for any possible non-Gaussian behavior of PGW in upcoming CMB missions will be more relevant than ever. These constraints will be further improved in next generation CMB missions. For example, LiteBIRD [14,15] targets to improve the constraint of the amplitude of tensor bispectra by three orders of magnitude [14,15]. The constraint on tensor-tensor-scalar mixed correlator will be improved by CMB-S4 [13]. Others have their own specific target. Keeping this in mind, in this article, we would like to explore all possible correlators, namely, the tensor-tensor-tensor (auto) correlator and tensortensor-scalar and tensor-scalar-scalar (both mixed) correlators and study the possible enhancement in signal strength in the light of TCC. The theoretical analysis is done using the model independent framework of Effective Field Theory (EFT) of inflation [39] to calculate the tensor-tensor-tensor and scalar-tensor-tensor correlators. To calculate the tensorscalar-scalar correlator we have used the EFT of inflation with broken space-time diffeomorphism [32] as this type of correlator can reflect a unique feature of broken spatial diffeomorphism. Once we have the correlators in our hand, we would move on to proposing templates for the corresponding nonlinearity parameters f N L for different shapes of interest and would also investigate if any of them could be the point of interest for future CMB missions by finding out the possible upper limits for each one allowed by the parameters of the theory under consideration. We would like to reiterate that because of the background EFT of inflation, our calculations and results are more or less generic and model-independent.

Swampland, TCC and Non Bunch Davies states
As already mentioned, it is the non Bunch Davies states that can give rise to a viable inflationary scenario respecting the SC and TCC. Let us begin our discussion with a brief review on the role of NBD states and how that can help bypass the SC and relax the bound on r coming from TCC in the same vein of [4,20,21]. This will also help us develop the rest of the article from a consistent theoretical setup. However it should be noted that in our analysis we would be primarily focusing on TCC to study the primordial correlators. To this end, we will first discuss the scenario with the minimal model, i.e., single field inflation with canonical kinetic term, and will subsequently move on to discussing the non-canonical inflationary models.

Canonical models of inflation
For a canonical inflation scenario the power spectra of scalar and tensor modes for NBD states are given by and, where ζ and γ represent scalar and tensor modes respectively; α and β's are Bogolyubov coefficients and the indices s and t signify scalar and tensor respectively as we have considered different vacua for scalar and tensor modes. Setting α s/t k = 1 and β s/t k = 0 one readily gets back the BD states. We also, introduce two quantities, The tensor-to-scalar ratio with NBD states can be written as, where is the first slow roll parameter and The first theoretical constraint on the Bogolyubov coefficients come via the Wronskian condition using which the Bogolyubov coefficients can further be parameterized as represents the number of NBD particles in BD state and θ (s/t) α (k) and θ (s/t) β (k) are the phase factors. Consequently, one can define the relative phase as This parametrisation would help us to reduce the number of parameters since the individual phase factors are degenerate. Also in our analysis we consider θ (s/t) (k) to be constant. Further, there is a theoretical constraint on the parameters β (s/t) k coming from backreaction condition [23,40]. For this one can model them as following [23] β s/t where η 0 describes the time when modes are below cut-off scale M (s/t) . One can readily check that for k > M (s/t) , β s/t → 0. Now, in order a theory of inflation to be valid within the regime of EFT, it has to satisfy M (s/t) > H . With the requirement that excited initial state does not spoil the single field slow roll inflationary conditions one can write the backreaction condition as [23]: From now we will represent the second slow roll parameter with η to avoid confusion with conformal time.
As mentioned earlier, the SC forces the first slow roll and as a result the single field consistency relation now violates the observational bound on r . From (2.5) and (2.6) one can see that for < 1 the observational constraint on r can still be satisfied even with a large coming from SC. Now there are two ways to make < 1. The first one is to have NBD states for scalar fluctuations and BD state for tensor fluctuation. One can choose a model where s = |α s k + β s k | 2 1 and suppress the value of r , without breaking the backreaction condition (2.12) as shown in [23]. But as explained in [4] this choice of s can ruin the constraint on scalar three point function. On the other hand, if one starts with tensor modes in NBD state and scalar modes in BD state one can find a region in the parameter space that satisfy the PGW bispectrum constraints and allow c ∼ 0.8; larger value of c ∼ 0.9 can also be achievable in this region of parameter space [20] thus satisfying the SC.
Secondly, TCC puts a very tight constraint on r < 10 −30 and hence from single field consistency relation for BD vac-uum, the first slow roll parameter < 10 −31 . This bound on the tensor to scalar ratio basically states that in future if any PGW gets detected that can not be due to inflation but can be produced from other sources during inflation [34][35][36][37]. So, apparently, PGW loses its 'Holy Grail' status in the light of TCC [41]. However, this theoretical bound on r placed by TCC can be relaxed significantly if one uses NBD states for both tensor and scalar modes [21]. Consequently, the consistency relation now modifies to (2.5), where the in (2.5) factor is re-written in terms of s/t defined earlier in this section. In this scenario TCC places the bound not on r alone but on the combined parameter / s < 10 −31 . So if one can enhance the t then the bound on r can be relaxed significantly. To do so we take a look at the backreaction condition (2.12). If one uses (2.1) to properly replace the M pl factor of (2.12) one arrives at 23 with M (t) = 10H we can have r ≤ 0.001 which is still within the detectable range of next generation CMB missions. Here one important thing to note is that β (s) 0 1 does not spoil the scalar non-Gaussianity bound as TCC requires a very small value of that guarantees that the bispectrum amplitude remains within the observational bound as it is evident from [21,42,43].
where f loc N L is the amplitude of local configuration of scalar bispectrum and k s and k l are respectively the shortest and longest wavelength probed by any CMB mission under consideration (e.g., Planck mission as of today).

Non-canonical models of inflation
Until now we have only discussed the scenario with canonical kinetic term. However, EFT of inflation can take into account inflationary models with non-canonical terms as well. This will modify the consistency relation [44] within the single field framework as well as the backreaction condition [27]. Effectively, the modifications come via the scalar and tensor sound speeds c s/t , including which (2.1) and (2.2) now look Consequently" the consistency relation now takes the form r = 16 t s c s c t (2.18) and the backreaction condition modifies to It is straightforward to check that for c s/t = 1 one gets back the canonical scenario. From (2.18) we can see in the context of SC, apart from the t / s term, there is another factor c s /c t which can also act as a suppression factor for a small value for c s . However, c s can not be arbitrarily small as there is another constraint on it coming from the bounds on scalar non-Gaussianities. In the equilateral limit, the amplitude of scalar bispectrum reads [45,46], Hence, together with the constraint from scalar non-Gaussianity f eq N L = −4 ± 40 [50], the value of factor c coming from SC becomes c ∼ 0.37 with c s > 0.067.
Clearly, the above value of c is not compatible with SC, as discussed earlier. So, in order to bypass the SC, one need take shelter of the NBD factor t s [20]. However, for TCC this non-trivial sound speed can help in relaxing the bound on r . From (2.19) one can write Together with (2.18) this equation tells us that in order to relax the bound on r , one needs to have a small value for the sound speed of tensor fluctuation such that the factor c s c t acts as an enhancing factor.
Combining all the above constraints on the parameters, we have c s / (s) < 10 −31 and r < 0.064. Figure 1 shows the region between c t and N (t) 0 that is allowed considering both the theoretical (TCC bound) and latest observational constraints (Planck 2018). The figure show that for particular combinations of c t and N (t) 0 one can have a bound on r which can match the upper bound set by Planck 2018, thereby bringing back its relevance as the smoking gun for inflationary cosmology.

The second order action from EFT
In this section we briefly summaries the major equations staring from the EFT of inflation [39], a model independent framework developed to analyze primordial fluctuations, that 0 and c t using Planck constraint on the amplitude of two point function will help us calculate the 3-point function for PGW with NBD states. As pointed out earlier, the prospects of EFT of inflation in exploring tensor non-Gaussianities have been investigated to some extent by the present authors in couple of articles [33,34]. In this article, our analysis takes under consideration NBD states as well as is consistent with theoretical constraint of TCC.
A consistent EFT of inflation can be given starting from the fact that time diffeomorphism is broken spontaneously after inflation and as a result Goldstone Boson is produced. In unitary gauge the Lagrangian can be written as, Here δ K is the fluctuation in extrinsic curvature. M i and M i denotes parameters that measures the strength of higher order fluctuations. Different combination of M i andM i pro-duce results for different kind of theory. Note that the scalar fluctuation is now represented by the Goldstone Boson π produced by broken time diffeomorphism and it is related to curvature perturbation ζ as ζ = −H π .
Another aspect to note here is that the sound speed of scalar and tensor fluctuations get modified due to the presence of higher order fluctuation term. From (3.1) the second order Lagrangian for scalar fluctuation can be written as, where scalar sound speed is given by Further, from (3.1) the second order tensor fluctuation can be written as, and the corresponding sound speed of tensor perturbation is given by As we have seen in the last section, these modified sound speeds of perturbations can significantly influence the backreaction condition and hence the parameter space to relax TCC bound. In the rest of the article, we will make use of the above actions as well as the non-trivial sound speeds for scalar and tensor perturbations to analyse the 3-point statistics for tensor modes using NBD initial state.

The tensor-tensor-tensor (γ γ γ ) correlator
Let us begin with calculating the auto-correlation of tensor modes that gives rise to the auto-bispectrum, i, e the tensortensor-tensor correlator. To calculate the bispectrum we have to look at the third order Lagrangian of tensor fluctuation and the most general Lagrangian can be written as [33,51], In the above action, the term proportional to M pl is the sole contribution of the so-called Einstein or R part and the term proportional toM 9 is the contribution of higher order gravitational fluctuations. As it will be revealed below, together they play a crucial role in determining the strength of tensor bispectra.
In order to calculate the 3-point correlation function of tensor fluctuations, we make use of the IN-IN formalism. With the NBD states the γ γ γ 3-point function for any general interaction Hamiltonian H I (t) can be written as, where s i = {+, ×} are the polarization indices. An essential point to note here is that the lower limit of the integration is different from the BD case which can give rise to a non-trivial oscillatory term proportional to e i K η 0 , where K can be any possible combination of (k 1 ± k 2 ± k 3 ). We will work with sub-horizon modes where |k i η 0 | 1. So, with NBD states the contribution of the R part to the γ γ γ correlator is given by with the following notations: α i = α k i , β i = β k i and so on, and the momentum conserving delta function is dropped here. Here on the right hand side the s i indices will determine the sign before the momentum hands k i (whenever the s i indices appear before k i ) which follows from the fact that the polarization tensors satisfy s i j (k) = s i j (−k) where if s is '+' thens is '×' and vice versa [47]. The g osc,R is the oscillatory part of the bispectrum which is explicitly written in 1. Also, a ' * ' stands for the complex conjugate of α i and β i , as the case may be.
The term F(x, y, z) appears due to the contraction of polarization tensors. Written explicitly, it reads, with x, y, z are, in this particular scenario, given by the combination of s i and k i as given in (4.3). For rapid oscillations i,e if K η 0 1, the contribution of g osc,R can be neglected as in this limit the average of the exponential will be zero. But in some configurations where the oscillation is not rapid for relevant modes, g osc,R can induce oscillatory features on the bispectrum [48,49]. From (A.1) one can check that for flattened limit configuration where, k 2 + k 3 = k 1 or similar configurations, K = k 2 + k 3 − k 1 vanishes and the threepoint function can get an enhancement by appropriate powers of η 0 as seen for scalar fluctuations [23], but the presence of the function F(x, y, z) in (4.3) and (A.1) ensures that the tensor bispectrum due to Einstein term always vanishes in flattened limit as one can see from (4.4) that F(x, y, z) contains the term (k i + k j − k l ) for which F(x, y, z) vanishes in flattened limit. In other momentum configurations the g osc,R will induce an oscillatory feature and can have an enhancement proportional to |kη 0 |, but projecting this oscillatory part onto different observational templates of bispectrum one power of |kη 0 | gets reduced [48] and observationally there is no enhancement in the signal due to this factor. Also these oscillatory features induced by g osc,R are hard to detect observationally, so for any kind of enhancement in the signal we focus solely on the non-oscillatory part.
To calculate the amplitude of bispectrum we will use the parameterisation (2.8) and (2.9). N (s/t) k can be modeled as Further, if M (s) lies within the observable mode, i.e if the exponential factor is nonnegligible, then the amplitude of bispectrum also gets suppressed. So, we consider N Following [19,38]. let us define the amplitude of tensor non-Gaussianity, i.e., the tensor nonlinearity parameter f N L as, In the limit N (s/t) 0 1, one finds that (s/t) ∼ N (s/t) 0 . It can also be shown that in the limit N (s/t) 0 1, Eq. (4.3) can be estimated as, where f (k i ) encodes all the momentum dependence. The (N (t) 0 ) 2 comes from different combinations of α t i and β t i in Eq (4.3).
Consequently, from (4.5) and (2.16) the nonlinearity parameter for γ γ γ correlation can be estimated as, Because of the above definition of f N L , it will always be proportional to r 2 , hence, even if we manage to produce a large value for N (t) 0 we can not generate a large bispectrum from the R part of the Einstein term alone.
However, we have in our hand contribution from another operator, namely, theM 9 operator, that arises from the higher order gravitational fluctuations in EFT. Let us now discuss the contribution ofM 9 operator to the 3-point tensor correlator. With this term playing the central role, the 3-point function for the EFT term can be written as, where, g osc,E FT is the oscillatory part. The TCC constraint < 10 −31 gives the maximum value for H as H ∼ 3 3/2 × 10 −20 M pl ; Consequently, with the current observational constraints mentioned earlier, the estimates for different limits of f N L for the operatorM 9 from the non-oscillatory part of the bispectrum are as: for equilateral limit f eq N L ∼ 10 −23 whereas for squeezed limit: f sq N L ∼ 10 −21 . These estimations are done under the assumption thatM 9 ∼ M pl and c t = 1 to estimate the maximum contribution of this operator. The details of f N L parameter is given in 1. It is also important to note that after projecting onto different templates g osc,E FT can have an effective enhancement by a factor of |kη 0 | (as one can see from (A.4)) but as the signal strength is extremely small this enhancement will be rather irrelevant from the observational perspective. This indicates that the contribution of this operator is highly suppressed for TCC as compared to the result reported in [33]. From this above discussion we can see that even with NBD states the signal strength for γ γ γ correlator arising from the combination of Einstein R and the EFT term is too feeble to be detected, at least in next generation CMB missions.

The tensor-tensor-scalar (γ γ ζ ) correlator
As the auto bispectrum, the mixed correlation between scalar and tensor modes can also help us to estimate the parameters of PGW. In [52] it is shown thatM 3 δ K μ ν δ K ν μ operator can contribute to γ γ ζ correlation in decoupling limit and it is very sensitive to the sound speed of tensor fluctuation. This is because (3.5) reveals that the tensor sound speed is modified due to the presence ofM 3 . So it is evident that a non zerō M 3 i, e c t = 1 will lead to non zero amplitude of γ γ ζ correlator.
In decoupling limit, the part of the action that is contributing to the above mixed correlator arises fromM 3 δ K μ ν δ K ν μ and can be written as, From the above action (5.1), we can calculate the γ γ ζ correlator that reads Here, as before, the s i are the polarization indices and α i and β i follow the definition as given earlier. The s i j (k)'s are the polarization tensors (not to be confused with the first slow roll parameter used in the analysis). Note that these tensors did not appear in the earlier correlators as they were completely contracted. Eq. (5.2) is calculated considering nontrivial sound speeds of scalar and tensor fluctuations.
However, current observations do not put any constraint on this correlator. So, in what follows we will mostly be interested to find out how much this signal can be enhanced with the theoretical constraints from TCC, and whether or not that can be of interest for upcoming CMB missions.
First let us analyze the contribution of the non-oscillatory part. From Eq. (5.2) it transpires that both β (s) k and β (t) k contribute to the mixed correlator. As a consequence, both N (s/t) 0 and θ (s/t) would eventually contribute to the signal strength of the corresponding mixed bispectrum. As we have argued earlier that to avoid TCC we need to have N Also (3.5) givesM 3 = (1 − c −2 t )M pl and using (2.16) one can estimate the correlator as, Therefore, using the same definition of f N L from (4.5) one can write, Here we make use of the fact that a large N 0 and as a consequence t > s , (5.6) suggests that for NBD state it is possible to have (r t / s > 1) and the mixed correlator gets an enhancement in the signal.
The signal strength of the above nonlinearity parameter for the γ γ ζ mixed correlator is enhanced due to the presence of the parameters of EFT and NBD states. However, before that, let us explore another nontrivial aspect of the correlator under consideration. We can see from (5.2) that the shape of the signal is changed with respect to the BD case and now it gets peaked at either of the following configurations: c t (k 2 + k 3 ) = c s k 1 , c s k 1 + c t k 3 = c t k 2 , c s k 1 + c t k 2 = c t k 3 . These configurations are not strictly the flattened limit where, k i + k j = k l but is nontrivially modified by the sound speed of perturbations. However, with c s/t = 1, they get back to the flattened limit. It is important to note here is that in these configurations the g osc,γ γ ζ affects the signal non-trivially. As one can see from (A.5) the dominant contribution of g osc,γ γ ζ in c t (k 2 + k 3 ) = c s k 1 configuration can be written as, Replacing c s k 1 with c t (k 2 +k 3 ) one can see that there is an enhancement of |c t kη 0 | 2 which after angular averaging one power of η 0 will drop [23] and it will enhance the signal by |c t kη 0 |. So in these configurations the amplitude of the correlator will get an additional enhancement factor of |c t kη 0 | on top of t / s . It is however important to note that in [52] the authors used a different route to show that the mixed correlator for c t = 1 can be larger than mixed correlator generated due to shift and lapse function for c t = 1. For BD case one can show that In the light of TCC we can see that, the mixed correlator for c t = 1 is now significantly larger than the c t = 1 case even with BD state. So, in principle, the chance of detection of this correlator should be much higher if c t = 1. As a result, the signal strength is significantly enhanced for a nontrivial tensor sound speed that arises from the higher curvature terms in EFT of inflation. However, as mentioned earlier, absence of any observational constraint on this correlator from present observation forbids us to do a numerical estimate of the corresponding signal strength. Even in absence of that, one can claim from the above analysis that these enhanced signals might of interest for future CMB missions and the estimates might be more relevant once the constraints come in.

The tensor-scalar-scalar (γ ζ ζ ) correlator
Till now we have analyzed the correlators with NBD states using EFT of inflation where only consideration was that the time diffeomorphism gets broken after the end of inflation. In order to calculate the tensor-scalar-scalar mixed correlator, we consider a scenario where spatial diffeomorphism is also broken along with the usual time diffeomorphism breaking. This formalism is more rich in structure and has enormous potential to give rise to novel features, e.g., the blue tilt of tensor power spectrum. Corresponding EFT for the scenario considering the breaking of space-time diffeomorphism can be written as follows [32]. The expectation value of the symmetry breaking fields look The fieldsφ 0 (t) andφ x i are identified as clock and ruler respectively, during inflation [32] and the parameter α measures the breaking of spatial diffeomorphism. In this framework, the general Lagrangian can be written as, where the parameters X, Y, Z defining F(X, Y, Z ), which can be any function of the above quantities respecting the symmetry, are given by First notice that due to the presence of α in the action (6.2), the expression for power spectrum of scalar fluctuations is non trivially modified in this scenario. Written explicitly, Here, the bar on the quantities represents their values as evaluated in the background and the subscripts represents derivative with respect to X , Y or Z . Note that the quan-tityF 2 y /2a 2 is slowly varying and can be treated as more or less a constant for our evaluation of observable parameters. The parameter α is also responsible for a tensor mass m 2 γ = α 2 (F z + α 2F zz )/a 2 ; but in the limit α 1 this can be ignored and the power spectrum for tensor modes closely resembles (2.17) so far as its numerical values are concerned.
The second nontrivial effect of the action (6.2) appears in the γ ζ ζ correlator via the interaction term F Y /(2a 2 ) (γ i j ∂ i π∂ j π/(2a 2 )). The γ ζ ζ correlator describes how much local quadruple will affect the power spectrum of scalar perturbations when a long wavelength tensor mode correlates with small wavelength scalar mode. The importance of this term in the action is that the behaviour of γ ζ ζ generated from this term does not satisfy the consistency relation of γ ζ ζ correlator and can be a distinct signature of spatial diffeomorphism breaking during inflation. With NBD state the correlator can be written as, The momentum dependence of the oscillatory term g osc,γ ζ ζ is exactly same as g osc,R and as we are only interested about the squeezed limit configuration of this correlator we will only focus on the non-oscillatory term to study the enhancement in the signal. In squeezed limit k 1 k 2 , k 3 if one chooses θ (s) = π and θ (t) = 0 one finds that, The interesting point to note here is that the second term in (6.8) can be produced with a BD state, but the first term is a unique contribution of NBD state which gets an enhancement in squeezed limit of O( k 2 k 1 ). Also it is important to note that for θ (s) = π choice the (N (s) 0 terms in the correlator vanish. For the first term of (6.8) the expression for the mixed correlator in squeezed limit can be rewritten using the expressions for scalar and tensor power spectra in NBD state, as follows It appears from the above expression that in order to obtain an enhancement of signal for f N L for this correlator, the following condition has to be satisfied: 1/ t (k 2 /k 1 ) > 1.
However, as we have mentioned before, TCC requires a large t ∼ 10 29 in order to have a detectable signal of PGW and thus suppressing this effect. On the other hand if we are not strict about the theoretical effect on θ (s) and set it any value other than π we can have, From (6.10) it is evident that in order to have a large signal for this mixed bispectrum, the parameters need to satisfy the following condition: s / t (k 2 /k 1 ) > 1 which is possible in some region of parameter space even for large s and t . Thus, a considerably strong signal for this correlator is achievable with NBD initial state for both scalar and tensor modes along with non-trivial sound speeds, satisfying TCC. The non-linearity parameter can be defined for this correlator as, It is easy to see from (6.9) and (6.10) that f N L | γ ζ ζ ∝ 1/ t (k 2 /k 1 ) for θ (s) = π and f N L | γ ζ ζ ∝ s / t (k 2 /k 1 ) for θ (s) = π . It is worth mentioning that this enhancement in the mixed correlator is reflected on the quadrupole moment and might be of interest for future CMB missions. This makes the present investigation for this correlator relevant from the point of view of observations.

Summary and outlook
Recently proposed TCC put stringent theoretical constraints on the amplitude of 2-point correlation function of PGW the so-called tensor-to-scalar ratio r . TCC bound reveals that with the usual BD vacuum, a detectable PGW signal can not be produced from vacuum fluctuation during inflation. However, subsequently, it was found that if the initial states are NBD both for scalar and tensor modes, one can relax the bound on r sufficiently bringing it back to the observational limit of next generation CMB missions with the help of the NBD parameters s and t (functions of scalar and tensor Bogolyubov coefficients, respectively). In this article we have explored the 3-point statistics of tensor modes respecting the constraints coming from TCC, thereby considering NBD states for both scalar and tensor modes to study any possible enhancement in the signal strength. Using a modelindependent framework of EFT of inflation that helps us keep the analysis and results more or less generic, we have calculated all possible correlators (auto and mixed) related to tensor non-Gaussianities, followed by proposing the possible templates for the non-linearity parameters f N L for different relevant shapes. We have also tried to explore if at all any of the bispectra could be observationally relevant for future CMB missions by simply calculating the possible upper bounds for each by tuning the parameters under consideration.
Our analysis reveals that the auto correlator γ γ γ generated from Einstein term R does not get significantly enhanced even with a choice of large Bogolyubov coefficient β s/t , since the definition of the amplitude of bispectrum always keeps it proportional to r 2 . The amplitude of the auto correlator due to the higher derivative EFT operator too is highly suppressed when one considers TCC. Thus, the prospects of detecting the tensor auto corrector are almost nil, at least in next generation CMB missions. We have also analyzed the mixed γ γ ζ correlator which is a good probe of sound speed of PGW. We have found that the shape of the correlator is modified due to the NBD state and gets peaked at flattened limit modified by the sound speed of perturbation. Also the amplitude of this correlator from the non-oscillatory part can be estimated as r t / s which can be large in the large s limit and in the flattened limit modified by sound speed the signal strength gets an additional enhancement from η 0 . Thus, this correlator can be relevant for future CMB missions. The behaviour of the other mixed correlator, namely, the γ ζ ζ generated due to space-time diffeomorphism breaking is also explored. We have used the approach of EFT of space-time diffeomorphism breaking for this scenario. This type of correlator can be a probe of spatial diffeomorphism breaking during inflation. We have found that using NBD state there will be an extra contribution on top of the BD contribution. In squeezed limit this contribution is proportional to k s /k l , where k s and k l are respectively the short and the long wavelength that can be probed by any CMB mission under consideration. The amplitude of this signal is highly dependent on the phase factor given by the θ (s) parameter. For any choice of θ (s) = π , the amplitude of this correlator can be significantly enhanced, thereby making it relevant to future CMB missions.
In a nutshell, if one begins with the theoretical constraints in the light of TCC for slow roll, single field inflation, the nonlinearity parameter f N L for tensor auto correlation would still be undetectable whereas the corresponding parameters for mixed correlations might be a point of interest for future CMB missions. Our results are more or less generic as we investigated for he parameters from a model-independent framework of EFT.
We must, however, admit that our analysis in this article is mostly theoretical and in order to find out the prospects of tensor non-Gaussianities in CMB, we made use of the simplistic numerical calculations based on current bounds given by Planck 2018 as well as the future goals of next generation missions like CMB-S4, LiteBIRD and COrE. A detailed forecast with particular mission is required to analyze the finer details of the results as well as to estimate the parameters with error in an observation-first approach. Such an analysis is beyond the scope of the present paper. We hope to be back with that in near future.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

A The oscillatory part of the three-point functions
The oscillatory part of the three-point functions are presented here. For the auto correlators the oscillatory part of the R-part can be written as,