⚛️ general relativity

Saturation Equations of State in Critical Gravitational Collapse: The Primordial Black Hole Threshold

Using general-relativistic simulations of spherically symmetric collapse with a lattice gas equation of state, this study demonstrates that pressure stiffening near saturation density increases the primordial black hole formation threshold by approximately 0.50% while leaving the critical mass-scaling exponent unchanged, thereby proving that saturation effects can stabilize gravitational collapse without altering universal scaling laws.

Original authors: Benaoumeur Bakhti

Published 2026-06-30

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Benaoumeur Bakhti

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The "Goldilocks" Moment of the Universe

Imagine the very early universe as a giant, chaotic soup of energy and particles. Sometimes, a small clump of this soup gets squeezed so tightly by its own gravity that it collapses into a Primordial Black Hole (PBH).

For decades, scientists have studied exactly how much squeezing is needed to make this happen. They found a "tipping point" (a threshold). If a clump is just a tiny bit too light, it bounces back and disperses. If it's just a tiny bit too heavy, it collapses into a black hole.

This paper asks a simple question: What happens if the "soup" gets harder to squeeze as it gets denser?

Most previous studies assumed the soup was like a gas that gets easier to compress the more you push it (or at least, stays the same). But in reality, matter often gets "stiff" or "hard" when it gets very crowded, like a crowded room where people can't move anymore. This paper tests what happens to the black hole tipping point when the matter gets "stiff" near its maximum density limit.

The Experiment: A "Lattice Gas" Toy Model

Since we can't go back in time to test the early universe, the author built a mathematical "toy model" to simulate it.

The Analogy: Imagine a grid of parking spots (a lattice).
The Rule: Each spot can hold only one car.
The Physics: As you try to park more cars, the spots fill up. When the lot is almost full, it becomes incredibly difficult to squeeze in one more car. The "pressure" to stop you from adding more cars shoots up.

The author used this "single-occupancy parking lot" model to represent the early universe's matter. They ran super-computer simulations to see how this "stiffening" changed the rules of gravitational collapse.

The Main Findings

The paper discovered two main things:

1. The "Stiff" Wall Pushes the Tipping Point Up

The Result: Because the matter gets harder to squeeze as it gets denser, it fights back against gravity more effectively.
The Analogy: Imagine trying to crush a soda can. If the can is empty (like the old models), it crushes easily. But if the can is filled with a material that gets harder and harder to compress the more you squeeze it (like our parking lot), you need to push much harder to crush it.
The Numbers: The study found that this "stiffening" raises the threshold for forming a black hole by about 0.5%.

Why does this matter? In the world of black holes, a tiny 0.5% change in the threshold is huge. It's like moving a goalpost by an inch in a game where the ball is microscopic. This small shift could change the predicted number of black holes in the universe by tens of percent.

2. The "Recipe" for Black Hole Size Stays the Same

The Result: Even though it takes more effort to make a black hole, the way the black hole grows (the scaling law) remains exactly the same as it was in the old, simpler models.
The Analogy: Imagine two different recipes for baking a cake. One uses a standard oven (radiation fluid), and the other uses an oven that gets slightly hotter as the cake rises (the stiff lattice gas).

The "stiff" oven makes it harder to get the cake to rise (higher threshold).
However, once the cake does rise, the relationship between how much batter you used and the final size of the cake is identical in both ovens.
The Numbers: The "scaling exponent" (a number that describes how mass scales with the collapse) remained 0.357 in both cases.

Why Did the Size Rule Stay the Same?

You might wonder: If the rules changed, why didn't the size rule change?

The author explains this with a "local vs. global" argument.

The Local View: The critical moment where a black hole is born happens at densities that are still relatively "low" compared to the maximum limit of the parking lot. At these lower densities, the "stiff" parking lot behaves almost exactly like the simple gas model.
The Analogy: It's like driving a car. If you drive on a road that is mostly smooth but has a few potholes at the very end, your driving style (the "exponent") for the smooth part of the road doesn't change, even if the potholes force you to start the trip from a different spot (the "threshold").

What This Paper Does NOT Say

It is important to stick to what the paper actually claims:

It is not a prediction for the real universe: The "lattice gas" is a mathematical toy, not a description of real nuclear matter. The author explicitly states this is a "proof of principle."
It does not claim universality: The author warns that the "size rule" (exponent) can change if the matter is very different. It only stayed the same here because the toy model was very similar to the standard model in the specific range where black holes form.
It does not predict a specific number of black holes: The 0.5% shift is a demonstration of a mechanism, not a forecast for what we will see in the sky.

Summary

Think of this paper as a controlled experiment in a physics lab. The scientist asked: "If we make the matter in the early universe slightly 'stiffer' as it gets crowded, does it change how black holes form?"

