NLO observables for QCD-like theories and application to pion dark matter

This paper derives next-to-leading order chiral perturbation theory expressions for QCD-like theories with non-degenerate fermion masses, applies them to extract low-energy constants from $Sp(4)$ lattice data, and demonstrates the critical role of these corrections in refining the viable parameter space for pion dark matter scenarios.

The Gist

Imagine the universe is built from a giant, invisible Lego set. For decades, physicists have been trying to understand the rules of how these Lego bricks snap together to form everything we see, including the mysterious "Dark Matter" that holds galaxies together.

This paper is like a new, highly detailed instruction manual for a specific, exotic type of Lego set that isn't used in our everyday world (the Standard Model) but might exist in the hidden sectors of the universe.

Here is the story of what the authors did, explained simply:

1. The Problem: The "Too Heavy" Lego Bricks

In the standard world, the forces that hold particles together are like a spring. When you pull them apart, they snap back. Physicists have a great way to describe these springs when they are light and easy to stretch (called "Leading Order" or LO).

However, in some theories about Dark Matter, these "springs" are very stiff and heavy. When you try to use the simple instructions (LO) to predict how these heavy bricks bounce off each other, the math breaks down. It's like trying to predict the flight of a bowling ball using the same simple rules you use for a ping-pong ball. You need a more complex set of rules that accounts for the extra weight and stiffness. This is what the authors call Next-to-Leading Order (NLO) corrections.

2. The Goal: Writing the "Advanced" Manual

The authors wanted to write these advanced rules for two specific types of exotic Lego sets:

• The "Pseudoreal" Set (Sp(4)): A complex, twisted arrangement of bricks.
• The "Real" Set (SO(4)): A slightly different, mirrored arrangement.

They calculated the exact formulas for how heavy these "Dark Pions" (the Lego bricks) would be, how they decay, and most importantly, how they crash into each other.

3. The Detective Work: Using a "Simulation" to Find the Constants

Here is the tricky part: The advanced manual has several "magic numbers" (called Low-Energy Constants or LECs) that the math can't predict on its own. These numbers depend on the specific material of the Lego bricks.

To find these numbers, the authors didn't build a physical model. Instead, they used supercomputer simulations (called Lattice QCD) that act like a virtual laboratory.

• They took data from other scientists who had already simulated these exotic Lego sets on a computer grid.
• They treated the computer data like a puzzle. They plugged the data into their new, complex formulas.
• By adjusting the "magic numbers" until the formulas matched the computer simulation perfectly, they successfully calibrated their manual.

4. The Big Discovery: The "Crash Test" Results

Once they had their calibrated manual, they ran a "crash test" to see how these Dark Matter particles would interact with each other in the real universe.

• The Old View (Simple Rules): If you used the simple rules, you might think Dark Matter could be a certain size and still fit the observations of our universe.
• The New View (Complex Rules): When they applied their new, advanced rules, the results changed significantly. The "crash test" showed that the particles interact much more strongly than previously thought.

The Analogy: Imagine you are trying to park a car in a tight spot.

• Simple Rules: You think, "I can squeeze in if I turn the wheel a little."
• Advanced Rules: You realize, "Oh, the car is actually much wider than I thought, and the ground is slippery. If I turn the wheel that much, I'll crash into the wall."

The authors found that for many theories of Dark Matter (specifically the "SIMP" scenario), the "crash" happens much sooner than expected. This means the "safe parking spots" (the viable parameter space) where Dark Matter could exist are much smaller and more restricted than we thought.

5. Why This Matters

The paper concludes that if we want to understand Dark Matter, we can't rely on the "back-of-the-envelope" calculations anymore. We need the full, complex math.

• For the "Pseudoreal" Set: They successfully calibrated the rules and showed that the "crash" limits are tighter.
• For the "Real" Set: They provided the formulas, but noted that we don't have enough computer simulation data yet to fully calibrate the "magic numbers" for this specific set.

In short: The authors built a more accurate map for a hidden part of the universe. They found that the terrain is rougher and the boundaries are tighter than the old maps suggested, forcing us to rethink where Dark Matter can actually live.

Technical Summary

Problem Statement
QCD-like theories with fermions in real or pseudoreal representations of the gauge group are central to various Beyond-the-Standard-Model (BSM) scenarios, including composite Higgs models and pion dark matter (specifically Strongly Interacting Massive Particles, or SIMPs). While the low-energy dynamics of these theories are described by Chiral Effective Field Theory (ChEFT), existing literature largely focuses on the mass-degenerate quark limit or leading-order (LO) approximations. However, for phenomenological applications such as SIMP dark matter, the relevant parameter space often involves relatively large pion masses where LO approximations are insufficient. Furthermore, viable dark matter scenarios may require non-degenerate quark masses to enlarge the parameter space or satisfy observational constraints. A critical gap exists in the precise determination of Next-to-Leading Order (NLO) corrections for theories with non-degenerate masses in real and pseudoreal representations, particularly regarding the low-energy constants (LECs) which are unknown in BSM contexts and must be extracted from lattice simulations.

