Primordial non-Gaussianity and the CMB bispectrum

We present a new formalism, together with efficient numerical methods, to directly calculate the CMB bispectrum today from a given primordial bispectrum using the full linear radiation transfer functions. Unlike previous analyses which have assumed simple separable ansatze for the bispectrum, this work applies to a primordial bispectrum of almost arbitrary functional form, for which there may have been both horizon-crossing and superhorizon contributions. We employ adaptive methods on a hierarchical triangular grid and we establish their accuracy by direct comparison with an exact analytic solution, valid on large angular scales. We demonstrate that we can calculate the full CMB bispectrum to greater than 1% precision out to multipoles l<1800 on reasonable computational timescales. We plot the bispectrum for both the superhorizon ('local') and horizon-crossing ('equilateral') asymptotic limits, illustrating its oscillatory nature which is analogous to the CMB power spectrum

Reduced bispectrum seeded by helical primordial magnetic fields

In this paper, we investigate the effects of helical primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) reduced bispectrum. We derive the full three-point statistics of helical magnetic fields and numerically calculate the even contribution in the collinear configuration. We then numerically compute the CMB reduced bispectrum induced by passive and compensated PMF modes on large angular scales. There is a negative signal on the bispectrum due to the helical terms of the fields and we also observe that the biggest contribution to the bispectrum comes from the non-zero IR cut-off for causal fields, unlike the two-point correlation case. For negative spectral indices, the reduced bispectrum is enhanced by the passive modes. This gives a lower value of the upper limit for the mean amplitude of the magnetic field on a given characteristic scale. However, high values of IR cut-off in the bispectrum, and the helical terms of the magnetic field relaxes this bound. This demonstrates the importance of the IR cut-off and helicity in the study of the nature of PMFs from CMB observations.Comment: 39 pages, figures improved, typos corrected. Version accepted for publication in JCA

Non-Gaussianity of the Cosmic Infrared Background anisotropies II : Predictions of the bispectrum and constraints forecast

Using a full analytical computation of the bispectrum based on the halo model together with the halo occupation number, we derive the bispectrum of the cos- mic infrared background (CIB) anisotropies that trace the clustering of dusty-star- forming galaxies. We focus our analysis on wavelengths in the far-infrared and the sub-millimeter typical of the Planck/HFI and Herschel/SPIRE instruments, 350, 550, 850, and 1380 um. We explore the bispectrum behaviour as a function of several models of evolution of galaxies and show that it is strongly sensitive to that ingredient. Contrary to the power spectrum, the bispectrum, at the four wavelengths, seems dominated by low redshift galaxies. Such a contribution can be hardly limited by applying low flux cuts. We also discuss the contributions of halo mass as a function of the redshift and the wavelength, recovering that each term is sensitive to a different mass range. Furthermore, we show that the CIB bispectrum is a strong contaminant of the Cosmic Microwave Background bispectrum at 850 um and higher. Finally, a Fisher analysis of the power spectrum, bispectrum alone and of the combination of both shows that degeneracies on the HOD parameters are broken by including the bispectrum information, leading to tight constraints even when including foreground residuals.Comment: 16 pages, Accepted by MNRA

Gridded and direct Epoch of Reionisation bispectrum estimates using the Murchison Widefield Array

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.Comment: 19 pages, 10 figures, accepted for publication in PAS

Observed parity-odd CMB temperature bispectrum

Parity-odd non-Gaussianities create a variety of temperature bispectra in the cosmic microwave background (CMB), defined in the domain: $\ell_1 + \ell_2 + \ell_3 = {\rm odd}$. These models are yet unconstrained in the literature, that so far focused exclusively on the more common parity-even scenarios. In this work, we provide the first experimental constraints on parity-odd bispectrum signals in WMAP 9-year temperature data, using a separable modal parity-odd estimator. Comparing theoretical bispectrum templates to the observed bispectrum, we place constraints on the so-called nonlineality parameters of parity-odd tensor non-Gaussianities predicted by several Early Universe models. Our technique also generates a model-independent, smoothed reconstruction of the bispectrum of the data for parity-odd configurations.Comment: 17 pages, 4 figures, 1 table. Accepted for publication in JCA

