We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $dotsigma^
4$, $dotsigma^2 (partialsigma^2)$, and $(partialsigma)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{dotsigma^4}$ and $g_{NL}^{(partialsigma)^4}$, in addition to the coefficient $g_{NL}^{rm loc}$ of $zeta^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{rm loc}, f_{NL}^{rm equil}, f_{NL}^{rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{rm loc} = (-3.80 pm 2.19) times 10^5$, $g_{NL}^{dotsigma^4} = (-3.20 pm 3.09) times 10^6$, and $g_{NL}^{(partialsigma)^4} = (-10.8 pm 6.33) times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.
The recent BICEP2 measurement of primordial gravity waves (r = 0.2^{+0.07}_{-0.05}) appears to be in tension with the upper limit from WMAP (r<0.13 at 95% CL) and Planck (r<0.11 at 95% CL). We carefully quantify the level of tension and show that it
is very significant (around 0.1% unlikely) when the observed deficit of large-scale temperature power is taken into account. We show that measurements of TE and EE power spectra in the near future will discriminate between the hypotheses that this tension is either a statistical fluke, or a sign of new physics. We also discuss extensions of the standard cosmological model that relieve the tension, and some novel ways to constrain them.
Cosmic microwave background observations are most commonly analyzed by estimating the power spectrum. In the limit where the CMB statistics are perfectly Gaussian, this extracts all the information, but the CMB also contains detectable non-Gaussian c
ontributions from secondary, and possibly primordial, sources. We review possible sources of CMB non-Gaussianity and describe statistical techniques which are optimized for measuring them, complementing the power spectrum analysis. The machinery of $N$-point correlation functions provides a unifying framework for optimal estimation of primordial non-Gaussian signals or gravitational lensing. We review recent results from applying these estimators to data from the WMAP satellite mission.
Models of inflation in which non-Gaussianity is generated outside the horizon, such as curvaton models, generate distinctive higher-order correlation functions in the CMB and other cosmological observables. Testing for violation of the Suyama-Yamaguc
hi inequality tauNL >= (6/5 fNL)^2, where fNL and tauNL denote the amplitude of the three-point and four-point functions in certain limits, has been proposed as a way to distinguish qualitative classes of models. This inequality has been proved for a wide range of models, but only weak
A wide range of multifield inflationary models generate non-Gaussian initial conditions in which the initial adiabatic fluctuation is of the form (zeta_G + g_{NL} zeta_G^3). We study halo clustering in these models using two different analytic method
s: the peak-background split framework, and brute force calculation in a barrier crossing model, obtaining agreement between the two. We find a simple, theoretically motivated expression for halo bias which agrees with N-body simulations and can be used to constrain g_{NL} from observations. We discuss practical caveats to constraining g_{NL} using only observable properties of a tracer population, and argue that constraints obtained from populations whose observed bias is <~ 2.5 are generally not robust to uncertainties in modeling the halo occupation distribution of the population.
Large-scale clustering of highly biased tracers of large-scale structure has emerged as one of the best observational probes of primordial non-Gaussianity of the local type (i.e. f_{NL}^{local}). This type of non-Gaussianity can be generated in multi
field models of inflation such as the curvaton model. Recently, Tseliakhovich, Hirata, and Slosar showed that the clustering statistics depend qualitatively on the ratio of inflaton to curvaton power xi after reheating, a free parameter of the model. If xi is significantly different from zero, so that the inflaton makes a non-negligible contribution to the primordial adiabatic curvature, then the peak-background split ansatz predicts that the halo bias will be stochastic on large scales. In this paper, we test this prediction in N-body simulations. We find that large-scale stochasticity is generated, in qualitative agreement with the prediction, but that the level of stochasticity is overpredicted by ~30%. Other predictions, such as xi independence of the halo bias, are confirmed by the simulations. Surprisingly, even in the Gaussian case we do not find that halo model predictions for stochasticity agree consistently with simulations, suggesting that semi-analytic modeling of stochasticity is generally more difficult than modeling halo bias.
