Mixed perturbations from an Inflaton and Two Curvatons - towards   understanding non-gaussianity in more complicated curvaton scenarios
  We analyse the situation where the primordial curvature perturbations are produced by the joint effects of an inflaton field and two curvaton fields. We present general equations which allow the reader to obtain f_NL for several different scenarios which differ in the order in which fields decay into radiation after inflation. In order to investigate the physics of these equations we analyse some simplified situations where the fields are harmonic and both curvatons are frozen at the same expectation value during inflation. We find quite complex behaviour - for a given situation where the inflaton contributes a fixed amount to the total curvature perturbation there are situations where f_NL is maximised if both curvatons share equally in contributing the rest and situations where f_NL is maximised if only one of the two curvatons contributes the rest. We are unable therefore to make any completely general extrapolations about the expected non-gaussianity from N curvatons. We find that as the curvaton contribution to the overall perturbation is gradually increased, f_NL rises to a maximum before falling again and that a given f_NL can correspond to many different parameter sets for the two curvatons. 
Maria Sueiro and Malcolm Fairbairn
astro-ph.CO