Loops produced by Abelian-Higgs cosmic string networks in lattice field-theory simulations decay in a time proportional to their initial length, tlife=(0.14±0.04) ℓinitt_{\rm life}=(0.14\pm0.04)\,\ell_{\rm init}, far faster than gravitational radiation allows, while carefully tuned oscillating loops have lifetimes scaling as ℓinit2\ell_{\rm init}^2 and Nambu-Goto-like velocities. The rapid decay is attributed to excitation of internal massive (dilatational) degrees of freedom on the string, with energy per unit length apparently saturated at e≃0.25 μ0e\simeq0.25\,\mu_0, but the instability that excites these modes in evolving networks has not been identified, and the fraction fNGf_{\rm NG} of loops surviving as Nambu-Goto-like gravitational radiators is bounded only crudely (fNG<0.1f_{\rm NG}<0.1 at 95% confidence, from 31 loops). Since essentially all cosmic-string gravitational-wave forecasts rest on Nambu-Goto loop dynamics, this gap directly limits predictions for pulsar timing arrays, LISA, and ground-based detectors, as well as particle and dark-matter forecasts. The problem is to establish the physical origin and amplitude control of internal-mode excitation in evolving field-theory string networks, to determine the Nambu-Goto-like loop fraction with adequate statistics, and to convert this into a calibrated decay model and field-theory-based multi-messenger forecasts. Later work has measured loop-level emission, demonstrated parametric instabilities of artificially excited strings, and built fNGf_{\rm NG}-parametrized constraint frameworks, but the network mechanism and fNGf_{\rm NG} itself remain undetermined.

Background

Cosmic strings are one-dimensional topological defects that can form at symmetry-breaking phase transitions in the early universe; if they exist, networks of them would source a stochastic gravitational-wave background (SGWB) within reach of pulsar timing arrays (PTAs), LISA, and ground-based detectors. The simplest field-theoretic realization is the Abelian-Higgs model, a complex scalar field coupled to a U(1)U(1) gauge field, whose strings are Nielsen-Olesen vortices: flux tubes with a finite core width set by the scalar and gauge masses msm_s and mvm_v. The standard analytical treatment replaces them by Nambu-Goto (NG) strings — infinitely thin objects whose motion extremizes the relativistic string worldsheet area. In the NG approximation, closed loops oscillate essentially without loss, radiating only gravitational waves, so their lifetime is set by GμG\mu (the string tension in gravitational units) and is astronomically long; all widely used SGWB forecasts for cosmic strings (built on the VOS, "velocity-dependent one-scale", and BOS, Blanco-Pillado-Olum-Shlaer, loop-distribution models used by PTA and LISA collaborations) are built on this assumption.

Lattice simulations of the full Abelian-Higgs field theory tell a different story. Network simulations reaching scaling were pioneered by Hindmarsh, Stuckey and Bevis (Phys. Rev. D 79, 123504 (2009)), who found that loops chopped off evolving networks disappear quickly rather than persisting as NG oscillators. Matsunami, Pogosian, Saurabh and Vachaspati (Phys. Rev. Lett. 122, 201301 (2019)) showed that artificially constructed loops decay by particle radiation with lifetime growing roughly as the square of their length. The central study of this problem, Hindmarsh, Lizarraga, Urio and Urrestilla, Phys. Rev. D 104, 043519 (2021), quantified the dichotomy: loops formed by intersections in a network have tlife=(0.14±0.04) ℓinitt_{\rm life}=(0.14\pm0.04)\,\ell_{\rm init} (linear regime, no trend with length, up to ℓinit≃6000\ell_{\rm init}\simeq6000 inverse mass units), while tuned, initially stationary, non-self-intersecting loops constructed by colliding pairs of wrapping strings have tlife∝ℓinit2t_{\rm life}\propto\ell_{\rm init}^2 (quadratic regime) and mean-square velocity vˉ2=0.500±0.004\bar v^2=0.500\pm0.004, the NG value in flat spacetime; network loops instead have vˉ2≃0.40±0.04\bar v^2\simeq0.40\pm0.04. The proposed explanation is the excitation of extra degrees of freedom carried on the string — plausibly the dilatational ("sausage") bound mode of the Nielsen-Olesen vortex, whose existence and coupling to massive radiation were established by Goodband and Hindmarsh (Phys. Rev. D 52, 4621 (1995)). A phenomenological model assigns these modes an energy per unit length ee and pressure pp, so that the string's total tension is μ=μ0+e\mu=\mu_0+e and its transverse tension is T=μ0−pT=\mu_0-p, with μ0\mu_0 the tension of the unexcited string; the modes radiate with power Prad=−κe2/φ0P_{\rm rad}=-\kappa e^2/\varphi_0 per unit length, where κ\kappa is a dimensionless coupling and φ0\varphi_0 is the symmetry-breaking field expectation value setting the string mass scale. Weak excitation (e≪μ0e\ll\mu_0, with e0=ϵφ0/ℓinite_0=\epsilon\varphi_0/\ell_{\rm init} seeded by curvature) gives the quadratic lifetime regime, while saturation of a long-wavelength component at e∼μ0e\sim\mu_0 (a constant κ~ϵaϵb\tilde\kappa\epsilon_a\epsilon_b decay rate) gives the linear regime; the observed velocity deficit implies e/μ0≃0.25e/\mu_0\simeq0.25. That network loops are generically saturated points to an intrinsic instability in the background of an evolving string, whose nature the authors call "somewhat mysterious": it is shut off if loops are prevented from collapsing (by spinning or wrapping), and one candidate is parametric resonance of a massive mode propagating in the effective 1+11+1-dimensional spacetime of an oscillating string.

