Fermionic Bubble Loop in Cosmological Collider Revisited: Exact signals from spectral and Mellin-Barnes methods

This paper presents an exact analytical calculation of fermionic bubble loop contributions to cosmological collider signals using parallel spectral and Mellin-Barnes methods, revealing that Yukawa interactions with the inflaton yield a vanishing bispectrum due to a field redefinition of tree-level counterparts.

The Gist

Imagine the universe as a giant, expanding balloon. In the very first moments after the Big Bang, this balloon was inflating so fast and was so hot that it acted like a massive particle accelerator, far more powerful than anything we could build on Earth. Physicists call this the "Cosmological Collider."

Usually, when we look at the leftover radiation from the Big Bang (the Cosmic Microwave Background), we see a smooth, boring pattern. But if heavy, exotic particles existed back then, they would have left a tiny, rhythmic "fingerprint" or "echo" in that pattern. Finding these echoes is like listening for a specific instrument in a noisy orchestra to figure out what kind of band was playing.

For a long time, scientists could easily predict the fingerprints of heavy particles that act like "balls" (scalars) or "spinning tops" (vectors). But they struggled with particles that act like "spinning electrons" (fermions). Why? Because calculating the behavior of these fermions involves incredibly complex math, specifically "loop diagrams."

Think of a loop diagram like a detour. Instead of a particle going straight from point A to point B, it briefly splits into two particles that travel in a circle before rejoining. Calculating this circle is mathematically messy and usually requires making rough guesses (approximations) because the equations are too hard to solve exactly.

What this paper does:
The authors, a team of physicists, decided to stop guessing. They used two completely different, high-powered mathematical "flashlights" to shine on the fermion loop problem and solve it exactly for the first time.

1. The "Spectral Decomposition" Method: Imagine you have a complex, tangled knot of string (the fermion loop). This method says, "Let's untangle it by realizing this knot is actually just a stack of many simple, straight strings (tree-level diagrams) of different lengths." They broke the complex loop down into an infinite sum of simpler, known pieces.
2. The "Mellin-Barnes" Method: This is like translating the problem into a different language (a mathematical space called "Mellin space"). In this new language, the complicated curves and waves turn into simple building blocks (Gamma functions). Once translated, the math becomes easy to solve, and then they translate the answer back.

The Big Surprise:
After doing all this heavy lifting and getting two different answers that perfectly matched each other, they tested their new formula on a very common scenario: Yukawa coupling.

In physics, Yukawa coupling is like a standard handshake between a heavy particle and the field that drove the Big Bang (the inflaton). It's the most basic, expected way these particles interact.

The authors expected to find a clear, rhythmic echo (a signal) in the data. Instead, they found nothing. The signal vanished completely.

Why did it vanish?
The paper explains this using a clever trick. Because the fermion loop is mathematically equivalent to a stack of simpler tree-level diagrams, they looked at those simpler diagrams. They found that for this specific type of interaction, the "echo" from one part of the stack perfectly cancels out the "echo" from another part. It's like two people shouting the same note but in opposite phases; the sound waves cancel each other out, leaving silence.

They also showed that this silence isn't a mistake; it's a fundamental property of the universe's geometry at that time. You can think of it as a "field redefinition"—a mathematical reshuffling of how we describe the particles—that proves the signal was never there to begin with.

The Takeaway:

• The Problem: Fermion loops were too hard to calculate exactly, so previous studies used approximations.
• The Solution: The authors solved the problem exactly using two different advanced math techniques that confirmed each other.
• The Result: When they applied their exact math to the most common type of interaction (Yukawa coupling), the predicted signal disappeared entirely.
• The Lesson: Previous studies that claimed to see these signals using approximations might have been seeing "ghosts" (artifacts of the math) rather than real physics. If you want to find fermion echoes in the universe, you can't look for them in this specific, simple setup; you'll need to look for more complex interactions or different conditions.

In short, the paper is a masterclass in doing the hard math correctly, only to discover that the universe is quieter than we thought in this specific scenario.

Technical Summary

Problem Statement
Fermionic degrees of freedom are fundamental to many Beyond the Standard Model (BSM) phenomenological scenarios and are essential ingredients in cosmological collider (CC) physics. However, their observational signatures have remained largely unexplored due to the computational difficulty of evaluating loop diagrams involving massive spinning fields. Traditional methods, such as the Schwinger-Keldysh formalism, require evaluating complicated nested time integrals over products of mode functions (typically Hankel or Whittaker functions), which often lack analytical solutions, especially when symmetries are broken. While significant theoretical progress has been made in computing loop-level correlators for scalar and integer-spin fields, fermionic fields have not been treated with the same level of analytical rigor. Furthermore, previous studies of fermionic signatures often relied on approximations, such as late-time expansions, whose precision was unclear. This work addresses the gap by providing an exact, analytical computation of cosmological collider signals arising from fermionic bubble loops with arbitrary couplings.

