Inflationary relics from an Ultra-Slow-Roll plateau

This paper investigates the formation of primordial black holes in ultra-slow-roll inflation scenarios with a sharp transition, revealing through numerical simulations that while both vacuum bubble and adiabatic perturbation channels contribute, the latter dominates the mass function by a factor of $10$to$100$ and is primarily driven by mean profiles.

The Gist

Imagine the early universe as a giant, rolling ball (the "inflaton field") trying to get down a very specific hill (the "potential"). Usually, this ball rolls smoothly and quickly. But in the scenario this paper explores, the hill has a strange, flat plateau in the middle.

Here is the story of what happens when the ball hits that flat spot, told through the lens of Primordial Black Holes (PBHs)—tiny black holes that formed right after the Big Bang and might make up the "Dark Matter" holding our galaxy together.

The Setup: The Flat Plateau

Think of the universe's expansion as a car driving down a road.

• Normal Inflation: The car is speeding down a steep hill. It moves fast and doesn't stop.
• Ultra-Slow-Roll (USR): Suddenly, the road flattens out into a long, perfectly level plateau. The car's engine (the inflaton) keeps running, but because the road is flat, the car slows down to a crawl. It almost stops.

This "crawling" phase is crucial. Because the car is moving so slowly, tiny bumps in the road (quantum fluctuations) get magnified. Instead of being small ripples, they become massive mountains.

The Two Ways Black Holes Form

The paper investigates two different ways these massive "mountains" can turn into black holes. Think of them as two different escape routes for the car.

1. The "Stuck Car" (Vacuum Bubbles)

Imagine the car is on the flat plateau. Most of the time, it slowly rolls forward. But sometimes, a random gust of wind (a quantum fluctuation) pushes the car backward.

• The Trap: If the wind pushes it back hard enough, the car gets stuck in a little dip on the plateau. It can't get over the edge to roll down the rest of the hill.
• The Result: This "stuck" patch of the universe keeps inflating forever, while the rest of the universe moves on. From the outside, this trapped patch looks like a bubble. Eventually, this bubble collapses into a black hole.
• The Paper's Finding: The authors found that while these "stuck bubbles" do form, they are actually quite rare compared to the other method. It's like trying to park a car in a tiny, specific spot in a massive parking lot; it happens, but it's not the main way cars get parked.

2. The "Crash" (Adiabatic Collapse)

Now, imagine the car isn't stuck, but the road itself is bumpy. Because the car was moving so slowly on the plateau, those bumps got huge.

• The Crash: When the car finally leaves the plateau and speeds up again, those huge bumps in the road become so dense that gravity takes over. The space itself collapses under its own weight, crushing everything into a black hole.
• The Paper's Finding: This is the dominant method. The paper calculates that this "crash" method produces about 10 to 100 times more black holes than the "stuck car" method. It's the main way the universe fills its dark matter inventory.

The "Shape" of the Problem

The authors didn't just look at the average scenario; they looked at the shape of the bumps.

• The Average Bump: They calculated what a "typical" large bump looks like.
• The Weird Bumps: They also looked at weird, jagged bumps that are slightly different from the average.
• The Surprise: They found that while weird bumps change the exact number of black holes slightly, the "average" bump is responsible for almost all of them. It's like saying that while some people are very tall and some are very short, the average height is what determines the size of the door you need to build.

The "Math Magic" (Templates)

To predict how many black holes form, scientists usually use a simple formula (a template) to translate "small wiggles" into "big black holes."

• The Old Formula: The old formula worked well for simple hills, but it broke down on this flat plateau. It was like trying to use a ruler to measure a curve; it just didn't fit.
• The New Formula: The authors created a new, more complex "template" (a generalized formula) that accounts for the weird physics of the flat plateau. They tested this new formula against super-computer simulations and found it was very accurate, except right at the very edge where things get chaotic.

The "Sponge" Effect (Quantum Diffusion)

What happens to the "stuck bubbles" after they form?

• Imagine a sponge. The bubble is the sponge, and the "trapped" part is the wet part.
• Over time, tiny bits of the sponge dry out (quantum diffusion). The bubble gets full of holes.
• The Result: Even though the bubble gets full of holes, the overall size of the sponge doesn't shrink. It still acts like a big object that will eventually collapse into a black hole. The "holes" just mean that inside the black hole, there are tiny, separate universes (baby universes) growing in the gaps.

