New gravitational-wave templates for metastable cosmic strings: Loop breaking versus network collapse

This paper introduces a new three-parameter gravitational-wave template for metastable cosmic strings that distinguishes between loop breaking and network collapse time scales, unifying standard and quasi-stable string models and providing an analytical expression to improve the analysis of future pulsar timing array data.

The Gist

Imagine the universe is a giant, invisible fabric. Sometimes, when this fabric gets stretched and cooled down in the very early universe, it doesn't smooth out perfectly. Instead, it gets knotted or torn, creating thin, string-like defects called cosmic strings. Think of these like incredibly strong, invisible rubber bands stretched across the cosmos.

For a long time, scientists thought these strings were permanent. They believed that once formed, they would just wiggle around forever, snapping off little loops that would vibrate and create a faint, constant hum of gravitational waves (ripples in space-time). This "hum" is what scientists call a Gravitational Wave Background (GWB).

Recently, astronomers using pulsar timing arrays (which act like giant, ultra-precise clocks made of spinning stars) heard a low-frequency rumble that looks a lot like this cosmic string hum. However, there was a problem: the "shape" of the sound predicted by the old, permanent-string theory didn't quite match the sound they actually heard.

The New Idea: Strings That "Unravel"
This paper proposes a new way to look at these strings. Instead of being permanent, these strings are metastable. This means they are like a house of cards or a sandcastle: they look stable for a while, but they are destined to collapse eventually.

The authors argue that we need to stop treating the collapse of these strings as a single event. Instead, they identify two different clocks ticking away at different speeds:

1. The "Loop Breaker" Clock ($t_{LB}$): This is the time it takes for the small, closed loops of string (the ones that are already detached and floating freely) to spontaneously snap and disappear. Imagine a rubber band that slowly develops a weak spot and eventually snaps.
2. The "Network Collapse" Clock ($t_{NC}$): This is the time it takes for the long strings (the ones still connected to the main network) to start falling apart. This happens because the long strings are attached to heavy "monopoles" (like heavy weights tied to the ends of the rubber bands). Eventually, these weights pull the strings into the observable universe, causing the whole network to crumble.

The Big Discovery: The Clocks Don't Have to Match
In previous models, scientists assumed these two clocks were synchronized. They thought the long strings would fall apart at the exact same moment the small loops started snapping.

This paper says: "Not necessarily!"

The authors show that these two clocks can run at very different speeds.

• Scenario A (The Old View): The clocks are synchronized. The network collapses just as the loops start breaking.
• Scenario B (The New View): The network collapses much earlier than the loops start breaking. Imagine the long rubber bands are cut and pulled apart by the weights (Network Collapse) while the small, detached loops are still sitting there, perfectly intact, for a long time afterward (Loop Breaking).

Why Does This Matter?
By separating these two clocks, the authors created a three-parameter model (String Tension, Loop Breaking Time, and Network Collapse Time). This allows for a much wider variety of "sounds" (spectral shapes) that the strings could produce.

They found that when the Network Collapse happens much faster than the Loop Breaking, the resulting gravitational wave signal changes shape. Specifically, the low-frequency part of the signal becomes steeper. This new shape fits the data from the 2023 pulsar timing array observations much better than the old models did.

The "Quasi-Stable" Connection
The paper also unifies two different ideas that scientists had been arguing about: "Metastable Strings" and "Quasi-Stable Strings."

• Think of "Metastable" as a string that is about to break.
• Think of "Quasi-Stable" as a string that is very close to breaking but hasn't yet.

The authors show that these aren't two different types of strings; they are just the same type of string viewed under different timing conditions. By adjusting the "Network Collapse" clock, you can smoothly slide from one description to the other.

The Bottom Line
The authors have created new "templates" (blueprints) for what the gravitational wave signal from these strings should look like. These new blueprints are more flexible and accurate. They suggest that the universe might contain strings that unravel in a specific way that perfectly explains the mysterious hum astronomers heard in 2023.

They also derived a simple mathematical formula that describes this signal when the network collapses quickly. This formula is now ready for other scientists to use when they analyze future data from gravitational wave detectors to see if it matches their new predictions.

Technical Summary

1. Problem Statement

Metastable cosmic strings are a compelling candidate to explain the stochastic gravitational-wave background (GWB) signal observed in 2023 by Pulsar Timing Array (PTA) collaborations (e.g., NANOGrav, EPTA). Unlike topologically stable strings, metastable strings can decay via the spontaneous nucleation of monopoles (Schwinger-like vacuum tunneling).

The Core Issue: Previous modeling of the GWB from metastable strings relied on a two-parameter model ($G\mu$ and a hierarchy parameter $p_\kappa$) based on the strict assumption that the time scale for loop breaking ($t_{LB}$) is identical to the time scale for network collapse ($t_{NC}$).

• $t_{LB}$: The time scale on which closed sub-Hubble loops break apart due to monopole nucleation.
• $t_{NC}$: The time scale on which the network of long strings collapses because finite segments attached to monopoles enter the Hubble horizon.

Recent theoretical work (Ref. [73]) suggested that $t_{NC}$ could be significantly smaller than $t_{LB}$ due to finite-temperature effects, thermal monopole production, or inflationary dynamics. The authors argue that treating $t_{NC}$ and $t_{LB}$ as distinct, independent parameters is necessary to accurately model the GWB spectrum and interpret PTA data.

2. Methodology

The authors develop a generalized framework to compute the GWB spectrum from metastable strings, moving from a two-parameter to a three-parameter model:

