The decay rate of metastable cosmic strings beyond the thin-string approximation

arXiv:2606.03008 · hep-ph, astro-ph.CO, hep-th · Submitted 2026-06-02 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "The decay rate of metastable cosmic strings beyond the thin-string approximation".

Jocelyn: In grand unified theories, cosmic strings are often metastable and their decay rate, set by monopole pair creation, is crucial for understanding their phenomenology and potential gravitational wave background.

Vera: First, who's behind it and why it matters.

Paper summary: Vera: So we're looking at this paper, "The decay rate of metastable cosmic strings beyond the thin-string approximation," which tackles how these cosmic strings actually die in grand unified theories. The core idea here seems to be that previous calculations relied too heavily on the thin string approximation, and they're using classical lattice simulations to get a more accurate picture of the decay action.

Jocelyn: Right, so if we can get a better handle on the decay rate, it directly impacts what we might see in gravitational wave backgrounds as these strings evolve. It sounds like the authors are focused on moving past those old estimates for how fast these strings break down.

Subrahmanyan: I agree with Jocelyn; understanding that decay rate is fundamental because it sets the timescale for when these structures disappear, which has massive implications for how we model cosmological evolution and detect them observationally <ref:2606.03008#pg1>.

Vera: Exactly, so the paper's thesis is that classical lattice simulations reveal a suppression of the bounce action compared to what was estimated before, which points toward a faster string decay than previously thought <ref:2606.03008#pg2>. It’s all about refining that calculation for the spontaneous creation of monopole pairs along the string, governed by Eq. (one point one) <ref:2606.03008#pg0>.

Jocelyn: And that suppression in action directly leads to an increase in the decay rate, which is a significant finding because it shifts things around in our understanding of these phenomena <ref:2606.03008#pg2>. How does this specific numerical approach help them achieve that result?

Subrahmanyan: The methodology involves moving beyond the thin string approximation by numerically evaluating the bounce action associated with spontaneous monopole formation on the string using gradient flow to find a more accurate expression for Eq. (one point one) <ref:2606.03008#pg3>. They are promoting field functions to depend on fictitious flow time tau and solving that gradient flow equation to find the minimum action S*(R) with a fixed monopole core radius R <ref:2606.03008#pg3>.

Vera: That sounds like a pretty involved numerical setup, especially dealing with those boundary conditions where they need Dirichlet conditions at the string core at rho=zero for certain fields to ensure regularity <ref:2606.03008#pg4>. It seems they had to be very careful about maintaining physical consistency in their simulation <ref:2606.03008#pg4>.

Jocelyn: Maintaining those constraints at the core must be tricky, especially when you're trying to capture a physics that goes beyond the simple thin string model <ref:2606.03008#pg1>. So, what are the main results they found after running these simulations?

Paper summary: Subrahmanyan: The simulation results show this significant suppression of the bounce action as we move from hierarchical to near degenerate symmetry breaking scales <ref:2606.03008#pg4>. They also noted that their continuum limit fits Eq. (four point one) exactly, which provides the exact solution for the bounce action for a benchmark case <ref:2606.03008#pg4>.

Vera: That is quite a strong result because it shows their method recovers the analytical estimates in the thin string limit while providing a corrected expression that resolves issues where previous solutions exceeded that limit for large m squared M/mu ratios <ref:2606.03008#pg4>. It confirms they found what they claim is "the minimal energy one" <ref:2606.03008#pg4>.

Jocelyn: If it is indeed the minimal energy solution, how does this translate into something tangible for our search for gravitational wave backgrounds? The paper discusses implications for gravitational wave searches in Section four point two <ref:2606.03008#pg1>.

Subrahmanyan: The calculated decay rate, which they relate to S B = pi kappa eff in Eq. (four point four), directly impacts the stochastic gravitational wave background spectrum <ref:2606.03008#pg1>. Specifically, the finite string width reduces the bounce action compared to the thin string limit, leading to a corresponding enhancement of the string decay rate <ref:2606.03008#pg1>.

Vera: So this shift in decay rate moves where we expect to see features in the gravitational wave spectrum—it shifts that turnover point to larger frequencies, which is very important for interpreting signals from pulsar timing arrays and LISA <ref:2606.03008#pg1>.

