Eigenvalue formulation of Stochastic Inflation and application to large perturbation generating inflationary features

This paper introduces a novel eigenvalue technique to solve the adjoint Fokker-Planck equation for the probability distribution of inflationary e-folds, revealing a previously overlooked power-law intermediate regime in quantum diffusion and characterizing how constant drift potentials qualitatively alter the distribution's peak and tail behavior in narrow- versus broad-well limits.

The Gist

The Big Picture: Predicting the "Unlikely"

Imagine the early universe as a giant, expanding balloon. Inside this balloon, there is a field (called the "inflaton") that drives the expansion. Usually, this field rolls down a gentle, smooth hill, creating a very predictable, calm universe. This is like a ball rolling slowly down a long, flat driveway.

However, sometimes this hill has a weird bump or a dip. When the field rolls over these features, it can get stuck or jitter around wildly. This jittering is caused by quantum mechanics—the universe's version of "static noise."

The authors of this paper are trying to answer a specific question: How likely is it for this field to get stuck in a weird spot for a really long time?

If the field gets stuck for a long time, it creates a massive burst of energy in that specific spot. When the universe cools down, these bursts can collapse into tiny, dense black holes called Primordial Black Holes (PBHs). These are the "dark matter" candidates the paper is interested in.

To find out how many of these black holes might exist, we need to know the probability of the field getting "stuck." This probability is described by a mathematical curve called a Probability Distribution Function (PDF).

The Problem: The Math is Too Hard

The paper explains that calculating this probability curve is incredibly difficult. It's like trying to predict exactly where a drunk person will end up after wandering through a maze for a long time. The math involved (Fokker-Planck equations) is usually solved using a mix of different tricks, but no one had found a single, self-contained "master key" (an eigenvalue technique) to solve it completely on its own.

The Solution: A New "Spectral" Key

The authors developed a new mathematical technique they call an eigenvalue formulation.

The Analogy: Tuning a Guitar
Imagine the universe's behavior is like a guitar string. When you pluck it, it doesn't just make one sound; it makes a complex chord made of many different notes (frequencies) vibrating at once.

• The notes are the "eigenvalues" (mathematical numbers that define the speed of decay).
• The shape of the vibration is the "eigenfunction."

The authors' new method breaks the complex problem of the field's movement down into these individual "notes." Instead of guessing the whole shape of the probability curve, they calculate each note individually and then stack them on top of each other to rebuild the full picture. This allows them to calculate the exact shape of the probability curve without needing to rely on other, less precise methods.

What They Found: Three Different "Zones"

Using this new method, they tested two scenarios: a field with no "drift" (just pure jittering) and a field with a constant "drift" (jittering while being pushed).

1. The Drift-Free Case (Pure Jittering)

Imagine a ball bouncing randomly in a box with no wind pushing it.

• The Peak: Most of the time, the ball exits the box quickly. The probability curve has a high peak here.
• The Middle (The Surprise): The authors found a hidden "middle zone" between the quick exit and the long wait. In this zone, the probability doesn't drop off smoothly; it follows a specific power law (it drops like $1/N^{1.5}$). They hadn't emphasized this "middle ground" in previous studies.
• The Tail: If the ball stays in the box for a very long time, the probability drops off exponentially (it becomes incredibly rare). This is the "tail" that determines how many black holes form.

2. The Constant-Drift Case (Jittering with a Push)

Now imagine the ball is in a box, but there is a gentle wind pushing it toward the exit.

• The Narrow Well (Small Box): If the box is small, the wind doesn't matter much. The ball still exits mostly by random bouncing. The probability curve looks almost the same as the drift-free case, just slightly tweaked.
• The Broad Well (Huge Box): If the box is massive, the wind becomes the dominant force.
  • The Peak: The ball exits much faster than random chance would suggest because the wind pushes it out. The peak of the probability curve is much higher and sharper.
  • The Tail: The "long tail" (the chance of the ball staying in for a huge amount of time) is strongly suppressed. The wind makes it almost impossible for the ball to stay stuck for a long time. This means fewer primordial black holes would form in this scenario compared to the drift-free case.

The "Piecewise" Puzzle

When dealing with the "Broad Well" (the huge box with strong wind), the math gets tricky. The authors realized that the "notes" (eigenvalues) behave differently depending on how high up the scale you go.

• For the first few notes, they behave one way.
• For the higher notes, they behave another way.

To solve this, they built a piecewise construction—like building a bridge where the first half is made of steel and the second half is made of wood, but they are joined perfectly so the bridge holds. They found that while this "patchwork" math works well for the tail of the curve, it creates "glitches" near the peak. To fix this, they used a different mathematical shortcut (involving special functions called Theta functions) that smoothed out the peak perfectly.

Summary of Results

1. New Tool: They created a self-contained mathematical method to calculate the probability of the universe's field getting "stuck."
2. Hidden Middle: They identified a specific "power-law" behavior in the middle of the probability curve that was previously overlooked.
3. Drift Matters:
   • If the field is just jittering (no drift), there is a moderate chance of forming black holes.
   • If the field is being pushed (drift) through a wide feature, the chance of it getting stuck long enough to form a black hole drops significantly.
4. Accuracy: Their method confirms previous results for simple cases but provides a much more detailed and accurate picture for complex scenarios involving "features" in the universe's potential.

