Vector dark matter production during inflation in the gradient-expansion formalism

This paper extends the gradient-expansion formalism to massive vector fields to model their non-linear production during inflation via inflaton couplings, demonstrating that while pure mass coupling enhances only longitudinal modes, the inclusion of kinetic coupling allows transverse modes to dominate or drives the system into a strong backreaction regime depending on the coupling function's behavior.

The Gist

Here is an explanation of the paper "Vector dark matter production during inflation in the gradient-expansion formalism," translated into everyday language with creative analogies.

The Big Picture: Hunting for the Invisible Ghost

Imagine the Universe is a giant, invisible ocean. We know there is a lot of "stuff" in it (Dark Matter) that holds galaxies together, but we can't see it, touch it, or smell it. It only interacts with normal matter through gravity, like a ghost that can only push you but never shake your hand.

Scientists have many ideas about what this ghost is. One popular idea is that it's a "Dark Photon." Think of a normal photon as a particle of light. A dark photon is like a "shadow twin" of light—it has mass and exists in a hidden sector, but it doesn't interact with our regular light.

The Problem: How do you create enough of these dark ghosts to fill the Universe?
If you try to make them using standard physics (like smashing particles together), they are too weak to be made in large numbers. The paper suggests a better way: Make them during the Big Bang's "Inflation" phase.

The Setup: The Great Cosmic Stretch

Imagine the early Universe as a balloon being blown up incredibly fast. This is Inflation.

• The Inflaton: This is the "air" blowing up the balloon. It's a field of energy that drives the expansion.
• The Dark Vector Field: This is the "dark ghost" we want to create.

Usually, if you just stretch a rubber band (the field), it snaps back to zero. But the authors show that if you tie the rubber band to the air pump (the Inflaton) in specific ways, the stretching actually creates more rubber bands out of thin air.

The Two Ways to Tie the Knot (Couplings)

The paper explores two main ways the "Inflaton" (the pump) can talk to the "Dark Vector" (the rubber band):

1. The Mass Coupling (The Heavy Anchor):
   Imagine the dark photon has a weight that changes depending on how fast the balloon is inflating. If the weight changes rapidly, it shakes the field, creating particles.

   • Result: This mostly creates Longitudinal particles. Think of these as "squishy" waves moving back and forth along the direction of travel.

2. The Kinetic Coupling (The Stretchy Rope):
   Imagine the stiffness of the rubber band itself changes as the balloon expands. If the band gets looser or tighter, it vibrates wildly.

   • Result: This creates Transverse particles. Think of these as "wiggly" waves moving side-to-side, like a guitar string.

The Big Challenge: The "Feedback Loop"

Here is where it gets tricky.

• Scenario A (Weak Production): If you make a few dark photons, they are like a few drops of water in a swimming pool. They don't change the pool. You can calculate them easily by looking at one drop at a time.
• Scenario B (Strong Production): If you make too many dark photons, they become a tsunami. They start pushing back against the balloon (the Inflaton), slowing down the expansion or changing how the balloon inflates. This is called Backreaction.

When this happens, you can't look at one drop at a time anymore because every drop is bumping into every other drop. The math becomes a messy, tangled knot that is nearly impossible to solve with standard methods.

The Solution: The "Gradient-Expansion Formalism" (GEF)

This is the paper's main innovation. The authors developed a new mathematical tool to solve this messy knot.

The Analogy: The Crowd vs. The Individual

• Old Method (Mode-by-Mode): Imagine trying to predict the weather by tracking every single raindrop individually. If there are a billion drops, you need a billion computers. If the drops start crashing into each other (non-linear), the math breaks.
• New Method (GEF): Instead of tracking drops, you track the density and flow of the rain. You ask, "How much water is in this bucket?" and "How fast is the water moving?"
  • The authors created a set of equations that track the average energy and twist of the dark field, rather than individual particles.
  • It's like switching from counting every grain of sand on a beach to measuring the shape of the dunes. It's much faster and handles the "tsunami" (strong backreaction) perfectly.

