stops. The linear analysis [@Kawasaki:2006zv] shows that the enhancement of the scalar fields occurs at the spacial scale $R_{\mathrm{peak}}\sim1/\left(0.35m_{\sigma}\right)$. Taking this into account, we set the box size of the simulation $V=L^3$ and the grid number $N_{\mathrm{grid}}$ as $L=100/m_{\sigma}$ and $N_{\mathrm{grid}}=256^3$. Then, $L>R_{\mathrm{peak}}$ and $L/N_{\mathrm{grid}}^{1/3}<R_{\mathrm{peak}}$ are satisfied even if we consider the growth of the scale factor. Thus, we can resolve the peak scale $R_{\mathrm{peak}}$ during the simulation. We choose the physical parameters as follows $$\label{eq:3:parameter}
{\mu} = { 1.38\times10^{-3}M_{{\rm pl}}},~
{ v} = { 3.45\times10^{-4}M_{{\rm pl}}},~
{ M} = { 0.545M_{{\rm pl}}},~
{ g} = { 2.0\times10^{-5}},~\mathrm{and}~
{ C_{N} = 0.04},$$ with which the spectral index $n_s-1=-0.035$ is within 1 $\sigma$ deviation from the central value of the Planck result [@Ade:2013zuv]. We set the initial amplitude of the fluctuations of the scalar fields by Rayleigh distribution [@Polarski:1995jg; @Khlebnikov:1996mc; @Khlebnikov:1996zt; @GarciaBellido:2002aj].

[c c]{}\

The inflation model we consider here has three scalar fields $\sigma,~\psi$ and $\phi$, where $\sigma$ is the inflaton of the hybrid inflation, $\psi$ is the waterfall field of the hybrid inflation and $\phi$ is the inflaton of the new inflation. We perform the simulation using these three real scalar fields.[^2] Using the lattice simulations, we calculate the 