⚛️ phenomenology

Gravitational Waves from Dimension-6 Assisted Peccei - Quinn Phase Transitions

This paper investigates how adding a dimension-6 operator to the minimal KSVZ axion model can trigger a first-order Peccei-Quinn phase transition, producing observable gravitational wave signals and specific dark matter phenomenology within a framework supported by ultraviolet completions.

Original authors: Nico Benincasa, Kristjan Müürsepp

Published 2026-08-17

📖 3 min read🧠 Deep dive

CC BY 4.0

Original authors: Nico Benincasa, Kristjan Müürsepp

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

Imagine the universe as a giant, invisible stage where the laws of physics play out like a cosmic drama. For decades, physicists have been puzzled by a specific glitch in this script: why does the strong nuclear force, which holds atomic nuclei together, seem to ignore a fundamental rule about time-reversal symmetry? This mystery, known as the "Strong CP problem," has led scientists to propose a clever fix involving a new, ghostly particle called the "axion." Think of the axion as a cosmic thermostat that automatically adjusts itself to zero out this glitch, keeping the universe stable. But here's the catch: we've never actually seen an axion. They are the "dark matter" of particle physics—everywhere, perhaps, but invisible. To find them, scientists look for clues in how the universe evolved, specifically during its hottest, earliest moments. When the universe cooled down, it underwent "phase transitions," much like water freezing into ice. Usually, these transitions are smooth, but if they happen suddenly and violently (a "first-order" transition), they can create ripples in spacetime itself, known as gravitational waves. Detecting these waves would be like hearing the echo of the universe's birth, potentially revealing the existence of these elusive axions.

This paper takes a fresh look at a specific model for these axions, called the KSVZ model, and asks a simple question: what if we tweak the rules just a little bit? The authors propose adding a tiny, extra ingredient—a "dimension-6 operator"—to the mathematical recipe that describes how the axion field behaves. In the standard version of this model, the universe's transition to create axions is usually a smooth, boring slide. However, the authors show that with this new ingredient, the transition can become a dramatic, explosive event. They ran detailed computer simulations to see if this explosion would be strong enough to create gravitational waves that our current and future detectors could actually hear.

The results are a mix of exciting possibilities and strict boundaries. The team found that for the transition to be explosive enough to create a detectable signal, the universe must have cooled in a very specific way, and the axion's "strength" (a value called the decay constant, $f_a$ ) must fall within a narrow range. If the axion is too weak or too strong, the transition remains too quiet to be heard. Interestingly, they discovered that for the axion to make up all the dark matter we see today, the transition would likely be too high-pitched for our current detectors to catch, requiring future, more sensitive instruments. However, if the axion only makes up a small part of the dark matter, the signal could be loud enough for today's detectors, like LIGO or Virgo, to potentially spot. The paper also rules out a scenario where the transition happens so slowly that the universe gets stuck in a "frozen" state, unable to complete the change; their calculations show that in those cases, the universe would remain trapped in a state that doesn't match what we observe today.

To make sure their idea isn't just a mathematical trick, the authors also built "ultraviolet completions"—real-world scenarios involving other heavy particles (like heavy neutrinos or extra axions) that could naturally create the extra ingredient they needed. They showed that these scenarios are physically possible and could even explain why we have the right amount of dark matter in the universe. In short, this paper suggests that by listening for a specific "crack" in the fabric of spacetime, we might finally catch a glimpse of the axion, solving one of physics' oldest mysteries and revealing a hidden layer of the universe's history.

Technical Summary: Gravitational Waves from Dimension-6 Assisted Peccei–Quinn Phase Transitions

Problem Statement
The strong CP problem remains a primary motivation for physics beyond the Standard Model (SM), typically addressed by the Peccei–Quinn (PQ) mechanism, which introduces a global $U(1)_{PQ}$ symmetry and a resulting axion field. In the minimal Kim–Shifman–Vainshtein–Zakharov (KSVZ) realization, the PQ phase transition (PT) is generally a second-order transition or a crossover, precluding the generation of a detectable stochastic gravitational wave (GW) background. Furthermore, while the QCD axion is a viable dark matter (DM) candidate, its relic abundance via the misalignment mechanism depends heavily on the PQ breaking scale ( $f_a$ ) and the thermal history. The authors seek to determine if extending the minimal KSVZ model with a non-renormalizable dimension-6 operator can induce a strong first-order phase transition (FOPT), thereby producing observable GWs while remaining consistent with DM constraints.

Methodology
The authors analyze an extension of the minimal KSVZ model where the scalar potential includes a dimension-6 operator invariant under $U(1)_{PQ}$ . The study proceeds through the following steps:

