Primordial Black Hole Dark Matter

• Primordial black holes are black holes formed from early Universe density fluctuations that act as non-baryonic dark matter candidates with unique mass distributions.
• Observational constraints from the CMB, microlensing, gravitational waves, and gamma-ray emissions narrow down the viable PBH mass windows for dark matter.
• Theoretical models predict PBH formation via inflationary mechanisms and phase transitions, with future surveys set to test and refine these dark matter scenarios.

Primordial black holes (PBHs) are black holes formed in the early Universe from the direct gravitational collapse of large density fluctuations, strongly out of equilibrium with the cosmic background. As a non-baryonic, collisionless, and dynamically cold relic, PBHs have been a compelling candidate for dark matter (DM), offering a purely gravitational explanation for the observed DM abundance without appealing to new elementary particles. Key advances in cosmic microwave background (CMB) observations, microlensing, gravitational-wave astronomy, and theoretical modeling have sharpened the mass windows, formation mechanisms, and observational signatures by which PBHs might constitute a significant or even dominant fraction of dark matter.

1. Origin, Formation Mechanisms, and Mass Function

PBHs can form when primordial curvature fluctuations re-enter the Hubble horizon during radiation domination and locally induce overdensities δ exceeding a threshold δ_c, set by the equation of state (typically δ_c ≈ 0.4–0.5 for a radiation fluid). The PBH mass at formation is typically of order the horizon mass,

$$M_{\rm PBH} \simeq \gamma\, M_H(t),\quad M_H(t) \sim 10^{15}\,{\rm g}\,(t/10^{-23}\,{\rm s}),$$

where γ ≈ 0.2–0.5 absorbs collapse details (Green, 2024, Carr, 2019, Carr et al., 2016).

Inflation models allow for sharply enhanced power spectra at specific small scales via mechanisms such as "ultra-slow-roll" plateaus, inflection points, spectator fields, phase transitions (e.g., QCD, electroweak), or tailored multi-field potentials (Correa et al., 2022, Frampton et al., 2010). Depending on the spectrum shape, PBHs may follow nearly monochromatic, power-law, or lognormal mass distributions:

$$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$

with $M_c$ the characteristic mass and $\sigma$ the width (Clesse et al., 2017). Nonlinear collapse, critical phenomena, and non-Gaussian perturbations induce further broadening and nontrivial low-mass tails (Carr et al., 2016, Green, 2024).

Alternate mechanisms include:

• Bubble collisions in first-order phase transitions
• Collapse of cosmic string loops or domain walls
• Early matter- or dust-dominated phases enabling rapid growth of subhorizon modes
• Accretion episodes, e.g., via magneto-rotational instability (MRI) in early positronium plasmas [(Curd et al., 2024) (abstract only)]

The initial fraction $\beta(M)$ of horizon patches collapsing into PBHs is exponentially sensitive to the local variance $\sigma(M)$:

$$\beta(M)\simeq \frac{\sigma(M)}{\delta_c\sqrt{2\pi}}\,\exp\left[-\frac{\delta_c^2}{2\sigma^2(M)}\right],$$

and relates to the present-day fraction $$f_{\rm PBH}(M) = \Omega_{\rm PBH}(M)/\Omega_{\rm DM}$$ through cosmological redshifting and the scaling difference between matter and radiation (Green, 2024, Carr et al., 2021).

2. Observational Constraints and Allowed Mass Windows

Astrophysical and cosmological datasets impose stringent constraints on PBH dark matter across many mass decades:

• Evaporation bounds: PBHs with $M < 5 \times 10^{14}$ g have lifetimes shorter than the Hubble time and are excluded as DM. Hawking radiation, producing photons, e$^\pm$, and neutrinos, severely limits $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$0 up to $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$1 g (Ray, 2022, Green, 2024).
• Gamma-ray excesses and 511 keV line: Galactic- and extragalactic gamma-ray backgrounds constrain PBHs with $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$2–$$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$3 g (Ray, 2022, Belotsky et al., 2014).
• Microlensing: MACHO, EROS, OGLE, Subaru/HSC, and Kepler surveys exclude $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$4 in the range $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$5–$$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$6 M$$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$7 for compact objects, though extended mass functions may avoid the strictest bounds (Carr et al., 2016, Clesse et al., 2017).
• Radio and sub-mm backgrounds: Radio transients from PBH explosions or tunneling (as in loop quantum gravity with $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$8) have additional—sometimes unique—signatures in burst statistics and $$\psi(M) \propto \frac{1}{M}\exp \left[ -\frac{(\ln M-\ln M_c)^2}{2\sigma^2} \right],$$9-$M_c$0 correlations (Vidotto, 2018).
• Stellar dynamics: Disruption of wide binaries or puffing-up of ultra-faint dwarfs constrains $M_c$1 for $M_c$2 M$M_c$3, with heating rates modeling via the Fokker-Planck equation indicating $M_c$4–$M_c$5 allowed if PBHs comprise all DM (Zhu et al., 2017).

