Dark Phase Transition: Testing the Model with Pulsar Timing Arrays
This work is motivated by recent hints from pulsar timing array (PTA) collaborations of a stochastic gravitational-wave background. Phase transitions in a dark Abelian Higgs sector—a minimal gauge theory with spontaneous symmetry breaking—are explored. Using dimensionally reduced high-temperature effective field theory, the impact of thermal resummation, higher-order matching corrections, and higher-dimensional operators on thermodynamics and the gravitational-wave signal is assessed. The parameter region favored by PTA lies near the boundary of the effective theory's validity, where contributions from higher operators grow; even within controlled limits the signal remains unfavorable compared to data, despite significant shifts from thermal corrections. Regimes of thermal and hydrodynamic coupling or decoupling between the dark and visible sectors are determined, and asymmetric dark matter freeze-out is consistent with the observed density and the coupling constants needed for strong transitions.
Context
Why it matters: Gravitational waves have opened a new window into the early Universe. While pulsars — spinning neutron stars, first discovered by Jocelyn Bell Burnell — revealed the existence of extreme objects, their ultra-precise 'ticks' now let us catch the trembling of spacetime from processes that occurred in the first moments after the Big Bang. Among these processes are phase transitions in dark sectors, where dark matter particles could be hiding. Moreover, these particles might influence the cooling of supernovae, imposing additional constraints.
Methods
How it was done: We constructed a three-dimensional effective field theory (EFT) for the dark abelian Higgs sector at finite temperature. Imagine studying a boiling cauldron, but instead of water, it's the Higgs field. To avoid counting each molecule, you use macroscopic parameters. Similarly, here we 'integrated out' fast thermal fluctuations, leaving only the slow modes that govern the dynamics of new-phase bubbles. This approach, developed in works on the Standard Model, allowed us to systematically include higher-order corrections and quantitatively assess uncertainties. Calculations went up to NNLO, including effects of a heavy fermion — a dark matter candidate whose mass suppresses its role at low energies. We also meticulously tracked entropy contributions from both sectors to properly compute the observable signal.
Results
What we found: We simulated the phase transition at temperatures from 1 to 100 MeV and computed the gravitational-wave spectrum in the nanohertz range — exactly where pulsar arrays are 'listening'. It turned out that for typical model parameters, the signal amplitude is orders of magnitude below the observed one. Even in the most 'violent' transitions, where bubbles of the new phase nucleate avalanche-like, predictions disagree with NANOGrav data at the 2σ level. Moreover, the parameter region favored by the data lies at the edge of our effective description's validity: there, higher-dimensional operators become important, and the high-temperature expansion begins to fray. We also found that the dark and visible sectors can be hydrodynamically decoupled: ordinary plasma barely participates in bubble motion, affecting the final loudness of the gravitational echo.
Implications
What this means for science: The simplest dark sector with U(1) gauge symmetry cannot explain the PTA signal without considerable strain. This strengthens the case for alternative sources, such as merging supermassive black holes or cosmic strings. However, our analysis does not close the book — it merely shows that more precise calculations and perhaps more complex models are needed, where the phase transition is stronger and predictions more robust. Meanwhile, the methodology based on dimensional reduction is becoming a benchmark for such studies, linking Standard Model physics and cosmology.
Future development
How this topic could develop: The next step is to include dimension-six operators in bubble nucleation and account for fluctuations around the critical bubble. Methods beyond the high-temperature expansion are already being developed, which will allow studying transitions with very strong supercooling. Additionally, lattice simulations could provide non-perturbative benchmarks, much like Vera Rubin once confirmed the existence of dark matter through galaxy observations.
Impact
Areas impacted: Cosmology of the early Universe and particle physics beyond the Standard Model. Gravitational-wave astronomy will also gain more precise templates for hunting exotic signals.
Next steps
Next steps: Build a perturbative theory with higher-dimensional operators and compare with non-perturbative lattice data; consider non-abelian gauge groups where the phase transition could be stronger.
Key open problems
Connection to unsolved problems in physics: The work directly touches on the nature of dark matter and the origin of baryon asymmetry. The asymmetric dark matter freeze-out scenario in this model naturally reproduces the observed density without contradicting a strong phase transition.
🎯 Fun fact: To detect gravitational waves from a phase transition at 10 MeV, pulsar timing arrays measure deviations in pulse arrival times with an accuracy down to tens of nanoseconds — that's like noticing a clock on Jupiter gains a billionth of a second per year. Curiously, similar temperatures are reached in the cores of [tag:supernova]supernovae[/tag].
Key numbers
- Phase transition temperature: 1–100 MeV
- Peak gravitational-wave frequency: ~ 1 nHz
- Typical transition duration (β/H): 1000–5000
- Dark matter mass: 10–100 GeV
- Gauge coupling constant g_d: 0.5–1