Self-consistent initial conditions for primordial black hole formation

For an arbitrarily strong, spherically symmetric super-horizon curvature perturbation, we present analytic solutions of the Einstein equations in terms of the asymptotic expansion over the ratio of the Hubble radius to the length-scale of the curvature perturbation to set initial conditions for numerical computations of primordial black hole formation. 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 Nth order terms in the asymptotic expansion with sources consisting of a non-linear combination of the lower order terms.