Primordial Black Holes

Primordial black holes connect rare early-universe fluctuations to nonlinear collapse, an evolving compact-object population, semiclassical evaporation, gravitational-wave backgrounds, and heterogeneous observations. This path is for advanced students and researchers who want to carry one declared population through that entire chain. It does not offer an undated exclusion plot; it teaches how to decide what a formation or constraint statement actually assumes.

Helpful background. Relativistic cosmology supplies homogeneous backgrounds, perturbations, thermal history, structure growth, light-cone observables, and cosmological inference. Black-hole dynamics and collapse supplies finite-amplitude collapse, trapped regions, critical phenomena, horizon growth, interiors, and nonlinear evolution; use Black-Hole Physics for a focused review. Probability and statistical inference supplies likelihoods, priors, posterior checks, uncertainty propagation, and selection effects. These are recommended rather than admission requirements, but develop any missing cosmology or probability foundation before calculating an abundance.

From primordial statistics to an observed population

By the end of the path, you should be able to:

• relate primordial perturbations and relativistic collapse to a formation fraction and mass distribution;
• model subsequent evolution, mass functions, and environmental effects;
• combine constraints while retaining selection, calibration, population, and model assumptions; and
• connect primordial black holes to scalar-induced gravitational waves and dark-matter questions without conflating formation channels.

The central dependency chain is

$$\begin{aligned} \text{primordial statistics}&\to\text{smoothed profile}\to\text{collapse}\\ &\to\text{initial population}\to\text{evolution}\\ &\to\text{observable likelihood}. \end{aligned}$$

A power spectrum alone does not fix a rare tail, a horizon mass is not automatically the black-hole mass, and a monochromatic limit is not automatically a limit on an extended or clustered population.

Check your preparation

1. How does a curvature perturbation become a gauge-controlled density or compaction variable suitable for a collapse criterion?
2. Why does a small change in the variance of a nearly Gaussian perturbation produce a much larger fractional change in a rare-tail abundance?
3. What information is lost when an upper-limit curve is detached from its likelihood, selection function, calibration, and population model?

Answers and review

1. Specify slicing, transfer through the background, smoothing window, and profile; then express collapse in a geometrically meaningful quantity appropriate to that prescription. Review cosmological perturbations through Relativistic Cosmology.
2. Collapse samples the exponentially suppressed tail, so its logarithm scales approximately as minus the threshold squared over twice the variance. The worked calculation makes this precise. Review Gaussian probabilities in probability and inference if needed.
3. One loses the conditions that map a population to detected or excluded events and therefore cannot consistently recast or combine the result. The curve alone is not a portable likelihood.

If answers 1 and 2 are secure, begin with the cosmological perturbation phases. Review either topic as needed before the formation calculation. An uncertain third answer is expected at entry but must be resolved before the observational-constraint and research-synthesis phases.

Route from formation to evidence synthesis

Phase 1 — Initial perturbations and collapse

Carry a single gauge, smoothing, profile, and statistical prescription from primordial production through the nonlinear formation calculation.

1. Establish the perturbative initial data with cosmological perturbations.
2. Study possible primordial origins in inflation and primordial fluctuations.
3. Follow nonlinear growth through relativistic structure formation.
4. Orient the full population problem with primordial black holes.
5. Develop collapse in primordial-black-hole formation from density perturbations.

Phase 2 — Population evolution and observational response

Evolve the population into the relevant epoch and pass it through each experiment’s response and likelihood before reading a constraint.

6. Add controlled mass loss through Hawking flux, backreaction, and evaporation dynamics.
7. Mark the extrapolation boundary with endpoint regimes and remnant questions.
8. Study primordial-black-hole observational constraints with each likelihood’s mass, spatial, and astrophysical assumptions attached.
9. Connect the same primordial statistics to scalar-induced gravitational waves.

Phase 3 — Evidence synthesis and source identification

Compare mutable constraints and then test whether the inferred population is distinguishable from astrophysical alternatives.

10. Compare mutable evidence in the primordial-black-hole constraint synthesis.
11. Finish by contrasting populations with identifying astrophysical black holes.

Worked reasoning: exponential abundance sensitivity

