Cosmic Strings: How a Precise Decay Calculation Reshapes the Gravitational-Wave Universe ✛ ✎ ✦

Original: "The decay rate of metastable cosmic strings beyond the thin-string approximation"
· Valerie Domcke, Yu Hamada
arXiv:2606.03008v1 · 2026-06-02 · arXiv non-exclusive · independent analysis · 4 min
New numerical simulations reveal that metastable cosmic strings decay faster than predicted by the thin-string approximation, demanding a revision of cosmological forecasts.

Context

One of the most intriguing predictions of cosmology, born from the Big Bang, is the existence of topological defects in spacetime, particularly cosmic strings, which can leave an imprint on the cosmic microwave background. In an expanding Universe, these objects, created during phase transitions, become powerful sources of gravitational waves, producing a stochastic background that detectors like LIGO—built with contributions from Rainer Weiss and Kip Thorne—aim to register. However, in grand unified theories, strings may be metastable: they decay by spontaneously nucleating magnetic monopole-antimonopole pairs. This process—quantum tunneling within quantum field theory—suppresses the gravitational-wave signal at low frequencies and depends critically on the exact value of the bounce action, which previously had only been calculated in the vanishingly thin string approximation.

Methods

The authors abandoned analytical simplifications and turned to numerical methods. They considered a minimal model with two stages of symmetry breaking, SU(2) → U(1) → nothing, where monopoles form first, then strings. To compute the bounce action—the key quantity determining the decay rate—they used the gradient flow method on a spacetime lattice. Instead of searching for the saddle point of the action functional directly, which is complicated by a negative mode, the simulations looked for the minimum action at a fixed monopole center position ($ρ_E$ = R) and then maximized the action over R. The field profiles of gauge bosons and Higgs scalars were described by a general symmetry-restricted ansatz without extra assumptions. Lattice simulations were performed with different lattice spacings to extrapolate the results to the continuum limit.

Results

Numerical experiments revealed a significant reduction in the bounce action compared to the classic formula $S_B$^(PV) = π $m_M^2$/μ. For a clear example with equal masses of Higgs and gauge particles in each of the two sectors, the bounce action drops below 60 (in units of $m_W^2$/($g^2$ m_γ^2)) at a mass ratio $m_W$/m_γ ≈ 3, whereas the thin-string limit gives about 60.5. Extrapolation to the continuum limit allowed the construction of a reliable fitting function that accurately describes $S_B$ across the parameter range. At large hierarchies ($m_W$/m_γ ≫ 1), the results smoothly converge to the thin-string limit, confirming the self-consistency of the method. This means the decay rate Γ ∝ exp(-$S_B$) increases by orders of magnitude, and cosmic strings live shorter than previously thought. Consequently, the characteristic frequency $f_{low}$, below which the gravitational-wave spectrum is cut off, shifts upward—from nanohertz to micro- and millihertz.

Implications

In practice, this amounts to a revision of the expected gravitational-wave background. Parameter regions previously considered safe for ground-based detectors may now fall within the LIGO/Virgo exclusion zone. On the other hand, the signal detected by pulsar timing array (PTA) collaborations like NANOGrav requires a more precise interpretation: to explain it with metastable strings, one must increase the ratio of symmetry breaking scales, which facilitates constructing consistent grand unified models with intermediate inflation. In other words, accounting for the finite thickness of strings makes physics at the grand unification scale more testable.

Future development

The proposed lattice simulation method can be generalized to arbitrary mass ratios of particles arising from symmetry breaking. The next step will be to include effects related to inhomogeneities along the strings as well as thermal production of monopoles. This will allow even more accurate predictions of the gravitational-wave spectrum and possibly explain the recent hint of a signal in PTA without conflicting with LIGO constraints.

Impact

The results will impact gravitational-wave astronomy, observational cosmology, and particle physics, since they connect specific experimental data to the parameters of grand unified theories.

Next steps

The immediate task is to apply this approach to more realistic gauge groups, such as SU(2)×U(1), and to include relativistic corrections to the string motion.

Key open problems

The work touches on fundamental questions at the intersection of quantum field theory and cosmology, including the problem of false vacuum decay in curved spacetime. Moreover, it provides a new tool for interpreting the enigmatic signal from pulsar timing arrays, based on observations of neutron stars—whether it is the first echo of string physics or has an astrophysical origin remains to be determined in the coming years.

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 examines the decay of cosmic strings via quantum nucleation of monopole-antimonopole pairs in an SU(2) theory with two stages of symmetry breaking. The authors performed lattice simulations and computed the exact bounce action, abandoning the infinitely thin string approximation. They show that the finite string thickness significantly accelerates the decay and shifts the cutoff frequency in the gravitational wave spectrum.

Where its strength is

A numerical calculation of the bounce action without analytical simplifications yielded a reliable dependence on the model parameters, eliminating a key uncertainty in predictions for the gravitational wave signal from metastable strings.

Where it could go next

Worth trying

  • Apply the obtained fitting function for the bounce action to compute the stochastic background and compare with pulsar timing array data.
  • Extend the lattice simulations to strings with currents to test how additional degrees of freedom affect the decay.
  • Investigate the effect of finite string thickness on the gravitational wave spectrum in scenarios with an intermediate expansion phase.

Why this may matter

Refining the lifetime of cosmic strings directly impacts predictions for gravitational wave experiments and could help connect observations with grand unified physics.

A cosmic string as thick as a billionth of a proton’s size can possess a tension equivalent to a million tons per centimeter of length — comparable to the mass of a skyscraper packed into a ruler.

S_B^{(\mathrm{PV})} = \frac{\pi m_M^2}{\mu}
S_B is the bounce action, m_M is the monopole mass, μ is the string tension

Key numbers

  • typical string tension (Gμ): 10-8
  • bounce action in the thin-string approximation: ≈ 60.5
  • interesting range of mass ratio m_W/m_γ: 3–9
  • cutoff frequency f_low: ~ nanohertz
  • action suppression due to finite thickness: up to 20%
Scientists
Albert EinsteinBarry BarishKip ThorneRichard FeynmanRussell HulseHenri Poincaré
Concepts
cosmic string cosmic string gravitational waves magnetic charge decay rate false vacuum decay damping rate magnetic monopole false vacuum stellar collapse mechanism lattice spacing bubble nucleation age of the universe inflationary cosmology planck time stellar stream relaxation time moment of inertia eternal inflation cosmological constant inflaton