Thermal evolution of dark matter and gravitational-wave production in the early universe from a symplectic glueball model

This paper investigates a dark matter model based on a symplectic gauge group, analyzing its thermodynamic properties near the confinement phase transition and exploring the resulting gravitational-wave production and relic abundances in the early universe.

The Gist

The Big Picture: A Hidden Universe of "Dark Glue"

Imagine the universe is like a giant, bustling city. We know a lot about the "visible" part of the city—the people, buildings, and cars (which represent normal matter like atoms and stars). But we also know there is a massive, invisible "ghost city" occupying about 85% of the space. This is Dark Matter. We can't see it, but we know it's there because its gravity holds galaxies together, like an invisible scaffolding.

For a long time, scientists have wondered: What is this ghost city made of?

This paper proposes a new theory. Instead of imagining dark matter as a single, mysterious particle (like a tiny, invisible marble), the authors suggest it might be made of clumps of "dark glue."

The Analogy: The "Dark Glue" Factory

To understand this, let's look at how our own visible world works.

• The Visible World: Inside an atom, there are particles called quarks held together by "glue" (particles called gluons). This glue is so strong that you can never pull a single piece of it apart. If you try to pull them, the energy creates new particles. The result is that you only ever see "clumps" of glue and quarks stuck together, called hadrons (like protons and neutrons).
• The Dark World: The authors suggest there is a parallel "Dark Sector" that works exactly the same way, but it has its own invisible glue and its own invisible particles. However, this dark glue doesn't stick to our visible atoms at all. It only talks to us through gravity.

In this model, the dark matter we see in the sky isn't a single particle; it's a glueball—a massive, heavy ball made entirely of this invisible dark glue.

The Experiment: Simulating a Cosmic Phase Change

The authors didn't just guess; they ran a massive computer simulation to see how this "Dark Glue" behaves when the universe is hot versus when it is cold.

Think of water.

• Hot Water (Steam): When water is hot, the molecules fly around freely. They are a gas.
• Cold Water (Ice): When it gets cold, the molecules lock together into a rigid crystal structure.

The universe went through a similar change. In the very early, hot universe, the "dark glue" was a hot, chaotic soup (like steam). As the universe expanded and cooled, it hit a critical temperature and suddenly "froze" into solid glueballs (like ice).

The authors used a technique called Lattice QCD (which is like building a giant 3D grid of pixels to simulate the laws of physics) to calculate exactly how this transition happens for their specific type of dark glue (based on a mathematical group called Sp(2)).

Key Findings from the Simulation

1. It's a Sudden Snap, Not a Slow Melt:
   When the dark glue cooled down, it didn't slowly turn into a solid. It happened all at once, like a sudden snap. In physics terms, this is a first-order phase transition.

   • The Analogy: Imagine a room full of people dancing wildly. Suddenly, a siren goes off, and everyone instantly freezes into a rigid pose. That sudden change releases a burst of energy.

2. The "Latent Heat" Explosion:
   Because the transition was so sudden, it released a huge amount of energy (called latent heat). The authors calculated exactly how much energy was released. This is important because that burst of energy didn't just disappear; it shook the fabric of space-time.

3. Ripples in Space-Time (Gravitational Waves):
   When that "snap" happened in the early universe, the sudden release of energy and the collision of the "freezing" bubbles created ripples in space-time. These are Gravitational Waves.

   • The Analogy: Imagine dropping a giant stone into a calm pond. The splash creates waves that travel outward. The authors calculated the "frequency" (pitch) of these waves. They found that these waves would have a pitch that future space-based detectors (like LISA) might be able to hear. It's like the universe is humming a specific note from its birth, and this model predicts what that note sounds like.

4. Why This Model is Special:
   Most previous studies looked at a different type of math (called SU(N)) for dark glue. This paper looks at a slightly different math (Sp(2)).

   • The Difference: In the "standard" dark glue models, there are some particles that are "odd" (like a left-handed glove). In this new Sp(2) model, all the particles are "even" (like a pair of matching socks). This changes how the dark matter might behave and how long it lasts. The authors found that despite this difference, the "freezing" process still happens in a very similar, explosive way.