As a diagnostic model, let a smoothed density variable $\delta$ be Gaussian with variance $\sigma^2$, and suppose collapse occurs when $\delta\gt\delta_c$. The formation fraction is

$$\beta =\int_{\delta_c}^{\infty}\frac{d\delta}{\sqrt{2\pi}\sigma} \exp\!\left(-\frac{\delta^2}{2\sigma^2}\right) =\frac12\operatorname{erfc}\!\left(\frac{\nu}{\sqrt2}\right), \qquad \nu=\frac{\delta_c}{\sigma}.$$

For $\nu\gg1$,

$$\beta\simeq\frac{e^{-\nu^2/2}}{\sqrt{2\pi}\,\nu}, \qquad \frac{d\log\beta}{d\log\sigma}\simeq\nu^2+1.$$

Thus a percent-level change in $\sigma$ can produce an order-$\nu^2$ percent-level change in the inferred abundance. This is why the window, transfer relation, collapse threshold, and non-Gaussian tail cannot be treated as cosmetic choices.

The calculation is a sensitivity model, not a universal formation law. In relativistic collapse, $\delta_c$ depends on the profile and equation of state; critical collapse spreads masses below the horizon scale; non-Gaussian statistics alter the tail; and the initial fraction must be redshifted and evolved before comparison with a present population. A numerical integral can accurately evaluate the wrong $\beta$ if these inputs are inconsistent.

Checkpoints

Checkpoint 1. Two studies quote different collapse thresholds. Is averaging the numbers a defensible way to estimate the abundance?

Conceptual hint

Ask whether their slicing, smoothing variable, window, profile, and equation of state match.

Setup hint

Make a comparison table containing slicing, collapse variable, smoothing window, profile family, and equation of state for each quoted threshold.

Next-step hint

Translate to one prescription where possible; otherwise propagate each matched threshold-and-variance pair and treat the spread as model dependence.

Solution

No. First translate both calculations to a common collapse variable and profile prescription. If they describe genuinely different profiles or backgrounds, propagate each threshold with its corresponding variance and treat the difference as model dependence; averaging mismatched definitions has no controlled meaning.

Checkpoint 2. May one constrain an extended mass function by requiring it to lie below every published monochromatic upper-limit curve point by point?

Conceptual hint

An experiment responds to an integral over the population, not independently to each plotted mass.

Setup hint

Write the expected observable as an integral of the evolved mass and spatial distribution against the experiment’s response and efficiency.

Next-step hint

Return that prediction to the original nuisance model and likelihood, then identify shared data or systematics before combining experiments.

Solution

Not generally. Recompute the expected observable by integrating the evolved mass and spatial distribution through that analysis’s response, efficiency, nuisance model, and likelihood. Only then combine conditionally independent data while accounting for shared systematics.

Transfer task

Choose a narrow feature in a primordial curvature spectrum. Carry it to both a primordial-black-hole population and a scalar-induced gravitational-wave spectrum. State the statistics beyond the two-point function, window and transfer conventions, collapse prescription, mass mapping, thermal history, and likelihood assumptions. Identify one change that shifts both signals and one formation scenario in which the gravitational-wave relation does not apply.

Capstone: one population from formation to constraint

Build a reproducible analysis for one declared primordial-black-hole scenario. It should:

• define the primordial random field, smoothing and profile prescription, collapse threshold, and mass mapping;
• compute an extended initial mass function with rare-tail sensitivity checks;
• evolve mass, number, and spatial statistics through the relevant thermal and astrophysical history;
• forward-model at least two observables with their calibration, selection, nuisance, and covariance assumptions; and
• report a conditional conclusion that distinguishes a dark-matter fraction, a formation-model constraint, and a generic black-hole-population statement.

Judge the work by formation consistency, population conservation, likelihood fidelity, and scope discipline. Exemplary work repeats the conclusion under a defensible alternative tail or collapse prescription and explains the change.

When you are ready to move on

You are ready to move on when you can derive the rare-tail sensitivity, solve both checkpoints without hints, complete the paired-signal transfer, and meet all four capstone criteria. Return to Relativistic Cosmology for broader early-universe dynamics, continue to Black-Hole Thermodynamics and Semiclassical Gravity for evaporation limits, or use the Research destination for dated observational assessments.