Analytic Approximations for Fermionic Preheating

This paper presents analytic approximations for non-perturbative fermion production during $\lambda\phi^{4}$ inflation, characterizing how the momentum spectrum and total number density depend on the coupling parameter $q$ and deriving mass constraints for fermionic dark matter.

The Gist

Imagine the early universe as a giant, vibrating drum. Right after the Big Bang, a field called the inflaton (the "drumhead") was oscillating wildly. Usually, we think of this drum just shaking and eventually settling down, releasing energy like a gentle breeze to create the particles we see today (like atoms and light). This gentle process is called "reheating."

But this paper explores a much more violent, explosive version of that process called Preheating. Instead of a gentle breeze, the drum shakes so hard it creates a shockwave that instantly blasts out new particles.

Here is the breakdown of what the authors found, using simple analogies:

1. The Players: The Drum and the Marbles

• The Inflaton (The Drum): A scalar field that oscillates back and forth. In this paper, they imagine it's vibrating with a specific, mathematically simple rhythm (like a perfect sine wave).
• The Fermions (The Marbles): These are the particles being created (like electrons or dark matter candidates). They are "fermions," which means they follow a strict rule: no two marbles can sit in the exact same seat (this is the Pauli Exclusion Principle).
• The Coupling ($q$): This is a knob that controls how strongly the drum shakes the marbles.
  • Small $q$: The drum shakes the marbles gently.
  • Large $q$: The drum shakes the marbles violently.

2. The Two Ways Marbles Get Created

The authors discovered that how the marbles fill up depends entirely on how hard the drum shakes (the value of $q$).

Scenario A: The "Resonance Peaks" (Small $q$)

Imagine the drum is shaking at a specific rhythm. If you push a child on a swing at just the right time, the swing goes higher and higher. This is resonance.

• What happens: The drum doesn't push all the marbles equally. Instead, it only pushes marbles that are moving at very specific speeds (momentums) that match the drum's rhythm.
• The Result: The marbles pile up in distinct, sharp "spikes" or resonance peaks. It's like the drum is only filling up specific seats in a stadium, leaving the rest empty.
• The Discovery: For weak shaking, the authors found a simple math formula to predict exactly which seats will be filled. They realized this is just a consequence of energy conservation: the drum gives energy to the marbles in chunks, and only marbles that can absorb that exact chunk get created.

Scenario B: The "Half-Filled Sphere" (Large $q$)

Now, imagine the drum is shaking so violently that it's a chaotic mess.

• What happens: The shaking is so fast and strong that the "no two marbles in the same seat" rule becomes the main driver. The marbles rush to fill up the lowest available seats as fast as possible.
• The Result: Instead of sharp spikes, you get a big, smooth "ball" of marbles. It looks like a half-filled sphere of marbles sitting at the bottom of a bowl. The "bulk" of the marbles are here, not in the spikes.
• The Discovery: When the shaking is this violent, the total number of marbles created follows a different mathematical rule (a power law) than when the shaking is gentle.

3. The "Magic Formula"

The authors did something clever. Usually, to figure out exactly how many marbles are created, you have to run a super-computer simulation that takes a long time.

• The Breakthrough: They derived a simple, semi-analytic formula (a "rule of thumb") that predicts exactly where those resonance spikes will happen for any strength of shaking.
• Why it matters: It's like having a map that tells you exactly where the treasure chests are buried, without needing to dig up the whole island first.

4. The Dark Matter Connection

Why do we care about these vibrating drums and marbles?

• Dark Matter: We know there is invisible "Dark Matter" holding galaxies together, but we don't know what it is made of.
• The Candidate: The authors suggest that these "marbles" (fermions) created during this violent preheating phase could be the Dark Matter.
• The Weight Limit: Because of the "no two marbles in the same seat" rule, if the marbles are too light, they would move too fast and fly apart, preventing galaxies from forming.
• The Conclusion: By calculating how many marbles were created, they set a lower weight limit on Dark Matter. If these particles are the Dark Matter, they must be at least a few thousand times heavier than an electron (a few keV). If they were lighter, the universe would look very different than it does today.

Summary

This paper is a study of how the universe's first "drum" shook to create the first particles.

1. Gentle shaking creates particles in specific, predictable "spikes" (Resonance).
2. Violent shaking creates a smooth "ball" of particles (Bulk).
3. The authors found a simple way to predict the spikes without complex math.
4. If these particles are the Dark Matter, they must be heavy enough to hold galaxies together.

It turns the complex, chaotic birth of the universe into a story about a drum, a swing, and a very strict rule about sitting in seats.

Technical Summary

1. Problem Statement

The paper addresses the non-perturbative production of fermions during the preheating phase of the early universe, specifically following $\lambda\phi^4$ inflation. While preheating is well-studied for bosons, fermionic production is more complex due to the Pauli exclusion principle and the specific dynamics of the Dirac equation in a time-dependent background.

Key challenges identified include:

• Momentum Spectrum Complexity: The momentum distribution of produced fermions is non-degenerate, featuring a "bulk region" (low momentum) and discrete "resonance peaks" (higher momentum).
• Dependence on Coupling ($q$): The relative contribution of the bulk region versus resonance peaks to the total number density depends heavily on the coupling parameter $q = h^2/\lambda$ (where $h$ is the Yukawa coupling and $\lambda$ is the inflaton quartic coupling).
• Lack of Analytic Predictions: Previous studies relied heavily on numerical solutions of the mode equation. There was a need for analytic approximations to predict resonance peak locations and total number densities across different regimes of $q$ without extensive numerical computation.
• Dark Matter Constraints: Understanding the precise momentum distribution is crucial for placing cosmological bounds on fermionic dark matter candidates produced via this mechanism.

