⚛️ phenomenology

A Tale of Two Gauges: Effective Field Theory for Relativistic Behavior of Cosmological Axions

This paper extends the Salehian, Namjoo & Kaiser effective field theory framework for relativistic axions to the synchronous gauge and to cases with non-zero anisotropic stress, while proposing a fluid interpretation that identifies the axion field as a perfect fluid at all times to facilitate its implementation in cosmological Boltzmann solvers.

Original authors: Hoang Nhan Luu, Chanda Prescod-Weinstein

Published 2026-09-11

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Hoang Nhan Luu, Chanda Prescod-Weinstein

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

The universe is filled with invisible substances that hold galaxies together and drive the expansion of space itself. Among the most promising candidates for these hidden components are axions, hypothetical particles that were originally proposed to solve a specific puzzle in particle physics but have since become leading contenders for dark matter. These particles are incredibly light, so light that if they exist in the vast quantities required to explain the cosmos, they would behave less like individual bullets and more like a vast, rippling wave stretching across the entire universe. This wave-like nature is what makes them so difficult to study. Unlike the steady, slow drift of ordinary matter, axions are predicted to vibrate with a frequency so rapid that they complete trillions of cycles in the time it takes for a single galaxy to rotate.

To understand the universe, cosmologists rely on computer simulations that track how matter clumps and moves over billions of years. However, trying to simulate a substance that vibrates trillions of times faster than the simulation's clock ticks is a computational nightmare. It would require a computer to take steps so small and numerous that the calculation would take longer than the age of the universe to finish. For decades, scientists have used a shortcut: they have smoothed out these rapid vibrations, treating the axion field as a calm, average fluid rather than a frantic, oscillating wave. This approach has worked well enough for many purposes, but it leaves a gap in our understanding. It fails to capture the subtle, high-frequency interactions that might occur when axions coexist with other cosmic ingredients like light, neutrinos, and ordinary matter.

In a new study, researchers Hoang Nhan Luu and Chanda Prescod-Weinstein have developed a more sophisticated way to model these particles, one that bridges the gap between the impossible task of tracking every single vibration and the oversimplified view of a smooth fluid. Their work focuses on a mathematical framework known as effective field theory, which allows them to separate the fast, jittery movements of the axion from its slow, large-scale evolution. By doing this, they can describe the axion as a perfect fluid at all times, capturing both its average behavior and the tiny, crucial corrections caused by its rapid oscillations. This method is significant because it works in the specific mathematical language used by the most powerful cosmological simulation tools currently available, allowing these tools to finally include axions with a level of precision that was previously out of reach.

The core challenge the authors address is the sheer speed of the axion's motion. In the standard model of cosmology, the universe expands over a timescale of billions of years. An axion with a mass typical of the "fuzzy dark matter" candidates vibrates on a timescale of months. To simulate the history of the universe with these particles, a computer would need to calculate the state of the axion field roughly one hundred billion times for every single step of cosmic expansion. This is not just slow; it is practically impossible. Previous attempts to solve this involved averaging the axion's position over one full vibration cycle, effectively blurring the rapid motion into a steady flow. While this "cycle-averaging" method is elegant and has been used in popular cosmology codes, it has limitations. It often assumes the axion is the only thing in the universe or that the universe is perfectly simple, ignoring the complex interplay that happens when multiple types of matter are present.

Luu and Prescod-Weinstein tackled this by refining the mathematical machinery used to describe the axion. They started by redefining the axion field not as a single number that changes over time, but as a "wavefunction" that contains both a slow-moving part and a fast-oscillating part. This is similar to how a sound wave can be described as a steady tone with a rapid vibration superimposed on it. The researchers then applied a systematic process to separate these two components. They treated the fast vibrations as a source of small, corrective forces that push and pull on the slow-moving background. By calculating these corrections step-by-step, they were able to write down a new set of equations that describe the axion's behavior without ever needing to resolve the individual, trillion-cycle vibrations.

A key innovation in this work is the choice of mathematical "gauge," or coordinate system, used to describe the universe. In cosmology, the same physical reality can be described using different coordinate systems, much like describing a landscape using a grid of latitude and longitude versus a grid of distance from a central point. One system, known as the Newtonian gauge, is intuitive and easy to visualize, but it can become unstable and difficult to use in complex computer simulations. The other, called the synchronous gauge, is less intuitive but is the standard language of the most widely used cosmological software. Previous studies had only developed this advanced averaging method for the intuitive system. Luu and Prescod-Weinstein extended the method to the synchronous gauge, ensuring that their results could be directly plugged into the tools that cosmologists use to predict what the universe should look like. They also generalized the math to handle scenarios where the universe is filled with a mix of different materials, not just axions.

