Dark Phase Transition: Testing the Model with Pulsar Timing Arrays ✛

Original: "A critical look at low-scale cosmological phase transitions in the PTA era"
· Simone Biondini, Philipp Schicho
arXiv:2607.02505v1 · 2026-07-02 · CC BY 4.0 · 3 min
Could a phase transition in a hidden sector at temperatures of tens of MeV generate the gravitational waves detected by pulsar timing arrays?

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.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The paper constructs a three-dimensional effective field theory for a dark Abelian Higgs sector at finite temperature, including higher-order corrections up to NNLO and the contribution of a heavy fermion. The authors computed the gravitational wave spectrum from the phase transition in the nanohertz range and compared it with pulsar timing array data, finding that the predictions are orders of magnitude below the observed signal. They also demonstrated the possibility of hydrodynamic decoupling between the dark and visible sectors.

Where its strength is

The systematic EFT-based approach allowed a quantitative estimate of uncertainties and showed that the simplest model cannot explain the PTA signal without fine-tuning, which is valuable for constraining dark sector parameters. Hydrodynamic decoupling is an important finding that could change the interpretation of the signals and should be considered in future work.

Where it could go next

Worth trying

  • Test how modifying the Higgs potential (e.g., adding a scalar singlet) would strengthen the transition and bring the spectrum closer to NANOGrav data.
  • Estimate the contribution of bubble collisions to the overall signal under hydrodynamic decoupling, using your EFT approach.
  • Compare your dark sector parameter constraints with those from supernova cooling to narrow down the allowed region.

Why this may matter

Refining the contribution of phase transitions to the nanohertz gravitational wave background is important for interpreting PTA data and may point to more complex dark sector physics or alternative signal sources.

Original abstract

as written by the authors arXiv:2607.02505v1

Motivated by the recent evidence for a stochastic gravitational-wave (GW) background reported by pulsar timing array (PTA) collaborations, we perform a precision study of low-scale phase transitions in a dark Abelian Higgs sector, a minimal gauge theory of spontaneous symmetry breaking relevant for cosmological phase transitions. Using dimensionally reduced high-temperature effective field theory, we quantify the impact of thermal resummation, higher-order matching corrections, and higher-dimensional operators on the phase-transition thermodynamics and the resulting GW signal. We find that the parameter region favored by current PTA observations lies close to the boundary of validity of the effective field theory, where higher-dimensional operators become increasingly important. Even within this controlled region, the predicted signal remains disfavored by the PTA data, despite the substantial shifts induced by higher-order thermal corrections. We further delineate parameter regions where the dark and visible sectors are thermally and hydrodynamically coupled or decoupled, and revisit the dark matter phenomenology, identifying asymmetric freeze-out as naturally compatible with both the observed relic abundance and the gauge couplings favored by strong phase transitions. Our results underscore the importance of systematically controlled finite-temperature calculations for reliable GW predictions from low-scale cosmological phase transitions.

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 supernovae.

m_{\text{eff}}^2(T) = -\mu^2 + \frac{T^2}{12}(4\lambda + 3g^2)
Square of the effective field mass at high temperature; the negative -μ² contribution can make the mass negative, triggering spontaneous symmetry breaking.
f_0 \sim \frac{T_* T_0}{M_{\text{Pl}}}
Today's peak frequency f₀ is proportional to the transition temperature T* and the current CMB temperature T₀, divided by the Planck mass Mₚₗ.

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
Scientists
Albert EinsteinBarry BarishKip ThorneRichard FeynmanRussell HulseVera Rubin
Concepts
cosmological phase transition geometric phase control gate voltage Effective field theory gauge field quantum phase transition