Curvaton preheating revisited

We study the thermalization process in the self-interacting curvaton preheating scenario. We solve the evolution of the system with classical lattice simulations with a recently released symplectic PyCOOL program during the resonance and the early thermalization periods and compare the results to the inflaton preheating. After this we calculate the generated non-gaussianity with the $\Delta N$ formalism and the separate universe approximation by running a large number of simulations with slightly different initial values. The results indicate a high level of non-gaussianity. We also use this paper to showcase the various post-processing functions included with the PyCOOL program that is available from https://github.com/jtksai/PyCOOL .