The researchers demonstrated that their new equations allow for a "fluid interpretation" of the axion field that remains valid at all times. In simpler terms, they showed that you can treat the axion as a fluid with a specific density, pressure, and flow, even while it is vibrating at its fastest. This is a departure from older methods where the fluid description was only an approximation that broke down under certain conditions. Their equations include "backreactions," which are the subtle effects the fast vibrations have on the slow evolution of the universe. For example, the rapid oscillation of the axion field slightly alters the expansion rate of the universe and the way gravity pulls on matter. By including these effects, the new model provides a more complete picture of how axions would influence the formation of galaxies and the cosmic microwave background radiation.

To prove their method works, the authors tested it against a known solution. They simulated the evolution of the axion field using both their new, simplified equations and the original, complex equations that track every vibration. The results showed that their new method could reconstruct the exact, rapid oscillations of the axion field with incredible accuracy, while simultaneously tracking the smooth, average behavior that drives cosmic evolution. The reconstructed values matched the exact solution almost perfectly, differing by less than one percent in some cases. This level of agreement confirms that the mathematical shortcuts they used do not lose important information. Instead, they successfully filter out the noise of the rapid vibrations while keeping the signal of the physical effects.

This work lays the theoretical foundation for a companion paper where the authors will implement these equations into a major cosmological simulation package. The goal is to create a tool that can accurately predict how axions would affect the large-scale structure of the universe, from the distribution of galaxies to the patterns seen in the cosmic microwave background. By providing a rigorous and flexible framework, this study opens the door to testing a wider range of axion models against real observational data. It moves the field beyond simple approximations and allows scientists to ask more precise questions about the nature of dark matter. If axions exist, they are likely the most abundant particles in the universe, and understanding their subtle, high-speed dance is essential to understanding the history and future of the cosmos itself.

Technical Summary: A Tale of Two Gauges: Effective Field Theory for Relativistic Behavior of Cosmological Axions

Problem Statement
Axions and axion-like particles (ALPs) are promising dark matter candidates, particularly those with ultralight masses ( $m \lesssim 10^{-18}$ eV). While their behavior on cosmological scales is of significant interest, modeling their relativistic evolution presents a computational challenge. The axion field undergoes "Compton" oscillations on a timescale of $m^{-1}$ , which is orders of magnitude shorter than the Hubble timescale ( $H^{-1}$ ) for typical axion masses (e.g., $m \sim 10^{-22}$ eV). Directly solving the exact equations of motion for these oscillations throughout cosmic history requires resolving $10^{11}$ to $10^{12}$ cycles, rendering standard numerical integration prohibitively expensive.

Previous approaches have attempted to mitigate this by "cycle-averaging" the field, treating axions as an effective fluid. However, these methods often rely on specific ansatzes (e.g., WKB) that assume the field can be decomposed into simple trigonometric functions, limiting their applicability to non-linear or complex axion theories. Furthermore, existing Effective Field Theory (EFT) frameworks for axions, such as the work by Salehian, Namjoo & Kaiser [1], were formulated specifically for the Newtonian gauge under the assumption of zero anisotropic stress ( $\Phi = \Psi$ ). This restricts their utility in realistic multi-component cosmological scenarios (e.g., $\Lambda$ CDM with photons and neutrinos) where anisotropic stress is non-negligible and where numerical solvers like CAMB and CLASS primarily operate in the synchronous gauge.

Methodology
The authors extend the EFT formalism to address these limitations through a systematic two-step expansion:

Time-Frequency Expansion (Mode Expansion): The authors decompose time-dependent variables into a Fourier series with frequencies $\nu m$ , separating "slow modes" ( $\nu=0$ ) from "fast modes" ( $\nu \neq 0$ ). A time-averaging operator is defined to isolate the slow modes, which evolve on the Hubble timescale, while the fast modes represent the rapid oscillations.
Perturbative Expansion in Small Parameters: The authors introduce a small parameter $\epsilon$ (representing relativistic corrections, e.g., $H/m$ , $k^2/m^2a^2$ ) and expand the dynamical variables in powers of $\epsilon$ . This allows for the "integration out" of high-frequency modes, replacing them with backreaction terms on the slow modes.

