Gravitational waves from axion inflation in the gradient expansion formalism. Part I. Pure axion inflation

Using the gradient expansion formalism, this paper presents a detailed parameter scan of gravitational wave production in pure axion inflation and finds that detectable signals require strong backreaction regimes that conflict with current constraints on dark radiation ($\Delta N_{\rm eff}$), thereby defining a critical target for future lattice studies.

The Gist

The Big Picture: Listening to the Baby Universe

Imagine the universe as a giant, expanding balloon. About 13.8 billion years ago, this balloon inflated incredibly fast in a period called inflation. Scientists believe that during this rapid expansion, the universe didn't just get bigger; it also got "noisy."

This noise comes in the form of Gravitational Waves (GWs)—ripples in the fabric of space-time itself. While we have detected waves from colliding black holes recently, this paper is about a specific type of "background hum" (a stochastic background) that might have been created in the very first split-second of the universe. The authors are trying to figure out if our future telescopes (like the Einstein Telescope or LISA) can hear this hum.

The Characters in the Story

1. The Axion (The Inflaton): Think of the axion as a giant, rolling ball (the inflaton field) that drives the expansion of the universe. It's rolling down a hill (its potential energy).
2. The Gauge Fields (The Magnetic Storm): As the axion rolls, it has a special, magical connection to invisible "magnetic" fields (gauge fields). It's like the axion is a conductor, and as it moves, it creates a massive, swirling storm of magnetic energy around it.
3. The "Pure" Model: This paper looks at a specific version of the story where the axion only talks to these magnetic fields and nothing else. There are no other particles (like electrons or protons) getting in the way. This is called Pure Axion Inflation (PAI).

The Mechanism: A Snowball Effect

Here is the core physics, explained simply:

• The Setup: As the axion rolls down its hill, it spins up the magnetic fields.
• The Feedback Loop: These magnetic fields get so strong that they push back on the axion. Imagine the axion is a car driving down a hill, but the magnetic fields act like a giant, invisible hand grabbing the car's bumper and trying to slow it down.
• The "Strong Backreaction": If the magnetic storm gets too strong, it creates a "friction" that drastically changes how the axion rolls. It might make the axion bounce or roll much slower than expected. This is the Strong Backreaction Regime.

The Discovery: The "Catch-22" of Detection

The authors ran massive computer simulations to see if we could detect the gravitational waves produced by this magnetic storm. They found a very frustrating, yet fascinating, rule:

You can only hear the signal if the storm is violent enough to break the rules.

Here is the dilemma they discovered:

1. The Weak Storm: If the magnetic fields are weak, the axion rolls smoothly. The gravitational waves produced are too faint. Future telescopes will hear nothing but silence.
2. The Strong Storm: To get a signal loud enough for our telescopes to hear, the magnetic fields must be incredibly strong. This triggers the "Strong Backreaction" (the friction mentioned above).
3. The Problem: When the storm is strong enough to be heard, it creates too much energy. It dumps so much extra radiation into the early universe that it breaks the "budget" allowed by the Big Bang theory.

The Analogy:
Imagine you are trying to hear a whisper from across a crowded room.

• If the person whispers normally (weak backreaction), you can't hear them.
• If the person screams (strong backreaction), you can definitely hear them.
• BUT, if they scream, the noise is so loud that it violates the "Quiet Zone" rules of the building (the $\Delta N_{eff}$ limit, which is a limit on how much extra energy/radiation the universe can have). If they scream, the building manager (the laws of physics/CMB data) kicks them out.

The Conclusion: A Dead End (For Now)

The authors conclude that, within the limits of their current computer model:

• Detectable signals are impossible without breaking the laws of physics regarding the early universe's energy budget.
• The "Strong Backreaction" is a double-edged sword: it makes the signal loud enough to find, but also makes it too loud to be allowed.

Why This Paper Matters

