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.
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.
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.
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.
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.
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.
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.
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.
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
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.
Nearby works in our archive
The closest work in the archive sits at 0.70 by meaning — the usual level for our corpus.
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.