2. Methodology

The authors employ a semi-analytic approach combining quantum field theory in curved spacetime with numerical validation.

• Model Setup: They assume a $\lambda\phi^4$ potential for the inflaton and a Yukawa coupling $h\phi\bar{\psi}\psi$ between the inflaton and fermions. The fermion bare mass $m_\psi$ is assumed negligible compared to the oscillation scale.
• Mode Equation: The fermion production is governed by a mode equation derived from the Dirac equation in a conformal background. The equation takes the form of an oscillator with a time-dependent frequency $\Omega_\kappa(\tau) = \sqrt{\kappa^2 + q f^2(\tau)}$, where $\kappa$ is the dimensionless comoving momentum and $f(\tau)$ is the Jacobi Elliptic Cosine function describing inflaton oscillations.
• Number Density Calculation:
  • They utilize Bogoliubov transformations to diagonalize the Hamiltonian and calculate the number density per mode, $n_\kappa(\tau) = |\beta|^2$.
  • For late times ($\tau \gg T$, where $T$ is the oscillation period), they approximate the rapidly oscillating $n_\kappa(\tau)$ with a smooth envelope function $\bar{n}_\kappa \approx F_\kappa/2$, where $F_\kappa$ is derived from the transfer matrix of the mode equation over one period.
• Regime Analysis: The total comoving number density is obtained by integrating the momentum distribution ($\int d\kappa \kappa^2 F_\kappa$). The authors divide the analysis into three regimes:
  1. Small $q$ ($q \lesssim 0.01$): Where resonance peaks are narrow and distinct.
  2. Large $q$ ($q \gtrsim 10$): Where the bulk region dominates and merges with the lowest resonance peak.
  3. Intermediate $q$: Where neither simple power law applies well.

3. Key Contributions

A. Semi-Analytic Prediction of Resonance Peaks

The authors derive a general semi-analytic relation to predict the momentum values ($\kappa_l$) of resonance peaks for any coupling $q$, without solving the differential equation numerically.

• Derivation: Based on energy conservation in the perturbative process $n\phi \to \Psi\bar{\Psi}$, they show that resonance occurs when the time-averaged frequency satisfies:
  $$\Omega_\kappa(q) = \frac{\pi}{T}(2l + 1), \quad l = 0, 1, 2, \dots$$
• Significance: This generalizes previous findings (which were limited to $q \ll 1$) to all $q$. It explains that only odd multiples of the fundamental frequency contribute due to the symmetry of the potential ($\phi \to -\phi$).

B. Analytic Power-Law Approximations for Number Density

The paper provides novel analytic approximations for the total number density of fermions as a function of $q$:

• Small $q$ Regime ($q \lesssim 0.01$):
  • The total number density is dominated by the first two resonance peaks ($l=0$ and $l=1$), contributing $\sim95\%$ of the total. The bulk region contributes $<5\%$.
  • Result: The total number density scales as $n \propto q^{1/2}$.
  • Specific fit: $I_S(q) \approx 0.38 \times q^{1/2}$.
• Large $q$ Regime ($q \gtrsim 10$):
  • The bulk region (driven by non-adiabaticity where effective mass crosses zero) dominates, merging with the lowest resonance peak.
  • Result: The total number density scales as $n \propto q^{3/4}$.
  • Specific fit: $I_L(q) \approx 0.13 \times q^{3/4}$.
• Intermediate Regime: No simple analytic fit is found, but numerical results show a transition between the two power laws.

C. Dark Matter Mass Bounds

The authors apply their results to constrain the mass of fermionic dark matter ($m_\psi$) produced via preheating, assuming they constitute all dark matter.

• They adapt bounds from structure formation (specifically the requirement that fermions become non-relativistic early enough to allow structure formation).
• Small $q$: The bound is strengthened. For $q=10^{-5}$, $m_\psi \gtrsim 10$ keV (compared to the standard $\sim 2$ keV bound).
• Large $q$: The momentum spectrum resembles a half-filled Fermi sphere, and the bound remains $m_\psi \gtrsim 2$ keV.

4. Key Results

1. Resonance Locations: The relation $\Omega_\kappa(q) = \frac{\pi}{T}(2l+1)$ accurately predicts peak locations across all $q$ values, verified numerically.
2. Dominance Shift:
   • For $q \lesssim 0.01$, the "bulk" (low momentum) is negligible; production is driven by discrete resonances.
   • For $q \gtrsim 10$, the bulk region (non-adiabatic production) dominates the total yield.
3. Scaling Laws: The discovery of distinct power-law scalings ($q^{1/2}$ vs. $q^{3/4}$) provides a robust tool for estimating fermion yields in cosmological models without full numerical simulation.
4. Generalizability: The authors argue these scaling laws and the resonance condition likely hold for other symmetric potentials (e.g., $m^2\phi^2$ inflation, hybrid inflation, axion oscillations), as verified for $m^2\phi^2$.

5. Significance

• Theoretical Efficiency: The derived analytic approximations significantly reduce the computational cost of modeling fermionic preheating in cosmological simulations, allowing for rapid estimation of dark matter abundances.
• Physical Insight: The work clarifies the transition between perturbative-like resonance production (small $q$) and non-adiabatic bulk production (large $q$), resolving ambiguities in previous literature regarding the momentum distribution.
• Phenomenological Impact: By providing tighter and more accurate mass bounds for fermionic dark matter, the paper refines the viability of preheating as a dark matter production mechanism. It highlights that the coupling strength $q$ is a critical parameter in determining whether preheating-produced fermions can serve as viable dark matter candidates.
• Dark Matter Candidates: It reinforces the potential of non-thermal fermions as dark matter, offering specific mass ranges ($2\text{--}10$ keV) depending on the coupling regime.