The open mass windows for PBH DM are currently (ranges are approximate and slightly model-dependent) (Carr et al., 2021, Green, 2024):

• $M_c$6–$M_c$7 g ("asteroid" to "sublunar" mass): Largely unconstrained as neither evaporation nor microlensing is sensitive in this regime.
• $M_c$8–$M_c$9 M$\sigma$0 ("stellar/intermediate" mass): Mildly allowed, with $\sigma$1 required by current dynamical, microlensing, and gravitational-wave bounds; broad lognormal spectra with peaks near $\sigma$2 may still saturate $\sigma$3 (Clesse et al., 2017).
• Planck-mass "relic" PBHs: If evaporation ends with a stable Planck-mass remnant, these evade all current non-gravitational constraints, but no observational signature is known (Carr et al., 2021).

3. Production of Dark Matter via PBH Evaporation and Exotic Channels

Beyond survival as gravitationally bound objects, PBHs can generate or transfer DM via Hawking evaporation, especially for light PBHs with $\sigma$4 g (Friedlander et al., 2023, Bernal et al., 2020, Bernal et al., 2020). In large extra dimensions (LEDs), quantum gravity with a lowered fundamental scale $\sigma$5 enhances micro black hole formation, evaporation, and possible Planckeon relics:

• Instantly Evaporating BHs (IEBHs): Promptly evaporate, emitting DM quanta, matching $\sigma$6 for $\sigma$7 and $\sigma$8 GeV (Friedlander et al., 2023).
• Delayed Evaporating BHs (DEBHs), Planckeon relics: Growth via accretion then eventual evaporation; stable relics can also contribute (Friedlander et al., 2023).
• Evaporation to dark sector: If DM self-interacts and thermalizes after PBH evaporation, the final yield can be "boosted" by entropy/energy conservation, allowing light (keV-scale) DM to evade structure-formation limits and lowering the PBH energy budget required (Bernal et al., 2020).
• Non-gravitational portal production: When freeze-out or freeze-in processes coexist, the PBH-induced entropy can dilute WIMP/FIMP yields, modifying allowed coupling ranges and re-opening (or tightening) parameter space (Bernal et al., 2020).

Quantum gravity scenarios introduce explosive decay or tunneling to white holes, with remnants potentially carrying DM (Vidotto, 2018). In models with a dark U(1) and heavy dark electrons, nearly extremal PBHs are stabilized against both Hawking and Schwinger processes, allowing stable, sub-lunar-mass DM objects from $\sigma$9 up to $\beta(M)$0 g (Bai et al., 2019).

4. Astrophysical Implications and Phenomenology

The presence of PBHs as a dark-matter component impacts cosmic structure and several critical observables:

• Gravitational waves: Binary mergers of PBHs, especially in the range $\beta(M)$1–$\beta(M)$2 M$\beta(M)$3, could account for some or all of the observed LIGO/Virgo merger rate. Forward projections for aLIGO design sensitivity show that $\beta(M)$4 can be excluded at $\beta(M)$5 for $\beta(M)$6–$\beta(M)$7 M$\beta(M)$8 if not detected as a distinct excess in 5–6 years (Kovetz, 2017, Clesse et al., 2017).
• Star cluster heating and core formation: Two-body relaxation between PBHs and stars in ultra-faint dwarf galaxies constrains the PBH mass to $\beta(M)$9–$\sigma(M)$0 if they dominate the DM halo, as more massive PBHs would overheat the stars and produce larger or hotter systems than observed (Zhu et al., 2017).
• Microlensing and stochastic backgrounds: PBH microlensing of multiply imaged quasars, bulge, or Andromeda stars manifests as distinctive, achromatic, symmetric variability, with Fourier spectra and lack of time-dilation consistent with a population of $\sigma(M)$1–$\sigma(M)$2 compact objects (Hawkins, 2011).
• Large-scale structure and seeding effects: Individual "seed" PBHs induce density perturbations $\sigma(M)$3, while Poisson noise from a PBH population collectively enhances small-scale power, potentially alleviating the "missing satellites" and small-scale structure problems of standard CDM (Carr, 2019).
• Clustered PBH signatures: Clusters or gravitationally bound pairs/mergers of PBHs may yield observable transient gamma-ray or neutrino signals (e.g., from rapid post-merger evaporation in the extremal scenario), subject to constraints from Fermi-LAT, HAWC, IceCube, and other high-energy observatories (Bai et al., 2019, Vidotto, 2018).

