Cosmic Strings: How a Precise Decay Calculation Reshapes the Gravitational-Wave Universe ✦
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.
🎯 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.
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%