Primordial black holes with an accurate QCD equation of state

Making use of definitive new lattice computations of the Standard Model thermodynamics during the quantum chromodynamic (QCD) phase transition, we calculate the enhancement in the mass distribution of primordial black holes (PBHs) due to the softening of the equation of state. We find that the enhancement peaks at approximately $0.7M_\odot$, with the formation rate increasing by at least two orders of magnitude due to the softening of the equation of state at this time, with a range of approximately $0.3M_\odot<M<1.4M_\odot$ at full width half-maximum. PBH formation is increased by a smaller amount for PBHs with masses spanning a large range, $10^{-3}M_\odot<M_{\rm PBH}<10^{3}M_\odot$, which includes the masses of the BHs that LIGO detected. The most significant source of uncertainty in the number of PBHs formed is now due to unknowns in the formation process, rather than from the phase transition. A near scale-invariant density power spectrum tuned to generate a population with mass and merger rate consistent with that detected by LIGO should also produce a much larger energy density of PBHs with solar mass. The existence of BHs below the Chandresekhar mass limit would be a smoking gun for a primordial origin and they could arguably constitute a significant fraction of the cold dark matter density. They also pose a challenge to inflationary model building which seek to produce the LIGO BHs without overproducing lighter PBHs.