⚛️ high-energy theory

Exact bispectra in strongly mixed multifield inflation

This paper utilizes the effective field theory of inflation with resummed curvature-isocurvature mixing to derive exact, closed-form expressions for all cubic primordial bispectra, revealing new strong-mixing shapes and universal cosmological collider signals that show a preference for strong mixing when confronted with Planck PR4 data.

Original authors: Lucas Pinol

Published 2026-09-30

📖 9 min read🧠 Deep dive

CC BY 4.0

Original authors: Lucas Pinol

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 first fraction of a second after the Big Bang, the universe underwent a period of explosive expansion known as inflation. This rapid stretching smoothed out the cosmos and stretched tiny quantum fluctuations into the seeds of all future galaxies. For decades, scientists have studied the afterglow of this event, the cosmic microwave background, to understand what happened during that brief, violent epoch. The standard story involves a single field driving this expansion, but many theories suggest that other heavy fields were present, interacting with the main driver. If these extra fields existed, they would have left a unique, oscillating signature in the distribution of matter, a kind of fossil record of the universe's infancy. Detecting this signal is the goal of the "cosmological collider" program, which treats the early universe as a particle accelerator far more powerful than anything we can build on Earth.

The challenge has always been that these signals are incredibly faint and difficult to calculate, especially when the interaction between the fields is strong. Previous methods relied on approximations that broke down when the mixing between fields was intense, leaving a vast region of theoretical possibility unexplored. A new study by Lucas Pinol at the École Normale Supérieure in Paris has changed this landscape. By developing a precise mathematical framework that handles these strong interactions exactly, the researcher has computed the exact shapes of the signals these fields would produce. This work does not just offer a better approximation; it provides a complete, exact map of the possible signals across the entire range of mixing strengths and field masses, revealing shapes that were previously invisible to standard analysis.

The core of this achievement lies in how the researcher treated the mathematics of the early universe. Instead of treating the interaction between the main inflation field and the heavy extra field as a small, manageable disturbance, the study placed the interaction at the very foundation of the calculation. This allowed for an exact solution to the equations governing the fields, even when they were strongly coupled. The result is a set of six distinct signal shapes, each corresponding to a different way the fields could have interacted. These shapes are not just theoretical curiosities; they are the specific patterns that would appear in the data if these heavy fields existed. The study reveals that when the interaction is strong, the signals take on entirely new forms that bear no resemblance to the standard patterns scientists have been looking for in the past. These new shapes are genuinely unique, with no correlation to the familiar templates used in current data analysis.

One of the most striking findings concerns the frequency of the oscillations in these signals. In the squeezed limit, where one part of the signal is much larger than the others, the oscillations carry a specific frequency determined by the mass of the heavy field. The study confirms that this frequency is universal, appearing at every level of interaction and for every type of exchange, regardless of how strong the mixing is. This frequency acts as a clock, ticking at a rate set by the mass of the particle that was produced during inflation. The study also uncovers a surprising behavior in the strength of these signals. When the mixing is weak, the signal is suppressed by a factor related to the mass, a phenomenon known as the Boltzmann factor. However, as the mixing becomes strong, this suppression changes. The signal does not simply vanish; instead, it follows a new exponential rule that depends on both the mixing strength and the mass. This rule connects all the previously studied regimes, showing that the weak-mixing suppression is just one edge of a much larger, continuous landscape.

To test these new ideas against reality, the researcher compared the six exact shapes against the latest data from the Planck satellite, which has mapped the cosmic microwave background with high precision. The analysis looked for any sign of these specific patterns in the data. While no definitive detection was made, the study found a mild preference for the strong-mixing regime. In this regime, the cosmological collider oscillations extend into configurations that are only slightly squeezed, making them potentially more visible than previously thought. The statistical significance of this preference reached a level of 2.8 sigma, which is a hint but not a discovery, especially before accounting for the fact that many different shapes were tested. Nevertheless, this result suggests that the strongest signals might be hiding in the data, waiting to be found with the right template.

The practical impact of this work is immediate and significant. Before this study, calculating these signals for strong mixing was an expensive computational task, often requiring hours of processing time for a single data point. The new framework reduces this to a matter of milliseconds. The researcher has released the code and the exact shapes for the entire range of mixing strengths and masses, allowing other scientists to instantly test any hypothesis against the data. This removes the computational bottleneck that previously limited the search for these signals. The study concludes that the strongly mixed regime, where the signal is largest, can now be searched systematically, just like the weakly mixed one. The door is now open to a comprehensive scan of the early universe's particle content, using the exact shapes of the cosmological collider to guide the search.

