Constraining the WMAP9 bispectrum and trispectrum with needlets

We develop a needlet approach to estimate the amplitude of general (including non-separable) bispectra and trispectra in the cosmic microwave background, and apply this to the WMAP 9-year data. We obtain estimates for the `orthogonal' bispectrum mode, yielding results which are consistent with the WMAP 7-year data. We do not observe the frequency-dependence suggested by the WMAP team's analysis of the 9-year data. We present 1-$\sigma$ constraints on the `local' trispectrum shape $\gnl/10^5= -4.1\pm 2.3$, the `$c1$' equilateral model $\gnl^{c_1}/10^6= -0.8\pm 2.9$, and the constant model $\gnl^{\rm{const}}/10^6= -0.2\pm 1.8$, together with a $95\%$ confidence-level upper bound on the multifield local parameter $\taunl<22000$. We estimate the bias on these parameters produced by point sources. The techniques developed in this paper should prove useful for other datasets such as Planck.