Beyond isolated curvature peaks: collective collapse and multiple Primordial Black Hole formation
Abstract
Primordial black hole (PBH) calculations usually treat rare curvature peaks as isolated collapsing regions. Using fully nonlinear 3+1 numerical-relativity simulations in a radiation-dominated Universe, we demonstrate that neighbouring primordial curvature perturbations need not map one-to-one onto PBHs: they may ultimately disperse, collapse collectively into a single PBH, or undergo distinct local collapses and form more than one PBH. In the family of two-component profiles studied here, the latter outcome is a pair of PBHs, identified by the coexistence of two disconnected apparent horizons on at least one time slice. We introduce a nonspherical quasi-local compaction diagnostic based on the Hawking mass and referenced to a round flat-FLRW sphere of equal area. It retains the angular structure of the curvature field and reduces to the standard Misner-Sharp compaction function in spherical symmetry. The global maximum of this diagnostic, 𝒦form, provides an empirical indicator of whether at least one PBH forms, with a transition near 𝒦form,c≈0.56 and modest profile-dependent scatter. For forming bimodal profiles, we supplement the nonlinear compactness with a signed linear surface strength evaluated on the same probing spheres. The viability of the weaker local branch and the competition between the local and common-enclosing branches improve the empirical discrimination between single- and double-PBH outcomes. Our results therefore show that the maximum curvature amplitude alone does not determine the collapse outcome, which also depends on the spatial extent, characteristic scales, and geometry of the surrounding curvature environment.
Related
- Simulating dust collapse: the impact of gauge choices on primordial black hole formation
- Primordial Black Holes from Vector-Induced Curvature Perturbations Sourced by Primordial Magnetic Fields
- Gravitational-wave signatures of primordial black hole clusters and the imprint of an early dense-core collapse
- How do late-time photons "record" primordial curvature perturbations?
- Reviving primordial black hole formation in slow first-order phase transitions
BibTeX
@misc{escriva2026beyond,
title = {{Beyond isolated curvature peaks: collective collapse and multiple Primordial Black Hole formation}},
author = {Albert Escrivà},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.03051},
eprint = {2609.03051},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.03051},
}
RIS
TY - UNPB
TI - Beyond isolated curvature peaks: collective collapse and multiple Primordial Black Hole formation
AU - Albert Escrivà
PY - 2026
DA - 2026/09/02
PB - arXiv
DO - 10.48550/arxiv.2609.03051
UR - https://arxiv.org/abs/2609.03051
AB - Primordial black hole (PBH) calculations usually treat rare curvature peaks as isolated collapsing regions. Using fully nonlinear 3+1 numerical-relativity simulations in a radiation-dominated Universe, we demonstrate that neighbouring primordial curvature perturbations need not map one-to-one onto PBHs: they may ultimately disperse, collapse collectively into a single PBH, or undergo distinct local collapses and form more than one PBH. In the family of two-component profiles studied here, the latter outcome is a pair of PBHs, identified by the coexistence of two disconnected apparent horizons on at least one time slice. We introduce a nonspherical quasi-local compaction diagnostic based on the Hawking mass and referenced to a round flat-FLRW sphere of equal area. It retains the angular structure of the curvature field and reduces to the standard Misner-Sharp compaction function in spherical symmetry. The global maximum of this diagnostic, 𝒦_form, provides an empirical indicator of whether at least one PBH forms, with a transition near 𝒦_form,c≈0.56 and modest profile-dependent scatter. For forming bimodal profiles, we supplement the nonlinear compactness with a signed linear surface strength evaluated on the same probing spheres. The viability of the weaker local branch and the competition between the local and common-enclosing branches improve the empirical discrimination between single- and double-PBH outcomes. Our results therefore show that the maximum curvature amplitude alone does not determine the collapse outcome, which also depends on the spatial extent, characteristic scales, and geometry of the surrounding curvature environment.
ER -