Jocelyn: That connection between the string physics and observable GW constraints is really compelling because it maps preferred regions in parameter space directly onto what we can actually measure with current or future instruments <ref:2606.03008#pg1>. It suggests the full solution offers a more accurate picture than relying solely on the thin string approximation for those searches.

Subrahmanyan: To put that bigger picture into context, this work demonstrates that going beyond the thin string approximation is necessary when we are in near-degenerate symmetry breaking regimes to accurately determine the decay rate of metastable cosmic strings <ref:2606.03008#pg1>. This finding simplifies concrete implementations in cosmological models because a given string lifetime corresponds to larger mass hierarchies between those symmetry breaking scales <ref:2606.03008#pg5>.

Vera: It really highlights how important it is to use these more detailed calculations when we are trying to constrain the various GUT models that predict these structures <ref:2606.03008#pg1>. This paper shows us where the observational constraints on gravitational waves can guide our theoretical modeling of the early universe.

Paper summary: Jocelyn: It sounds like this paper is a really valuable piece for anyone working on pulsar timing array data, showing exactly how the underlying string physics dictates what kind of background we should be looking for <ref:2606.03008#pg1>. So, if we look at the authors' conclusions regarding their title and the overall implications, what do you both think is the main point they are driving home?

Subrahmanyan: They are showing that classical lattice simulations successfully recover the thin string limit while simultaneously providing a corrected expression for the bounce action <ref:2606.03008#pg1>. This means we have a robust way to calculate string decay rates, which has significant implications for constraining GUT models and interpreting astrophysical gravitational wave observations <ref:2606.03008#pg5>.

Vera: I think the main point is that this improved calculation gives us a more reliable tool than what was previously available when dealing with near-degenerate breaking scales <ref:2606.03008#pg4>. It moves our understanding of these strings forward by providing a solution that is confirmed to be the minimal energy one <ref:2606.03008#pg4>.

Jocelyn: And for us, it means when we look at potential SGWB signals, we can use this refined framework to better map out what specific parameter space regions might correspond to those signals <ref:2606.03008#pg1>. It’s a practical tool for connecting theory and observation.

Subrahmanyan: Indeed, the paper proves that when you go beyond the thin string approximation, you get a more accurate picture of the physics governing how these strings decay <ref:2606.03008#pg1>. This refined understanding is critical for testing our cosmological models derived from grand unified theories <ref:2606.03008#pg5>.

Vera: So, to wrap up, this work confirms that when you look at the decay rate of metastable cosmic strings beyond the thin-string approximation, we get a suppression of the bounce action that leads to a faster decay and better constraints on gravitational wave phenomenology <ref:2606.03008#pg4>.

Jocelyn: And for us, it means we can interpret those potential signals from pulsar timing arrays with much more confidence because the underlying physics model is more precise now <ref:2606.03008#pg1>. We're really excited to see how this refined picture helps us guide our next observational surveys.

Subrahmanyan: That’s the essence of it; a better calculation means we can build more accurate constraints on the symmetry breaking scales in GUTs <ref:2606.03008#pg5>. The work provides a solid foundation for connecting high-energy theory to observable cosmological signals <ref:2606.03008#pg1>.

Vera: It’s been a really insightful discussion on how moving beyond the thin string approximation is essential for accurately determining these decay rates and understanding the implications for gravitational wave searches <ref:2606.03008#pg1>.

Conclusion: Vera: So, we’ve been diving deep into how these classical lattice simulations refine our understanding of cosmic string decay rates, and now we’re coming to the conclusion of this paper, "The decay rate of metastable cosmic strings beyond the thin-string approximation."

Jocelyn: Exactly; after seeing all that technical detail about the SU(two) theory and those complex numerical methods, we need to distill what this means for us on the pulsar timing array side.

Subrahmanyan: From a theoretical standpoint, these authors successfully moved past the limitations of the thin-string approximation by finding a more accurate expression for that bounce action.

Vera: That’s right; they're essentially showing us that when symmetry breaking scales are close together, we have to use this more detailed approach to get the actual decay rate correct.

Jocelyn: And what they’ve found is that this refined calculation directly affects the predicted stochastic gravitational wave background spectrum.

Subrahmanyan: They've shown that this suppression of action translates into a faster string decay, which in turn means the turnover point in the gravitational wave spectrum shifts to higher frequencies.

Vera: It’s fascinating because it connects these high-energy particle physics models right down to observable signals we’re trying to hunt with instruments like LISA.

