Primordial Black Holes: Early Universe Probes

• Primordial black holes (PBHs) are black holes formed in the early universe via density perturbations and topological defects, serving as probes for cosmology and high-energy physics.
• Their mass spectrum, defined by horizon mass at formation, can be sharply peaked or extended, with constraints imposed by Hawking radiation, gravitational lensing, and cosmic background observations.
• PBHs influence gravitational wave signals and structure formation while offering a unique window into quantum gravity and early-universe phase transitions.

Primordial black holes (PBHs) are black holes formed in the early Universe through mechanisms independent of stellar evolution. As non-baryonic compact objects with a range of possible masses, PBHs provide a unique probe of early-universe cosmology, inflationary dynamics, high-energy particle physics, and are leading non-particle candidates for cold dark matter. They are subject to a wide array of cosmological, astrophysical, and particle-physics constraints owing to their gravitational, evaporative, and lensing signatures across cosmic time.

1. Formation Mechanisms

1.1. Collapse of Large Density Perturbations

During the radiation-dominated epoch, a comoving region of size $R$ can form a PBH upon horizon entry (at $t \sim R/c$) if its smoothed density contrast

$\delta = \frac{\delta\rho}{\rho}$

exceeds a critical value, $\delta_c \sim 0.3$–$0.5$. The mass of such a PBH is typically set by the horizon mass at that time: $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$ (Green, 2014). The Press–Schechter formalism with Gaussian statistics gives the initial PBH formation fraction as

$$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$

with $\sigma^2(M)$ the variance of $\delta$ on mass scale $M$. Critical-collapse numerics show the PBH mass–overdensity relation is

$t \sim R/c$0

where for radiation, $t \sim R/c$1 and $t \sim R/c$2.

1.2. Topological Defect Collapse: Cosmic Strings and Bubble Collisions

• Cosmic-String Loops: Oscillating loops with tension $t \sim R/c$3 can collapse if they contract within their Schwarzschild radius, forming PBHs of mass $t \sim R/c$4 with a broad spectrum $t \sim R/c$5. The collapse probability is $t \sim R/c$6, and $t \sim R/c$7 is required by evaporation constraints (Green, 2014).
• Cusp Collapse on Strings: Cusps are generically present on cosmic string loops. A segment near a cusp collapses under the hoop conjecture if $t \sim R/c$8, with rest-mass $t \sim R/c$9 and universal dimensionless spin $\delta = \frac{\delta\rho}{\rho}$0. Constraints from PBH evaporation give nearly model-independent upper limits, $\delta = \frac{\delta\rho}{\rho}$1, and the PBH DM fraction is suppressed below $\delta = \frac{\delta\rho}{\rho}$2 (Jenkins et al., 2020).
• Bubble Collisions in First-Order Phase Transitions: Collisions of true-vacuum bubbles can trap regions of the false vacuum; rare multibubble collisions may lead to collapse for sufficiently slow phase transitions, producing PBHs with mass at the horizon scale at percolation (Green, 2014).

1.3. Scalar Field and Topological Soliton Mechanisms

A scalar condensate (e.g., along SUSY flat directions) can fragment into Q-balls; Poisson fluctuations in the Q-ball distribution can seed overdense regions that collapse to PBHs. For SUSY, $\delta = \frac{\delta\rho}{\rho}$3 g; a general charged scalar can yield black holes up to $\delta = \frac{\delta\rho}{\rho}$4—relevant for LIGO (Cotner et al., 2017).

Topological defect networks (defects with $\delta = \frac{\delta\rho}{\rho}$5, $\delta = \frac{\delta\rho}{\rho}$6) can also generate large fluctuations and PBHs through Poisson statistics.

2. Mass Spectrum and Abundance

The PBH mass function, $\delta = \frac{\delta\rho}{\rho}$7, for nearly instantaneous formation is sharply peaked,

$\delta = \frac{\delta\rho}{\rho}$8

(Green, 2014), or for critical collapse,

$\delta = \frac{\delta\rho}{\rho}$9

where $\delta_c \sim 0.3$0 is the PDF of the density field. Non-instantaneous or broad mechanisms (e.g., scalar fragmentation, clustered formation) can yield extended or log-normal mass functions appropriate for interpreting observational constraints (Dolgov, 2017).

PBH abundance constraints are parametrized by $\delta_c \sim 0.3$1, with

$\delta_c \sim 0.3$2

(Green, 2014).

3. Evolution and Evaporation Physics

PBHs lose mass via Hawking radiation at the rate

$\delta_c \sim 0.3$3

with evaporation time

$\delta_c \sim 0.3$4

(Green, 2014, MacGibbon et al., 2015). PBHs with $\delta_c \sim 0.3$5 g have fully evaporated by the current epoch.

Accretion in the expanding universe is subdominant for cosmological mean densities: only PBHs with $\delta_c \sim 0.3$6–$\delta_c \sim 0.3$7 g grow significantly, reaching the Eddington limit and potentially growing to SMBH scales only in special, overdense environments (Rice et al., 2017). Hence, PBHs in the relevant cosmological mass range are generally relics with unchanged mass after formation, aside from Hawking evaporation for $\delta_c \sim 0.3$8 g.

