⚛️ high-energy theory

Massive de Sitter Correlators as Finite Mellin-Barnes Integrals

This paper utilizes the Mellin-Barnes representation of de Sitter propagators to express tree-level and one-loop massive particle exchange diagrams as finite, convergent integrals, enabling the direct computation of cosmological collider signatures like the triple-exchange bispectrum and complete one-loop triangle without requiring analytic continuation.

Original authors: Amara McCune

Published 2026-10-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Amara McCune

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the earliest moments of our universe, a period known as inflation, space expanded faster than the speed of light. During this rapid growth, the cosmos was not empty; it was filled with heavy particles that existed at energies far beyond what we can create in any laboratory on Earth. These particles left a faint, rhythmic imprint on the distribution of matter we see today. Much like a tuning fork that rings when struck, these heavy particles caused the early universe to vibrate with specific frequencies. By listening to these vibrations in the cosmic microwave background and the large-scale structure of galaxies, scientists hope to identify the mass and spin of particles that have been invisible for billions of years. This field of study is often called cosmological collider physics, treating the entire early universe as a giant particle detector.

The challenge has always been that calculating the signals from these heavy particles is incredibly difficult. When physicists try to map out how these particles interacted, the mathematics becomes so complex that it often breaks down. Previous methods could only handle simple interactions or required approximations that failed when the particles were very heavy or when multiple particles exchanged energy in a single event. For more complicated scenarios, such as three particles exchanging energy or particles looping back on themselves, the existing mathematical tools could not provide a clear answer at the physical conditions that actually occurred in the universe. The calculations either diverged into infinity or required moving into a mathematical realm that did not correspond to reality, making it impossible to compare the theory directly with observational data.

A new approach presented by Amara McCune changes this landscape by introducing a unified way to calculate these interactions. The researcher developed a method that translates the complex time-dependent evolution of these particles into a single, finite integral. This technique uses a specific mathematical representation that keeps the time ordering of events intact, allowing the calculation to remain stable and convergent even at the exact physical conditions of the early universe. Instead of getting stuck in mathematical dead ends, this method allows the researcher to compute the signals for multiple particle exchanges and even for one-loop diagrams, where particles travel in a closed circle, without needing to force the numbers into a different mathematical shape.

The results of this work are immediate and practical. The researcher successfully calculated the signal for a triple-particle exchange, a scenario that had previously resisted a convergent solution. This calculation revealed that the oscillatory signals from these heavy particles do not vanish exponentially as some theories had suggested; instead, they remain strong enough to be detected. The method also provided the first complete analytical calculation for a one-loop triangle diagram involving particles of arbitrary mass. In terms of speed, the new method is remarkably efficient, taking only a fraction of a second to compute a specific configuration, whereas previous techniques required minutes or even hours. This speed makes it possible to generate the detailed templates needed to search for these signals in current and future astronomical data.

The study confirms that the universe acts as a sensitive detector for heavy particles, even when they are produced in complex chains of interactions. By showing that these signals are not suppressed by the extreme conditions of the early universe, the work suggests that we may be able to detect a wider range of particles than previously thought. The method applies to any mass and any spin, offering a robust tool for future searches. It bridges the gap between the theoretical predictions of heavy particle physics and the actual data collected by telescopes, turning the abstract mathematics of the early universe into a concrete guide for discovery. The findings indicate that the rhythmic patterns left by these heavy particles are within reach, waiting to be identified in the vast data sets of modern cosmology.

Technical Summary: Massive de Sitter Correlators as Finite Mellin–Barnes Integrals

Problem Statement
Massive fields present during inflation imprint oscillatory signatures on primordial correlation functions, serving as a probe for physics at energies near the Hubble scale ( $H$ ). While tree-level single-exchange diagrams are solvable in closed form, and double-exchange results exist for specific limits, higher-order tree diagrams (such as triple exchange) and one-loop diagrams (such as the triangle) have historically lacked convergent analytic representations at general kinematics and arbitrary masses. Existing solutions often rely on hypergeometric series that do not converge at physical kinematics or require analytic continuation, while numerical transport methods are computationally expensive (minutes per configuration). Furthermore, one-loop calculations for general masses have been restricted to soft limits or numerical evaluations at equal masses.