The paper derives these equations in two distinct gauges:

Synchronous Gauge: The authors derive the EFT equations for linear perturbations in the synchronous gauge for the first time. This gauge is preferred for numerical implementation in standard Boltzmann solvers.
Newtonian Gauge: The authors re-derive the EFT equations in the Newtonian gauge without assuming $\Phi = \Psi$ , thereby generalizing the framework to include non-zero anisotropic stress from other cosmological species.

The derivation involves redefining the axion field $\phi$ in terms of a complex wavefunction $\psi$ (Eq. 13) to factor out oscillatory terms. The authors then systematically compute corrections up to order $O(\epsilon^4)$ for background and perturbation variables, including the wavefunction, metric perturbations, and fluid quantities (density, pressure, velocity).

Key Contributions

Synchronous Gauge EFT: The paper provides the first derivation of axion EFT equations in the synchronous gauge. This includes explicit expressions for the slow-mode evolution of the wavefunction perturbation $\delta\psi_s$ , metric perturbations $h_s$ and $\eta_s$ , and axion fluid variables.
Generalized Newtonian Gauge EFT: The authors extend the Newtonian gauge formalism to scenarios with non-zero anisotropic stress. They demonstrate that the previous assumption of $\Phi = \Psi$ is not strictly necessary for the validity of the EFT up to the derived order, making the framework applicable to multi-component universes.
Fluid Interpretation: Unlike previous effective treatments that identified axions as imperfect fluids with viscosity, this work proposes a fluid interpretation where the axion field is identified as a perfect fluid at all times (both in exact and effective regimes). The authors distinguish between "slow-mode" fluid quantities (averaged properties of the exact field) and "effective" fluid quantities derived from imperfect fluid equations, noting that the differences are negligible at low orders but theoretically distinct.
Reconstruction Capability: A unique feature of this EFT formalism is the ability to reconstruct the exact oscillatory behavior of the axion field from the slow-mode solutions. By adding the calculated high-frequency corrections (backreactions) to the slow modes, one can recover the full time-dependent solution without solving stiff differential equations.

Results
The paper presents a complete set of coupled differential equations for the slow-mode variables in both gauges (summarized in Sec. III D and Sec. IV B).

Background Equations: The authors derive effective equations for the background wavefunction $\tilde{\psi}_s$ and Hubble parameter $\tilde{H}_s$ that include relativistic corrections up to $O(\epsilon^4)$ .
Perturbation Equations: The paper provides the evolution equations for metric perturbations ( $\eta_s, h_s$ in synchronous; $\Phi_s, \Psi_s$ in Newtonian) and axion perturbations ( $\delta\tilde{\psi}_s$ ). These equations explicitly include terms representing the backreaction of high-frequency modes on the slow evolution.
Validation: The authors validate the formalism by comparing the reconstructed EFT solutions against exact numerical solutions of the original Klein-Gordon equations. Figure 1 demonstrates that the EFT slow modes accurately capture the average behavior of the exact oscillatory variables, and the reconstructed curves match the exact solutions with subpercent agreement.
Gauge Consistency: The results in the Newtonian gauge are shown to reduce to the results of [1] in the limit where anisotropic stress vanishes ( $\Psi_s \to \Phi_s$ ), confirming consistency with previous work while offering a more general formulation.

Significance
The authors state that this work lays the theoretical foundation for a companion paper [2], which details the implementation of these EFT equations into a modified version of the CLASS Boltzmann solver. The significance of this work lies in:

Numerical Efficiency: It enables the precise modeling of axion perturbations in realistic cosmological settings without the computational cost of resolving rapid oscillations.
Generality: By formulating the theory in the synchronous gauge and allowing for non-zero anisotropic stress, the framework is directly applicable to standard cosmological codes and multi-component universes.
Systematic Precision: The EFT approach allows for arbitrary precision control through higher-order relativistic corrections, offering a systematic path to improve accuracy beyond simple cycle-averaging.
Theoretical Clarity: It clarifies the distinction between exact field dynamics and effective fluid descriptions, providing a rigorous "perfect fluid" interpretation for axions that remains valid even when coupled to other species.

The paper concludes that the EFT formalism is a robust tool for studying the phenomenology of ultralight axions, from fuzzy dark matter to evolving dark energy, and provides the necessary theoretical machinery for their inclusion in next-generation cosmological analyses.

More like this