The answer is: Yes, it makes it harder to form them (raising the threshold), but it doesn't change the fundamental rules of how they grow once they start. This gives scientists a new tool to understand how different types of matter might have influenced the birth of black holes in the early universe.

Technical Summary: Saturation Equations of State in Critical Gravitational Collapse

Problem Statement
The threshold and scaling laws governing gravitational critical collapse are known to depend sensitively on the matter equation of state (EOS). While primordial black hole (PBH) studies typically assume a simple radiation fluid ( $p = \omega\rho$ with constant $\omega=1/3$ ), realistic dense matter (such as nuclear or quark matter) exhibits "pressure stiffening" as it approaches a maximum density. This feature—absent in standard radiation models—implies that matter cannot be compressed without limit. The central problem addressed is how such a saturation-induced stiffening modifies the critical threshold ( $\delta_c$ ) for PBH formation and whether it alters the critical mass-scaling exponent ( $\gamma$ ).

Methodology
To isolate the gravitational consequences of a density cap without introducing complex microphysics, the authors employ a controlled toy model: a single-occupancy lattice gas.

Equation of State: The model utilizes the analytically tractable EOS $p = -T \ln(1-\rho)$ , where $\rho$ is the density fraction ( $0 < \rho < 1$ ) and $T$ is a dimensionless stiffness parameter. This EOS features a hard density cap at $\rho=1$ and a density-dependent sound speed $c_s^2 = T/(1-\rho)$ , which diverges as saturation is approached.
Numerical Framework: The authors perform general-relativistic simulations of spherically symmetric collapse using the Misner-Sharp-Hernandez formalism. The simulations are conducted on a flat Friedmann-Robertson-Walker (FRW) background in the radiation-dominated era ( $T=1/3$ ).
Numerical Techniques: To handle the diverging sound speed and maintain causality ( $c_s \leq 1$ ), the study employs adaptive spectral resolution (increasing Chebyshev points in high-density regions) and CFL-aware adaptive timestepping. A causal cap is applied to the sound speed, and a smooth density limiter prevents numerical overshoots beyond $\rho=1$ .
Analysis: The critical threshold $\delta_c$ is determined via bisection searches bracketing the boundary between black hole formation and perturbation dispersal. The critical exponent $\gamma$ is extracted from the power-law scaling of black hole mass near the threshold ( $M_{BH} \propto (\delta - \delta_c)^\gamma$ ).

Key Results

Threshold Shift: The saturation stiffening increases the critical threshold for PBH formation. Relative to the polytropic radiation reference, the lattice model yields a shift of:
$\Delta \delta_c = 0.00471 \pm 0.00004$
This represents an upward shift of approximately 0.50%. The shift is accumulated within the causal density window ( $\rho < 1-T$ ) and is attributed to the enhanced pressure resistance provided by the density-dependent stiffening.
Critical Exponent Invariance: Despite the shift in the threshold, the critical mass-scaling exponent remains unchanged within numerical precision:
$\gamma = 0.357 \pm 0.001$
This value is consistent with the standard radiation-fluid value ( $\gamma \approx 0.3558$ ).
Stiffness Dependence: The threshold shift exhibits an approximate power-law dependence on the stiffness parameter $T$ ( $\Delta \delta_c \propto T^{0.5}$ ). Increasing $T$ strengthens the pressure support at all densities, monotonically raising $\delta_c$ .
Profile Independence: The threshold offset is robust across various initial curvature profiles (Gaussian, narrow, broad, and broken power-spectrum), confirming the shift is a property of the EOS rather than the perturbation shape.

Interpretation and Significance
The paper argues that the agreement in the critical exponent $\gamma$ does not imply universality across all equations of state. Citing Maison's work, the authors note that $\gamma$ generally varies with the EOS parameter $k$ in $p=k\rho$ models. The observed agreement here arises because the lattice EOS coincides with the radiation fluid at low densities ( $c_s^2(0) = 1/3$ ) and the near-critical dynamics sample a density range where the deviation from the radiation sound speed is a "mild perturbation." Thus, the locally relevant sound speed remains close to $1/3$ , preserving the radiation exponent.

Claims and Limitations
The authors explicitly frame this work as a proof of principle rather than a prediction for specific cosmological scenarios.

Mechanism Demonstration: The primary contribution is demonstrating that saturation-induced stiffening can stabilize gravitational collapse, shifting the PBH threshold upward, while introducing a linear-response framework (via a stiffness integral) to assess such effects.
Model Limitations: The lattice gas is an idealization. The authors do not claim it describes relativistic degenerate matter or specific early-universe physics. Consequently, the numerical value of the 0.5% shift is illustrative of the mechanism's magnitude, not a forecast for realistic PBH abundance.
Future Directions: The paper suggests applying this pipeline to more realistic EOSs, particularly the QCD equation of state near the chiral crossover, where sound speed variations are known to affect PBH formation. It also outlines extensions to non-spherical collapse and modified gravity as future steps.