The Bottom Line

This paper is a detailed map of how the universe might have created a population of tiny black holes.

1. The Main Event: The most common way these black holes form is through the "crash" of huge density waves (Adiabatic channel).
2. The Side Event: A smaller number form from "stuck" patches of the universe (Bubble channel).
3. The Takeaway: If you want to find Dark Matter in the form of these tiny black holes, you should look for the signatures of the "crash" method, as it is the heavy lifter. The "stuck bubble" method is a fascinating side story, but it's not the main show.

In short: The universe had a flat spot where things slowed down, creating massive waves. Most of these waves crashed into black holes, while a few got stuck in bubbles. The crash was the winner.

Technical Summary

Here is a detailed technical summary of the paper "Inflationary relics from an Ultra-Slow-Roll plateau" by Escrivà, Garriga, and Pi.

1. Problem Statement

The paper investigates the formation of Primordial Black Holes (PBHs) in single-field inflationary models featuring an Ultra-Slow-Roll (USR) phase with a flat potential plateau followed by a sharp transition to standard slow-roll.

While USR phases are known to significantly enhance the power spectrum of curvature perturbations (leading to PBH formation via the collapse of large adiabatic fluctuations), the authors address a specific, often overlooked mechanism: relic vacuum bubbles. In scenarios where the inflaton field encounters a flat plateau, quantum fluctuations can trap regions of the field on the plateau, preventing them from rolling down. These trapped regions evolve into vacuum bubbles (baby universes connected by transient wormholes), which appear as PBHs to an external observer.

The core problem is to quantify and compare the abundance of PBHs produced via two coexisting channels:

1. Adiabatic Channel: Gravitational collapse of large, non-Gaussian curvature perturbations.
2. Bubble Channel: Formation of vacuum relics due to the inflaton getting trapped on the plateau.

The authors aim to determine which channel dominates, how the statistical distribution of initial conditions (specifically shape dispersion) affects the formation thresholds, and to derive accurate analytical templates for the non-linear curvature perturbation $\zeta$ in this specific USR context.

2. Methodology

The study employs a combination of analytical derivations, statistical analysis of initial conditions, and fully non-linear numerical simulations.

• Model Setup:

  • A single scalar field $\phi$ with a piecewise potential $V(\phi)$ featuring a flat plateau ($V \approx V_0$) connected to slopes via parabolic regions.
  • The background dynamics are solved using the Klein-Gordon equation in a FLRW background, characterized by the USR phase where the second Hubble flow parameter $\epsilon_2 = -6$.
  • The transition to slow-roll is sharp, implying a transient mass scale $M \gg H$.

• Statistics of Initial Conditions:

  • The authors analyze the probability distribution of initial field profiles $\delta\phi(r)$ and momentum profiles $\delta\pi(r)$ at the end of the USR phase ($N_*$).
  • They define a mean profile $\bar{\delta\phi}$ derived from the power spectrum and calculate shape dispersions ($\Delta\phi, \Delta\pi$) representing deviations from the mean.
  • Realizations are parameterized by integers $(n, m)$, representing deviations in standard deviations from the mean for the field and momentum, respectively.

• Numerical Simulations:

  • Bubble Formation: The non-linear Klein-Gordon equation is solved in spherical symmetry to determine the critical threshold amplitude $\mu_c^{bub}$ required for the field to become trapped on the plateau.
  • Adiabatic Collapse: For profiles that do not form bubbles, the authors use the $\delta N$ formalism to map the initial field perturbations to the non-linear curvature perturbation $\zeta$. They then use the relativistic hydrodynamic code SPriBHoS-II to simulate the gravitational collapse of these perturbations to find the critical threshold $\mu_c^{adi}$ for PBH formation.

• Analytical Template:

  • The authors derive a generalized mapping $\zeta[\zeta_G]$ (where $\zeta_G$ is the Gaussian curvature perturbation) that accounts for the specific dynamics of USR (where $\delta\pi$ decays as $a^{-2}$ rather than $a^{-3}$). This extends the standard logarithmic template used in barrier models.

3. Key Contributions

