The Challenge of Observing Patchy Reionization with CMB Optical-Depth Fluctuations

In the early history of our universe, a profound transformation occurred. For hundreds of thousands of years after the Big Bang, the cosmos was a foggy, neutral gas, opaque to light. Then, as the first stars and galaxies ignited, their intense ultraviolet radiation began to strip electrons from atoms, turning the vast spaces between galaxies into a transparent, ionized plasma. This era, known as cosmic reionization, is one of the final major phase transitions in the universe's story. While we know this event happened roughly between 13 and 12 billion years ago, the exact details of how it unfolded—whether it happened quickly or slowly, and whether the ionized regions grew like expanding bubbles or spread out more evenly—remain some of the biggest open questions in cosmology.

Scientists look for clues to this ancient history in the cosmic microwave background, the faint afterglow of the Big Bang that fills the sky. As light from the early universe traveled toward us, it passed through the newly freed electrons of the reionization era. These electrons acted like tiny mirrors, scattering the light and leaving a subtle imprint on the cosmic background radiation. By measuring the average amount of this scattering, astronomers can determine when reionization finished. However, the average value tells only part of the story. If the process was patchy, with some regions ionizing early and others late, it would create tiny fluctuations in the scattering across the sky. These fluctuations hold the key to understanding the shape and structure of the reionization era, offering a glimpse into how the first galaxies carved out their place in the dark.

A team of researchers has now taken a closer look at these fluctuations, using powerful computer simulations to untangle two different causes behind the signal. The team, led by Nick Takoudes, Hanjue Zhu, and Nickolay Gnedin, utilized a suite of advanced simulations called Cosmic Reionization on Computers, or CROC. These simulations model the complex dance of dark matter, gas, and radiation as the first galaxies form and the universe transitions from neutral to ionized. By constructing a virtual "light cone"—a way of looking back in time through the simulated universe—the researchers created detailed maps of how the electron density varied across the sky during this critical epoch.

The researchers faced a significant challenge in interpreting these maps. The fluctuations in the scattering signal arise from two distinct sources. The first is the patchy nature of reionization itself: the fact that some areas became ionized while others remained neutral, creating a mosaic of different ionization levels. The second source is simply the clumpiness of the gas. Even if the entire universe were ionized at the exact same moment, the gas would still be denser in some places and thinner in others, simply because matter naturally clumps together under gravity. This density variation would also create fluctuations in the scattering signal, even without any patchy ionization. To understand the true story of reionization, the team had to separate these two effects.

Using their simulations, the team broke down the electron density into these two components: the variations caused by the changing ionization state, and the variations caused by the underlying density of the gas. They then calculated the power spectrum of these fluctuations, a statistical measure that describes how strong the variations are at different angular scales on the sky. The results were surprising. Contrary to the hope that the signal would be dominated by the patchy nature of reionization, the researchers found that the fluctuations in the gas density were the primary driver of the signal. Across almost all the scales they examined, the clumpiness of the gas contributed more to the fluctuations than the patchy ionization did.

This finding held true even when the researchers accounted for the different ways the simulations were constructed and the different large-scale environments they modeled. In some simulations, they artificially enhanced the density of the universe to see how it affected the timing of reionization. They found that while the overall timing changed, the dominance of the density signal remained. In fact, the researchers discovered that the patchy and density components often worked against each other. On small scales, the denser regions of gas tended to recombine electrons more quickly and shield themselves from radiation, leading to a lower ionization fraction than the surrounding, less dense areas. This created a negative correlation, where the density signal and the patchy signal partially canceled each other out, further complicating the picture.

The study also looked beyond the reionization era itself. The simulations ended when the universe was about a billion years old, but the universe continued to evolve for billions of years after that. The researchers added an analytical estimate for the signal coming from the later, fully ionized universe. They found that this later contribution, driven entirely by the growth of cosmic structure in a fully ionized gas, was actually larger than the signal from the reionization era itself. When combined, the total signal is overwhelmingly dominated by the density of the gas, both during the reionization era and in the billions of years that followed.

These results suggest that interpreting the scattering fluctuations as a direct map of the patchy reionization process is misleading. The signal is not a clean window into the shape of the ionized bubbles; it is heavily contaminated by the natural clumpiness of the gas. To truly understand the morphology of reionization, future studies will need to carefully model and subtract the density contribution, rather than assuming the entire signal comes from the patchy ionization. The researchers conclude that while the fluctuations in the cosmic microwave background do contain information about reionization, extracting that information requires a much more nuanced approach than previously thought, one that acknowledges the dominant role of the underlying cosmic web.