The findings also clarify what happens to the signal as the mixing strength changes. In the weak-mixing limit, the results perfectly match the known perturbative calculations, confirming the accuracy of the new method. As the mixing grows, the signals evolve, and the different interaction channels begin to converge toward a single, universal shape. This convergence happens when the mixing is very strong, making the specific details of the interaction less important than the overall strength of the coupling. The study also identifies a specific "sweet spot" for the mixing strength where the signal is maximized for a given model. This optimal point occurs at the bottom of the strong-mixing regime, suggesting that the most promising signals are not at the extremes of weak or infinite mixing, but in a specific intermediate zone.

The research also addresses the behavior of the signals in the squeezed limit, where one momentum is much smaller than the others. Here, the study confirms that the oscillations are driven by a single frequency, set by the effective mass of the heavy field. This frequency is robust, appearing consistently across all interaction orders. The amplitude of these oscillations is governed by a simple exponential rule that interpolates between the known weak-mixing suppression and the strong-mixing enhancement. This rule provides a clear prediction for how the signal should behave as the mixing strength varies, offering a concrete target for future observations. The study also notes that the phase of the oscillation is not universal; it depends on the specific mixing strength and the type of interaction, adding another layer of detail that can help distinguish between different models.

In the confrontation with the Planck data, the study found that the strongest hints of these signals appear in the strong-mixing region, where the oscillations invade mildly squeezed configurations. The data shows a slight preference for these shapes, with the signal-to-noise ratio reaching nearly 2.8 sigma for the most promising channels. While this is not a definitive discovery, it is a significant step forward, demonstrating that the new templates are capable of finding signals that older methods missed. The study also provides upper bounds on the strength of the interactions, constraining the possible values of the coupling constants. These bounds are consistent with current limits but are derived from a more complete and accurate theoretical framework.

The work represents a major advance in the field of cosmological collider physics. By providing exact solutions for the bispectra generated by strong mixing, it removes a long-standing barrier to understanding the early universe. The release of the code and the exact shapes ensures that this progress is immediately available to the broader scientific community. The study shows that the strongly mixed regime is not a theoretical dead end but a rich landscape of new physics, full of unique signatures waiting to be discovered. As the data from future missions becomes available, these exact templates will be essential tools for unlocking the secrets of the universe's first moments. The ability to search the entire plane of mixing and mass with millisecond precision transforms the search for heavy fields from a slow, approximate process into a rapid, systematic exploration.

The study also highlights the importance of looking beyond the standard templates. The new shapes found in the strong-mixing regime are genuinely different from the equilateral and local shapes that have dominated the search so far. They are not just variations of the old patterns; they are distinct functions with their own unique properties. This suggests that the universe may have produced signals that we have been overlooking because we were looking for the wrong shapes. The mild preference for strong mixing in the Planck data serves as a reminder that the answer might lie in the regions we have not fully explored. The study urges a re-evaluation of the data with these new, exact shapes in mind, potentially opening the door to a discovery that has been just out of reach.

Ultimately, this paper provides a new lens through which to view the early universe. It replaces approximations with exactness, and limitations with possibilities. The cosmological collider is no longer just a theoretical idea; it is a concrete program with a clear path forward. The signals are there, the shapes are known, and the tools to find them are now in hand. The next step is to apply these tools to the data, to see if the universe has indeed left the fingerprints of heavy fields in the cosmic microwave background. The journey to understand the origin of the universe continues, and this study has provided a new, precise map for the road ahead.

Technical Summary: Exact Bispectra in Strongly Mixed Multifield Inflation

Problem Statement
The paper addresses a critical limitation in the "cosmological collider" program: the inability to analytically compute primordial bispectra in the regime of strong mixing between the Goldstone boson of broken time translations ( $\pi_c$ ) and a massive isocurvature scalar ( $\sigma$ ). While the perturbative regime (weak mixing, $\lambda \ll 1$ ) is well-understood, the strongly mixed regime ( $\lambda \gtrsim 1$ ) has historically required expensive numerical integration of ordinary differential equations (e.g., via CosmoFlow) to evolve correlators over the entire history of inflation. This computational cost has hindered systematic scans over the parameter space of mixing strength ( $\lambda = \rho/H$ ) and effective mass ( $\mu_{\text{eff}}$ ), preventing a comprehensive search for cosmological collider signals in current and future data. Furthermore, existing perturbative templates fail to capture the qualitative changes in bispectrum shapes and the enhancement of signals that occur at strong mixing.

Methodology
The author extends the Effective Field Theory (EFT) of inflationary fluctuations to include a massive scalar with unit sound speed. The core methodological innovation is the exact resummation of the quadratic mixing into the free Hamiltonian, rather than treating it as a perturbation.

