Primordial black holes (PBHs) are an important tool in cosmology to probe the primordial spectrum of small-scale curvature perturbations that reenter the cosmological horizon during radiation domination epoch. We numerically solve the evolution of sp
herically symmetric highly perturbed configurations to clarify the criteria of PBHs formation using a wide class of curvature profiles characterized by five parameters. It is shown that formation or non-formation of PBHs is determined essentialy by only two master parameters.
Primordial black holes (PBHs) are an important tool in cosmology to probe the primordial spectrum of small-scale curvature perturbations that reenter the cosmological horizon during radiation domination epoch. We numerically solve the evolution of sp
herically symmetric highly perturbed configurations to clarify the criteria of PBHs formation using an extremely wide class of curvature profiles characterized by five parameters, (in contrast to only two parameters used in all previous papers) which specify the curvature profiles not only at the central region but also at the outer boundary of configurations. It is shown that formation or non-formation of PBHs is determined entirely by only two master parameters one of which can be presented as an integral of curvature over initial configurations and the other is presented in terms of the position of the boundary and the edge of the core.
For an arbitrary strong, spherically symmetric super-horizon curvature perturbation, we present analytical solutions of the Einstein equations in terms of asymptotic expansion over the ratio of the Hubble radius to the length-scale of the curvature p
erturbation under consideration. To obtain this solution we develop a recursive method of quasi-linearization which reduces the problem to a system of coupled ordinary differential equations for the $N$-th order terms in the asymptotic expansion with sources consisting of a non-linear combination of the lower order terms. We use this solution for setting initial conditions for subsequent numerical computations. For an arbitrary precision requirement predetermined by the intended accuracy and stability of the computer code, our analytical solution yields optimal truncated asymptotic expansion which can be used to find the upper limit on the moment of time when the initial conditions expressed in terms of such truncated expansion should be set. Examples of how these truncated (up to eighth order) solutions provide initial conditions with given accuracy for different radial profiles of curvature perturbations are presented.