1. String Tension ($G\mu$): Controls the overall amplitude of the signal.
2. Loop Breaking Time Scale ($t_{LB}$): Parametrized by the hierarchy parameter $p_\kappa$, governing the exponential decay of loop number density.
3. Network Collapse Time Scale ($t_{NC}$): Treated as an independent free parameter, governing when the production of new loops ceases.

Key Technical Steps:

• Formalism: They utilize the Velocity-Dependent One-Scale (VOS) model to describe the evolution of the long-string network.
• Loop Number Density: They derive a modified expression for the loop number density $n_{meta}(\ell, t)$:
  $$n_{meta}(\ell, t) = n_{stable}(\ell, t) \cdot \Theta(t_{NC} - t_*) \cdot \exp\left[-\frac{A(\ell, t)}{t_{LB}^2}\right]$$
  where $t_*$ is the formation time of a loop, $\Theta$ is the Heaviside step function (halting loop production at $t_{NC}$), and the exponential term accounts for loop decay.
• Spectral Integration: They compute the GWB energy density spectrum $\Omega_{GW}(f)$ by integrating over the loop population, accounting for redshift and the specific break-up mechanisms.
• Regime Analysis: They analyze the parameter space defined by the relationship between $t_{LB}$, $t_{NC}$, and the age of the Universe ($t_0$).
• Analytical Derivation: For the regime where $t_{LB} \gg t_0$ (quasi-stable loops) but $t_{NC} < t_0$ (collapsing network), they derive a compact analytical expression for the GWB spectrum.

3. Key Contributions

• Generalization to Three Parameters: The paper establishes that $t_{NC}$ and $t_{LB}$ are logically distinct. This breaks the previous assumption that string breaking dynamics are uniform across all scales (super-Hubble vs. sub-Hubble).
• Unified Description: The authors demonstrate that their three-parameter model unifies the description of "standard metastable strings" (where $t_{NC} \approx t_{LB}$) and "quasi-stable strings" (where $t_{LB} \gg t_0$ but $t_{NC} < t_0$).
• New Spectral Shapes: They identify previously unexplored spectral behaviors, particularly in the low-frequency tail of the GWB spectrum.
• Analytical Template: They provide a fully analytical expression (Eq. 51) for the GWB spectrum in the limit of large $p_\kappa$ (quasi-stable loops), which matches numerical results with high precision.

4. Results

The authors identify three distinct regimes in the $p_\kappa$–$t_{NC}$ parameter space:

1. Breaking Loops, Collapsing Network ($t_{NC} < t_{LB} < t_0$):

   • Both loops and the network decay within the history of the Universe.
   • Result: The low-frequency spectral index changes from $f^2$ (standard metastable case) to $f^3$ when $t_{NC} \ll t_{LB}$. This steepening can significantly improve the fit to PTA data.

2. Quasi-Stable Loops, Quasi-Stable Network ($t_0 < t_{NC} < t_{LB}$):

   • Both time scales exceed the age of the Universe.
   • Result: The network behaves indistinguishably from stable strings. The spectrum is flat ($f^0$) in the relevant band, which is inconsistent with the rising PTA signal.

3. Quasi-Stable Loops, Collapsing Network ($t_{NC} < t_0 < t_{LB}$):

   • Loops are effectively stable (do not break), but the network collapses because long segments enter the horizon.
   • Result: This regime corresponds to "quasi-stable strings." The authors derive a compact analytical formula (Eq. 51) for this case.
   • Low-Frequency Behavior: The spectrum exhibits a $f^3$ power law at low frequencies, transitioning to a plateau.
   • Fit to Data: The analytical template shows excellent agreement with the NANOGrav 15-year (NG15) data, particularly for $G\mu \sim 10^{-7}$ and specific $t_{NC}$ values.

Numerical Findings:

• The paper presents 18 benchmark spectra covering these regimes.
• They show that fixing $G\mu$ and $t_{NC}$ does not uniquely determine the spectrum if $t_{LB}$ is varied; the spectral shape evolves as $t_{LB}$ increases, eventually converging to the quasi-stable analytical limit.
• They confirm that a string tension of $G\mu = 10^{-7}$ is marginally ruled out by LIGO-Virgo-KAGRA (LVK) O4a data, while $G\mu \sim 10^{-7.5}$ remains viable.

5. Significance

• Phenomenological Impact: The new templates offer a broader variety of spectral shapes, allowing for a better fit to the 2023 PTA signal than previous models. Specifically, the ability to tune $t_{NC}$ independently of $t_{LB}$ allows for steeper low-frequency slopes ($f^3$) that match the data more naturally.
• Theoretical Implications: The results relax constraints on Grand Unified Theory (GUT) model building. Previously, models required a small hierarchy parameter ($p_\kappa \lesssim 10$) to fit the data. The new analysis shows that models with large hierarchies ($p_\kappa \gtrsim 10$, corresponding to energy scale separations of factors $\sim 10\text{--}30$) are also viable, provided $t_{NC}$ is appropriately suppressed.
• Future Data Analysis: The derived analytical expressions provide ready-to-use templates for future PTA data analysis, enabling more robust constraints on the underlying particle physics parameters (monopole mass, string tension) without needing complex numerical simulations for every data point.

In conclusion, this paper resolves a critical ambiguity in metastable string modeling by decoupling network collapse from loop breaking, providing a more flexible and physically accurate framework for interpreting the emerging gravitational-wave background signal.