Methodology
The authors develop two parallel, independent analytical methods to compute the fermionic bubble loop contributions to cosmological correlators. Both methods are applied to a general seed integral involving two internal fermion propagators with potentially different masses ($m_1, m_2$) and arbitrary couplings to the inflaton.

1. Spectral Decomposition Method:

   • This approach is rooted in the Källén-Lehmann representation in de Sitter (dS) spacetime.
   • The authors utilize a mathematical identity relating the product of two hypergeometric functions to an infinite series of a single hypergeometric function.
   • This allows the fermionic bubble (a product of propagators) to be recast as an infinite sum of tree-level exchange signals with varying effective masses.
   • The method establishes a direct relationship between the fermionic bubble and the scalar bubble via time-differential operators, allowing results to be derived from known scalar spectral densities.
   • The final result is expressed as a spectral integral over a density function $\rho(\nu)$, which is then evaluated to yield non-local and local signal components.

2. Mellin-Barnes (MB) Transformation Method:

   • This method employs a partial Mellin-Barnes representation for the massive fermionic propagators (Hankel functions), converting them into products of Gamma functions.
   • This transformation renders both the momentum and time integrals essentially trivial.
   • The result is reconstructed by evaluating contour integrals and picking up residues from two distinct families of poles:
     • Spectrum poles: Originating from the fermionic mass parameters, these contribute to the oscillatory cosmological collider signals.
     • UV poles: Originating from the loop integration structure, these contribute to both signal and background terms.
   • The extraction of local signals requires the application of Barnes' lemma to handle nested sums.

Key Contributions and Results

• Exact Analytical Expressions: The paper provides the first exact analytical expressions for cosmological collider signals from fermionic bubble loops with arbitrary couplings, capturing the full non-analytic structure (both local and non-local) without relying on late-time approximations.
• Methodological Consistency: The two distinct methods (Spectral and MB) yield totally different mathematical expressions (a single spectral series vs. a four-fold power series). The authors perform detailed numerical checks and kinematic limit analyses (single and hierarchical squeezed limits) to verify that both methods produce identical results, providing a non-trivial cross-check.
• Spectral Density for Fermions: The work derives the explicit spectral density $\rho_{fermion}(\nu)$ for fermionic bubbles, showing it involves four Gamma functions related to the sum and difference of the fermion masses ($M = m_1 + m_2$ and $\delta = m_1 - m_2$).
• Vanishing Yukawa Signals: A central phenomenological result is the finding that the cosmological collider signals in the bispectrum (three-point function) generated by fermionic bubble loops with Yukawa couplings vanish identically.
  • This result holds for the bubble topology with pure de Sitter mode functions.
  • The vanishing is traced back to the spectral decomposition: the bubble contribution maps to a series of tree-level diagrams with quadratic mixing.
  • Through a field redefinition argument, the authors demonstrate that each of these tree-level contributions vanishes individually, leading to the vanishing of the full loop signal.
  • This contradicts previous results based on approximate methods, suggesting those earlier findings may be unreliable.
• Background Contributions: The paper also derives the background contributions (meromorphic parts) using the MB method, noting that these terms are UV divergent and require renormalization, unlike the finite signal parts.

Significance and Claims
The authors claim that this work provides a necessary and timely revisit of fermionic correlators using modern analytical tools. By establishing exact results, they clarify the behavior of fermionic loops, which are otherwise difficult to compute.

The most significant claim is the vanishing of CC signals for Yukawa-type couplings in the bubble topology. The authors emphasize that this is a specific result for the bubble diagram with pure de Sitter mode functions. They suggest that previous non-zero results obtained via approximation methods may be artifacts of those approximations.

The paper modestly outlines the limitations and future directions:

• The vanishing result relies on dS isometry; introducing a chemical potential (which breaks this symmetry) or considering different topologies (e.g., triangle loops) could yield non-vanishing signals.
• The analysis is currently restricted to spin-1/2 fermions; extending this to spin-3/2 fields (gravitinos) is a noted future direction.
• The background contributions, while derived, require further investigation regarding renormalization and alternative calculation methods.

In summary, the paper establishes a rigorous framework for computing fermionic loop signals, demonstrates the equivalence of spectral and MB techniques in this context, and reveals a surprising cancellation mechanism for Yukawa interactions that fundamentally alters the expected phenomenological landscape for this specific class of models.