New Bounds on Heavy QCD Axions from Big Bang Nucleosynthesis

This paper establishes robust upper bounds on heavy QCD axion lifetimes as low as 0.017 seconds by analyzing Big Bang Nucleosynthesis constraints on neutron-to-proton ratio modifications from hadronic decays, a limit that surpasses projected future CMB bounds and remains largely insensitive to hadronic uncertainties.

The Gist

Here is an explanation of the paper "New Bounds on Heavy QCD Axions from Big Bang Nucleosynthesis," translated into everyday language with creative analogies.

The Big Picture: A Cosmic Kitchen Disaster

Imagine the universe, just a few minutes after the Big Bang, as a giant, bubbling kitchen. In this kitchen, the chefs (protons and neutrons) are trying to bake a very specific cake: Helium-4.

In our standard recipe (the Standard Model of physics), the chefs work at a steady pace. They swap ingredients back and forth until the oven gets too hot, at which point they stop swapping and start baking. The result is a cake with a very precise amount of Helium. Astronomers have looked at the oldest stars in the universe and measured this cake; it matches the standard recipe almost perfectly.

The Problem: What if there was a secret ingredient added to the kitchen that we didn't know about? Specifically, what if a mysterious particle called a Heavy QCD Axion showed up, decayed, and started throwing extra "hadrons" (strong-force particles like pions and kaons) into the mix?

These extra particles are like a chaotic sous-chef running around the kitchen, violently swapping the chefs' ingredients. Because the "strong force" is incredibly powerful (about a trillion times stronger than the weak force that usually governs these swaps), even a tiny amount of this chaos would ruin the recipe, changing the amount of Helium baked.

The Detective Work: The "Neutron-to-Proton Ratio"

The scientists in this paper are cosmic detectives. They know the final amount of Helium is a fingerprint of what happened in those first few minutes.

• The Neutron/Proton Ratio: Think of neutrons and protons as two types of dough. Before the cake bakes, they are constantly turning into each other.
• The Freeze-Out: Eventually, the universe cools down, and this swapping stops. The ratio at that exact moment determines how much Helium gets baked.
• The Intruder: If a Heavy Axion decays and injects energetic particles, it acts like a high-speed blender, forcing neutrons to turn into protons (or vice versa) much faster than nature intended. This changes the final Helium count.

The New Discovery: Catching the Culprit Earlier

Previous studies looked for these Axions by checking if they messed up the temperature of the universe (via the Cosmic Microwave Background). But this paper says: "Wait a minute! We can catch them much earlier and much more easily by looking at the Helium cake."

The team found that if these Heavy Axions exist and live for a certain amount of time, they would have caused a massive "kitchen disaster" that would have left the universe with the wrong amount of Helium. Since the universe does have the right amount of Helium, these Axions cannot exist with those specific properties.

