Clustering of primordial black holes from quantum diffusion during inflation

Abstract

We study how large fluctuations are spatially correlated in the presence of quantum diffusion during inflation. This is done by computing real-space correlation functions in the stochastic-δ N formalism. We first derive an exact description of physical distances as measured by a local observer at the end of inflation, improving on previous works. Our approach is based on recursive algorithmic methods that consistently include volume-weighting effects. We then propose a “large-volume” approximation under which calculations can be done using first-passage time analysis only, and from which a new formula for the power spectrum in stochastic inflation is derived. We then study the full two-point statistics of the curvature perturbation. Due to the presence of exponential tails, we find that the joint distribution of large fluctuations is of the form P(ζR1, ζR2) = F(R1,R2, r) P(ζR1)P( ζR2), where ζR1 and ζR2 denote the curvature perturbation coarse-grained at radii R1 and R2, around two spatial points distant by r. This implies that, on the tail, the reduced correlation function, defined as P(ζR1 > ζc, ζR2 > ζc)/[P(ζR1 > ζc) P(ζR2 > ζc)]-1, is independent of the threshold value ζc. This contrasts with Gaussian statistics where the same quantity strongly decays with ζc, and shows the existence of a universal clustering profile for all structures forming in the exponential tails. Structures forming in the intermediate (i.e. not yet exponential) tails may feature different, model-dependent behaviours.