The Conclusion: A Viable Candidate

The paper concludes that this "Dark Glueball" model is a very strong candidate for what Dark Matter actually is.

• It explains why dark matter is heavy and clumpy.
• It explains why it doesn't interact with light (it's made of invisible glue).
• It predicts a specific "sound" (gravitational wave signature) that we might be able to detect in the near future.

The authors admit that while they have calculated the "thermodynamics" (the heat and pressure) perfectly using their supercomputer, some details about how these glueballs might eventually decay or interact with our world are still a bit fuzzy. However, the core finding is solid: If dark matter is made of this specific type of "dark glue," the early universe would have made a loud, detectable noise that we might finally be able to hear.

Technical Summary

Technical Summary: Thermal Evolution of Dark Matter and Gravitational-Wave Production in the Early Universe from a Symplectic Glueball Model

Problem and Motivation
The nature of dark matter remains a fundamental open problem in physics. While observational evidence strongly supports its existence, its particle identity is unknown. Theoretical models proposing dark matter as a composite state of a strongly coupled, non-Abelian gauge theory (analogous to Quantum Chromodynamics, QCD) are appealing because they naturally generate mass scales without fine-tuning and can provide stable, neutral, and self-interacting candidates. Previous literature has extensively studied such "dark glueball" models based on $SU(N)$ gauge groups. However, the behavior of theories based on other gauge groups, specifically symplectic groups $Sp(N)$, remains less explored. This work investigates the viability of a dark matter model based on the compact symplectic group $Sp(2)$ (the real form of $Sp(4, \mathbb{C})$ intersected with $U(4)$). A primary motivation is to determine if this theory undergoes a first-order confinement/deconfinement phase transition, which could produce a detectable stochastic gravitational-wave background, and to characterize its equation of state (EoS) non-perturbatively.

Methodology
The authors perform a non-perturbative study of the $Sp(2)$ Yang-Mills theory using lattice regularization on a Euclidean, isotropic, hypercubic lattice.

• Simulation Setup: The theory is defined by the Wilson action with link variables $U_\mu(x)$ in the defining representation of $Sp(2)$. Simulations are conducted at finite temperature $T = 1/(a_{lat} N_\tau)$, where $a_{lat}$ is the lattice spacing and $N_\tau$ is the extent in the Euclidean time direction.
• Algorithm: The authors utilize a Markov-chain Monte Carlo (MCMC) algorithm combining local heat-bath and over-relaxation updates. To ensure ergodicity and reduce autocorrelation, updates are applied to embedded $SU(2)$ subgroups within the $Sp(2)$ matrices.
• Scale Setting: The lattice scale is set non-perturbatively by determining the critical coupling $\beta_c$ for the deconfinement transition for various $N_\tau$ values. This is achieved by analyzing the fourth-order Binder cumulant of the Polyakov loop and performing finite-size scaling. The resulting $\beta$ vs. $a_{lat}$ relation is interpolated and matched with a two-loop perturbative expression in an improved lattice scheme to cover a wider range of couplings.
• Equation of State (EoS): The thermodynamic quantities (pressure $p$, energy density $\epsilon$, entropy density $s$, and trace of the energy-momentum tensor $\Delta$) are computed using the "integral method." This involves integrating the difference between the average plaquette trace at finite temperature and at zero temperature with respect to the coupling $\beta$.
• Continuum Extrapolation: Results are obtained at several lattice spacings ($N_\tau = 5, 6, 7$) and extrapolated to the continuum limit ($a_{lat} \to 0$) by fitting data as a function of $1/N_\tau^2$. Finite-volume effects are suppressed by using large spatial volumes ($N_s/N_\tau \gtrsim 4$).

Key Contributions and Results

