⚛️ phenomenology

Saxion Decay and Dark Radiation in the KSVZ Axion Model

This paper calculates the dark radiation contribution ( $\Delta N_{\rm eff}$ ) from axions produced by the decay of a curvature-displaced Peccei-Quinn condensate in a quartic KSVZ model, finding that the dominant axion channel yields $\Delta N_{\rm eff}$ values of approximately 0.033 and 0.12 while demonstrating that extending the curvature release duration significantly suppresses axion radiation relative to visible radiation without altering primordial isocurvature amplitudes.

Original authors: Sidhartha Samtani, Amartya Sengupta

Published 2026-10-23

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Sidhartha Samtani, Amartya Sengupta

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 threads that hold the fabric of reality together, and among the most persistent mysteries is why the laws of physics seem to treat matter and antimatter with such perfect symmetry, except for one stubborn exception. This exception, known as the strong CP problem, suggests that the universe should behave differently than it does, a discrepancy that has puzzled physicists for decades. To solve this, scientists proposed a mechanism that introduces a new, ghostly particle called the axion. This particle acts like a cosmic dial, constantly adjusting itself to cancel out the unwanted asymmetry, keeping the universe in balance. While the axion is a leading candidate for dark matter—the invisible mass that holds galaxies together—its existence also implies the presence of a heavier partner, a field that stretches and shrinks as the universe evolves. Understanding how these fields behave in the chaotic moments after the Big Bang is crucial, because their decay could have left behind a faint, invisible heat that we can still detect today.

In a new study, researchers have traced the life cycle of this heavier partner, known as the saxion, to see how it influences the universe's thermal history. They focused on a specific scenario where the universe began with a large, smooth field that was pushed away from its natural resting place by the intense curvature of space during the earliest moments of inflation. As the universe expanded and this curvature faded, the field was forced to roll back toward its equilibrium, a journey that caused it to oscillate wildly. The team calculated how these oscillations would behave as the field settled down, and more importantly, how the energy stored in these movements would eventually be released. They found that the saxion does not simply vanish; instead, it decays into two types of particles: the invisible axions and standard particles that make up the visible world. The ratio of these two outcomes determines how much "dark radiation"—a form of energy that behaves like light but cannot be seen—was added to the early universe.

The researchers discovered that the speed at which the universe's curvature faded played a decisive role in the final outcome. If the curvature relaxed quickly, the field was thrown into a violent, chaotic motion, retaining a large amount of energy that was later converted into a significant amount of dark radiation. However, if the curvature faded more gradually, the field could track its path smoothly, losing much of its energy to friction and damping before it ever had a chance to decay. In their simulations, the team tested two specific histories. In one, the field settled relatively quickly, resulting in a dark radiation contribution that is just below the current limits of what astronomers can observe. In the other, a slightly more energetic scenario, the contribution was higher, pushing the predicted value beyond what current measurements allow. This suggests that the universe likely experienced a smoother, more gradual release of curvature, or that the specific properties of the axion field are tuned to keep this extra radiation low.

A key finding of the work is that even a tiny amount of this saxion energy, representing less than one percent of the total energy in the early universe, can produce a measurable effect on the number of effective neutrino species, a standard way physicists count the types of light particles in the cosmos. The study shows that by changing the duration of the curvature release by a small amount, the amount of dark radiation produced can drop by a factor of nearly thirty, while the initial conditions that set the stage for the universe's structure remain unchanged. This sensitivity means that future observations of the cosmic microwave background, the afterglow of the Big Bang, could potentially rule out entire classes of theories about how the axion field behaved. The researchers also confirmed that the decay process is dominated by the production of axions rather than visible particles, meaning that the invisible sector receives the lion's share of the energy, further constraining how much heat this process adds to the cosmic soup.

