The shape of primordial non-Gaussianity and the CMB bispectrum

We present a set of formalisms for comparing, evolving and constraining primordial non-Gaussian models through the CMB bispectrum. We describe improved methods for efficient computation of the full CMB bispectrum for any general (non-separable) primordial bispectrum, incorporating a flat sky approximation and a new cubic interpolation. We review all the primordial non-Gaussian models in the present literature and calculate the CMB bispectrum up to l <2000 for each different model. This allows us to determine the observational independence of these models by calculating the cross-correlation of their CMB bispectra. We are able to identify several distinct classes of primordial shapes - including equilateral, local, warm, flat and feature (non-scale invariant) - which should be distinguishable given a significant detection of CMB non-Gaussianity. We demonstrate that a simple shape correlator provides a fast and reliable method for determining whether or not CMB shapes are well correlated. We use an eigenmode decomposition of the primordial shape to characterise and understand model independence. Finally, we advocate a standardised normalisation method for $f_{NL}$ based on the shape autocorrelator, so that observational limits and errors can be consistently compared for different models.Comment: 32 pages, 20 figure

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

General CMB and Primordial Bispectrum Estimation I: Mode Expansion, Map-Making and Measures of f_NL

We present a detailed implementation of two bispectrum estimation methods which can be applied to general non-separable primordial and CMB bispectra. The method exploits bispectrum mode decompositions on the domain of allowed wavenumber or multipole values. Concrete mode examples constructed from symmetrised tetrahedral polynomials are given, demonstrating rapid convergence for known bispectra. We use these modes to generate simulated CMB maps of high resolution (l > 2000) given an arbitrary primordial power spectrum and bispectrum or an arbitrary late-time CMB angular power spectrum and bispectrum. By extracting coefficients for the same separable basis functions from an observational map, we are able to present an efficient and general f_NL estimator for a given theoretical model. The estimator has two versions comparing theoretical and observed coefficients at either primordial or late times, thus encompassing a wider range of models, including secondary anisotropies, lensing and cosmic strings. We provide examples and validation of both f_NL estimation methods by direct comparison with simulations in a WMAP-realistic context. In addition, we show how the full bispectrum can be extracted from observational maps using these mode expansions, irrespective of the theoretical model under study. We also propose a universal definition of the bispectrum parameter F_NL for more consistent comparison between theoretical models. We obtain WMAP5 estimates of f_NL for the equilateral model from both our primordial and late-time estimators which are consistent with each other, as well as with results already published in the literature. These general bispectrum estimation methods should prove useful for the analysis of nonGaussianity in the Planck satellite data, as well as in other contexts.Comment: 41 pages, 17 figure

Primordial non-Gaussianity and Bispectrum Measurements in the Cosmic Microwave Background and Large-Scale Structure

The most direct probe of non-Gaussian initial conditions has come from bispectrum measurements of temperature fluctuations in the Cosmic Microwave Background and of the matter and galaxy distribution at large scales. Such bispectrum estimators are expected to continue to provide the best constraints on the non-Gaussian parameters in future observations. We review and compare the theoretical and observational problems, current results and future prospects for the detection of a non-vanishing primordial component in the bispectrum of the Cosmic Microwave Background and large-scale structure, and the relation to specific predictions from different inflationary models.Comment: 82 pages, 23 figures; Invited Review for the special issue "Testing the Gaussianity and Statistical Isotropy of the Universe" for Advances in Astronom

General CMB and Primordial Trispectrum Estimation

We present trispectrum estimation methods which can be applied to general non-separable primordial and CMB trispectra. We present a general optimal estimator for the connected part of the trispectrum, for which we derive a quadratic term to incorporate the effects of inhomogeneous noise and masking. We describe a general algorithm for creating simulated maps with given arbitrary (and independent) power spectra, bispectra and trispectra. We propose a universal definition of the trispectrum parameter $T_{NL}$, so that the integrated bispectrum on the observational domain can be consistently compared between theoretical models. We define a shape function for the primordial trispectrum, together with a shape correlator and a useful parametrisation for visualizing the trispectrum. We derive separable analytic CMB solutions in the large-angle limit for constant and local models. We present separable mode decompositions which can be used to describe any primordial or CMB bispectra on their respective wavenumber or multipole domains. By extracting coefficients of these separable basis functions from an observational map, we are able to present an efficient estimator for any given theoretical model with a nonseparable trispectrum. The estimator has two manifestations, comparing the theoretical and observed coefficients at either primordial or late times. These mode decomposition methods are numerically tractable with order $l^5$ operations for the CMB estimator and approximately order $l^6$ for the general primordial estimator (reducing to order $l^3$ in both cases for a special class of models). We also demonstrate how the trispectrum can be reconstructed from observational maps using these methods.Comment: 38 pages, 9 figures. In v2 Figures 4-7 are altered slightly and some extra references are included in the bibliography. v3 matches version submitted to journal. Includes discussion of special case

Rapid Separable Analysis of Higher Order Correlators in Large Scale Structure

We present an efficient separable approach to the estimation and reconstruction of the bispectrum and the trispectrum from observational (or simulated) large scale structure data. This is developed from general CMB (poly-)spectra methods which exploit the fact that the bispectrum and trispectrum in the literature can be represented by a separable mode expansion which converges rapidly (with $n_\textrm{max}={\cal{O}}(30)$ terms). With an effective grid resolution $l_\textrm{max}$ (number of particles/grid points $N=l_\textrm{max}^3$), we present a bispectrum estimator which requires only ${\cal O}(n_\textrm{max} \times l_\textrm{max}^3)$ operations, along with a corresponding method for direct bispectrum reconstruction. This method is extended to the trispectrum revealing an estimator which requires only ${\cal O}(n_\textrm{max}^{4/3} \times l_\textrm{max}^3)$ operations. The complexity in calculating the trispectrum in this method is now involved in the original decomposition and orthogonalisation process which need only be performed once for each model. However, for non-diagonal trispectra these processes present little extra difficulty and may be performed in ${\cal O}(l_\textrm{max}^4)$ operations. A discussion of how the methodology may be applied to the quadspectrum is also given. An efficient algorithm for the generation of arbitrary nonGaussian initial conditions for use in N-body codes using this separable approach is described. This prescription allows for the production of nonGaussian initial conditions for arbitrary bispectra and trispectra. A brief outline of the key issues involved in parameter estimation, particularly in the non-linear regime, is also given