1. Generalized $\zeta[\zeta_G]$ Template: The paper derives a new analytical template for the non-linear curvature perturbation in USR scenarios with sharp transitions. Unlike previous models relying on a potential barrier, this model has no barrier. The derived formula involves solving a cubic equation (via Cardano's formula) and accurately captures the non-Gaussianity, including the logarithmic divergence that signals the transition to the bubble channel.
2. Characterization of Shape Dispersion: The study systematically quantifies how deviations from the mean initial profile (parameterized by $n$ and $m$) affect the critical thresholds for both bubble formation and adiabatic collapse. It demonstrates that while dispersion shifts thresholds, the mean profile dominates the statistical weight.
3. Type-II Fluctuations: The authors confirm that vacuum bubbles formed in this USR scenario are generically surrounded by Type-II curvature fluctuations (characterized by a "neck" where $1 + r\zeta'(r) = 0$), similar to bubble formation in barrier models.
4. Quantum Diffusion Analysis: The paper analyzes the effect of quantum diffusion on the trapped bubbles. It shows that while diffusion creates a "sponge-like" fractal structure (perforating the bubble with slow-roll regions), the overall comoving size of the relic remains fixed, ensuring the eventual formation of a PBH upon horizon re-entry.

4. Key Results

• Dominance of the Adiabatic Channel:

  • The adiabatic channel (standard collapse of curvature perturbations) dominates the PBH abundance.
  • The adiabatic contribution is larger than the bubble-induced contribution by a factor of **$\sim 10$ to $100$** (specifically, the ratio$f_{adi}/f_{bub} \approx 50$ for the mean profile).
  • Consequently, vacuum bubbles contribute only $\sim 2\%$ to the total PBH dark matter abundance in this model.

• Critical Thresholds:

  • Bubble Threshold ($\mu_c^{bub}$): Determined numerically. For the mean profile, $\mu_c^{bub} \approx 2.61 \times 10^{-5}$.
  • Adiabatic Threshold ($\mu_c^{adi}$): Determined numerically. For the mean profile, $\mu_c^{adi} \approx 2.46 \times 10^{-5}$.
  • The difference between the two thresholds is roughly one quantum step ($\delta\phi \sim H/2\pi$).
  • The analytical template $\zeta[\zeta_G]$ agrees well with numerical $\delta N$ results, deviating only near the logarithmic divergence point.

• Mass Functions:

  • The PBH mass functions for both channels are computed. The adiabatic channel produces a mass distribution concentrated near the scale set by the USR enhancement ($k_{peak} \approx 10^{13} \text{Mpc}^{-1}$).
  • The bubble channel produces a sub-dominant mass function with a similar scale but lower amplitude.

• Impact of Dispersion:

  • Positive deviations ($n, m > 0$) increase the formation thresholds and are statistically suppressed.
  • Negative deviations ($n, m < 0$) lower the thresholds but are also statistically suppressed relative to the mean.
  • Conclusion: The mean profile $(n,m)=(0,0)$ provides the overwhelming majority of the PBH abundance for both channels.

5. Significance

• Robustness of USR PBH Production: The study confirms that USR plateaus can produce PBHs even without explicit potential barriers, validating the mechanism for a broad class of inflationary models.
• Clarification of Competing Channels: It resolves the question of whether vacuum bubbles (often cited as a dominant source in barrier models with large non-Gaussianity $\beta \gtrsim 3.1$) compete with adiabatic collapse in USR scenarios. The result is that for the specific case of a sharp USR transition ($\beta=3$), the adiabatic channel remains dominant.
• Observational Predictions: The paper notes that the induced gravitational waves from the enhanced curvature perturbations would peak in the mHz range (detectable by LISA, Taiji, TianQin) if these PBHs constitute dark matter. The sub-dominant bubble channel does not significantly alter this prediction but adds a minor feature to the mass function.
• Methodological Advance: The derivation of the generalized $\zeta[\zeta_G]$ template and the rigorous treatment of shape dispersion provide a robust framework for analyzing PBH formation in non-standard inflationary dynamics where the separate universe approximation breaks down.

In summary, the paper establishes that while vacuum bubbles are a generic feature of USR plateaus, the standard adiabatic collapse mechanism is the primary source of PBHs in this specific regime, with the mean initial profile driving the majority of the production.