They extended this method (which was previously used for massless fields) to handle massive fields, which is much harder because massive fields have that extra "squishy" (longitudinal) movement.

What They Found (The Results)

The team ran simulations using a simple model of the early Universe (a quadratic hill) and tested different "knots" (couplings):

1. If you only use the Mass Coupling:

   • You mostly get "squishy" (longitudinal) waves.
   • If the coupling is strong, these waves grow so big they push back on the inflation, extending the life of the Big Bang's expansion by a few extra moments.

2. If you use Kinetic Coupling (and it's decreasing):

   • You get mostly "wiggly" (transverse) waves. These are the dominant energy source.

3. If you use Kinetic Coupling (and it's increasing):

   • The "squishy" (longitudinal) waves take over. They grow explosively and quickly push the system into the "strong backreaction" zone, where the dark matter fights back against the inflation.

The Takeaway

This paper is a mathematical toolkit upgrade.

• Before: Scientists could only calculate dark matter production when it was weak and quiet.
• Now: They have a method (GEF) to calculate what happens when the dark matter goes wild, creates a tsunami, and changes the history of the Universe's expansion.

They proved that by tweaking how the dark field connects to the inflation field, you can create different types of dark matter (squishy vs. wiggly) and potentially explain why the Universe has the amount of dark matter it does today. It's a step forward in understanding the "ghost" that holds our universe together.

Technical Summary

Here is a detailed technical summary of the paper "Vector dark matter production during inflation in the gradient-expansion formalism" by Lysenko, Sobol, and Vilchinskii.

1. Problem Statement

Dark matter (DM) remains one of the most significant open problems in cosmology and particle physics. While massive vector fields (dark photons) are promising candidates, their production during inflation faces specific challenges:

• Conformal Invariance: Standard massless gauge fields are conformally invariant in a Friedmann–Lemaître–Robertson–Walker (FLRW) background, preventing their production from vacuum fluctuations.
• Massive Vector Fields: While a mass term breaks conformal invariance, allowing for production, it typically only efficiently generates longitudinal modes. Generating significant transverse modes usually requires non-minimal couplings to curvature or the inflaton.
• Strong Backreaction: When vector fields are produced in large quantities, they can significantly alter the background inflationary dynamics (backreaction). Standard "mode-by-mode" Fourier analysis becomes inadequate in this nonlinear regime because the evolution of each mode depends on all others, making the system of equations computationally intractable.

Goal: To develop a robust framework for calculating the production of massive vector dark matter during inflation, specifically addressing the nonlinear dynamics and strong backreaction regimes, and to analyze the interplay between transverse and longitudinal polarization components.

2. Methodology: Gradient-Expansion Formalism (GEF)

The authors extend the Gradient-Expansion Formalism (GEF), previously used for massless gauge fields, to include massive vector fields.

• Core Concept: Instead of evolving individual Fourier modes (which becomes impossible in the nonlinear regime), the GEF works in coordinate space. It tracks the evolution of a set of bilinear functions (correlators) of the vector field, such as $\langle E^2 \rangle$, $\langle B^2 \rangle$, and $\langle E \cdot B \rangle$.
• Extension to Massive Fields:
  • The authors decompose the vector field into transverse ($\perp$) and longitudinal ($\parallel$) components.
  • They derive coupled systems of first-order differential equations for bilinear quantities for both sectors.
  • Boundary Terms: A crucial innovation is the inclusion of "boundary terms" on the right-hand side of the evolution equations. These terms act as quantum sources, pumping energy from vacuum modes as they cross the Hubble horizon, counteracting the redshift due to cosmic expansion.
• Model Setup:
  • Action: A massive Abelian vector field $A_\mu$ coupled to the inflaton $\phi$ via:
    1. Kinetic (Dilatonic) Coupling: $I_1(\phi) F_{\mu\nu}F^{\mu\nu}$
    2. Axial Coupling: $I_2(\phi) F_{\mu\nu}\tilde{F}^{\mu\nu}$
    3. Field-Dependent Mass: $m(\phi) A_\mu A^\mu$
  • Benchmark Model: A quadratic inflaton potential ($V \propto \phi^2$) with Ratra-type exponential coupling functions ($I_1, m \propto e^{\beta \phi/M_P}$).
  • Approximation: The study focuses on the "low-mass" limit ($\mu = m/H\sqrt{I_1} \ll 1$), which is typical for dark photon models, but retains the time-dependence of the mass ($M = \dot{m}/mH$).

3. Key Contributions

1. Formalism Extension: Successfully generalized the GEF to massive vector fields, deriving the complete system of equations for both transverse and longitudinal polarizations, including the necessary boundary terms derived from Whittaker and Hankel functions at horizon crossing.
2. Nonlinear Dynamics: Provided a computationally efficient method to study the strong backreaction regime, where the energy density of the vector field becomes comparable to the inflaton, a scenario where standard Fourier methods fail.
3. Polarization Hierarchy Analysis: Systematically investigated how different coupling types (mass-only vs. kinetic+mass) dictate the relative abundance of transverse vs. longitudinal modes.
4. Validation: Validated the GEF against the standard mode-by-mode Fourier solution in the weak backreaction regime, demonstrating high accuracy (errors $\sim 1\%$ for dominant components).

4. Key Results

The numerical analysis of the benchmark model yielded several critical findings:

• Pure Mass Coupling ($I_1 = \text{const}$):

  • Only longitudinal modes are significantly enhanced.
  • Transverse modes remain at vacuum levels.
  • For strong couplings, the longitudinal energy density grows large enough to cause strong backreaction, prolonging inflation by a few e-folds and inducing damped oscillations in the energy density.

• Kinetic + Mass Coupling ($I_1(\phi)$ varying):

  • Decreasing Coupling Function ($\beta > 0$ in $e^{\beta \phi}$): As the coupling decreases during inflation, transverse modes dominate the energy density. The longitudinal component remains subdominant.
  • Increasing Coupling Function ($\beta < 0$): As the coupling increases, the longitudinal component grows rapidly and becomes the dominant contributor to the energy density. This scenario drives the system quickly into the strong backreaction regime.
  • Hierarchy Flip: The sign of the coupling constant $\beta$ effectively flips the hierarchy of energy densities. For negative $\beta$, the longitudinal magnetic-like energy ($\rho_F$) dominates over the transverse electric energy ($\rho_E$).

• Backreaction Effects:

  • In the strong backreaction regime, the vector field energy density stops growing and exhibits damped oscillations.
  • The energy transfer from the inflaton to the vector field maintains the gauge-field energy density almost constant, creating a "vacuum-like" state that extends the duration of inflation.

5. Significance and Conclusion

• Theoretical Advancement: This work resolves the technical bottleneck of studying massive vector DM production in the nonlinear regime. The GEF is shown to be a powerful, faster, and computationally cheaper alternative to Fourier mode summation when backreaction is strong.
• Phenomenological Insight: The study clarifies that the nature of the coupling (kinetic vs. mass) is the deciding factor in whether transverse or longitudinal modes dominate the dark matter relic abundance. This has direct implications for observational signatures, as transverse and longitudinal modes interact differently with standard matter and gravitational waves.
• Future Outlook: The authors note that while kinetic mixing with the Standard Model photon was neglected here to isolate the production mechanism, the formalism is ready to incorporate it. This paves the way for future studies connecting inflationary production to experimental detection prospects for dark photons.

In summary, the paper establishes that mass coupling alone favors longitudinal modes, while kinetic coupling can enhance transverse modes, and the interplay between these couplings determines whether the system remains in a perturbative regime or enters a strong backreaction phase that alters the inflationary history.