The Verdict: They have set a new, stricter rule. If these Heavy Axions exist, they must decay extremely quickly—in less than 0.02 seconds (that's 1/50th of a second). If they lived any longer, they would have ruined the Helium recipe, and we wouldn't be here to talk about it.

The "Waterfall" Analogy: How the Limit Works

The paper describes a fascinating mechanism called the "Waterfall."

1. The Dam: Imagine a dam holding back a massive lake of neutrons and protons. In the normal universe, the water level (the ratio) changes slowly.
2. The Axion Injection: If an Axion decays, it's like someone blowing a hole in the dam. Suddenly, water rushes out, and the level changes violently and instantly.
3. The Waterfall: The ratio of neutrons to protons "falls" rapidly to a new, chaotic level.
4. The Lock: Eventually, the Axions run out (they all decay), and the hole in the dam closes. The water level tries to go back to the normal recipe, but it's too late. The "baking" has already started with the wrong ratio.

The scientists calculated exactly how big the hole could be before the cake was ruined. They found that even a tiny hole (a very short-lived Axion) is enough to spoil the Helium.

Why This Paper is Special (The "Upgrade")

The authors didn't just do a quick calculation; they upgraded the entire simulation engine. Here's what they improved:

• The "Ghost" Particles ($K_L$): They paid special attention to a particle called the Long-Lived Neutral Kaon ($K_L$). It's like a ghost in the kitchen; it doesn't interact much, so it flies around at high speeds. Previous studies ignored it or treated it simply. This team realized these "ghosts" bounce around and change the chaos in subtle but important ways. They mapped out exactly how these ghosts move and collide.
• The "Secondary" Chaos: They realized that when a particle hits another, it doesn't just stop; it creates new particles (secondary hadrons). It's like a billiard ball hitting a rack of balls, which then scatter and hit other balls. They tracked this entire chain reaction, which previous studies missed.
• Robustness: They tested their theory against all kinds of uncertainties. They asked, "What if our math on particle collisions is off by 50%?" They found that even with huge errors in the details, the main conclusion (the 0.02-second limit) stays rock solid. It's like saying, "Even if we're wrong about the exact size of the flour bag, the cake will still be ruined if you add a cup of salt."

The Bottom Line

This paper is a major victory for the "Big Bang Nucleosynthesis" method. It proves that looking at the ancient Helium in the universe is a more powerful tool than looking at the afterglow of the Big Bang (CMB) for finding these specific heavy particles.

In simple terms: The universe's oldest cake is the ultimate lie detector. It tells us that if these Heavy Axions exist, they must be incredibly short-lived ghosts that vanish before they can mess up the recipe. If they lingered even a tiny bit longer, the universe would taste very different today.

Technical Summary

Here is a detailed technical summary of the paper "New Bounds on Heavy QCD Axions from Big Bang Nucleosynthesis" by Jung et al.

1. Problem Statement

The paper addresses the constraints on heavy QCD axions (with masses $m_a \gtrsim 300$ MeV) using Big Bang Nucleosynthesis (BBN).

• Context: The QCD axion solves the Strong CP problem. While standard axion models predict very light masses, "heavy" axion models exist where the axion mass is generated by additional mechanisms (e.g., extra confinement groups or mirror QCD).
• Physics Mechanism: Heavy QCD axions couple dominantly to gluons ($\mathcal{L} \supset \frac{\alpha_s}{8\pi} \frac{a}{f_a} G\tilde{G}$). If kinematically allowed ($m_a > 2m_\pi$), they decay predominantly into hadrons (pions, kaons, nucleons).
• The Conflict: If these axions decay during or shortly after the neutron freeze-out epoch ($t \sim 1$ s), the injected hadrons interact via the strong force with the primordial plasma. These strong interactions ($\sigma \sim \text{mb}$) are orders of magnitude larger than the weak interactions governing the standard neutron-to-proton ($n/p$) ratio. Even a small abundance of decaying axions can drastically alter the $n/p$ ratio, thereby changing the primordial Helium-4 ($^4$He) abundance, which is precisely measured.
• Gap: Previous studies often focused on axion-like particles (ALPs) decaying to photons or neglected specific hadronic channels (like $K_L$) and secondary scattering effects. This paper provides the first comprehensive BBN analysis specifically for heavy QCD axions.

2. Methodology

The authors employ a rigorous numerical framework to calculate the modification of the $n/p$ ratio ($X_n$) caused by axion decays.

A. Thermal Production and Abundance

• Initial Conditions: They assume a high reheating temperature allowing axions to thermalize with the Standard Model (SM) sector.
• Freeze-out: Axions decouple while relativistic. The initial abundance $Y_a = n_a/s$ is calculated across the QCD deconfinement-confinement transition ($T \sim 1-2$ GeV).
• Uncertainty Handling: They define a range $[Y_a^{\text{min}}, Y_a^{\text{max}}]$ to account for uncertainties in the freeze-out temperature across the phase transition, finding the bounds are largely insensitive to the precise value of $Y_a$.

B. Axion Decay and Hadron Yields

• Mass Regimes:
  • $m_a < 2$ GeV: Uses a data-driven method based on Chiral Perturbation Theory and Vector Meson Dominance, constrained by experimental data (Ref. [75]).
  • $m_a > 2$ GeV: Uses Perturbative QCD (pQCD) at NNLO for the total width. For branching fractions, they utilize hadronization generators PYTHIA 8.306 and Herwig 7.3.0.
• Key Observation: While PYTHIA and Herwig show discrepancies in baryon yields (up to $\mathcal{O}(1)$), the final lifetime bound is robust because it is logarithmically dependent on the number of hadrons.

C. Hadronic Interactions (Major Improvements)

The paper introduces critical updates to the treatment of hadronic injections, which are essential for accuracy:

1. $K_L$ Dynamics: Unlike charged pions/kaons, $K_L$ are neutral and do not thermalize electromagnetically. The authors track the momentum distribution of $K_L$ from axion decay and calculate elastic scattering ($N K_L \to N K_L$) to "reshape" the spectrum before they decay or undergo inelastic scattering.
2. Secondary Hadrons: They include secondary hadrons produced from the scattering of primary injected hadrons (e.g., $n K^- \to \pi^- \Lambda \to p \pi^- \pi^-$). They demonstrate that neglecting these leads to incorrect net conversion rates (e.g., a process that looks like $n \to p$ might actually result in a net $p \to n$ due to secondary pions).
3. Cross Sections: They perform a partial-wave analysis and simultaneous fitting of isospin-related scattering amplitudes for $N\pi$, $NK$, and $N\bar{N}$ processes, updating cross sections for $K^\pm$ and $K_L$ which were previously poorly constrained or oversimplified.

D. Numerical Evolution

• Boltzmann Equations: They solve a coupled system of equations for:
  • Background cosmology (Hubble rate, photon/neutrino temperatures, axion energy density).
  • Axion number density (including decay).
  • Hadron number densities ($\pi^\pm, K^\pm, K_L, N, \bar{N}$).
  • Neutron fraction $X_n$.
• Modified Cosmology: They account for the possibility of an axion-dominated era (if axions are massive and abundant) which alters the Hubble expansion rate and delays neutron freeze-out.
• Constraint: The constraint is derived by requiring the deviation in the primordial $^4$He mass fraction ($Y_p$) to be within the $2\sigma$observational limit ($|\delta Y_p|/Y_p < 2.45%$).

3. Key Contributions

1. First BBN Bounds on Heavy QCD Axions: Establishes the first robust constraints specifically for heavy QCD axions that decay hadronically, distinct from generic ALP constraints.
2. Methodological Advances in Hadronic Injection:
   • Properly treats the non-thermalized $K_L$ spectrum and its elastic scattering.
   • Includes secondary hadron production from scattering, showing it is essential for calculating the net $n/p$ conversion.
   • Updates hadronic cross sections using isospin relations and modern fitting techniques.
3. Robustness Analysis: Demonstrates that the derived lifetime bound is insensitive to:
   • Uncertainties in hadronic cross sections (factor of 2 changes result in only $\sim 1\%$ shift in the bound).
   • Uncertainties in branching fractions (differences between PYTHIA and Herwig).
   • Initial axion abundance ($Y_a$).
   • This insensitivity arises because the bound depends logarithmically on these parameters.

4. Results

• Lifetime Bound: The authors derive a robust upper bound on the axion lifetime:
  $$\tau_a \lesssim 0.017 \text{ s} \quad (\text{for } m_a > 300 \text{ MeV})$$
  This bound is nearly flat across the mass range $0.3 \text{ GeV} < m_a < 100 \text{ GeV}$.
• Comparison with Other Probes:
  • Colliders: Limited to $\tau_a \lesssim 10^{-9}$ s (detector geometry).
  • Beam Dumps: Limited to $\tau_a \lesssim 10^{-6}$ s.
  • CMB ($N_{\text{eff}}$): Probes $\tau_a \gtrsim 0.1$ s.
  • Significance: The BBN bound fills the critical gap between $10^{-2}$s and$0.1$ s, a region previously unconstrained for hadronically decaying particles.
• Parameter Space ($m_a$ vs. $f_a$):
  • The lifetime bound translates to constraints on the decay constant $f_a$.
  • For $m_a \sim 1$ GeV, the bound excludes $f_a$ values roughly between $10^7$and$10^{11}$ GeV (depending on the specific model assumptions for the axion mass generation).
  • The BBN constraint is stronger than projected future CMB constraints (e.g., from CMB-S4) for masses where hadronic decays are kinematically allowed.

5. Significance

• Exclusion of Heavy Axion Models: The results rule out a significant portion of the parameter space for heavy QCD axion models that solve the Strong CP problem without introducing a "quality problem" (where global symmetries are broken by quantum gravity).
• General Applicability: The improved methodology for treating hadronic injections (specifically $K_L$ dynamics and secondary scattering) is applicable to other long-lived particle (LLP) scenarios, such as gravitinos, dark photons, or heavy neutral leptons, potentially tightening constraints on those models as well.
• Theoretical Insight: The paper highlights that the "waterfall" effect (where the $n/p$ ratio rapidly relaxes from a quasi-stable value set by hadronic interactions back to the SM value) is the dominant mechanism setting the bound, making the result robust against model details.

In summary, this paper establishes that Big Bang Nucleosynthesis is the most powerful probe for heavy QCD axions with lifetimes in the range of $0.01$to$0.1$ seconds, excluding models that would otherwise be viable based on collider or CMB data alone.