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.