Methodology
The authors construct the NLO chiral Lagrangian for $N_F=2$ fermions in both pseudoreal ($SU(4)/Sp(4)$) and real ($SU(4)/SO(4)$) representations of the gauge group, explicitly accounting for non-degenerate quark masses ($m_u \neq m_d$).

1. Theoretical Framework: Using the CCWZ formalism, the authors derive NLO expressions ($O(p^4)$) for physical observables: pion masses, quark condensates, decay constants, and $2 \to 2$ meson scattering amplitudes. The power counting scheme assigns chiral dimension 2 to fermion masses, allowing for a consistent expansion.
2. Lattice Data Integration: To fix the unknown NLO LECs for the $Sp(N_c=4)$ theory with fundamental fermions, the authors utilize recent lattice data from Refs. [33] and [34].
   • Mass spectrum and decay constant data (non-degenerate case) from Ref. [33] are used, with infinite-volume extrapolations performed where feasible.
   • Scattering length data (degenerate case) from Ref. [34] (Lüscher method) are incorporated.
3. Fitting Procedure: A Bayesian Markov Chain Monte Carlo (MCMC) fit is performed to determine linear combinations of the renormalized LECs ($\tilde{L}_i$). The fit accounts for uncertainties in both dependent and independent variables (PCAC mass ratios and dimensionless mass ratios) and assesses discretization errors by analyzing data at different lattice spacings ($\beta$ values).
4. Phenomenological Application: The fitted LECs are applied to calculate the non-relativistic $2 \to 2$ self-scattering cross-sections for pion dark matter in the $SU(4)/Sp(4)$ and $SU(4)/SO(4)$ theories. These results are compared against astrophysical constraints, specifically the Bullet Cluster bound on dark matter self-interactions.

Key Contributions

• Analytical Formulas: The paper provides the complete set of NLO formulas for masses, condensates, decay constants, and scattering amplitudes for $N_F=2$ theories with non-degenerate quark masses in both real and pseudoreal representations. This extends previous work which was limited to degenerate masses or LO.
• Extraction of LECs: The authors present the first determination of NLO LECs for the $Sp(N_c=4)$ gauge theory using lattice data that includes mass splitting. They identify four linear combinations of LECs ($\tilde{L}_1, \tilde{L}_2, \tilde{L}_3, \tilde{L}_4$) that can be reliably constrained by current data.
• Mass Splitting Analysis: The study explicitly demonstrates how NLO corrections generate mass splitting among pion multiplets in the non-degenerate case, a feature absent at LO. For the $Sp(4)$ case, the fit robustly predicts the singlet pion to be lighter than the four-plet.
• Refined Phenomenology: By applying the fitted LECs to SIMP dark matter scenarios, the authors quantify the impact of NLO corrections on the self-interaction cross-section.

Results

• LEC Values: The fit yields specific values for the linear combinations of LECs at $\beta=7.2$ (closest to the continuum limit). The results indicate that $\tilde{L}_4$ is negative, implying the singlet pion is the lightest state in the $Sp(4)$ theory.
• Scattering Cross-Sections: The inclusion of NLO corrections significantly alters the predicted $2 \to 2$ scattering cross-sections. For typical SIMP parameters ($M_\pi \approx 0.2$ GeV), the NLO cross-section deviates from the LO result by approximately 20% at $M_\pi/F_\pi \approx 1.6$ and up to 100% at higher ratios.
• Parameter Space Constraints: When applied to the SIMP scenario, the NLO corrections tighten the constraints on the viable parameter space. The authors find that for the $Sp(4)$ theory with $N_c=4$, the region satisfying both the relic abundance (via $3 \to 2$ freeze-out) and the Bullet Cluster self-interaction bounds is severely restricted or potentially eliminated when NLO corrections are included, compared to the LO analysis.
• Real Representation ($SO(4)$): For the $SU(4)/SO(4)$ theory, where no specific lattice data for non-degenerate masses was available, the authors modeled the LECs based on the $Sp(4)$ results. They found that the LO mass hierarchy (where a doublet of pions is lightest) generally holds, though NLO corrections can induce hierarchy flips in specific regions of parameter space.

Significance and Claims
The paper claims to provide a more accurate and reliable description of QCD-like theories with (pseudo)real representations by moving beyond the LO approximation and incorporating non-degenerate mass effects. The authors emphasize that NLO corrections are not merely perturbative refinements but are crucial for determining the viable parameter space of dark matter models like SIMPs. They argue that previous LO-based phenomenological studies may have overestimated the allowed parameter space.

The work serves as a bridge between lattice QCD simulations and BSM phenomenology, demonstrating that lattice data can be used to fix EFT parameters for theories relevant to dark matter and composite Higgs models. The authors modestly note that their analysis relies on controlled approximations due to limited lattice data (e.g., inability to fit all individual LECs, lack of continuum extrapolation for some datasets, and restriction to $N_F=2$). They conclude that their results motivate further lattice calculations, particularly for real representations and higher-order (NNLO) analyses, to further refine the phenomenological viability of these dark matter scenarios.