Non-Gaussianity after many-field reheating

We numerically investigate reheating after quadratic inflation with up to 65 fields, focusing on the production of non-Gaussianity. We consider several sets of initial conditions, masses, and decay rates. As expected, we find that the reheating phase can have a significant effect on the non-Gaussian signal, but that for this number of fields a detectable level of non-Gaussianity requires the initial conditions, mass range, and decay rates to be ordered in a particular way. We speculate on whether this might change in the N-flation limit.