The observational stake is encapsulated in the parameter fNGf_{\rm NG}, the fraction of loops in a network that follow NG dynamics and survive long enough to radiate gravitationally — operationally, loops surviving beyond 0.25 ℓinit0.25\,\ell_{\rm init}. From the fact that none of 31 large network loops survived beyond 0.25 ℓinit0.25\,\ell_{\rm init}, the 2021 paper estimates fNG<0.1f_{\rm NG}<0.1 at 95% confidence (a "rule of three" bound that includes zero: with zero NG-like survivors among n=31n=31 sampled loops, the 95% upper limit on the survival fraction is 3/n≃0.13/n\simeq0.1); conversely, if the NANOGrav PTA signal were due to hidden-sector strings, consistency with cosmic-microwave-background bounds would require fNG≳10−3f_{\rm NG}\gtrsim10^{-3}. Subsequent work has sharpened parts of this picture: Blanco-Pillado et al. (JCAP 05 (2023) 035) compared field-theory loops trajectory-by-trajectory with NG reconstructions; Baeza-Ballesteros, Copeland, Figueroa and Lizarraga (Phys. Rev. D 112, 043540 (2025)) measured the simultaneous particle and gravitational-wave emission of loops on the lattice; and the effective-action program of Aurrekoetxea et al. (arXiv:2603.05243) derived how massive worldsheet bound states couple to the string's Goldstone modes and to worldsheet curvature. What none of these delivers is the network-level mechanism, statistics, or forecast — the open problem stated here.

Problem Statement

Establish a validated, field-theory-based account of how loops decay in Abelian-Higgs cosmic string networks, sufficient to replace Nambu-Goto-based observational modeling. Concretely, in the formulation of Hindmarsh, Lizarraga, Urio and Urrestilla (2021), the following linked questions must be answered:

  1. Identify the physical origin and detailed nature of the instability that excites the internal massive degrees of freedom (plausibly the dilatational bound mode) on network-produced strings: which mechanism — parametric resonance of a massive mode in the effective 1+11+1-dimensional worldsheet spacetime, nonlinear mode coupling, reconnection- or curvature-induced excitation, or another process — operates in evolving networks, and under what conditions it activates or is suppressed (the 2021 simulations show it shutting off when loops are prevented from collapsing, e.g. by spinning or wrapping, and being avoidable only by careful tuning of initial conditions).

  2. Determine the excitation pathway and amplitude distribution: how initial conditions and network interactions populate long- and short-wavelength excitations, and what controls whether the excitation energy per unit length ee stays small (e≪μ0e\ll\mu_0, giving the quadratic lifetime regime tlife∝ℓinit2t_{\rm life}\propto\ell_{\rm init}^2) or saturates (e∼μ0e\sim\mu_0, giving the linear regime tlife=(0.14±0.04)ℓinitt_{\rm life}=(0.14\pm0.04)\ell_{\rm init} with the observed velocity deficit vˉ2≃0.40\bar v^2\simeq0.40).

  3. Determine the fraction fNGf_{\rm NG} of loops in realistic cosmological string networks that evolve as Nambu-Goto-like oscillators, replacing the current bound (fNG<0.1f_{\rm NG}<0.1 at 95% confidence, from 31 network loops) with a measurement from much larger-statistics simulations exploring cosmologically relevant parameter ranges and including the effects of cosmic expansion.

  4. On this basis, recalculate the gravitational-wave and particle-emission signatures of field-theory string networks — the stochastic gravitational-wave background and the high-energy-particle (e.g. gamma-ray and dark-matter) yields — so as to replace Nambu-Goto-based forecasts with robust field-theory-based ones, with fNGf_{\rm NG} or its mechanistic successor as a determined input rather than a free parameter.

  5. Validate and refine the phenomenological decay model (parameters κ,ϵ,κ~\kappa,\epsilon,\tilde\kappa and its equation-of-state assumptions), testing whether it predicts the distribution of loop lifetimes, the energy spectra of the emitted massive radiation, and the conditions for the transition between the quadratic and linear lifetime regimes.