4. Observational Constraints

4.1. Hawking Radiation, $\delta_c \sim 0.3$9-Rays, and BBN

• Diffuse $0.5$0-Ray Background: For $0.5$1, non-detection of extragalactic $0.5$2-rays constrains $0.5$3 (Green, 2014).
• BBN Constraints: Evaporation products during $0.5$4 affect light-element abundances; $0.5$5 (Green, 2014).
• CMB Spectral Distortions and Cosmic Ray Antiprotons: Place comparable or slightly weaker bounds.

4.2. Gravitational Lensing and Dynamical Effects

• Femtolensing: $0.5$6: $0.5$7.
• Microlensing: $0.5$8: $0.5$9 (EROS, MACHO, Kepler).
• Wide Binaries, Halo Friction, Disk-Heating:
  • Disruption of wide binaries: $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$0, $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$1.
  • MW halo friction/disk heating provide further exclusion at larger masses.

Clustering can relax microlensing bounds if PBHs form in clusters of $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$2; the microlensing event rate scales as $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$3, and the CMB bounds from accretion are also weakened due to enhanced PBH velocities in dense clusters. This opens a viable window $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$4–$$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$5 for clustered PBH DM (García-Bellido et al., 2017).

4.3. Dark Matter Window

Synthesizing all constraints, only the mass window

$$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$6

($$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$7) remains where PBHs could still make up all of the dark matter (Green, 2014, Liang et al., 17 Jan 2025). For some broad mass functions and special formation scenarios, $$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$8–$$M_H(t) \simeq \frac{c^3 t}{G} \sim 10^{15}\Bigl(\frac{t}{10^{-23}\,\mathrm{s}}\Bigr)\,\mathrm{g}$$9 windows can survive, particularly if clustering is significant (García-Bellido et al., 2017, Dolgov, 2017).

5. PBHs as Probes of Cosmology and Early-Universe Physics

Formation of PBHs probes the primordial power spectrum at scales far smaller than the CMB, $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$0. For PBH formation to be significant,

$$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$1

at these scales, several orders of magnitude above the $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$2 measured by Planck at CMB scales (Green, 2014). Stringent upper limits on $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$3 translate directly into bounds on small-scale power. Inflation models that predict "spikes" or "blue" spectra often violate these bounds. Non-Gaussian statistics can enhance PBH formation, further tightening joint constraints with non-Gaussianity parameters such as $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$4.

PBHs from first-order phase transitions, magneto-hydrodynamic turbulence, Affleck–Dine condensate fragmentation, or string-defect collapse offer direct mapping from particle physics parameters (e.g., transition scale, string tension) onto PBH observables, establishing synergy between early-universe microphysics and late-time cosmological signatures (Dent et al., 6 Oct 2025, Liang et al., 17 Jan 2025, Cotner et al., 2017, Jenkins et al., 2020).

6. Role in Gravitational Wave and Structure Formation

PBHs in the stellar-mass range have garnered significant attention as progenitors of the binary BH merger events detected by LIGO/Virgo, capable of yielding high merger rates for $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$5–$$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$6, depending on their mass spectrum and spatial distribution (Sasaki et al., 2018). PBH binaries formed in the early universe, through spatial clustering or gravitational decoupling, yield characteristic eccentricity and spin distributions distinct from those of stellar-origin BHs.

High-mass PBHs ($$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$7) with substantial abundance can accelerate the collapse of minihalos and advance Population III star formation, but are constrained by observations of early galaxies and 21-cm signals (Koulen et al., 6 Jun 2025). Lower-mass PBHs can delay or suppress star formation through tidal disruption of gas clouds.

Stochastic gravitational wave backgrounds are induced by both the PBH merger history and second-order effects from large scalar perturbations at PBH-forming scales, providing a target for future GW detectors (Inomata et al., 2021).

7. Quantum Gravity, Particle Physics, and Future Constraints

PBHs with $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$8 g evaporate via Hawking radiation, emitting a spectrum of photons, leptons, and hadrons. Stringent observational limits from the extragalactic $$\beta(M) = \frac{1}{2} \mathrm{erfc}\left(\frac{\delta_c}{\sqrt{2}\,\sigma(M)}\right)$$9-ray background, cosmic-ray antiprotons, and BBN restrict the PBH abundance at low masses (Carr, 2014, MacGibbon et al., 2015). Detection of final-stage PBH bursts would furnish unique evidence for quantum black hole thermodynamics.

The possibility of Planck-mass relics or exotic evaporation signatures constrains scenarios with new light degrees of freedom or quantum gravity effects. Fisher-matrix analyses combining GW, CMB spectral distortions, lensing statistics, and 21-cm cosmology are expected to further delimit PBH parameter space across multiple orders of magnitude in mass.

Across all scenarios, PBHs remain a multi-probe for the small-scale cosmological perturbations, high-energy particle physics, and the microphysical properties of the early universe. Ongoing and forthcoming observational programs targeting gravitational waves, CMB distortions, microlensing, and transient $\sigma^2(M)$0-ray bursts will empirically test the PBH hypothesis in the coming decade.