Simulation of primordial black holes with large negative non-Gaussianity

Abstract

AbstractIn this work, we have performed numerical simulations of primordial black hole (PBH) formation in the Friedman-Lemaître-Robertson-Walker universe filled by radiation fluid, introducing the local-type non-Gaussianity to the primordial curvature fluctuation. We have compared the numerical results from simulations with previous analytical estimations on the threshold value for PBH formation done in the previous paper [1], particularly for negative values of the non-linearity parameterfNL. Our numerical results show the existence of PBH formation of (the so-called) type I also in the casefNL≲ -0.336, which was not found in the previous analytical expectations using the critical averaged compaction function. In particular, although the universal value for the averaged critical compaction function 𝒞̅𝒸= 2/5 found previously in the literature is not satisfied for all the profiles considered in this work, an alternative direct analytical estimate has been found to be roughly accurate to estimate the thresholds, which gives the value of the critical averaged density with a few % deviation from the numerical one forfNL≳ -1.