  We investigate non-Gaussianities in self-interacting curvaton models treating
both renormalizable and non-renormalizable polynomial interactions. We scan the
parameter space systematically and compute numerically the non-linearity
parameters f_NL and g_NL. We find that even in the interaction dominated regime
there are large regions consistent with current observable bounds. Whenever the
interactions dominate, we discover significant deviations from the relations
f_NL ~ 1/r_decay and g_NL ~ 1/r_decay valid for quadratic curvaton potentials,
where r_decay measures the curvaton contribution to the total energy density at
the time of its decay. Even if r_decay << 1, there always exists regions with
f_NL ~ 0 since the sign of f_NL oscillates as a function of the parameters.
While g_NL can also change sign, typically g_NL is non-zero in the low-f_NL
regions. Hence, for some parameters the non-Gaussian statistics is dominated by
g_NL rather than by f_NL. Due to self-interactions, both the relative signs of
f_NL and g_NL and the functional relation between them is typically modified
from the quadratic case, offering a possible experimental test of the curvaton
interactions.

Cosmology and Nongalactic Astrophysics