5. Theoretical Developments, Model-Dependent Loopholes, and Open Issues

Several theoretical considerations can either broaden or close the allowed parameter space for PBH DM:

• Non-Gaussian statistics and profile dependence: The high-δ tail relevant for PBH formation may depart strongly from Gaussian expectation due to inflationary quantum diffusion or nontrivial potential structure, leading to vast modifications in $\sigma(M)$4 predictions (Green, 2024, Carr et al., 2016).
• Critical collapse and threshold corrections: Realistic mapping from the primordial power spectrum to the PBH mass function requires inputs from numerical relativity—critical scaling, non-sphericity, and compaction function mapping are all important (Carr et al., 2016, Green, 2024).
• Extended mass functions: Accurately confronting model predictions with data demands integral constraints of the form

$\sigma(M)$5

where $\sigma(M)$6 encodes experimental bounds for each channel (Carr et al., 2016).

• Planck-mass relics: If evaporation leaves stable relics of $\sigma(M)$7, these evade all astrophysical constraints except those deriving from their (possibly undetectable) gravitation alone (Carr et al., 2021, Carr et al., 2016).
• Clustering and minihalo effects: PBHs might cluster subgalactically, altering both merger rates and local constraints (potentially relaxing CMB or dynamical bounds) (Green, 2024, Clesse et al., 2017).
• Cosmological phase impacts: Early matter domination (e.g., co-decaying dark sectors), modified thermal histories, and nonstandard reheating enable PBH production in otherwise excluded mass regimes (Georg et al., 2019).

6. Prospects for Detection and Future Probes

The PBH dark matter hypothesis remains testable across a broad array of astrophysical probes:

• Gravitational-wave observatories: Next-generation networks (Einstein Telescope, Cosmic Explorer, LISA/BBO) will push merger constraints and resolve the sub-Chandrasekhar tail of PBH mergers, potentially providing a smoking-gun signal (Clesse et al., 2017, Kovetz, 2017).
• MeV and sub-mm gamma-ray surveys: Upcoming AMEGO, GECCO, and survey missions are designed to close the "evaporation gap" and test the ultra-light PBH window $\sigma(M)$8–$\sigma(M)$9 g by looking for Hawking-radiated photons (Ray, 2022).
• 21 cm tomography and CMB spectroscopy: SKA and HERA (21 cm), as well as spectral distortion imprints in PIXIE/PRISM, will tighten bounds on PBH accretion at cosmic dawn and on dissipation from small-scale power (Green, 2024).
• Microlensing and astrometric campaigns: LSST, Euclid, WFIRST, and Gaia will constrain the low-mass end of the PBH mass distribution, including lunar-mass scales and isolated PBHs in the Milky Way halo (Green, 2024).
• Stochastic/background measurements: Stochastic gravitational-wave backgrounds from PBH binaries or evaporation (including quantum-gravity-induced bursts) act as integral tests for aggregate PBH populations (Vidotto, 2018, Clesse et al., 2017).

Active theoretical efforts target improved calculations of collapse statistics, cluster evolution, feedback effects, and CMB/GW signatures, with the field converging toward either a definitive detection or the closure of the remaining PBH DM mass windows.

---

Key References:

(Green, 2024, Carr et al., 2016, Ray, 2022, Clesse et al., 2017, Zhu et al., 2017, Bernal et al., 2020, Banks et al., 2020, Carr, 2019, Georg et al., 2019, Correa et al., 2022, Vidotto, 2018, Frampton et al., 2010, Belotsky et al., 2014, Hawkins, 2011, Bernal et al., 2020, Friedlander et al., 2023, Bai et al., 2019, Curd et al., 2024).