Effective Potential Construction: The effective potential $V_{\text{eff}}$ is constructed including tree-level terms, zero-temperature 1-loop Coleman-Weinberg corrections, finite counterterms (to fix the vacuum expectation value and curvature at the tree-level minimum), and thermal corrections. The tree-level potential is defined as $V_{\text{tree}}(\phi) = \frac{m^2}{2}\phi^2 + \frac{\lambda}{4}\phi^4 + \frac{1}{8\Lambda^2}\phi^6$ , where $\Lambda$ is the cutoff scale.
Phase Transition Dynamics: The authors compute the critical temperatures, nucleation temperatures ( $T_n$ ), and percolation temperatures ( $T_p$ ). They distinguish between thermal tunneling (O(3) symmetry) and zero-temperature quantum tunneling (O(4) symmetry). Special attention is paid to scenarios involving Symmetry Non-Restoration (SNR) and multi-step transitions.
Gravitational Wave Calculation: The GW power spectrum is calculated based on three sources: bubble collisions, sound waves in the plasma, and magnetohydrodynamic turbulence. The authors utilize state-of-the-art numerical fits recommended by the LISA Cosmology working group. They evaluate the signal against the sensitivity curves of current (aLIGO, aVirgo, KAGRA) and future (LISA, ET, CE, etc.) detectors, as well as constraints from the effective number of relativistic species ( $\Delta N_{\text{eff}}$ ).
Parameter Space Scan: A numerical scan is performed over the parameters $m/\Lambda$ , $\lambda$ , and the Yukawa coupling $y_F$ for cutoff scales $\Lambda \in [10^{10}, 10^{12}]$ GeV. The scan is constrained by the requirements that the potential possesses two minima separated by a barrier at $T=0$ and that the symmetry-breaking minimum is the global one.
UV Completions: The paper presents three renormalizable UV completions that generate the effective dimension-6 operator: a heavy real scalar singlet, a Type-I seesaw extension involving a pseudo-Majoron, and a two-axion system with a hidden confining gauge group.

Key Contributions and Results

Induction of First-Order Transitions: The inclusion of the dimension-6 operator allows for a tree-level barrier between the symmetric ( $\phi=0$ ) and broken ( $\phi=f_a$ ) vacua even at zero temperature. This enables strong FOPTs in parameter regions where the minimal KSVZ model would only exhibit a crossover.
Analytical Constraints: The authors derive analytical bounds on the Yukawa coupling $y_F$ required to maintain a local minimum at the origin and a global minimum at $f_a$ at zero temperature. These bounds are visualized in the parameter space of $\lambda$ vs. $m/\Lambda$ .
Non-Standard Thermal Histories: The study identifies scenarios of Symmetry Non-Restoration (SNR) where the PQ symmetry remains broken at high temperatures due to a negative quartic coupling. This can lead to a two-step PT: a transition from the broken phase to the symmetric phase, followed by a return to the broken phase at lower temperatures. If both steps are first-order, this could produce a distinctive double-peaked GW signal.
GW Observability:
- For lower cutoff scales ( $\Lambda \sim 10^{10}$ GeV), the peak GW frequency falls within the sensitivity range of current ground-based detectors (aLIGO, aVirgo, KAGRA), provided the transition is sufficiently strong.
- For higher scales ( $\Lambda \sim 10^{12}$ GeV), the peak frequency shifts to higher values, potentially requiring future high-frequency detectors.
- The paper finds that for the parameter space where the PQ transition completes successfully, the ratio $f_a/\Lambda$ is tightly constrained to $0.42 \lesssim f_a/\Lambda \lesssim 2/3$ .
Dark Matter Phenomenology:
- In the post-inflationary scenario, the axion relic abundance is determined by the misalignment angle and $f_a$ .
- The observed DM abundance ( $\Omega_{\text{DM}}h^2 \approx 0.12$ ) is reproduced for $\Lambda \simeq (3-5) \times 10^{11}$ GeV.
- For $\Lambda < 3 \times 10^{11}$ GeV, the axion constitutes only a fraction of the DM; for $\Lambda > 5 \times 10^{11}$ GeV, the axion overproduces DM (unless topological defects are considered, which would lower the saturation scale).
- The UV completions (specifically the Majoron and two-axion models) offer alternative DM candidates that can account for the full abundance in regions where the axion alone underproduces it.
Quantum Tunneling Failure: The authors demonstrate that for parameter points where thermal percolation fails, zero-temperature O(4) quantum tunneling also fails to complete the transition. The bounce action $S_4$ remains significantly above the critical threshold ( $S_4 \gtrsim 296$ vs. $S_{\text{crit}} \sim 70-110$ ), implying these regions are cosmologically excluded as the universe would remain trapped in the false vacuum.

Significance
The paper claims that higher-dimensional operators provide a viable pathway to strong FOPTs within the KSVZ framework, a feature absent in the minimal model. This creates a "multimessenger" perspective where GW detectors can complement traditional axion search experiments. The work motivates the use of Effective Field Theories (EFTs) to trigger strong PQ transitions and provides specific UV completions that link the dimension-6 operator to neutrino physics and multi-axion scenarios. The authors emphasize that their model-independent approach is justified by the distinct phenomenological signatures of these UV completions, particularly the potential for multi-peaked GW signals from SNR scenarios and the ability to explain DM abundance in conjunction with GW observability.

More like this