Decay of a scalar condensate in two different approaches

This paper establishes the equivalence between the parametric resonance and coherent state $S$-matrix approaches to scalar condensate decay by modifying the latter to eliminate unwanted diagrams and explicitly demonstrating their agreement at lower orders in a double expansion.

The Gist

The Big Picture: The "Cosmic Melting Ice Cube"

Imagine the very early Universe as a giant, super-cold block of ice. In physics, this "ice" is called a scalar condensate. It's not just one particle; it's a massive, coherent wave of energy, like a giant ocean wave that hasn't broken yet.

As the Universe expands and cools, this giant wave starts to wobble and eventually breaks apart. When it breaks, it releases tiny droplets of water (new particles) and heats up the surrounding air. In cosmology, this process is called reheating, and it's what allowed the Universe to become the hot, particle-filled place we live in today.

The big question this paper asks is: How fast does this giant wave break apart?

The authors found that physicists have been trying to answer this question using two completely different toolkits. One toolkit is like watching a wave crash (mathematical waves), and the other is like counting the droplets using a particle counting machine (Feynman diagrams).

The paper's main achievement is proving that both toolkits give the exact same answer, even though they look like they are speaking different languages.

---

The Two Approaches: Two Ways to Watch the Wave

1. The "Parametric Resonance" Approach (The Wave Watcher)

The Analogy: Imagine you are pushing a child on a swing. If you push at just the right rhythm, the swing goes higher and higher. This is "resonance."

In this approach, physicists look at the giant wave (the condensate) as a background rhythm. They ask: "If I have a tiny ripple (a daughter particle) sitting in this rhythm, will it grow?"

• They solve a complex equation (the Mathieu equation) to see if the ripple grows exponentially.
• If it grows, that growth rate tells them how fast the energy is leaking out of the giant wave.
• The Vibe: This is like watching a video of the ocean and measuring how fast the waves are getting bigger. It's very visual and relies on solving differential equations.

2. The "Feynman-Diagrammatic" Approach (The Particle Counter)

The Analogy: Imagine you are a detective trying to figure out how a bank vault was robbed. Instead of watching the whole vault, you look at every single possible way a thief could have entered. You draw a map of every path (a diagram) the thief could take.

In this approach, physicists treat the giant wave not as a smooth wave, but as a giant crowd of particles all moving together. They use Feynman diagrams (stick-figure drawings of particle interactions) to calculate the probability of the wave breaking apart.

• They draw circles and lines representing particles bumping into each other.
• They sum up all the possible "paths" the energy could take to escape.
• The Vibe: This is like counting every possible way a crowd of people could exit a stadium through different doors. It's very mechanical and relies on counting and summing up probabilities.

The Problem: They Look Totally Different

For a long time, it wasn't clear if these two methods were actually measuring the same thing.

• The Wave Watcher says: "The energy leaks out because the rhythm makes the ripples grow."
• The Particle Counter says: "The energy leaks out because the particles are bumping into each other in specific patterns."

It's like trying to measure the speed of a car. One person measures the time it takes to pass a tree (Wave Watcher), and the other person counts the number of engine rotations (Particle Counter). They use totally different math, so it's hard to believe they will get the exact same speed.

The Solution: The "Double Expansion" Bridge

The authors of this paper built a bridge between these two worlds. They realized that to compare them fairly, they had to zoom in very closely.

They introduced a "Double Expansion," which is like looking at the problem through two different zoom lenses simultaneously:

1. Zoom Lens A (Amplitude): How big is the wave's wobble? (Small wobble vs. big wobble).
2. Zoom Lens B (Velocity): How fast are the new particles moving? (Slow vs. fast).

By calculating the answer step-by-step for small wobbles and slow particles, they could compare the results of the Wave Watcher and the Particle Counter side-by-side.

The Result:
When they did the math for the first few steps (the "lower orders"), the numbers matched perfectly.

• The "Wave Watcher" calculated a specific growth rate.
• The "Particle Counter" calculated the exact same rate by summing up all the particle diagrams.

This proves that the two approaches are just two different ways of describing the exact same physical reality. One is solving a wave equation; the other is doing a particle count. They are two sides of the same coin.

Why Does This Matter?