1. Phase Transition Nature: The study confirms that the $Sp(2)$ Yang-Mills theory undergoes a first-order deconfinement phase transition at a finite critical temperature $T_c$. This is evidenced by the crossing of Binder cumulant curves for different volumes and the presence of a non-zero latent heat.
2. Thermodynamic Quantities: The authors provide the first continuum-extrapolated results for the EoS of $Sp(2)$ theory in the range $0.7 \le T/T_c \le 2.75$.
   • The trace of the energy-momentum tensor $\Delta/T^4$ exhibits a peak near $T_c$, characteristic of a first-order transition.
   • The pressure, energy density, and entropy density are computed and found to approach the Stefan-Boltzmann limit slowly, indicating significant interactions even in the deconfined phase.
   • The speed of sound $c_s^2$ is derived, showing behavior qualitatively and quantitatively similar to $SU(3)$ Yang-Mills theory.
3. Latent Heat: A finite latent heat $L_h$ is determined. The continuum-extrapolated value is found to be $L_h/T_c^4 \simeq 1.69(12)$. This value is slightly larger than that of $SU(3)$($\simeq 1.4$) but consistent with scaling expectations based on the dimension of the gauge group (10 for $Sp(2)$ vs. 8 for $SU(3)$).
4. Confining Phase Modeling: In the low-temperature confining phase, the EoS is well-described by a gas of non-interacting massive glueballs. The lightest glueball is a scalar ($J^{PC}=0^{++}$) with a mass $m \approx 3.58 \sqrt{\sigma}$. The results show an exponential growth in the density of states near $T_c$, consistent with a Hagedorn-like spectrum.
5. Gravitational-Wave Production: The authors estimate the gravitational-wave (GW) spectrum generated by the first-order phase transition.
   • Using the calculated latent heat and an estimated interface tension, they derive the transition strength parameter $\alpha \simeq 0.62(4)$, classifying the transition as "intermediate" to "strong."
   • The GW spectrum is modeled considering contributions from bubble collisions, sound waves, and hydrodynamic turbulence.
   • The dominant contribution is predicted to come from sound waves, with a peak frequency in the milli-Hertz (mHz) range. This places the signal potentially within the sensitivity reach of future space-based detectors like LISA, DECIGO, and BBO, assuming the phase transition occurs at temperatures around hundreds of GeV.

Dark Matter Interpretation and Significance
The paper interprets the lightest scalar glueball ($0^{++}$) as a dark matter candidate.

• Stability: Unlike $SU(N)$ theories with $N \ge 3$, the $Sp(2)$ group does not support states with negative charge-conjugation ($C=-1$). Consequently, the mechanism proposed in some $SU(N)$ models where a $C=-1$ state is stabilized by discrete symmetry (allowing it to be dark matter while $C=+1$ states decay) is not applicable here. Stability in the $Sp(2)$ model relies on the absence of lighter states to decay into and potentially suppressed decay channels via higher-dimensional operators (e.g., to gravitons).
• Relic Abundance: The authors discuss the thermal evolution of the dark sector, including the "cannibalism" scenario (3-to-2 self-annihilation) which can alter the relic density. They provide order-of-magnitude estimates for the dark matter abundance, noting that the results are sensitive to the reheating temperature and the ratio of dark-to-visible sector temperatures.
• Significance: The paper claims that the $Sp(2)$ theory is a viable, "economic" dark matter model depending on a single dimensionful parameter. Its primary significance lies in providing the first rigorous, non-perturbative lattice determination of the EoS and latent heat for a symplectic gauge theory. This establishes that such theories can support strong first-order transitions capable of generating observable gravitational waves. The authors emphasize that while their lattice results for the EoS are robust, predictions for GW spectra and dark matter relic densities rely on semi-quantitative estimates for non-equilibrium parameters (like bubble wall velocity and nucleation rates) which are not directly calculable in their current equilibrium lattice setup.

Limitations and Future Work
The authors explicitly state that their GW and dark matter abundance predictions are benchmark examples rather than first-principle quantitative predictions, as they rely on estimates for dynamical quantities (bubble wall velocity, interface tension, nucleation rates) that require real-time simulations or effective field theory approaches beyond the scope of the current equilibrium study. Future work could involve refining the scale setting, extending the temperature range, and performing dedicated non-perturbative calculations for interface tension and bubble nucleation rates. Additionally, the absence of $C=-1$ states in $Sp(2)$ precludes certain stability mechanisms available in $SU(N)$ models, a feature that distinguishes this model's phenomenology.