Tensor Non-Gaussianity Search: Current Status and Future Prospects
- Department of General Education, National Institute of Technology, Kagawa College, Takamatsu, Japan
Primordial gravitational waves (GWs) are said to be a smoking gun in cosmic inflation, while, even if they are detected, the specification of their origins are still required for establishing a true inflationary model. Testing non-Gaussianity in the tensor-mode anisotropies of the cosmic microwave background (CMB) is one of the most powerful ways to identify sources of GW signals. In this paper, we review studies searching for tensor non-Gaussianities employing the CMB bispectrum and forecast future developments. No significant signal has so far been found from temperature and E-mode polarization data, while orders-of-magnitude improvements in detection limits can be achieved by adding the information of B-mode polarization. There is already an established methodology for bispectrum estimation, which encourages a follow-up investigation with next-decadal CMB B-mode surveys.
Recent cosmic-variance-limited-level measurements of the scalar sector in primordial fluctuations using the temperature and E-mode polarization fields of the cosmic microwave background (CMB), indicate that the Universe experienced an inflationary expansion at very early stages (Aghanim et al., 2018; Akrami et al., 2018). Now, compelling evidence is expected to lie in the tensor sector; namely, the primordial gravitational wave (GW). Even without any particular source, primordial GWs arise naturally from inflationary vacuum fluctuations. Since their amplitude directly reflects the inflationary energy scale, various observational projects aim at hunting them. There, the measurements of the large-scale CMB B-mode polarization have attracted everyone's attention as it is a distinctive observable of the tensor mode (see e.g., Guzzetti et al., 2016; Kamionkowski and Kovetz, 2016 and references therein).
In order to investigate the primordial tensor sector using the CMB anisotropies, the power spectrum (2-point correlation) is primarily employed. In single-field inflation with Einstein gravity, GWs are nearly Gaussian and hence the higher-order correlations vanish (Acquaviva et al., 2003; Maldacena, 2003). In this case, only the power spectrum becomes an informative observable. This paper, however, will address the potential of the bispectrum (3-point correlation) as a probe of inflationary physics.
One possible way, making it more informative, is the addition of some extra source fields. In the presence of the gauge field, the anisotropic stress fluctuations are formed, generating GWs. Because this is a non-linear process, the resulting GWs are non-Gaussian (NG) and therefore a large tensor bispectrum can be realized. In particular, the production of NG GWs becomes more efficient when the gauge field is coupled to the axion (Namba et al., 2016; Agrawal et al., 2018a). Moreover, production due to higher-spinning particles has been studied intensively in recent years (Goon et al., 2018; Dimastrogiovanni et al., 2019). There is also the possibility of a post-inflationary generation of NG GWs due to magnetic fields (Shiraishi et al., 2011c). Another possible way to produce large tensor NGs comes from non-trivial non-linear gravitational interactions predicted in some modified theories of gravity (e.g., Gao et al., 2011, 2013; Maldacena and Pimentel, 2011; Akita and Kobayashi, 2016; Bartolo and Orlando, 2017; Domenech et al., 2017; Naskar and Pal, 2018; Anninos et al., 2019; Ozsoy et al., 2019). Then, the induced bispectra can have distinctive shapes compared with that from the usual Einstein term.
In this sense, the CMB tensor-mode bispectrum is a key indicator of the inflationary particle content and/or high energy gravity. Since a general formalism was established in Shiraishi et al. (2010, 2011b), various signatures of various theoretical models in the CMB bispectra have been studied [see (Shiraishi et al., 2011c, 2012; Shiraishi, 2012, 2013) for primordial magnetic field (PMF) models (Shiraishi et al., 2013, 2016; Namba et al., 2016), for the axion inflation ones (Shiraishi et al., 2011a,b; Domenech et al., 2017; Tahara and Yokoyama, 2018; Bartolo et al., 2019), for modified gravity ones and (Meerburg et al., 2016; Kothari and Nandi, 2019) for other motivations]. One of the most interesting phenomena discovered there (unseen in the scalar bispectrum analysis) is that the tensor bispectrum can yield the non-vanishing signal in not only even but also odd ℓ1 + ℓ2 + ℓ3 multipoles. The GW bispectra generate the CMB auto and cross bispectra including the B-mode polarization. In the usual theories where the GW bispectra are parity-invariant, the odd ℓ1 + ℓ2 + ℓ3 signal can arise in TTB, TEB, EEB, and BBB (Shiraishi, 2013; Meerburg et al., 2016; Domenech et al., 2017; Tahara and Yokoyama, 2018). Moreover, parity-odd GW bispectra motivated by parity-breaking theories can also induce the non-zero odd ℓ1 + ℓ2 + ℓ3 signal in TTT, TTE, TEE, TBB, EEE, and EBB (Shiraishi et al., 2011a, 2013, 2016; Shiraishi, 2012, 2016;Bartolo et al., 2019).
Now, any type of theoretical bispectrum template is testable with the temperature and E/B-mode polarization data by use of a general bispectrum estimator for the even (Fergusson et al., 2010, 2012; Fergusson, 2014; Shiraishi et al., 2019) and odd ℓ1 + ℓ2 + ℓ3 domain (Shiraishi et al., 2014, 2019). The magnitudes of some representative tensor NG templates have already been estimated from the temperature map of WMAP (Shiraishi and Sekiguchi, 2014; Shiraishi et al., 2015, 2018) and the temperature and E-mode polarization from Planck (Ade et al., 2016a,b; Akrami et al., 2019). Absence of NG reported there, constrains the various inflationary models mentioned above.
In this paper, besides reviewing previous studies, we discuss future prospects of the tensor NG search assuming the detection of B-mode polarization in next-generation CMB experiments. We found that, in terms of observational limits on the size of some scale-invariant templates, an improvement up to three orders of magnitude is expected assuming LiteBIRD-level sensitivities (Hazumi et al., 2012; Matsumura et al., 2013; Matsumura et al., 2016).
This paper is organized as follows. In the next section, we summarize theoretical scenarios predicting characteristic primordial tensor bispectra. In section 3, we compute the CMB bispectra and show their behavior. In section 4, we present the current observational constraints on tensor NGs and discuss future prospects assuming B-mode polarization surveys. The final section is devoted to the conclusions of this paper.
2. Theoretical Motivations
In this section, we briefly review the primordial tensor bispectra predicted in some theoretical scenarios. For convenience, let us work with the helicity basis; thus, the primordial perturbation is represented as (Shiraishi et al., 2011b). Here, λ denotes the helicity and takes 0, ±1, and ±2 for the scalar, vector, and tensor modes, respectively. In practice, we identify and with the curvature perturbation ζk and the GW , respectively.
Assuming statistical homogeneity, the primordial bispectrum takes the form:
Then, let us split this in two:
In what follows, how the shape, represented by , and the angular structure, represented by , changes depending on the model is argued.
2.1. Extra Sources
Efficient tensor NG production is realized by adding extra source fields to the theory. In what follows, we discuss the vector field as a source.
2.1.1. Inflationary Axion-Gauge Coupling
A characteristic tensor NG is realized in inflationary models involving the coupling between the axion ϕ and a gauge field like , with F and the gauge field strength tensor and its dual, respectively (Cook and Sorbo, 2013; Maleknejad, 2016; Namba et al., 2016; Dimastrogiovanni et al., 2017; Agrawal et al., 2018a,b). In this case, the chirality of the gauge field is transferred into the GW sector and hence only the plus mode of GW, h(+2), survives. The resulting production is efficient at specific scales, and the tensor-tensor-tensor bispectrum peaks for equilateral configurations, k1 ~ k2 ~ k3 (Cook and Sorbo, 2013).
The detailed spectral feature varies with the shape of the coupling f(ϕ) or the axion potential (Namba et al., 2016; Agrawal et al., 2018b). For example, adopting f(ϕ) ∝ ϕ and a nearly flat potential, one can obtain a scale-invariant shape as (Shiraishi et al., 2013).
where is the polarization tensor defined in the helicity basis obeying , , and (Shiraishi et al., 2011b), and is the usual scalar equilateral bispectrum template, given by:
Here is the amplitude of the curvature power spectrum. A newly introduced parameter quantifies the relative size of the λ = +2 bispectrum to and satisfies1
The expected size of the GW bispectrum, of course, depends on the case. If the axion is identified with the inflaton field, the amount of scalar production exceeds the GW one (Barnaby et al., 2011). Then, the scalar bispectrum measurement gives much tighter constraints on the model (Ade et al., 2016b; Akrami et al., 2019). In contrast, in multifield models where the inflaton field and the axion coexist, more effective GW production occurs. In the model where the axion is coupled to the U(1) gauge field, a characteristic bump appears in the GW bispectrum, and, depending on its location, it is detectable by the CMB BBB bispectrum measurement (Namba et al., 2016; Shiraishi et al., 2016). However, at the same time, a similar bump in the GW power spectrum is measured with a higher signal-to-noise ratio from the CMB BB power spectrum. The strongest constraints are obtained through the GW power spectrum measurement, and the GW bispectrum provides complementary information. In the model including a SU(2)-gauge coupling, the GW bispectrum can dominate the scalar one, and more interestingly, in some regions of the parameter space, the GW bispectrum has high detectability compared to the GW power spectrum (Agrawal et al., 2018a,b). The GW bispectrum is nearly scale-invariant and therefore parameterized by , which is related to the tensor-to-scalar ratio r and the energy density fraction of the SU(2) gauge field ΩA according to Agrawal et al. (2018a).
This indicates that can take detectably large values, i.e., , when ΩA is smaller than r2.
2.1.2. Primordial Magnetic Fields
If PMFs, Bi, are created via some non-trivial mechanism before inflation, they are stretched beyond the horizon by the inflationary expansion and there remain relic fields at large scales. Such a scenario has been widely argued to explain galactic or extragalactic magnetic fields observed at present (for reviews see Widrow, 2002; Kulsrud and Zweibel, 2008). In this case, PMFs form the anisotropic stress fluctuations, and they source the GW at super-horizon scales until neutrino decoupling. The induced GW is called a passive mode and behaves as the initial condition of the CMB fluctuations after reentering the horizon (Shaw and Lewis, 2010).
Assuming Gaussianity of PMFs, the induced GW, , becomes a chi-square random field (Brown and Crittenden, 2005). When the PMF power spectrum has a scale-invariant shape, in analogy with the usual local type NG, the tensor-tensor-tensor bispectrum has the squeezed shape (Shiraishi et al., 2011c, 2012):
where the usual scalar local bispectrum template reads:
and is given in the same form as Equation (4). A newly introduced parameter represents the relative size of the λ = +2 bispectrum to and satisfies:
As the GW bispectrum has a sextic dependence on PMFs, takes the form:
where B1Mpc is the PMF strength smoothed on 1Mpc. This means that is realized if nano-Gauss-level PMFs exist. Note that, in the same mechanism, the scalar local type bispectrum is generated, however, it is subdominant compared with the GW one (Shiraishi, 2013).
This is the result under the absence of a helical (parity-odd) term in the PMF power spectrum, however, the bispectrum shape is modified if it exists (Shiraishi, 2012).
2.2. Modified Gravity
The above scenarios rely on some extra source fields. Let us discuss another NG GW production through the modification of Einstein gravity.
2.2.1. Weyl Gravity
In the context of the Weyl gravity, there may exist the Chern-Simons term , where W is the Weyl tensor and is its dual. Like the axion inflation models, this term also sources chiral GWs (Lue et al., 1999; Alexander and Martin, 2005). The tensor-tensor-scalar bispectrum, whose angular structure reads:
surpasses the tensor-tensor-tensor and tensor-scalar-scalar ones in amplitude (Bartolo and Orlando, 2017). Due to the parity-violating nature in the GW sector, the helical bispectrum obeys:
One can also construct the cubic actions f(ϕ)W3 and (Maldacena and Pimentel, 2011; Shiraishi et al., 2011a; Soda et al., 2011). The tensor-tensor-tensor bispectra generated from both terms have identical angular dependence as Equation (4). In contrast, the helicity dependence is completely different (Shiraishi et al., 2011a): the former gives rise to the non-helical (parity-even) contribution:
while the latter induces the helical (parity-odd) one:
2.2.2. Massive Gravity
In single-field slow-roll inflation with Einstein gravity, the tensor-scalar-scalar bispectrum is affected by slow-roll suppression (Maldacena, 2003). However, in a massive gravity model (Domenech et al., 2017), the size of the bispectrum is controlled by the graviton mass. The bispectrum, taking the form:
with kt ≡ k1 + k2 + k3, is amplified at long-wavelength tensor and short-wavelength scalar configurations (k1 ≪ k2 ~ k3). A newly introduced parameter denotes the relative size of the λ = +2 bispectrum to , obeying:
With an upper bound on the tensor-to-scalar ratio r ≲ 0.1, the prediction in Einstein gravity, , indicates . In contrast, massive gravity model modifies this as:
and can then be realized depending on the strength parameter of a non-linear interaction λsst. Note that the spectral index of the GW power spectrum is also a possible observable of this model (Domenech et al., 2017) but has been unconstrained so far.
3. CMB Bispectra From Tensor Non-gaussianities
Next, let us discuss the CMB bispectra generated from GW NGs.
3.1. General Formalism for CMB Angular Bispectrum
The CMB field is characterized by temperature and two linear polarizations called E and B modes. All of these are spin-0 fields, and the temperature and E-mode fields have even parity, while the B mode is parity-odd. Since temperature (X = T), E-mode (X = E), and B-mode (X = B) fields are distributed on the 2D sphere, they can be expanded using the spherical harmonic basis according to
where we have assumed the spatial flatness of the Universe. The spherical harmonic coefficients are given by the sum of the scalar, vector, and tensor modes as:
Note that the scalar mode cannot generate a B-mode signal due to zero helicity. The coefficients of the scalar (z = s), vector (z = v), and tensor (z = t) modes have the structure:
where x changes depending on the parity of the CMB field according to x = 0 for X = T, E and x = 1 for X = B, and recall that λ represents helicity as λ = 0 for z = s, λ = ±1 for z = v, and λ = ±2 for z = t. The linear transfer function represents the time evolution of the CMB fluctuations originating from the primordial perturbation . Using Equation (22), one can form the CMB bispectrum as:
The primordial curvature perturbation and the primordial GW act as initial conditions of the scalar and tensor CMB anisotropies, respectively; thus, and . Employing the harmonic expansion:
the above formula is rewritten as:
This is further simplified once the explicit formula of is given. The computation procedure is as follows:
1. Expand using the (spin-weighted) spherical harmonics in terms of , , and .
2. Perform the angular integrals, , , and , of the resultant spherical harmonics and convert them into products of Wigner symbols.
3. Simplify the resultant products of the Wigner symbols by adding angular momenta.
If the primordial bispectrum respects rotational invariance, via the above computation, the Wigner 3j symbol is singled out in and hence the resulting CMB bispectrum takes the form:
In this case, the non-zero signal is confined to the tetrahedral domain:
3.2. Practical Examples
In the following, we demonstrate the CMB bispectrum computation in the context of the practical cases discussed in section 2.
3.2.1. Three Tensors
First, we discuss the CMB bispectra sourced by primordial tensors where we examine the axion model template (3) (Shiraishi et al., 2013), the W3 model (14) (Shiraishi et al., 2011a), the model (15) (Shiraishi et al., 2011a), and the PMF model (8) (Shiraishi et al., 2011c). Since takes the identical form, can be computed in a similar way.
For the first step, in order to simplify the R.H.S. of Equation (26), we expand and the Dirac delta function:
For the second step, the angular integrals of the resultant spherical harmonics are performed, such that:
As a final step, the summation of the Wigner 3j symbols appearing above over angular momenta is computed according to:
Consequently, we obtain:
Here, the prefactor is due to rotational invariance of the primordial bispectrum. Inserting this into Equation (25) and simplifying the resultant equation, we find that the CMB bispectra take the rotational-invariant form (27) with:
The ℓ-space domain containing the non-zero signal is determined by the summation over the helicities λ1, λ2, and λ3. In the model, performing the summation of Equation (15) leads to the non-vanishing condition: (Shiraishi et al., 2011a)
In the axion model, however, the helicity dependence in Equation (3) does not yield any restriction, allowing a non-vanishing signal in the whole tetrahedral domain (Shiraishi et al., 2013).
The left and center panels of Figure 1 show the intensity distributions of the temperature bispectra from the axion model template (3) and the PMF model template (8), respectively. As expected, it is confirmed that the signal comes mostly from the equilateral (ℓ1 ~ ℓ2 ~ ℓ3) and squeezed (ℓ1 ≪ ℓ2 ~ ℓ3, ℓ2 ≪ ℓ3 ~ ℓ1 and ℓ3 ≪ ℓ1 ~ ℓ2) configurations, respectively. The decaying nature for ℓ≳100 due to the lack of the tensor-mode integrated Sachs-Wolfe contribution (Shiraishi et al., 2011c, 2013) is also visually apparent.
Figure 1. Intensity distributions of the CMB temperature bispectra from the axion model template (3), the PMF model one (8), and the massive gravity model one (16) in the ℓ-space tetrahedral domain where the axes correspond to ℓ1, ℓ2, and ℓ3, respectively. Here, the bispectra are rescaled using a constant Sachs-Wolfe template (Fergusson et al., 2010). Dense red (blue) color expresses a larger positive (negative) signal. The current and future observational limits on , , and are discussed in sections 4.2 and 4.3, respectively.
3.2.2. Two Tensors and One Scalar
Here, we compute the CMB tensor-tensor-scalar bispectrum from the model (13) (Bartolo et al., 2019). The angular dependence in (12) is decomposed according to:
By means of the above methodology, we obtain (Bartolo et al., 2019):
Performing the summation over λ1 and λ2 with Equation (13), we find that the non-vanishing signal obeys (Bartolo et al., 2019):
3.2.3. One Tensor and Two Scalars
We here compute the CMB bispectrum from the massive gravity model template (16) (Domenech et al., 2017). The spherical harmonic expansion of (17) then reads:
The right panel of Figure 1 describes the intensity distribution of the resultant temperature bispectrum showing explicitly that it is amplified in the squeezed limit. Differently from the case, the squeezed signal survives even for ℓ≳100 because of the presence of the scalar-mode acoustic oscillation.
4. Measurements of CMB Tensor-Mode Bispectra
In this section, we examine the present and future constraints on GW NGs through the CMB bispectrum measurements.
4.1. Optimal Tensor-Mode Bispectrum Estimation
Let us begin by explaining how to estimate the overall magnitude of the primordial bispectrum, dubbed as fNL, from the CMB data.
An optimal fNL estimator takes the form (Komatsu et al., 2005):
where Cℓ is the CMB angular power spectrum, is the theoretical template of the angle-averaged CMB bispectrum one wants to measure [corresponding to in Equation (27) for fNL = 1],
is the bispectrum reconstructed from the observed aℓm, and
is the Fisher matrix. The bestfit value and the error bar on fNL is obtained computing from the observed data and the simulation maps, respectively. The straightforward computation of Equation (44) is not a wise manner since the numerical operation required here is quite time-consuming and practically unfeasible for ℓmax≳100.
To reduce the computational cost, factorizing the summation as is a tested approach2. This idea was implemented in Komatsu et al. (2005) for the first time, enabling the analysis for originally separable templates such as local, equilateral and orthogonal NGs. Later, it was generalized by developing the modal decomposition technique of any originally non-separable (Fergusson et al., 2010, 2012). The original version of the modal methodology is for the analysis of even ℓ1 + ℓ2 + ℓ3 domain, while the spin-weighted version developed recently (Shiraishi et al., 2014, 2015) can also cover odd ℓ1 + ℓ2 + ℓ3 domain. As seen in section 3, the tensor-mode CMB bispectrum is generally non-separable and does not always vanish in the odd ℓ1 + ℓ2 + ℓ3 domain. In this sense, the modal methodology is an indispensable tool for the tensor NG search. In fact, except for the results obtained brute forcing the ℓ-space summation with ℓmax = 100 (Shiraishi and Sekiguchi, 2014), all observational constraints reported so far have been obtained by it Shiraishi et al. (2015, 2018); Ade et al. (2016a,b); Akrami et al. (2019). In what follows, for simplicity, we explain the modal methodology for the auto bispectra (TTT, EEE, and BBB) based on Fergusson et al. (2010, 2012), Shiraishi et al. (2014, 2015), while one can deal with the cross bispectra (TTE, TTB, TEE, TEB, TBB, EEB, and EBB) in the very similar way (Fergusson, 2014; Shiraishi et al., 2019).
Let us introduce the reduced bispectrum according to:
Here and hereinafter, we follow the notation for the symmetric operation The spined weight (x, y, z) is fixed to be, e.g., (0, 0, 0) and (1, 1, −2) for even and odd ℓ1 + ℓ2 + ℓ3 analysis, respectively. The reduced bispectrum of the theoretical template is decomposed in the odd/even ℓ1 + ℓ2 + ℓ3 domain separately, according to:
where the separable modal basis Qijk is composed of the products of the eigenfunctions qi(ℓ) ∈ ℝ as:
The σ factor, defined by:
makes the L.H.S of Equation (49) real, so the modal coefficients are always real. The convergence speed of the modal decomposition relies on the choice of the vℓ weighting and qi(ℓ). Employing the Qijk templates, we define a matrix according to:
where the triples ijk in Qijk are labeled by means of a single index n, and the bracket 〈⋯〉o/e denotes the summation in terms of the odd/even ℓ1 + ℓ2 + ℓ3 domain as:
Then, the modal coefficients are computed according to the inverse operation:
Using Equation (49) and the identity:
the sum over ℓ1, ℓ2, and ℓ3 in the estimator (44) is rewritten into the sum over finite eigenmodes as:
where the coefficients, reading:
are computed from the filtered maps:
The most time-consuming task is the computation of 〈⋯〉o/e in , however, it requires, at most, operations (assuming that the modal decomposition converges within a reasonable time), and hence the separable estimator (56) makes the tensor bispectrum analysis feasible.
4.2. Current Observational Limits
Table 1 summarizes the current CMB limits on , , and obtained from the WMAP and Planck maps3. Regarding and , the consistency between the WMAP T only results and Planck T only ones is confirmed there. Note that the unimproved constraints are expected because, as inferred from Figure 1, decays rapidly for ℓ≳100 and hence the estimator sum (44) saturates well below the WMAP resolution (ℓmax = 500). In contrast, for the current constraint from WMAP is expected to be much improved by Planck since the squeezed-limit signal survives even for high ℓ (Shiraishi et al., 2011b; Domenech et al., 2017). We also confirm the improvement of the constraint on by analyzing the temperature and E-mode polarization data jointly.
As seen in Table 1, no significant deviation from Gaussianity has been found so far.
4.3. Future Prospects
Finally, we discuss future prospects of detecting these tensor NGs by adding B-mode polarization to the data analysis. Here, we evaluate expected 1σ errors and from BBB, and from BTT, through the computation of the Fisher matrix. The covariance matrices [corresponding to the denominator of Equation (46)] are given by and , respectively. Note that these are estimated under the diagonal covariance matrix approximation. The B-mode power spectrum is computed by summing up the primordial contribution parameterized by the tensor-to-scalar ratio r, the lensing B-mode one and the experimental noise spectrum as . Then, let us examine three cleanliness levels of the B-mode data: a perfectly-delensed and noiseless full-sky case (i.e., ), a non-delensed and noiseless full-sky one (i.e., Nℓ = 0), and a non-delensed realistic case where experimental uncertainties due to beam, noise, mask and residual foreground in LiteBIRD (Hazumi et al., 2012; Matsumura et al., 2013, 2016) are considered4. The first example provides the theoretical limits. In this case, the size of the covariance matrix is determined by r alone and, accordingly, the errors simply scale like:
Comparing this with Equations (7) or (19), one can see that the error on a model parameter of the axion model, , or that of the massive gravity model, λsst, is proportional to r−1/2.
The left panel of Figure 2 describes the ℓmax dependence of ΔfNL in the perfectly-delensed and non-delensed cases with r = 10−3 and Nℓ = 0. From the former results (corresponding to the dot-dashed lines), the sensitivity improvement by increasing ℓmax is confirmed. From the latter results (corresponding to the solid lines), it is visually apparent that and immediately saturate because dominates over for . In contrast, is free from such a degradation since the dominant signal in the Fisher matrix comes from large-scale B modes, more precisely, long-wavelength B and short-wavelength T squeezed configurations () (Meerburg et al., 2016; Domenech et al., 2017). However, of course, would also saturate at very small scales, i.e., ℓT≳3000, where the scalar-mode lensing contamination dominates . Higher-order lensing contributions introduce non-vanishing off-diagonal components in the covariance matrix; thus, the above simple Fisher matrix computation would no longer be credible. In the left panel of Figure 2, the results of the LiteBIRD-like experiment are not shown, however, they have similar ℓmax scalings and the slightly larger overall sizes in comparison with the no-delensed and noiseless full-sky results.
Figure 2. Expected 1σ errors: (red lines) and (blue lines) from BBB, and (green lines) from BTT, as a function of the maximum multipole number ℓmax (left panel) and the tensor-to-scalar ratio r (right panel). The left panel is depicted adopting r = 10−3, and in the right panel is computed with ℓmax = 2000. The linestyle discriminates the cleanliness level of the B-mode data: a perfectly-delensed and noiseless full-sky case (dot-dashed lines in left panel), a non-delensed and noiseless full-sky one (solid lines in both panels), and a LiteBIRD-like realistic one (dashed lines in right panel).
In the right panel of Figure 2, the r dependence of ΔfNL for the no-delensed and noiseless case (corresponding to the solid lines) and a LiteBIRD-like experiment (corresponding to the dashed lines) is presented. For large r, dominates over within a wide range of ℓ and hence the r dependence reaches the ideal case (59). On the other hand, for decreasing r, or Nℓ starts to dominate and errors converge. We find that LiteBIRD could measure an signal of , , or .5 If detected, as discussed in section 2, we will have compelling evidence for a deviation from Einstein gravity or the existence of extra source fields, which will help establish an improved picture of the early Universe.
Finally, it is worth mentioning the impact of late-time secondary bispectra. Via gravitational lensing, primordial temperature and E-mode polarization fields induce secondary BTT signals (Hu, 2000; Lewis et al., 2011). Fortunately, this has a very weak correlation with the bispectrum template (Domenech et al., 2017). Similarly, secondary BBB is also induced, while this is subdominant compared to primordial BBB at interesting or (Shiraishi et al., 2016). Besides of these leading-order contributions, higher-order ones exist due to post-Born lensing (Marozzi et al., 2016; Pratten and Lewis, 2016). Secondary polarized bispectra can also be generated via the Sunyaev-Zel'dovich effect (Sunyaev and Zeldovich, 1972) and extragalactic foreground (Jung et al., 2018; Coulton and Spergel, 2019). These impacts are still non-trivial and should therefore be solved for more precise discussions.
The search for tensor NG is theoretically well-motivated, and the methodology for testing with the CMB temperature and E/B-mode polarization data has already been established. Observational constraints on some templates from the WMAP temperature and the Planck temperature and E-mode polarization data do exist. Any significant signal has not yet been discovered, implying the smallness of tensor NGs. There is still a chance of approaching such a small signal by including B-mode polarization in the data analysis. As we have found, an signal of , , or would be detectable in LiteBIRD; current constraints are expected to be updated in the next decade.
Throughout this paper, we have focused mainly on scale-invariant cases, while specific scale dependences can be generated for tensor NGs depending on the shapes of the inflationary potentials and non-linear interactions. The detectability is then enhanced at specific multipoles (Shiraishi et al., 2016) or outside the CMB scales. In the latter case, the information from other probes such as an intereferometric GW survey (e.g., Bartolo et al., 2018; Tsuneto et al., 2019), the galaxy one (e.g., statistical anisotropy in the galaxy power spectrum Jeong and Kamionkowski, 2012) and the 21-cm line one also becomes indispensable. A more comprehensive analysis based on multi-wavelength observations remains an interesting and important future issue.
The author confirms being the sole contributor of this work and has approved it for publication.
This work was supported by the JSPS Grant-in-Aid for Research Activity Start-up Grant Number 17H07319 and the JSPS Grant-in-Aid for Early-Career Scientists Grant Number 19K14718. Numerical computations were in part carried out on Cray XC50 at the Center for Computational Astrophysics, National Astronomical Observatory of Japan.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
I thank the editors of Status and Prospects of Cosmic Microwave Background Analysis and especially, Maurizio Tomasi and Michele Liguori for their invitation to contribute to this special issue and for all their support. I also appreciate the careful reviews of the manuscript by the referees.
4. ^In the third case, we compute Nℓ by taking into account the contamination due to residual foregrounds in the same manner as Shiraishi et al. (2016); namely, we assume that foregrounds due to galactic dust emission and synchrotron radiation are subtracted using 9 channels (corresponding to 40-89 GHz and 280-402 GHz) and therefore reduce to the 2% level in CMB maps. Any higher-order contribution is not considered for simplicity
Ade, P. A. R., Aghanim, N., Arnaud, M., Arroja, F., Ashdown, M., Aumont, J., et al. (2016a). Planck 2015 results. XIX. Constraints on primordial magnetic fields. Astron. Astrophys. 594:A19. doi: 10.1051/0004-6361/201525821
Ade, P. A. R., Aghanim, N., Arnaud, M., Arroja, F., Ashdown, M., Aumont, J., et al. (2016b). Planck 2015 results. XVII. Constraints on primordial non-Gaussianity. Astron. Astrophys. 594:A17. doi: 10.1051/0004-6361/201525836
Bartolo, N., Domcke, V., Figueroa, D. G., García-Bellido, J., Peloso, M., Pieroni, M., et al. (2018). Probing non-gaussian stochastic gravitational wave backgrounds with LISA. J. Cosmol. Astropart. Phys. 1811:034. doi: 10.1088/1475-7516/2018/11/034
Bartolo, N., and Orlando, G. (2017). Parity breaking signatures from a Chern-Simons coupling during inflation: the case of non-Gaussian gravitational waves. J. Cosmol. Astropart. Phys. 1707:034. doi: 10.1088/1475-7516/2017/07/034
Bartolo, N., Orlando, G., and Shiraishi, M. (2019). Measuring chiral gravitational waves in Chern-Simons gravity with CMB bispectra. J. Cosmol. Astropart. Phys. 1901:050. doi: 10.1088/1475-7516/2019/01/050
Bucher, M., Racine, B., and van Tent, B. (2016). The binned bispectrum estimator: template-based and non-parametric CMB non-Gaussianity searches. J. Cosmol. Astropart. Phys. 1605:055. doi: 10.1088/1475-7516/2016/05/055
Delabrouille, J., de Bernardis, P., Bouchet, F. R., Achúcarro, A., Ade, P. A. R., Allison, R., et al. (2018). Exploring cosmic origins with CORE: survey requirements and mission design. J. Cosmol. Astropart. Phys. 1804:014. doi: 10.1088/1475-7516/2018/04/014
Dimastrogiovanni, E., Fasiello, M., Tasinato, G., and Wands, D. (2019). Tensor non-gaussianities from non-minimal coupling to the inflaton. J. Cosmol. Astropart. Phys. 1902:008. doi: 10.1088/1475-7516/2019/02/008
Domenech, G., Hiramatsu, T., Lin, C., Sasaki, M., Shiraishi, M., and Wang, Y. (2017). CMB scale dependent non-Gaussianity from massive gravity during inflation. J. Cosmol. Astropart. Phys. 1705:034. doi: 10.1088/1475-7516/2017/05/034
Fergusson, J. R., Liguori, M., and Shellard, E. P. S. (2010). General CMB and primordial bispectrum estimation I: mode expansion, map-making and measures of fNL. Phys. Rev. D82:023502. doi: 10.1103/PhysRevD.82.023502
Gao, X., Kobayashi, T., Shiraishi, M., Yamaguchi, M., Yokoyama, J., and Yokoyama, S. (2013). Full bispectra from primordial scalar and tensor perturbations in the most general single-field inflation model. Progress Theoret. Exp. Phys. 2013:053E03. doi: 10.1093/ptep/ptt031
Gao, X., Kobayashi, T., Yamaguchi, M., and Yokoyama, J. (2011). Primordial non-Gaussianities of gravitational waves in the most general single-field inflation model. Phys. Rev. Lett. 107:211301. doi: 10.1103/PhysRevLett.107.211301
Hazumi, M., Borrill, J., Chinone, Y., Dobbs, M. A., Fuke, H., Ghribi, A., et al. (2012). LiteBIRD: a small satellite for the study of B-mode polarization and inflation from cosmic background radiation detection. Proc. SPIE Int. Soc. Opt. Eng. 8442:844219. doi: 10.1117/12.926743
Jung, G., Racine, B., and van Tent, B. (2018). The bispectra of galactic CMB foregrounds and their impact on primordial non-Gaussianity estimation. J. Cosmol. Astropart. Phys. 1811:047. doi: 10.1088/1475-7516/2018/11/047
Matsumura, T., Akiba, Y., Arnold, K., Borrill, J., Chendra, R., Chinone, Y., et al. (2016). LiteBIRD: mission overview and focal plane layout. J. Low Temperat. Phys. 184, 824–831. doi: 10.1007/s10909-016-1542-8
Namba, R., Peloso, M., Shiraishi, M., Sorbo, L., and Unal, C. (2016). Scale-dependent gravitational waves from a rolling axion. J. Cosmol. Astropart. Phys. 1601:041. doi: 10.1088/1475-7516/2016/01/041
Shiraishi, M., Hikage, C., Namba, R., Namikawa, T., and Hazumi, M. (2016). Testing statistics of the CMB B -mode polarization toward unambiguously establishing quantum fluctuation of the vacuum. Phys. Rev. D94:043506. doi: 10.1103/PhysRevD.94.043506
Shiraishi, M., Nitta, D., Yokoyama, S., and Ichiki, K. (2012). Optimal limits on primordial magnetic fields from CMB temperature bispectrum of passive modes. J. Cosmol. Astropart. Phys. 1203:041. doi: 10.1088/1475-7516/2012/03/041
Shiraishi, M., Nitta, D., Yokoyama, S., Ichiki, K., and Takahashi, K. (2011b). CMB bispectrum from primordial scalar, vector and tensor non-gaussianities. Prog. Theor. Phys. 125, 795–813. doi: 10.1143/PTP.125.795
Shiraishi, M., Nitta, D., Yokoyama, S., Ichiki, K., and Takahashi, K. (2011c). Cosmic microwave background bispectrum of tensor passive modes induced from primordial magnetic fields. Phys. Rev. D83:123003. doi: 10.1103/PhysRevD.83.123003
Shiraishi, M., and Sekiguchi, T. (2014). First observational constraints on tensor non-Gaussianity sourced by primordial magnetic fields from cosmic microwave background. Phys. Rev. D90:103002. doi: 10.1103/PhysRevD.90.103002
Shiraishi, M., Yokoyama, S., Nitta, D., Ichiki, K., and Takahashi, K. (2010). Analytic formulae of the CMB bispectra generated from non-Gaussianity in the tensor and vector perturbations. Phys. Rev. D82:103505. doi: 10.1103/PhysRevD.82.103505
Tsuneto, M., Ito, A., Noumi, T., and Soda, J. (2019). Searching for bispectrum of stochastic gravitational waves with pulsar timing arrays. J. Cosmol. Astropart. Phys. 1903:032. doi: 10.1088/1475-7516/2019/03/032
Keywords: CMB, non-Gaussianity, gravitational wave, inflation, observational constraints, B-mode
Citation: Shiraishi M (2019) Tensor Non-Gaussianity Search: Current Status and Future Prospects. Front. Astron. Space Sci. 6:49. doi: 10.3389/fspas.2019.00049
Received: 24 March 2019; Accepted: 28 June 2019;
Published: 23 July 2019.
Edited by:Maurizio Tomasi, University of Milan, Italy
Reviewed by:Kazuya Koyama, University of Portsmouth, United Kingdom
Daan Meerburg, University of Groningen, Netherlands
Copyright © 2019 Shiraishi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Maresuke Shiraishi, email@example.com