1. Validation: It gives physicists confidence. If two totally different math methods give the same answer, we know the answer is likely correct.
2. Flexibility: Sometimes the "Wave Watcher" method is easier to use; sometimes the "Particle Counter" method is easier. Now that we know they are equivalent, physicists can switch between them depending on which tool is better for the specific problem they are solving.
3. Cosmology: Understanding exactly how the early Universe heated up helps us understand why the Universe is the way it is today. If the math is wrong, our models of the Big Bang are wrong.

Summary

Think of the early Universe as a giant, vibrating drum.

• Method A listens to the sound of the drum vibrating to see how fast it loses energy.
• Method B counts the tiny air molecules hitting the drum to see how fast they carry energy away.

This paper says: "Don't worry, both methods tell you the drum is losing energy at the exact same speed." They have proven that the physics of the early Universe is consistent, no matter which mathematical lens you use to look at it.

Technical Summary

1. Problem Statement

The paper addresses the theoretical description of the decay (energy dissipation) of a coherently oscillating scalar condensate (e.g., an inflaton field) into daughter particles. This process is crucial in cosmology for understanding reheating and preheating in the early Universe.

Two distinct quantum field theoretic (QFT) approaches exist in the literature to calculate this decay rate:

1. Parametric-Resonance Approach: Solves the equation of motion for the daughter particle's mode function in the presence of the time-dependent background field. It identifies exponentially growing modes (instability bands) and calculates the growth rate.
2. Feynman-Diagrammatic Approach: Treats the scalar condensate as a coherent state and computes the decay rate using the $S$-matrix and perturbation theory (Feynman diagrams).

The Core Issue: While these methods are known to yield similar results in specific limits (e.g., narrow resonance where daughter mass $\approx$ half the parent mass), their equivalence was not rigorously established for general amplitudes and velocities. Furthermore, the diagrammatic approach in previous literature (Ref. [7]) included unwanted diagrams and lacked a clear manifest connection to the vacuum-to-vacuum transition amplitude. The authors aim to:

• Refine the diagrammatic approach to isolate relevant physical processes.
• Demonstrate the mathematical equivalence of both approaches.
• Explicitly compute the decay rate using a double expansion in the oscillation amplitude ($\theta$) and the daughter particle velocity ($\beta$).

2. Methodology

The authors consider a toy model with a real scalar field $\phi$ (condensate) and a real scalar field $\chi$ (daughter), interacting via a term $\frac{1}{2}\mu\phi\chi^2$.

A. Refined Feynman-Diagrammatic Approach

The authors modify the standard coherent state $S$-matrix formalism:

• Coherent State: The condensate is described as a coherent state $|\phi\rangle$, allowing the replacement of field operators $\hat{\phi}$ with classical background values in normal-ordered products.
• Vacuum-to-Vacuum Transition: The decay rate $\Gamma$ is extracted from the imaginary part of the effective action, specifically the vacuum-to-vacuum transition amplitude $\langle \phi | \hat{S} | \phi \rangle \propto e^{-\Gamma L^3 T / 2}$.
• Pinch Technique: To handle apparent divergences arising when multiple propagators in a loop share the same momentum (a common issue in bubble diagrams with coherent backgrounds), the authors employ a pinch technique. They replace identical propagator masses $m_\chi$ with distinct parameters $\xi_i$, perform the calculation, and take the limit $\xi_i \to m_\chi$ after applying cutting rules. This regularizes the divergences and isolates the finite physical contribution.
• Diagram Selection: They focus on "bubble diagrams" where the background field $\phi$ is inserted as external lines. Only diagrams that can be cut into two pieces (contributing to the imaginary part) are retained.

B. Parametric-Resonance Approach

• Mathieu Equation: The equation of motion for the daughter particle $\chi$ in the background $\phi(t) = A_\phi \cos(m_\phi t)$ is reduced to a Mathieu equation.
• Floquet Theory: The solution is sought in the form $u(z) = e^{\mu z}P(z)$, where $\mu$ is the Floquet exponent. The real part of $\mu$ corresponds to the growth rate $\lambda$.
• Tridiagonal Matrix: The authors solve the recurrence relations for the Fourier coefficients of the periodic function $P(z)$ by constructing a tridiagonal matrix. They derive analytical expressions for the matrix elements to higher orders in the amplitude parameter $\theta$.

C. Double Expansion Framework

To compare the two methods, the authors expand the decay rate $\Gamma$ in two small parameters:

1. $\theta$ (Amplitude): Proportional to the coupling and oscillation amplitude ($\theta \propto \mu A_\phi / m_\phi^2$).
2. $\beta$ (Velocity): The velocity of the daughter particle in the absence of the background ($\beta \propto \sqrt{1 - 4m_\chi^2/n_\phi^2 m_\phi^2}$).

They define the expansion as $\Gamma = \sum \Gamma^{(p,q)} \theta^{2p} \beta^{1-2q}$, where $p$ represents the order in amplitude and $q$ represents the order in velocity.

3. Key Contributions

1. Formal Equivalence: The paper establishes that the parametric-resonance approach (solving the Schrödinger-like equation for mode functions) and the Feynman-diagrammatic approach (path integral/S-matrix) compute the exact same physical quantity: the vacuum-to-vacuum transition amplitude of the daughter field in the presence of the background.
2. Refinement of Diagrammatic Method: The authors clarify how to compute the $S$-matrix for a coherent state without including unphysical disconnected diagrams, explicitly linking the imaginary part of the effective action to the decay rate.
3. Resolution of Divergences: They provide a rigorous treatment of the "pinch" singularities in bubble diagrams using the $\xi$-parameter regularization, showing how divergent terms cancel out to yield finite physical results.
4. Explicit Higher-Order Calculations: The paper provides explicit analytical results for the decay rate up to Next-to-Leading Order (NLO) and Next-to-Next-to-Leading Order (NNLO) in the double expansion for processes involving $n_\phi = 1$ and $n_\phi = 2$ parent particles converting to two daughter particles ($n_\phi \phi \to 2\chi$).

4. Results

The authors computed the decay rate $\Gamma$ for $n_\phi = 1$ and $n_\phi = 2$ using both methods and found perfect agreement:

• Parametric-Resonance Results:

  • Derived the growth factor $\lambda$ by solving the determinant of the tridiagonal matrix.
  • Obtained $\Gamma$ as an integral over the instability band width, expanded in powers of $\theta$ and $\beta$.
  • Found that for a fixed order in $\theta$, the method yields a specific range of $\beta$ orders (diagonal lines in the $(p,q)$ plane).

• Feynman-Diagrammatic Results:

  • Evaluated bubble diagrams ($p=1, 2, 3$) using the pinch technique and cutting rules.
  • Calculated the imaginary parts of diagrams involving $2p$ insertions of the background field.
  • Found that for a fixed order in $\theta$ (fixed $p$), the method yields the full series in $\beta$ (vertical lines in the $(p,q)$ plane).

• Verification of Equality:

  • $n_\phi = 1$: The Leading Order (LO), NLO, and NNLO terms derived from the Mathieu equation matched exactly with the sums of the corresponding Feynman diagrams (e.g., Eq. 5.1 vs. Eq. 4.8, 4.13, 4.16).
  • $n_\phi = 2$: Similarly, the LO and NLO results matched (Eq. 5.4, 5.5 vs. Eq. 4.14, 4.17).
  • The paper demonstrates that the "narrow resonance" limit in the parametric approach corresponds to the "small amplitude" (perturbative) limit in the diagrammatic approach.

5. Significance

• Unification of Frameworks: This work bridges the gap between non-perturbative mode-function methods (often used in preheating simulations) and standard perturbative QFT methods. It validates the use of Feynman diagrams for condensate decay, provided the coherent state formalism and pinch technique are correctly applied.
• Computational Efficiency: It clarifies which method is more efficient for specific regimes. The diagrammatic approach is powerful for obtaining all velocity orders at a fixed amplitude, while the resonance approach is natural for finding the instability bands.
• Cosmological Implications: By establishing the equivalence and providing higher-order corrections, the paper improves the precision of reheating temperature calculations in inflationary models. It confirms that energy dissipation via parametric resonance can be consistently described within standard QFT perturbation theory in the narrow-resonance regime.
• Future Directions: The authors note that while the equivalence holds for small $T$ (oscillation periods) and narrow resonance, the behavior at large $T$ (where back-reaction and broad resonance dominate) requires non-equilibrium QFT or lattice simulations, which are beyond the scope of this perturbative study.

In summary, the paper successfully demonstrates that two seemingly different approaches to scalar condensate decay are mathematically equivalent, providing a robust theoretical foundation for calculating reheating rates in the early Universe.