Methodology
The paper introduces a unified representation for both tree-level and one-loop massive de Sitter correlators using the Mellin–Barnes (MB) formalism. The approach builds on two established components:

Propagator Representation: Utilizing the MB representation of the de Sitter propagator (due to Sleight and Taronna), where the time ordering of each line is encoded in a pair of Gamma functions, $\Gamma(\delta)\Gamma(1-\delta)$ .
Vertex Integration: Applying the partial MB method to perform time integrals at each vertex in closed form, yielding Gamma functions of the vertex energies.

The core innovation lies in the contour analysis. The author demonstrates that for tree diagrams with positive vertex powers ( $r_v > 0$ ), there exist straight integration contours in the complex plane that ensure absolute convergence at physical kinematics. Crucially, these contours are mass-independent, allowing for direct evaluation without analytic continuation. For one-loop diagrams, the loop momentum enters only through a flat-space momentum integral with complex propagator powers. The author shows that by commuting the momentum integral with the MB contours, one can isolate the physical signal, handling ultraviolet divergences (in the bubble) and secular logarithms (in the triangle) via residue analysis of the pole at $S = 3/2$ .

Key Contributions and Results

Unified Tree-Level Representation: The author derives a single integral formula (Eq. 4) for arbitrary tree diagrams of heavy exchange. This formula expresses the correlator as a $2L$ -fold integral (where $L$ is the number of internal lines) of Gamma functions over Mellin variables.
Triple-Exchange Bispectrum: The paper provides the first convergent analytic representation for the triple-exchange bispectrum at arbitrary heavy masses and general kinematics.
- Performance: Evaluation takes approximately 0.1 seconds per leg assignment, significantly faster than transport methods (minutes).
- Validation: The results agree pointwise with direct Schwinger–Keldysh integration and the CosmoFlow transport code to a median precision of 0.1% across 48 configurations, including scalene triangles.
- Squeezed Limit: The author computes the squeezed limit coefficient to all orders in the Boltzmann factor, confirming that multi-exchange signals are not exponentially suppressed relative to single exchange.
One-Loop Triangle and Bubble:
- Bubble: The method reproduces known spectral decompositions for the bubble diagram at equal masses, matching term-by-term to $10^{-6}$ .
- Triangle: The paper computes the complete one-loop triangle at general mass. For conformally coupled legs, the result confirms the conformal symmetry prediction of a logarithmic form $-4\pi D_0 \ln k_T + b$ . For $\dot{\pi}$ legs (relevant for the inflaton bispectrum), the logarithmic divergence is absent, and the integral converges on straight contours. The results match brute-force integration to a median of $5.5 \times 10^{-6}$ .
Oscillatory Signatures: The formalism explicitly captures the oscillatory signatures of heavy particles (the "cosmological collider" signal) in the squeezed limit, characterized by exponents $(k_{soft}/k_{hard})^{1/2 \pm i\mu}$ for single exchange and $(k_{soft}/k_{hard})^{2 \pm i(\mu_B + \mu_C)}$ for the triangle pair oscillation.

Significance
The paper claims that this representation covers both trees and loops within a single framework, utilizing contours that do not depend on the heavy masses. This allows for the direct evaluation of correlators at the physical kinematics where cosmological data (such as Planck and BOSS) are analyzed. The computational efficiency makes the representation suitable for generating templates for data analysis. Furthermore, the work establishes a physical picture where, in squeezed limits, each soft line pays a Boltzmann factor while spectator lines contribute only $O(1)$ factors, clarifying the scaling of multi-exchange signals. The author notes that this approach extends the analytic control of massive diagrams in de Sitter space, bearing on the analytic structure and loop consistency of quantum field theory during inflation.

More like this