Numerical simulations of primordial black hole formation via delayed first-order phase transitions

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Ning, Zhuan, Zeng, Xiang-Xi, Cai, Rong-Gen, Wang, Shao-Jiang
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866911719184400384
author Ning, Zhuan
Zeng, Xiang-Xi
Cai, Rong-Gen
Wang, Shao-Jiang
author_facet Ning, Zhuan
Zeng, Xiang-Xi
Cai, Rong-Gen
Wang, Shao-Jiang
contents We perform fully nonlinear, spherically symmetric numerical simulations of superhorizon false-vacuum-domain (FVD) collapse in a coupled gravity-scalar-fluid system to study primordial black hole (PBH) formation during delayed first-order phase transitions (FOPTs). Using adaptive mesh refinement to resolve the bubble wall, we identify three dynamical outcomes: type B (supercritical) PBHs with an interior baby universe and a bifurcating trapping horizon, type A (subcritical) PBHs with an apparent horizon formed by direct wall collapse, and dispersal with no PBH formation. To separate these three cases, we evaluate two commonly used PBH-formation criteria: the time scale ratio $t_\mathrm{H}/t_\mathrm{V}$ (horizon crossing time versus vacuum-energy domination time) and the local density contrast $δ(t_\mathrm{H})$ at horizon crossing. For the parameter space explored, we find that $t_\mathrm{H}/t_\mathrm{V}$ is a more robust predictor of outcome: type B PBHs form when $t_\mathrm{H}/t_\mathrm{V} \gtrsim 1$ (critical range $\sim 1.1 - 1.6$ in our survey), type A PBHs arise when $t_\mathrm{H}/t_\mathrm{V}$ is below this threshold but remains above a lower bound (typical range $\sim 0.35 - 0.7$), and no-PBH dispersal occurs when $t_\mathrm{H}/t_\mathrm{V}$ falls below this lower bound. When a clear thin-wall FVD boundary exists, $δ(t_\mathrm{H})$ can correspondingly distinguish different outcomes (roughly $δ_c(t_\mathrm{H}) \sim 1 - 1.7$ for type B and $δ_c(t_\mathrm{H}) \sim 0.35 - 0.5$ for type A), but is highly sensitive to wall structure and model details and thus less universal. These results offer new insights into the dynamics of FVD collapse, quantify practical PBH-formation thresholds, and pave the way for precise predictions of PBH abundance from delayed FOPTs.
format Preprint
id arxiv_https___arxiv_org_abs_2601_21878
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Numerical simulations of primordial black hole formation via delayed first-order phase transitions
Ning, Zhuan
Zeng, Xiang-Xi
Cai, Rong-Gen
Wang, Shao-Jiang
General Relativity and Quantum Cosmology
Cosmology and Nongalactic Astrophysics
High Energy Physics - Phenomenology
We perform fully nonlinear, spherically symmetric numerical simulations of superhorizon false-vacuum-domain (FVD) collapse in a coupled gravity-scalar-fluid system to study primordial black hole (PBH) formation during delayed first-order phase transitions (FOPTs). Using adaptive mesh refinement to resolve the bubble wall, we identify three dynamical outcomes: type B (supercritical) PBHs with an interior baby universe and a bifurcating trapping horizon, type A (subcritical) PBHs with an apparent horizon formed by direct wall collapse, and dispersal with no PBH formation. To separate these three cases, we evaluate two commonly used PBH-formation criteria: the time scale ratio $t_\mathrm{H}/t_\mathrm{V}$ (horizon crossing time versus vacuum-energy domination time) and the local density contrast $δ(t_\mathrm{H})$ at horizon crossing. For the parameter space explored, we find that $t_\mathrm{H}/t_\mathrm{V}$ is a more robust predictor of outcome: type B PBHs form when $t_\mathrm{H}/t_\mathrm{V} \gtrsim 1$ (critical range $\sim 1.1 - 1.6$ in our survey), type A PBHs arise when $t_\mathrm{H}/t_\mathrm{V}$ is below this threshold but remains above a lower bound (typical range $\sim 0.35 - 0.7$), and no-PBH dispersal occurs when $t_\mathrm{H}/t_\mathrm{V}$ falls below this lower bound. When a clear thin-wall FVD boundary exists, $δ(t_\mathrm{H})$ can correspondingly distinguish different outcomes (roughly $δ_c(t_\mathrm{H}) \sim 1 - 1.7$ for type B and $δ_c(t_\mathrm{H}) \sim 0.35 - 0.5$ for type A), but is highly sensitive to wall structure and model details and thus less universal. These results offer new insights into the dynamics of FVD collapse, quantify practical PBH-formation thresholds, and pave the way for precise predictions of PBH abundance from delayed FOPTs.
title Numerical simulations of primordial black hole formation via delayed first-order phase transitions
topic General Relativity and Quantum Cosmology
Cosmology and Nongalactic Astrophysics
High Energy Physics - Phenomenology
url https://arxiv.org/abs/2601.21878