We study the primordial bispectrum of curvature perturbation in the uniform- density slicing generated by the interaction between the inflaton and isotropic background gauge fields. We derive the action up to cubic order in perturbation and take into account all the relevant effects in the leading order of slow-roll expansion. We first treat the quadratic vertices perturbatively and confirm the results of past studies, while identifying their regime of validity. We then extend the analysis to include the effect of the quadratic vertices to all orders by introducing exact linear mode functions, allowing us to make accurate predictions long after horizon crossing where the features of both the power spectrum and the bispectrum are drastically different. It is shown that the spPPectra become constant and scale-invariant in the limit of large e-folding. As a result, we are able to impose reliable constraints on the parameters of our theory using the recent observational data coming from Planck.

We propose a new method to numerically calculate higher-order correlation functions of primordial fluctuations generated from any early-universe scenario. Our key-starting point is the realization that the tree-level In-In formalism is intrinsically separable. This enables us to use modal techniques to efficiently calculate and represent non-Gaussian shapes in a separable form well suited to data analysis. We prove the feasibility and the accuracy of our method by applying it to simple single-field inflationary models in which analytical results are available, and we perform non-trivial consistency checks like the verification of the single field consistency relation. We also point out that the i epsilon prescription is automatically taken into account in our method, preventing the need for ad-hoc tricks to implement it numerically.