The shape of the CMB lensing bispectrum

Lensing of the CMB generates a significant bispectrum, which should be detected by the Planck satellite at the 5-sigma level and is potentially a non-negligible source of bias for f_NL estimators of local non-Gaussianity. We extend current understanding of the lensing bispectrum in several directions: (1) we perform a non-perturbative calculation of the lensing bispectrum which is ~10% more accurate than previous, first-order calculations; (2) we demonstrate how to incorporate the signal variance of the lensing bispectrum into estimates of its amplitude, providing a good analytical explanation for previous Monte-Carlo results; and (3) we discover the existence of a significant lensing bispectrum in polarization, due to a previously-unnoticed correlation between the lensing potential and E-polarization as large as 30% at low multipoles. We use this improved understanding of the lensing bispectra to re-evaluate Fisher-matrix predictions, both for Planck and cosmic variance limited data. We confirm that the non-negligible lensing-induced bias for estimation of local non-Gaussianity should be robustly treatable, and will only inflate f_NL error bars by a few percent over predictions where lensing effects are completely ignored (but note that lensing must still be accounted for to obtain unbiased constraints). We also show that the detection significance for the lensing bispectrum itself is ultimately limited to 9 sigma by cosmic variance. The tools that we develop for non-perturbative calculation of the lensing bispectrum are directly relevant to other calculations, and we give an explicit construction of a simple non-perturbative quadratic estimator for the lensing potential and relate its cross-correlation power spectrum to the bispectrum. Our numerical codes are publicly available as part of CAMB and LensPix.Comment: 32 pages, 10 figures; minor changes to match JCAP-accepted version. CMB lensing and primordial local bispectrum codes available as part of CAMB (http://camb.info/

The CMB Bispectrum

We use a separable mode expansion estimator with WMAP data to estimate the bispectrum for all the primary families of non-Gaussian models. We review the late-time mode expansion estimator methodology which can be applied to any non-separable primordial and CMB bispectrum model, and we demonstrate how the method can be used to reconstruct the CMB bispectrum from an observational map. We extend the previous validation of the general estimator using local map simulations. We apply the estimator to the coadded WMAP 5-year data, reconstructing the WMAP bispectrum using $l<500$ multipoles and $n=31$ orthonormal 3D eigenmodes. We constrain all popular nearly scale-invariant models, ensuring that the theoretical bispectrum is well-described by a convergent mode expansion. Constraints from the local model \fnl=54.4\pm 29.4 and the equilateral model \fnl=143.5\pm 151.2 (\Fnl = 25.1\pm 26.4) are consistent with previously published results. (Here, we use a nonlinearity parameter \Fnl normalised to the local case, to allow more direct comparison between different models.) Notable new constraints from our method include those for the constant model \Fnl = 35.1 \pm 27.4 , the flattened model \Fnl = 35.4\pm 29.2, and warm inflation \Fnl = 10.3\pm 27.2. We investigate feature models surveying a wide parameter range in both the scale and phase, and we find no significant evidence of non-Gaussianity in the models surveyed. We propose a measure \barFnl for the total integrated bispectrum and find that the measured value is consistent with the null hypothesis that CMB anisotropies obey Gaussian statistics. We argue that this general bispectrum survey with the WMAP data represents the best evidence for Gaussianity to date and we discuss future prospects, notably from the Planck satellite

Position-dependent power spectrum of the large-scale structure: a novel method to measure the squeezed-limit bispectrum

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called "squeezed configurations," in which one wavenumber, say $k_3$, is much smaller than the other two, i.e., $k_3\ll k_1\approx k_2$. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the "position-dependent power spectrum," i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological $N$-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.Comment: 23 pages, 6 figures; dependence on cosmological parameters added; JCAP accepte