Jocelyn: That shift in frequency is a crucial piece of information for us, as it tells us exactly where we should focus our observational efforts when looking at PTA data.

Subrahmanyan: This work suggests that the lifetime of these strings is linked to larger mass hierarchies between the symmetry breaking scales, which simplifies how we can implement these ideas in cosmological models.

Vera: It really helps ground the abstract theory in concrete constraints, showing us where those theoretical predictions actually intersect with what’s possible to measure astrophysically.

Jocelyn: So, if we take this result as a starting point for our future searches, what specific regions of parameter space are we now better equipped to explore with more confidence?

Theoretical Physics Department, CERN · Department of Physics, The University of Osaka

hep-ph, astro-ph.CO, hep-th

Submitted: 2026-06-02

Updated: 2026-10-07

Comments: v2: published version, 29 pages, 10 figs

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 90/100

The gist: In grand unified theories, cosmic strings are often metastable and their decay rate, set by monopole pair creation, is crucial for understanding their phenomenology and potential gravitational wave

Key concepts

Cosmic String Decay
Metastable cosmic strings can lose energy by spontaneously creating monopole-antimonopole pairs along their length. The rate of this process determines how long the string survives before decaying, which affects observable phenomena like gravitational waves.
Bounce Action
This is a numerical quantity calculated using gradient flow to find the minimum action required for spontaneous monopole formation on the string. It serves as a key parameter in Equation (1.1) that dictates the string's decay rate.
Thin String Approximation
This is a simplified model used previously, treating cosmic strings as infinitely thin lines. The simulation showed that this approximation fails in near-degenerate regimes, meaning the full, wider string solution provides a more accurate description of the physics.

Terminology

Summary

In grand unified theories, cosmic strings are often metastable and their decay rate, set by monopole pair creation, is crucial for understanding their phenomenology and potential gravitational wave background. The core finding of this work is that classical lattice simulations reveal a suppression of the bounce action compared to previous estimates, indicating a faster string decay than previously thought.

Minimal Setup and Theory

The framework begins with an SU(2) gauge theory featuring two successive stages of symmetry breaking controlled by an adjoint scalar field and a doublet scalar field. The first stage, occurring at scale V, produces 't Hooft-Polyakov monopoles due to the non-trivial first homotopy class of the vacuum manifold. The second stage, at scale v, breaks the remaining U(1) symmetry and gives rise to cosmic strings. The decay mechanism is modeled as the spontaneous creation of monopole antimonopole pairs along the string with a rate given by Eq. (1.1):

[6–8]

Γ ∼ µ2π e − πm2M/µ,

per unit length of string, where mM is the monopole mass and µ the string tension.

Numerical Calculation of Bounce Solution

The paper moves beyond the thin string approximation by numerically evaluating the bounce action associated with spontaneous monopole formation on the string using gradient flow. This method aims to find a more accurate expression for Eq. (1.1) that overcomes shortcomings in previous treatments, particularly in the near degenerate case where m2M ≳ µ. The simulation involves promoting field functions to depend on fictitious flow time τ and solving the gradient flow equation:

∂τX(τ, ρ, ρE) = −δS/δX(τ, ρ, ρE)

where X represents the set of field functions. The bounce solution is then obtained by maximizing S∗(R), where S∗(R) is the minimum action with the constraint that the monopole worldline has a constant radius R.

Boundary Conditions and Constraints

To ensure regularity and physical consistency in the numerical simulation, specific boundary conditions are imposed at different spatial boundaries. For instance, at the core of the string (ρ = 0), Dirichlet conditions are required for certain scalar fields (a = b = h = d = u3 = u4 = u5 = u6) to ensure regularity, while Neumann boundary conditions are imposed on other functions (e.g., ∂ρc, ∂ρf, etc.). These constraints lead to specific requirements at ρ=0, such as the cancellation of a divergence in the gauge-invariant term F aµνF µν a.

Simulation Results and Implications

The simulation results show a significant suppression of the bounce action (and hence an increase in the decay rate) as we move from hierarchical to near degenerate symmetry breaking scales. The continuum limit of these simulations is well-fitted by Eq. (4.1), which provides the exact solution for the bounce action for a benchmark case, and its asymptotic value in the thin string limit is consistent with analytical estimates. This improved calculation resolves issues encountered in previous work where solutions exceeded the thin string limit for large mW/mγ ratios, confirming that the found solution is indeed the minimal energy one.

