Calculating the power spectrum in stochastic inflation by Monte Carlo simulation and least squares curve fitting

Abstract

The stochastic-δ N formalism is widely used to study inflation models in which the quantum diffusion of inflatons dominates the background dynamics, leading to interesting phenomena such as the production of primordial black holes. Among numerical approaches to calculate the curvature perturbation spectrum Pζ(k) in this formalism, the Monte Carlo simulation-based approach has been proposed as a promising choice, especially in multifield cases. In this approach, we generate many paths of inflatons from the initial points to the end of inflation, obtain statistics of δ N from the paths, and then estimate Pζ(k). However, this method involves a nested Monte Carlo simulation, which requires generating many branch paths from each trunk path at the point corresponding to the scale k of interest, resulting in a high computational cost. In this paper, we propose a new Monte Carlo-based approach that utilizes least squares fitting, introducing two novel features for reducing computational cost. First, we devise a simple estimator of a key statistic δNX2, the variance of δ N conditioned on the branching point, to avoid nesting path generation. Second, via least squares fitting of a parametric function to the sampled values of the estimator, we obtain not just an estimate of Pζ(k) for a single value of k but an approximating function of Pζ(k) over a range of k of interest. We also conduct numerical demonstrations for concrete inflation models, which show the usefulness of our method.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…