Exact Linear Theory: The coupled linear system for $\pi_c$ and $\sigma$ is solved exactly for any mixing strength $\lambda$ using confluent hypergeometric functions. The solution is recast as a single integral over a dimensionless weight function $\omega_a(u)$ , representing the mode as a superposition of plane waves with "dressed" frequencies.
Diagrammatic Rules: A new set of Schwinger-Keldysh diagrammatic rules is established for this resummed theory. In this framework:
- Propagators are "mixed" (no undressed $\sigma$ lines).
- Exchange diagrams (single, double, triple exchange of the massive field) are reinterpreted as contact diagrams where the massive field legs are attached to a single vertex.
- The interaction picture fields obey the exact linear equations, allowing the cubic interaction vertices to be evaluated directly.
Master Formulas: The bispectrum for every cubic interaction of dimension four or less is reduced to a single one-dimensional integral (Schwinger parameter $\xi$ ) over products of "leg kernels." These kernels are themselves integrals over the weight $\omega_a(u)$ .
Analytic Limits: The paper derives closed-form expressions for the squeezed limit (cosmological collider clock) in both weak and strong mixing regimes, utilizing hypergeometric identities and asymptotic expansions.

Key Contributions and Results

Exact Bispectrum Shapes: The paper provides five master formulas (one for each cubic operator: $\dot{\pi}_c^3$ , $\dot{\pi}_c^2\sigma$ , $(\partial_i\pi_c)^2\sigma$ , $\dot{\pi}_c\sigma^2$ , and $\sigma^3$ ) that are exact in the mixing strength $\lambda$ . These formulas reduce the computation of a bispectrum shape to a few milliseconds on a single CPU core, a reduction of computational cost by roughly five orders of magnitude compared to previous numerical flow methods.
Universal Strong-Mixing Shapes: At strong mixing ( $\lambda \gg 1$ ), the paper unveils genuinely new bispectrum shapes that are uncorrelated with standard templates (equilateral, local, orthogonal). These shapes exhibit sign changes within the momentum triangle and peak in folded configurations.
The Collider Clock Amplitude: A central result is the derivation of the exact scaling of the cosmological collider oscillation amplitude across the entire $(\lambda, \mu_{\text{eff}})$ $(λ, μ_{eff})$ plane. The amplitude scales as:
$\mathcal{A} \propto e^{\frac{\pi}{2}(\lambda - 2\mu_{\text{eff}})}$
This single exponent interpolates between:
- The perturbative Boltzmann suppression $e^{-\pi\mu_{\text{eff}}}$ at weak mixing.
- A suppression of $e^{-\pi\mu_{\text{eff}}/2}$ along lines of fixed bare mass (where $\lambda \approx \mu_{\text{eff}}$ ).
- An exponential growth $e^{\pi\lambda/2}$ at fixed effective mass.
Closed-Form Coefficients:
- Weak Mixing: The paper recovers known perturbative results for single and double exchange in closed form. Crucially, it provides the first closed-form expression for the triple-exchange leading coefficient at arbitrary mass, expressed in terms of hypergeometric functions.
- Strong Mixing: The soft leg coefficients and hard leg factors are derived, showing that the frequency of the clock is universally set by the dressed mass $\mu_{\text{eff}}$ at every exchange order, never $2\mu_{\text{eff}}$ or $3\mu_{\text{eff}}$ (the "1/4 theorem").
Data Confrontation: The six exact shapes were confronted with the Planck PR4 binned bispectrum. The analysis reveals a mild preference for strong mixing ( $\lambda \sim 2\text{--}4$ ), where collider oscillations invade mildly squeezed configurations. The best-fit signals reach up to $2.8\sigma$ (before look-elsewhere correction) for the no-exchange, single-exchange, and triple-exchange channels. Upper bounds on the dimensionless cubic couplings are established for the weak-mixing regime.

Significance and Claims
The paper claims to have formalized a framework that makes the computation of multifield inflation bispectra "essentially free" in terms of computational cost, enabling systematic scans over the full parameter space of mixing and mass.

Theoretical Impact: It resolves the long-standing technical barrier of the "mixed linear system" by providing exact linear solutions and a diagrammatic reduction to Schwinger integrals. It demonstrates that strong mixing fundamentally alters the shape of non-Gaussianity, creating new templates that cannot be captured by perturbative expansions.
Observational Impact: By releasing the code and an interactive explorer, the author enables the community to search for these specific strong-mixing signatures in CMB data without the prohibitive computational cost previously required. The paper suggests that the regime where the cosmological collider signal is largest (strong mixing) is now accessible to data analysis.
Modesty: The paper acknowledges that the hard-leg factor in the strong-mixing amplitude is currently measured numerically rather than derived analytically. It also notes that the loop expansion parameter for the resummed theory is not yet established, and the confrontation with Planck data is a "first" that requires further refinement (e.g., look-elsewhere corrections and joint coupling treatments) in future work. The results are presented as exact within the tree-level, dimension-four EFT framework provided.

Technical Summary: Exact Bispectra in Strongly Mixed Multifield Inflation

More like this