In this Letter we show that there is a unique nonminimal derivative coupling of the standard model Higgs boson to gravity such that it propagates no more degrees of freedom than general relativity sourced by a scalar field, reproduces a successful inflating background within the standard model Higgs parameters, and finally does not suffer from dangerous quantum corrections.

In this Letter, combining peak theory and the numerical analysis of gravitational collapse in the radiation dominated era, we show that the abundance of primordial blacks holes, generated by an enhancement in the inflationary power spectrum, is extremely dependent on the shape of the peak. Given the amplitude of the power spectrum, we show that the density of primordial black holes generated from a narrow peak is exponentially smaller than in the case of a broad peak. Specifically, for a top-hat profile of the power spectrum in Fourier space, we find that to have primordial black holes comprising all of the dark matter, one would only need a power spectrum amplitude an order of magnitude smaller than suggested previously, whereas in the case of a narrow peak, one would instead need a much larger power spectrum amplitude, which in many cases would invalidate the perturbative analysis of cosmological perturbations. Finally, we show that, although critical collapse gives a broad mass spectrum, the density of primordial black holes formed is dominated by masses roughly equal to the cosmological horizon mass measured at horizon crossing.

In this paper, we argue and show numerically that the threshold to form primordial black holes from an initial spherically symmetric perturbation is, to an excellent approximation, universal, whenever given in terms of the compaction function averaged over a sphere of radius $r_m$, where $r_m$ is the scale on which the compaction function is maximum. This can be understood as the requirement that, for a black hole to form, each shell of the averaged compaction function should have an amplitude exceeding the so-called Harada-Yoo-Kohri limit. For a radiation dominated universe we argued, supported by the numerical simulations, that this limit is ${\delta }_c=0.40$, which is slightly below the one quoted in the literature. Additionally, we show that the profile dependence of the threshold for the compaction function is only sensitive to its curvature at the maximum. We use these results to provide an analytic formula for the threshold amplitude of the compaction function at its maximum in terms of the normalized compaction function curvature at $r_m$.