Warm inflation beyond the Gaussian limit: primordial black holes and induced gravitational waves

In the earliest moments of the universe, a fraction of a second after the Big Bang, space itself underwent a period of rapid expansion known as inflation. During this time, tiny quantum fluctuations were stretched into vast cosmic structures, eventually seeding the galaxies we see today. Most theories assume these fluctuations were perfectly random, following a bell-curve distribution where extreme events are vanishingly rare. However, if the universe was filled with a hot, dense bath of radiation during this expansion, the story changes. In this "warm inflation" scenario, the random jolts from the radiation bath could create rare, massive spikes in the density of the universe. When these spikes re-entered the cosmic horizon, they could collapse directly into primordial black holes. These ancient black holes, if they exist, would leave a dual signature: they would form the black holes themselves, and the violent collision of space-time required to create them would generate a background hum of gravitational waves.

A team of researchers has now investigated how to distinguish between two different ways these random jolts might behave. They focused on two specific properties of the radiation bath: how long the bath "remembers" a jolt before forgetting it, and how the strength of individual jolts is distributed. The first property shapes the overall pattern of the gravitational waves, while the second controls the frequency of the most extreme, rare events that create black holes. The researchers built a mathematical model to simulate a universe where these two properties could be adjusted independently. They found that while the gravitational waves are largely blind to how the jolts are distributed, the formation of black holes is incredibly sensitive to it. This means that two universes could produce the exact same gravitational wave signal but generate vastly different numbers of black holes, depending on the statistical nature of the random kicks.

The study centers on a model where the random force driving inflation is treated as a stream of independent kicks, similar to a series of distinct taps rather than a continuous flow. The researchers examined how the "memory" of the bath, which determines how long a kick influences the system, affects the power spectrum—the map of how strong the fluctuations are at different scales. They found that this memory reshapes the spectrum and directly influences the strength of the induced gravitational waves. However, the distribution of the kick strengths, which determines how often a massive, rare kick occurs, does not change the gravitational wave signal at the leading order. Instead, this distribution dictates the "tail" of the probability curve, the region where the most extreme fluctuations live. It is from this tail that primordial black holes are born.

By running simulations with different statistical distributions for these kicks, the team discovered that a modest change in the distribution's shape could dramatically alter the number of black holes formed without changing the gravitational wave signal. Specifically, they looked at a measure of how "heavy" the tails of the distribution are, known as kurtosis. They found that increasing the kurtosis from a value of 0.1 to 0.3 could boost the rate of black hole formation by factors ranging from ten to a million, even while keeping the gravitational wave amplitude exactly the same. In some cases, a non-Gaussian bath could produce the same number of black holes as a standard Gaussian bath but with gravitational waves that were up to five times weaker. This suggests that the presence or absence of primordial black holes could serve as a powerful diagnostic tool for understanding the statistics of the early universe, independent of the gravitational wave signal.

The researchers tested the robustness of these findings by varying several aspects of their model, such as the shape of the response kernel, which describes how the universe reacts to the kicks, and the method used to count potential black holes. They found that while these changes could slightly adjust the numbers, they never eliminated the massive enhancement in black hole formation caused by the non-Gaussian statistics. Even when they accounted for the fact that black holes form at the peaks of the fluctuations rather than at every point, or when they broadened the spectrum of fluctuations, the enhancement remained significant. The study confirms that the gravitational wave signal is primarily a measure of the covariance, or the overall correlation, of the random force, while the black hole abundance is a direct probe of the higher-order statistics, or the specific shape of the rare events.

This distinction offers a new way to interpret future observations. If a gravitational wave background is detected by instruments like the Laser Interferometer Space Antenna, its amplitude alone cannot tell us how many primordial black holes were formed. The number of black holes depends on the hidden statistics of the source, which the gravitational waves do not reveal. To make a correct prediction, scientists would need to know the specific statistical properties of the radiation bath. The study provides a concrete example where a joint measurement of the gravitational wave amplitude and the black hole abundance could pin down these statistics. For instance, if the gravitational wave signal is measured to be at a certain level, the number of black holes could vary by orders of magnitude depending on whether the source follows a standard Gaussian distribution or a more complex, heavy-tailed one.

The authors also explored how this model relates to other theories of the early universe. Unlike some other models where the distribution of fluctuations is skewed, their model assumes a symmetric distribution where the average kick is zero, meaning the primary deviation from the standard bell curve is in the kurtosis rather than the skewness. This makes the induced gravitational waves, which are blind to the kurtosis at the leading order, a unique tool for isolating the black hole channel. The research suggests that the "memory" of the bath and the "statistics" of the kicks are two distinct handles on the early universe. Changing the memory shifts both the gravitational waves and the black hole abundance together, while changing the statistics shifts only the black holes. This separation allows scientists to disentangle the physical properties of the early universe's radiation bath.

Ultimately, the work highlights that the search for primordial black holes and induced gravitational waves is not just a hunt for two separate signals, but a way to reconstruct the physical conditions of the universe's infancy. The researchers conclude that if a signal is detected, the challenge will not be finding it, but interpreting its amplitude correctly. A signal that is too strong or too weak for the number of black holes observed could indicate that the early universe was governed by a non-Gaussian radiation bath. The study sets the stage for future work that will require a more detailed microscopic calculation of how these kicks arise from the physics of the radiation bath, moving from a general model to a specific prediction for warm inflation. The tools developed in this paper provide a clear path to understanding how the random jolts of the early universe shaped the cosmos we see today.