Even though the news is "bad" (we probably won't see this signal), the paper is a huge success for science because:

1. It sets a target: They used a clever math trick (the Gradient Expansion Formalism) to scan millions of possibilities quickly. They found the "danger zone" where the signal is too loud.
2. It guides future research: They are telling other scientists: "Hey, if you want to find this signal, you need to look at the 'Strong Backreaction' zone. But be careful! Our simple model says it's forbidden. Maybe you need a more complex model (like a full 3D lattice simulation) that accounts for the axion field getting messy and uneven, which might soften the blow."

Summary in One Sentence

The paper finds that the only way to hear the gravitational waves from this specific type of early-universe inflation is to crank the volume up so high that it violates the universe's energy limits, suggesting that either we won't hear it, or our current understanding of the "volume knob" needs a more complex model to fix.

Technical Summary

1. Problem Statement

The paper investigates the production of Stochastic Gravitational Wave Backgrounds (SGWB) during Pure Axion Inflation (PAI). In this model, a pseudoscalar inflaton field ($\phi$) couples to a pure Abelian gauge field via a Chern-Simons interaction ($\phi F\tilde{F}$).

• The Mechanism: The rolling inflaton breaks parity spontaneously, leading to the tachyonic amplification of one helicity of the gauge field. This produces abundant helical gauge fields, which in turn source gravitational waves (GWs) through anisotropic stress.
• The Challenge: While strong gauge-field production can generate observable GW signals, it also introduces a "friction" term on the inflaton (backreaction). If this backreaction becomes too strong, it alters the inflationary dynamics, potentially leading to an overproduction of GWs that violates cosmological bounds on dark radiation ($\Delta N_{\text{eff}}$).
• The Gap: Previous studies have relied on lattice simulations (computationally expensive, limited parameter range) or analytical approximations valid only in the weak-backreaction regime. There was a lack of a systematic, high-resolution parameter scan covering the transition from weak to strong backreaction to determine if observable GW signals are compatible with cosmological constraints.

2. Methodology

The authors utilize the Gradient Expansion Formalism (GEF), a powerful numerical technique that captures non-linear dynamics in the limit of vanishing axion gradients (homogeneous inflaton field).

• Model Setup:
  • Inflaton Potential: Assumed to be quadratic ($V(\phi) = \frac{1}{2}m^2\phi^2$) near the end of inflation, which is the standard benchmark for strong backreaction studies.
  • Coupling: Linear axion-vector coupling $I(\phi) = (\beta/M_P)\phi$.
  • Parameters: The scan focuses on the axion coupling constant $\beta$ and the inflaton mass $m$.
• Numerical Approach (GEF):
  • Instead of solving the full lattice equations for inhomogeneous fields, the GEF solves evolution equations for bilinear expectation values of the gauge fields (e.g., $\langle E \cdot \text{rot}^n E \rangle$).
  • This allows for a computationally efficient exploration of a vast parameter space.
  • The system is evolved from the slow-roll attractor until the end of inflation (defined by $\ddot{a} < 0$).
• GW Calculation:
  • The induced tensor power spectrum ($P_T^{\text{ind}}$) is computed using the mode functions of the gauge fields derived from the GEF background evolution.
  • The GW energy density spectrum $\Omega_{\text{GW}}(f)$ is calculated, accounting for cosmological redshift and transfer functions.
• Constraints Applied:
  • $\Delta N_{\text{eff}}$: The integrated energy density of GWs must not exceed the bound on extra relativistic degrees of freedom ($\Delta N_{\text{eff}} \lesssim 0.5$).
  • CMB Constraints: Limits on the tensor-to-scalar ratio ($r$) and non-Gaussianities constrain the inflaton mass $m$.
  • Detector Sensitivity: Signal-to-Noise Ratios (SNR) are calculated for future detectors: Einstein Telescope (ET), LISA, and current HLVO3 (LIGO-Virgo-KAGRA).

3. Key Contributions

1. First Detailed Parameter Scan: The paper presents the first extensive scan of the $(\beta, m)$ parameter space for Abelian PAI specifically targeting the regime close to the onset of strong backreaction.
2. GEF Benchmark: It establishes a rigorous benchmark using the Gradient Expansion Formalism to define the "homogeneous" limit of strong backreaction, providing a target for future, more expensive lattice simulations that include inhomogeneities.
3. Identification of a "No-Go" Region: The study identifies a sharp transition in parameter space where the system moves from a weak-backreaction regime (undetectable GWs) to a strong-backreaction regime (observable GWs but cosmologically forbidden).

4. Key Results

• The Trade-off: The authors find a strict dichotomy: Observable GW signals (SNR > 1) for ET, LISA, or HLVO3 can only be realized in parameter regions that trigger strong backreaction.
• Conflict with $\Delta N_{\text{eff}}$: The same strong backreaction required to produce an observable signal leads to an overproduction of gravitational radiation. This results in $\Delta N_{\text{eff}} > 0.5$, violating Big Bang Nucleosynthesis (BBN) and CMB constraints.
• Sharp Transition: There is a very narrow boundary (relative mass variation $\delta m/m \sim 0.02$) separating the two regimes.
  • Below the threshold: Backreaction is negligible; GW signals are too weak to be detected.
  • Above the threshold: Backreaction becomes strong, extending inflation and amplifying GWs to observable levels, but simultaneously violating $\Delta N_{\text{eff}}$ bounds.
• Phenomenological Fit: The authors derive a boundary equation for the excluded region:
  $$\log_{10}\left(\frac{m}{M_P}\right) \lesssim -0.387\beta - 1.34 \log_{10}\left(\frac{\beta}{10}\right) + 0.51$$
  Any parameter point above this line (in the direction of stronger coupling/lower mass) leads to $\Delta N_{\text{eff}}$ violation.
• Comparison with Lattice Data: The GEF results are consistent with existing lattice simulations (e.g., Refs. [44, 46]) in the weak-backreaction regime. However, the GEF predicts a harder exclusion boundary. The authors note that lattice simulations including axion inhomogeneities (gradients) might dampen the GW production slightly, potentially allowing a "window" of detectability that the homogeneous GEF misses.

5. Significance and Outlook

• Theoretical Implication: Within the context of the homogeneous GEF approximation, Pure Axion Inflation coupled to a pure Abelian gauge sector cannot produce a detectable SGWB signal without violating cosmological bounds. The "sweet spot" for detection is effectively closed by the $\Delta N_{\text{eff}}$ constraint.
• Guidance for Future Work: The paper defines a clear "target region" for future lattice studies. Since GEF assumes a homogeneous inflaton, the authors suggest that including axion gradients (which are known to dampen backreaction effects in lattice simulations) might smooth the transition between regimes. This could potentially open a narrow window where GWs are observable but $\Delta N_{\text{eff}}$ is satisfied.
• Distinction from Fermionic Models: The paper distinguishes PAI from Fermionic Axion Inflation (FAI), where the Schwinger effect generates charged fermions that dampen gauge-field production. The companion paper (Ref. [84]) suggests FAI might avoid these constraints, making the distinction between pure and fermionic models crucial for GW phenomenology.

Conclusion: The study concludes that for homogeneous pure axion inflation, the detection of gauge-field-induced gravitational waves is inextricably linked to the violation of dark radiation bounds. This result serves as a critical benchmark, urging future lattice simulations to investigate whether inhomogeneous effects can reconcile observability with cosmological constraints.