In short, the authors built a better calculator to predict how often the early universe might have created tiny black holes, revealing that the "wind" (drift) in the universe's landscape plays a crucial role in whether these black holes can form.

Technical Summary

Technical Summary: Eigenvalue Formulation of Stochastic Inflation and Application to Large Perturbation Generating Inflationary Features

Problem Statement
Primordial Black Holes (PBHs) are a leading dark matter candidate and potential progenitors for supermassive black holes. Their formation relies on the gravitational collapse of rare, large-amplitude inflationary density fluctuations. Accurately predicting the abundance of PBHs requires a precise determination of the non-Gaussian tail of the probability distribution function (PDF) of these perturbations. While stochastic inflation is a powerful framework for calculating these tails, existing methods often rely on a combination of techniques, such as characteristic functions, to derive the full PDF. There is a lack of a closed-form, self-contained eigenvalue solution for the PDF of the stochastic number of e-folds, $N$, within the literature.

Methodology
The authors develop a new, self-contained eigenvalue technique (spectral method) to solve the adjoint Fokker-Planck Equation (FPE) governing the PDF, $P(N)$, of the stochastic number of e-folds. The methodology proceeds as follows:

1. Formalism: The stochastic dynamics of the coarse-grained inflaton field $\Phi$ are described by an adjoint FPE. The authors express the general solution as a spectral decomposition: $P(N; \Phi) = \sum_n c_n \Psi_n(\Phi) e^{-\Lambda_n N}$, where $\Psi_n$ are eigenfunctions and $\Lambda_n$ are eigenvalues of the adjoint Fokker-Planck operator.
2. Boundary Conditions: The problem is defined on a field interval $[\phi_A, \phi_R]$ with an absorbing boundary at $\phi_A$ (marking the end of the diffusion regime) and a reflecting boundary at $\phi_R$ (preventing diffusion toward the slow-roll regime).
3. Eigenvalue Determination: The eigenfunctions satisfy a Sturm-Liouville problem with a specific weight function $w(\Phi)$. The authors derive the coefficients $c_n$ and the normalization constant $B$ using Fourier's trick and the properties of the Dirac delta function derivative as an initial condition.
4. Application to Models: The formalism is applied to two specific regimes relevant to PBH formation:
   • Drift-free quantum diffusion: Quantum diffusion along a flat potential with no classical drift.
   • Constant-drift inflation: Inflation with a constant classical drift term, analyzed in both "narrow-well" (diffusion-dominated) and "broad-well" (drift-dominated) limits.

Key Contributions and Results

• Self-Contained Eigenvalue Solution: The paper provides the first closed-form eigenvalue solution for the PDF in stochastic inflation that does not rely on partially incorporating other techniques (like characteristic functions).
• Drift-Free Diffusion:
  • The method recovers the known exact PDF for drift-free diffusion, confirming the validity of the technique.
  • Intermediate Power-Law Regime: The authors identify an intermediate regime between the peak and the exponential tail of the PDF where the behavior follows a power law, $P(N) \propto N^{-3/2}$. This regime, previously unemphasized, arises from the collective contribution of higher modes in the spectral sum.
  • The tail of the PDF exhibits an exponential decay, $P(N) \propto e^{-\pi^2 N / (8\varepsilon)}$, where $\varepsilon$ is the dimensionless diffusion width.
• Constant-Drift Inflation:
  • Narrow-Well Limit ($\alpha \ll 1$): In this regime, where diffusion dominates over drift, the PDF is qualitatively similar to the drift-free case but with a mildly suppressed tail. The eigenvalues receive small corrections proportional to the drift parameter $\alpha$.
  • Broad-Well Limit ($\alpha \gg 1$): In this regime, classical drift dominates. The authors find that determining the full set of eigenvalues requires a piecewise construction of the spectrum because the behavior of the eigenvalues changes depending on the mode number $n$ relative to $\alpha$.
  • PDF Characteristics: In the broad-well limit, the PDF exhibits a qualitatively different shape: an enhanced peak and a strongly suppressed exponential tail compared to the drift-free case.
  • Average E-folds: In the narrow-well limit, the stochastic average number of e-folds $\langle N \rangle$ is significantly smaller than the classical prediction $N_{cl}$. Conversely, in the broad-well limit, $\langle N \rangle \approx N_{cl}$, reflecting the dominance of classical drift over stochastic transport.

Significance and Claims
The paper claims that the developed eigenvalue formalism is a robust, self-contained tool for determining the full PDF of large inflationary perturbations. Its significance lies in:

1. Accuracy: It provides a rigorous method to calculate the non-perturbative tails of the PDF, which are critical for estimating PBH abundance.
2. Novelty: The identification of the $N^{-3/2}$ power-law behavior in the intermediate regime offers new insight into the structure of the PDF, potentially affecting the abundance calculations for ultra-compact mini-halos.
3. Versatility: The formalism is general enough to be applied to other inflationary scenarios, including hilltop features and spectator fields, though the authors note that applying it to constant-$\eta_H$ inflation (where $\eta_H \neq 0$) requires more complex spectral analysis reserved for future work.

The authors emphasize that while the piecewise construction of the spectrum in the broad-well limit is analytically useful for the tail, it introduces artifacts near the peak of the PDF. They demonstrate that a large-$\alpha$ analytic solution provides a more accurate approximation for the peak than the piecewise construction, preserving the coherence of the spectral structure.