The study relied on detailed computer simulations that followed the field from its inflationary state through its relaxation and eventual decay, accounting for the complex interactions between the field and the hot plasma of the early universe. They found that while the field can generate large fluctuations that might threaten to disrupt the smoothness of the cosmos, these fluctuations remain small enough in most scenarios to preserve the universe's large-scale structure. However, in the most extreme cases, the fluctuations grow so large that the simple calculations break down, hinting that a more complex, chaotic evolution might be required to fully understand the most energetic histories. The work does not claim to have solved the mystery of dark radiation, but it provides a precise map of how different cosmic histories lead to different observable outcomes. By narrowing down the range of possible scenarios, the study offers a clearer path for astronomers to test the existence of the axion and its heavier partner, turning a theoretical possibility into a concrete target for observation.

Technical Summary: Saxion Decay and Dark Radiation in the KSVZ Axion Model

Problem Statement
The Peccei–Quinn (PQ) mechanism solves the strong CP problem by introducing the QCD axion. If the PQ symmetry breaks before inflation and remains unbroken, the observable universe inherits a homogeneous axion angle. However, quantum fluctuations during inflation perturb this angle, generating cold-dark-matter isocurvature modes that are tightly constrained by Cosmic Microwave Background (CMB) observations. A known mechanism to suppress these fluctuations involves a positive coupling between the PQ field and spacetime curvature, which displaces the PQ condensate to a large radius during inflation. As curvature decreases post-inflation, the radial minimum contracts toward its vacuum value, potentially exciting large radial oscillations (saxions). These oscillations can amplify axion fluctuations via tachyonic growth and parametric resonance, potentially violating isocurvature bounds. Furthermore, the energy remaining in these radial oscillations eventually decays, producing relativistic axions that contribute to dark radiation, quantified as $\Delta N_{\text{eff}}$ . The interplay between the suppression of primordial isocurvature and the generation of dark radiation depends critically on the duration of the curvature release and the subsequent relaxation of the radial field.

Methodology
The authors perform a numerical and analytical study within the Kim–Shifman–Vainshtein–Zakharov (KSVZ) model featuring one heavy colored Dirac fermion. The analysis proceeds through several stages:

Model Setup: The Lagrangian includes a complex scalar $\Phi$ with a quartic potential, a non-minimal curvature coupling $\xi R |\Phi|^2$ , and the KSVZ Yukawa interaction. The radial field is parameterized as $\rho = f_a + s$ , where $s$ is the saxion.
Background Evolution: The authors model the transition from inflation to radiation domination using a smooth interpolation function controlled by a release width parameter $\Delta$ (in e-folds). They fix the inflationary Hubble rate ( $H_*$ ) and the scalar power normalization to satisfy Planck constraints, varying only the quartic coupling $\lambda_\Phi$ and the release duration $\Delta$ .
Field Evolution: The homogeneous radial trajectory is evolved from the inflationary minimum through the curvature release. Transverse fluctuations (axion modes) are evolved on this background to calculate the growth of the variance $V_{\text{band}}$ . The evolution stops if the variance exceeds a cutoff ( $10^{-2}$ ), indicating a breakdown of the linear approximation.
Action Matching: Once the field settles into small harmonic oscillations near the vacuum, the surviving condensate abundance is determined by calculating the oscillation action $J$ . This action is used to set the initial conditions for the late-time fluid evolution.
Decay and Thermal Effects: The saxion decays into axions and KSVZ fermions. The authors calculate the decay widths, including thermal corrections and plasma response (damping and noise) to assess their impact on the relaxation. The energy transfer to the visible bath and the relativistic axion population is tracked until the decay epoch.
Dark Radiation Calculation: The final contribution to the effective number of neutrino species, $\Delta N_{\text{eff}}$ , is computed based on the ratio of axion radiation energy to visible radiation energy after the decay is complete.

Key Contributions

