Calculating the squeezed bispectrum in stochastic inflation by Monte Carlo simulation

Abstract

The bispectrum of cosmological perturbations is an important probe of inflation and the underlying particle physics. However, calculating it becomes challenging in inflation models where inflaton dynamics is dominated by quantum diffusion rather than slow-roll, especially in the multi-field case. In this paper, employing the stochastic-$δ\mathcal{N}$ formalism, we propose a Monte Carlo-based method to calculate the bispectrum in the squeezed limit, as an extension of the method for the power spectrum that we previously proposed. Our method involves generating paths of inflatons' time evolution that branch only several times, avoiding nested path generation, which incurs a prohibitive computational cost. As numerical demonstrations, we apply the proposed method to single- and double-field chaotic inflation for validation, and to hybrid inflation with mild waterfall, for which, to the best of our knowledge, the bispectrum is calculated for the first time in this work.