We explore the possibility that Dark Matter (DM) may be explained by a nonuniform background of approximately stellar mass clusters of Primordial Black Holes (PBHs) by simulating the evolution from recombination to the present with over 5000 realisations using a Newtonian \({ N }\)-body code. We compute the cluster rate of evaporation and extract the binary and merged sub-populations along with their parent and merger tree histories, lifetimes and formation rates, the dynamical and orbital parameter profiles, the degree of mass segregation and dynamical friction and power spectrum of close encounters. Overall, we find that PBHs can constitute a viable DM candidate, and that their clustering presents a rich phenomenology throughout the history of the Universe. We show that binary systems constitute about 9.5% of all PBHs at present, with mass ratios of \({ \bar{q}_{\rm B} = 0.154 }\), and total masses of \({ \bar{m}_{\rm T,\,B} = 303\,M_\odot}\). Merged PBHs are rare, about 0.0023% of all PBHs at present, with mass ratios of \({ \bar{q}_{\rm B}= 0.965 }\) with total and chirp masses of \({ \bar{m}_{\rm T,\,B}= 1670\,M_\odot}\) and \({ \bar{m}_{c,{\rm M}} = 642\,M_\odot }\), respectively. We find that cluster puffing up and evaporation leads to bubbles of these PBHs of order 1 kpc containing at present times about 36% of objects and mass, with one-hundred pc-sized cores. We also find that these PBH sub-haloes are distributed in wider PBH haloes of order hundreds of kpc, containing about 63% of objects and mass, coinciding with the sizes of galactic halos. We find at last high rates of close encounters of massive Black Holes (\({ M \sim 1000\,M_\odot}\)), with \({ \Gamma^{\mathrm{S}} = (1.2_{-0.9}^{+5.9}) \times 10^{7}~{\rm yr^{-1}~\rm Gpc^{-3}}}\) and mergers with \({\Gamma^{\mathrm{M}} = 1337 \pm 41~{\rm yr^{-1}~\rm Gpc^{-3}} }\).

中文翻译：

---

N体模拟中原始黑洞的聚类动力学

我们通过使用牛顿\（{N}来模拟从重组到现在的进化过程，并利用5000多个实现，探索了暗物质（DM）可能由原始恒星黑洞（PBH）近似恒星质量簇的非均匀背景所解释的可能性。 \）-主体代码。我们计算蒸发的群集速率，并提取二元和合并的子种群及其父级和合并树的历史记录，寿命和形成率，动力学和轨道参数概况，质量偏析程度以及动摩擦和闭合功率谱遇到。总体而言，我们发现PBH可以构成可行的DM候选者，并且它们的聚集在整个宇宙历史上都呈现出丰富的现象学。我们表明，二进制系统目前占所有PBH的9.5％，质量比为\（{\ bar {q} _ {\ rm B} = 0.154} \），总质量为\（{\ bar {m} _ {\ rm T，\，B} = 303 \， M_ \ odot} \）。合并的PBH很少见，目前约占所有PBH的0.0023％，质量比为\（{\ bar {q} _ {\ rm B} = 0.965} \），总质量和线性质量为\ {{\ bar {m } _ {\ rm T，\，B} = 1670 \，M_ \ odot} \）和\（{\ bar {m} _ {c，{\ rm M}} = 642 \，M_ \ odot} \） ， 分别。我们发现，簇状膨胀和蒸发导致这些PBH的气泡出现，数量为1 kpc，目前包含约36％的物体和质量，具有100个pc尺寸的芯。我们还发现，这些PBH子光环分布在数百kpc数量级的更宽的PBH光环中，包含约63％的物体和质量，与银河光环的大小相吻合。我们发现大量黑洞（\（{M \ sim 1000 \，M_ \ odot} \）的近距离接触率很高，