Quantitative Evolution of Curvature-Displaced Condensates: The paper provides a detailed numerical evolution of the PQ condensate through the curvature release, moving beyond analytic estimates to track the specific survival of the condensate and the amplification of fluctuations.
Two-Channel Decay Analysis: The authors rigorously treat the decay of the saxion into both axions and KSVZ fermions. They demonstrate that even with an optimized Yukawa coupling to maximize fermion decay (minimizing the axion branching fraction), the axion channel dominates the decay width due to kinematic constraints and the specific structure of the KSVZ model.
Thermal and Backreaction Assessment: The study includes a systematic check of thermal forces, vacuum corrections, and collisional damping. It establishes that for the benchmark parameters, these effects are subdominant and do not prevent the condensate from reaching the fluctuation cutoff in certain regimes.
Isocurvature vs. Dark Radiation Trade-off: The work explicitly connects the suppression of primordial isocurvature (via large inflationary radii) with the subsequent production of dark radiation, showing how the duration of curvature release acts as a critical parameter linking these two observables.

Results

Benchmark Histories: Two reference histories (A and B) are presented with quartic couplings $\lambda_\Phi = 3 \times 10^{-11}$ $λ_{Φ} = 3 \times 1 0^{- 11}$ and $10^{-11}$ $1 0^{- 11}$ , respectively, and a release width $\Delta = 0.20$ $Δ = 0.20$ .
- Benchmark A: Predicts $\Delta N_{\text{eff}} \simeq 0.033$ . The saxion energy fraction peaks at $\sim 0.87\%$ .
- Benchmark B: Predicts $\Delta N_{\text{eff}} \simeq 0.12$ . The saxion energy fraction peaks at $\sim 3.1\%$ .
- Benchmark C, with a smaller coupling ( $\lambda_\Phi = 3 \times 10^{-12}$ ), reaches the fluctuation cutoff ( $V_{\text{band}} = 10^{-2}$ ) during relaxation, implying a breakdown of the linear approximation and an undetermined late abundance.
Dependence on Release Duration: A primary finding is the strong sensitivity of $\Delta N_{\text{eff}}$ to the duration of curvature release. At a fixed $\lambda_\Phi = 10^{-11}$ , increasing the release width $\Delta$ from 0.20 to 0.50 reduces the predicted $\Delta N_{\text{eff}}$ from $\sim 0.12$ to $\sim 0.0044$ . This represents a reduction factor of approximately 28.
Mechanism of Reduction: The reduction in dark radiation is attributed to the field tracking the contracting minimum more closely during a gradual release. This results in smaller radial excitations and less coherent energy available for late decay, while the primordial angular power (and thus the isocurvature constraint) remains fixed.
Observational Constraints: The calculated $\Delta N_{\text{eff}}$ for Benchmark A ($0.033$) yields a total $N_{\text{eff}} \approx 3.077$ , which lies just below the 95% upper limit reported by recent ACT/Planck analyses ( $N_{\text{eff}} < 3.08$ ). Benchmark B ( $N_{\text{eff}} \approx 3.165$ ) exceeds this bound and is disfavored under the assumptions of that analysis.

Significance and Claims
The paper claims to demonstrate that a subdominant saxion population (contributing less than 4% of the total energy density) can generate an appreciable dark radiation signal. The study highlights that the "curvature release" mechanism, often invoked to suppress isocurvature, simultaneously controls the amplitude of dark radiation. The authors assert that the duration of this release is a decisive factor: a more gradual release suppresses both the transverse fluctuation growth (preserving the isocurvature bound) and the coherent energy available for decay (lowering $\Delta N_{\text{eff}}$ ).

The authors remain modest regarding the limitations of their calculation. They note that for histories where fluctuations grow large (like Benchmark B), nonlinear effects and backreaction could alter the surviving condensate abundance, a regime that would require lattice simulations to resolve. Furthermore, they acknowledge that their decay estimates rely on a weak-Yukawa regime and a negligible Higgs portal; stronger couplings or additional portals could introduce significant damping or new decay channels, potentially changing the parameter space allowed by observational bounds. The work concludes that while the current benchmarks offer testable predictions, a full cosmological assessment requires a joint treatment of the physical prior $\Delta N_{\text{eff}} \ge 0$ and the isocurvature constraints, alongside a more complete treatment of nonlinear fluctuation evolution.

More like this