Implications for Gravitational Wave Searches

The calculated decay rate, governed by SB = πκeff (Eq. 4.4), directly impacts the stochastic gravitational wave background (SGWB) spectrum. The finite string width reduces the bounce action compared to the thin string limit, leading to a corresponding enhancement of the string decay rate. This shift moves the turnover point of the GW spectrum to larger frequencies, which is crucial for interpreting tentative signals from pulsar timing arrays (PTAs), LISA, and LIGO/Virgo/KAGRA. The results map preferred regions in parameter space onto observable GW constraints, showing that the full solution provides a more accurate picture than the thin string approximation alone.

Conclusion

The work demonstrates that going beyond the thin string approximation is necessary for accurately determining the decay rate of metastable cosmic strings, especially in near-degenerate symmetry breaking regimes. The numerical method successfully recovers the thin string limit while providing a corrected expression for the bounce action, which has significant implications for constraining GUT models and interpreting astrophysical gravitational wave observations. The results suggest that a given string lifetime corresponds to larger mass hierarchies between symmetry breaking scales, simplifying concrete implementations in cosmological models.

The gist: Classical lattice simulations reveal a suppression of the bounce action compared to previous estimates, indicating a faster string decay than previously thought.

How it works

  1. The setup involves an SU(2) gauge theory with two successive symmetry breaking stages driven by adjoint and doublet scalar fields, leading to 't Hooft-Polyakov monopoles and cosmic strings.

  2. The decay rate is governed by the bounce action SB, calculated via gradient flow to find the saddle point of the action associated with monopole formation on the string.

Improvements for AI systems

As a fastidious and diligent researcher, I have thoroughly reviewed this paper, The decay rate of metastable cosmic strings beyond the thin-string approximation, by Domcke and Hamada. This work provides a crucial theoretical bridge between simplified models (thin string approximation) and more realistic numerical simulations (gradient flow) for calculating the decay rate of cosmic strings.

Here are the specific improvements that can be made to AI systems, along with what those improved systems could achieve:

)

  1. AI Systems can perform highly accurate, non-perturbative calculations of topological defect dynamics in Grand Unified Theories (GUTs).

  2. AI Systems can predict the precise spectral characteristics of Stochastic Gravitational Wave Backgrounds (SGWB) from metastable cosmic strings, including frequency cutoffs determined by string decay rates.

  3. AI Systems can map complex model parameter spaces (e.g., mass ratios and couplings) to observable gravitational wave signatures in real-time, allowing for rapid hypothesis testing of high-energy physics scenarios.

Specific Capabilities of the Improved AI System:

  1. [Calculation of Bounce Action and Decay Rate] The AI system can numerically solve the non-linear field equations (using advanced techniques like gradient flow simulations described in Section 3) to determine the exact bounce action for monopole nucleation on a string, moving beyond approximations like the thin-string limit.

  2. [Phenomenological Parameter Space Mapping] Given a set of fundamental model parameters (like symmetry breaking scales, gauge couplings, and mass ratios), the AI can rapidly compute the effective decay rate and map this onto predicted gravitational wave spectra (Section 4.2). This allows researchers to instantly determine if a given GUT model is viable for detection by LIGO/Virgo/KAGRA or Pulsar Timing Arrays (PTA) based on the calculated SGWB amplitude and frequency dependence.

  3. [Constraint Inference from Observational Data] The system can ingest observational constraints (like the tentative NANOGrav signal or LVK upper bounds shown in Figure 10) and use Bayesian inference to constrain the underlying microphysical parameters of the cosmic string model (e.g., inferring preferred values for the effective parameter κeff).

  4. [Model Selection for Inflationary Scenarios] The AI can identify specific mass hierarchies (e.g., mild hierarchy between non-abelian and abelian sectors, as highlighted in Section 5) that maximize the string lifetime while maintaining a decay rate consistent with observed SGWB constraints, thus guiding the construction of viable hybrid inflation models.

  5. [Distinguishing Physical Mechanisms] The system can perform comparative analysis between different decay mechanisms (e.g., monopole creation vs. thermal production, as discussed in Section 5) by analyzing how different physical assumptions affect the resulting GW spectral tilt and amplitude predictions, helping to isolate which physical process dominates the observed signal in astrophysical data.

Sources

Related papers