Limits on non-Gaussianities from WMAP data

We develop a method for constraining the level of non-Gaussianity of density perturbations when the three-point function is of the 'equilateral' type. Departures from Gaussianity of this form are produced by single field models such as ghost or Dirac Born Infeld inflation and in general by the presence of higher order derivative operators in the effective Lagrangian of the inflaton. We show that the induced shape of the three-point function can be very well approximated by a factorizable form, making the analysis practical. We also show that, unless one has a full sky map with uniform noise, in order to saturate the Cramer Rao bound for the error on the amplitude of the three-point function, the estimator must contain a piece that is linear in the data. We apply our technique to the WMAP data obtaining a constraint on the amplitude f_NL^equil of 'equilateral' non-Gaussianity: -366 < f_NL^equil < 238 at 95% C.L. We also apply our technique to constrain the so-called 'local' shape, which is predicted for example by the curvaton and variable decay width models. We show that the inclusion of the linear piece in the estimator improves the constraint over those obtained by the WMAP team, to -27 < f_NL^local < 121 at 95% C.L.