A PBH with an initial mass larger than \(\sim 10^{15}\)  g (i.e., heavy PBH) can survive in today's universe [1]}. In this case, the relic density of PBHs is obtained as [2]}, [1]}
\(\Omega _{\rm PBH} &=& 1.42\times 10^{17} \left(\frac{\gamma }{0.2}\right)^{1/2}\left(\frac{0.67}{h}\right)^{2} \left(\frac{106.75}{g_\ast }\right)^{1/4}\nonumber \\&& \times \left(\frac{10^{15} {\rm g}}{M_{\rm PBH}}\right)^{1/2} \beta ,\) 