These five components jointly serve one objective — a predictive account of network-loop decay linking microscopic field dynamics to observables.

References

  1. M. Hindmarsh, J. Lizarraga, A. Urio, J. Urrestilla, "Loop decay in Abelian-Higgs string networks", Phys. Rev. D 104, 043519 (2021). DOI: 10.1103/PhysRevD.104.043519; arXiv:2103.16248; https://arxiv.org/abs/2103.16248
  2. J. J. Blanco-Pillado, D. Jiménez-Aguilar, J. Lizarraga, A. Lopez-Eiguren, K. D. Olum, A. Urio, J. Urrestilla, "Nambu-Goto dynamics of field theory cosmic string loops", JCAP 2023 (05) 035. DOI: 10.1088/1475-7516/2023/05/035; arXiv:2302.03717; https://arxiv.org/abs/2302.03717
  3. J. Baeza-Ballesteros, E. J. Copeland, D. G. Figueroa, J. Lizarraga, "Particle and gravitational wave emission by local string loops: Lattice calculation", Phys. Rev. D 112, 043540 (2025). DOI: 10.1103/ym9p-scmw; arXiv:2408.02364; https://arxiv.org/abs/2408.02364
  4. J. Baeza-Ballesteros, E. J. Copeland, D. G. Figueroa, J. Lizarraga, "Gravitational wave emission from a cosmic string loop: Global case", Phys. Rev. D 110, 043522 (2024). DOI: 10.1103/PhysRevD.110.043522; arXiv:2308.08456; https://arxiv.org/abs/2308.08456
  5. M. Hindmarsh, J. Kume, "Multi-messenger constraints on Abelian-Higgs cosmic string networks", JCAP 2023 (04) 045. DOI: 10.1088/1475-7516/2023/04/045; arXiv:2210.06178; https://arxiv.org/abs/2210.06178
  6. J. Kume, M. Hindmarsh, "Revised bounds on local cosmic strings from NANOGrav observations", JCAP 2024 (12) 001. DOI: 10.1088/1475-7516/2024/12/001; arXiv:2404.02705; https://arxiv.org/abs/2404.02705
  7. J. J. Blanco-Pillado, D. Jiménez-Aguilar, J. M. Queiruga, J. Urrestilla, "Parametric Resonances in Axionic Cosmic Strings", JCAP 2023 (04) 043. DOI: 10.1088/1475-7516/2023/04/043; arXiv:2212.06194; https://arxiv.org/abs/2212.06194
  8. J. C. Aurrekoetxea, J. J. Blanco-Pillado, A. García Martín-Caro, J. M. Queiruga, "On curvature corrections for field theory cosmic strings" (2026). arXiv:2603.05243; https://arxiv.org/abs/2603.05243
  9. Z. Zhao, L. Bian, J. Shu, "Testing the Nambu-Goto approximation of cosmic string by lattice field theory simulations", Phys. Rev. D 114, 023539 (2026). DOI: 10.1103/52nd-mg9m; arXiv:2507.00685; https://arxiv.org/abs/2507.00685
  10. A. Alonso-Izquierdo, J. Blanco-Pillado, D. Miguélez-Caballero, S. Navarro-Obregón, J. Queiruga, "Excited Abelian-Higgs vortices: Decay rate and radiation emission", Phys. Rev. D 110, 065009 (2024). DOI: 10.1103/PhysRevD.110.065009; arXiv:2405.06030; https://arxiv.org/abs/2405.06030
  11. M. Hindmarsh, S. Stuckey, N. Bevis, "Abelian Higgs cosmic strings: Small-scale structure and loops", Phys. Rev. D 79, 123504 (2009). DOI: 10.1103/PhysRevD.79.123504; arXiv:0812.1929; https://arxiv.org/abs/0812.1929
  12. D. Matsunami, L. Pogosian, A. Saurabh, T. Vachaspati, "Decay of Cosmic String Loops due to Particle Radiation", Phys. Rev. Lett. 122, 201301 (2019). DOI: 10.1103/PhysRevLett.122.201301; arXiv:1903.05102; https://arxiv.org/abs/1903.05102
  13. M. Goodband, M. Hindmarsh, "Bound states and instabilities of vortices", Phys. Rev. D 52, 4621 (1995). DOI: 10.1103/PhysRevD.52.4621; arXiv:hep-ph/9503457; https://arxiv.org/abs/hep-ph/9503457
  14. A. Dimitriou, D. G. Figueroa, P. Simakachorn, B. Zaldívar, "Cosmic string gravitational wave backgrounds at LISA: I. Signal survey, template reconstruction, and model comparison", JCAP 2026 (05) 037. DOI: 10.1088/1475-7516/2026/05/037; arXiv:2508.05395; https://arxiv.org/abs/2508.05395