Pith. sign in

REVIEW 6 minor 52 references

The paper claims that the Abelian–Higgs cosmic string carries a localized massive gauge-field mode — the vector mode — that remains bound for every value of the scalar self-coupling λ, decays by a calculable power law, and resonantly transf

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-01 21:16 UTC pith:IAT3Y73V

load-bearing objection Solid, honest paper that completes the vector-mode story for Abelian-Higgs strings; the CCM is the weak link but the field-theory evidence stands.

arxiv 2607.16120 v1 pith:IAT3Y73V submitted 2026-07-17 hep-th astro-ph.COhep-ph

Localized Vector Modes on Local Cosmic Strings

classification hep-th astro-ph.COhep-ph PACS 11.27.+d11.15.-q
keywords cosmic stringsAbelian-Higgs modelvector modesinternal modes of solitonsparametric instabilityquasinormal modesshape modecollective coordinate model
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper aims to establish the vector mode — a fluctuation of the gauge field along the string axis — as a genuine, long-lived internal degree of freedom of Abelian–Higgs cosmic strings. The authors argue that, unlike the scalar shape mode, this vector mode never merges into the bulk continuum; it remains a bound state for arbitrarily large λ, so for λ > 1.5 it is the only massive discrete internal mode of the string. They derive and verify a characteristic power-law decay of the mode's amplitude through nonlinear emission of radiation, and show that when the scalar radiation channel closes, quasinormal modes take over as Feshbach resonances. They then demonstrate, analytically and in 3+1-dimensional field-theory simulations, that the vector mode triggers a parametric instability in the translational zero modes at resonance k₀ = ω_z/2, exciting both transverse directions simultaneously — a behavior their effective model explains via a cubic coupling between the two zero modes and the vector mode.

Core claim

The central discovery is that the longitudinal gauge-field fluctuation A_z = C(t)ψ(r) is a solution of the single-channel Schrödinger-type problem −ψ'' − ψ'/r + f(r)²ψ = ω_z²ψ, with f(r) the Higgs profile of the vortex. Because the Higgs condensate vanishes on the string axis, the effective gauge mass is suppressed inside the core, making the potential f(r)² attractive; a classical theorem for two-dimensional Schrödinger operators guarantees at least one bound state for any λ. The paper's key non-perturbative claim is that this bound state persists for arbitrarily large λ, becoming an increasingly shallow, broad mode as the core narrows, while the shape mode disappears at λ ≈ 1.5. On the dyn

What carries the argument

The load-bearing object is the vector mode profile ψ(r), defined by the radial eigenvalue equation −ψ'' − ψ'/r + f(r)²ψ = ω_z²ψ; the effective potential f(r)² is the 'soft waveguide' that traps the gauge field where the Higgs condensate is suppressed. Around it the paper builds a four-degree-of-freedom collective-coordinate model with amplitudes X (x-zero mode), Y (y-zero mode), C (vector mode), and S (shape mode). The crucial new mechanism is the cubic term I_XYC X Y C in the effective Lagrangian, which arises only when the y-modulation is taken with a relative phase sin(k₀z) rather than cos(k₀z); this term converts the naive Mathieu resonance at k₀ = ω_z into the observed one at k₀ = ω_z/2

Load-bearing premise

The explanatory model of the instability assumes a four-mode truncation in which the y-zero mode carries a sin(k₀z) modulation, phase-shifted by π/2 relative to the x-zero mode, so that a cubic X Y C term exists and produces the k₀ = ω_z/2 resonance; if this phase choice or the truncation is not robust, the mechanism explaining the observed instability would change, although the instability itself is a simulation result.

What would settle it

In a 3+1-dimensional simulation with periodic z-boundaries, scan k₀ across ω_z/2 for fixed λ, amplitude C₀, and initial X₀; the instability band should be centered at k₀ = ω_z/2 with width set by C₀ I_XYC/(2 I_dX²). Observing growth at k₀ = ω_z instead, or no y-zero-mode component with sin(k₀z) modulation appearing, would falsify the cubic-coupling mechanism. Alternatively, measuring the zero-mode growth rate for λ ≫ 1 and checking that the resonance frequency tracks ω_z(λ)/2 for all λ would test whether the instability is generic or an artifact of the BPS degeneracy.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • For λ > 1.5, after the shape mode merges into the continuum, the vector mode is the only massive discrete internal excitation of Abelian–Higgs strings, making it the dominant channel for internal energy storage in type-II regimes.
  • Excited vector modes relax by a universal power-law amplitude decay Ĉ ~ 1/√(Γt), with a decay rate that peaks near λ ≈ 3.8 and then falls off as quasinormal modes take over.
  • A vector-excited string is parametrically unstable at wavenumber k₀ = ω_z/2, transferring energy to both transverse Goldstone modes; this provides a microscopic mechanism for converting internal oscillations into string wiggles and gravitational-wave-carrying motion.
  • In the BPS case the instability can cascade to the shape mode through a second resonance, giving a vector→zero→shape energy-transfer chain; away from BPS this cascade is generically absent while the zero-mode instability persists.
  • The persistence of the vector mode at all λ means any complete low-energy effective theory of Abelian–Higgs strings must include a massive vector field on the worldsheet, with tensorial couplings distinct from the shape mode's scalar coupling.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the vector mode is as long-lived as the shape mode and survives at large λ, then cosmic-string network simulations that ignore internal modes may systematically miss energy stored in these excitations; quantifying this requires evolving loops with vector-mode initial data, which the paper does not do.
  • The resonance condition k₀ = ω_z(λ)/2 is a parameter-free prediction that could be tested on closed string loops, where discrete wavenumbers sweep as the loop shrinks; one would expect bursts of transverse oscillation whenever a loop's mode crosses the half-frequency condition.
  • The 1+1-dimensional decay law and the 3+1-dimensional instability are computed in flat, straight strings; on curved or oscillating strings the vector mode may mix with worldsheet curvature, possibly modifying the decay rate or triggering the instability sooner — an extension left implicit.
  • A direct observational consequence, if internal modes survive in networks, is a reshuffling of the gravitational-wave spectrum: more power in transverse motion could alter loop chopping and gravitational-wave emission rates relative to the Nambu–Goto picture.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 6 minor

Summary. This paper studies the spectrum and nonlinear dynamics of localized vector excitations (longitudinal A_z fluctuations) on Abelian-Higgs cosmic strings. The authors compute the vector-mode eigenfrequency and profile from the radial equation (2.25), and show that, unlike the scalar shape mode, the mode never merges into the continuum for any λ. The nonlinear decay of the mode is analyzed by expanding around the vortex plus radiation fields; the resulting power-law envelope law, Ĉ(t)=1/√(Γ(λ)t+Ĉ(0)^{-2}) (Eq. 3.21), is derived from a semianalytical treatment of the inhomogeneous radial equations (3.9)-(3.10) and verified against 1+1D field-theory simulations (Fig. 7). For λ≳3.8 the scalar radiation channel closes and the decay proceeds through a quasinormal/Feshbach mode; the paper shows that keeping only the vector channel gives a good approximation (Fig. 6). In 3+1D, the paper reports a parametric instability between the vector mode and the transverse zero modes, with resonance at k0=ω_z/2, exciting both x and y zero modes. A 2-DoF collective-coordinate model predicts the wrong resonance (k0=ω_z), while a 4-DoF model including the second zero mode and the shape mode, with a sin(k0z) modulation for the y mode, reproduces the observed resonance and the early-time dynamics. The paper candidly reports the limitations of its effective models.

Significance. If correct, this is a significant contribution. The vector mode provides a massive discrete internal degree of freedom of local strings that persists for all λ, and it can mediate a new, qualitatively different energy-transfer channel between the internal and translational sectors. The main quantitative outputs—Eq. (3.21) and the resonance condition k0=ω_z/2—are concrete, falsifiable predictions, and the paper supports them with an independent semianalytical calculation and dedicated field-theory simulations. Particularly commendable is the honest reporting of the failure of the initial 2-DoF CCM, and the explicit identification of the cubic XYC mechanism in the 4-DoF model. The numerical methods (Appendix C) are described in sufficient detail for reproducibility. The stress-test concern about the relative phase δ in the 4-DoF CCM does not, on close reading, undermine the central claims: the phase is motivated by the structure of the source term 2iA_z∂_zφ in the scalar equation, and the parametric instability itself is a direct field-theory simulation result, independent of the CCM.

minor comments (6)
  1. [Sec. 5, Conclusions] The statement that for λ>1.5 the vector mode 'constitutes the only discrete internal excitation supported by the string' is broader than the evidence presented. The spectral analysis in Sec. 2 covers the longitudinal vector sector (Eq. 2.25) and the scalar–angular shape-mode sector (Eq. 2.13), but not all possible angular-momentum or coupled channels. This overclaim is not needed for the paper's main quantitative results, and I recommend qualifying it, e.g., 'among the sectors studied here' or 'the only known discrete massive mode.'
  2. [Sec. 4, final paragraph] The claim that the vector–zero-mode instability is 'a general feature of the model, independent of the chosen λ' is supported only by simulation at λ=1 (BPS). The CCM argument shows the resonance condition k0=ω_z(λ)/2, but the coupling I_XYC in Eq. (D.21) could in principle vanish or change sign with λ; no non-BPS simulation or explicit check of I_XYC(λ) is provided. I suggest softening the wording or adding a non-BPS representative run.
  3. [Sec. 2.1, Eqs. (2.21)-(2.22)] The notation in Eqs. (2.21)-(2.22) is confusing: the same symbol f is used for the background profile and for the equation written after the perturbation is introduced, and the factor (1-a²) in Eq. (2.21) should presumably be (1-a)² to match the vortex equation (2.5). Please clarify that these are the background equations augmented by the A_z²φ term, not a coupled perturbation system.
  4. [Eq. (3.12)] In the boundary-condition line for R_θ, the wavenumber inside k_ϕ = √(4ω_z²−1) should be k_θ. This is a typographical slip but could confuse readers.
  5. [Appendix B] The Sturm-comparison argument for an arbitrarily large number of bound states as λ→0 is somewhat sketchy: the comparison potential V(r) and the constant c are not fully specified, and the matching procedure is only outlined. Since this is an auxiliary point, a brief clarification would suffice.
  6. [Sec. 3, quasinormal regime] In the discussion following Eq. (3.22), the paper states that the approximation (3.23) assumes the energy to excite the scalar component of the quasinormal mode is negligible. This is supported by Fig. 6, but it would be useful to state explicitly that this is an assumption for λ≳3.8 and to comment on its expected range of validity.

Circularity Check

0 steps flagged

No significant circularity: central quantitative claims are derived from the vortex fluctuation equations and checked against independent field-theory simulations.

full rationale

I walked the derivation chain. The vector-mode spectrum (Sec. 2.1, Eq. (2.25)) is a direct eigenvalue problem in the vortex background; the persistence of the bound state is supported by an external theorem (Ref. [39]) and by the Sturm-comparison argument of Appendix B, not by a self-referential fit. The decay law of Sec. 3 is derived from the sourced linearized fluctuation equations (3.9)-(3.10) together with energy conservation, yielding the power-law (3.21); the coefficients Aphi and Atheta are obtained by solving those inhomogeneous equations, and the resulting Gamma(lambda) is then compared with independent 1+1D field-theory simulations (Figs. 6-7). This is a genuine semianalytic prediction, not a fitted parameter renamed as a prediction. The parametric instability of Sec. 4 is established directly by 3+1D field-theory simulations (Figs. 9-13). The paper candidly reports that its initial 2-DoF collective-coordinate model fails, predicting k0 = omega_z while the simulation shows k0 = omega_z/2 (Fig. 9). The subsequent 4-DoF extension is constructed after observing that the y-zero mode acquires a sin(k0 z) modulation; importantly, the paper derives this modulation from the field equations via the 2 i A_z partial_z phi term rather than fitting it to the final resonance condition, and it verifies the expectation in snapshots. The extended model then makes new, testable predictions, such as the chain resonance vector -> zero -> shape, which are confirmed in Fig. 12. This is transparent model revision, not circularity: the central instability result does not depend on the CCM, and the CCM phase choice is not used to manufacture the main quantitative claims. Self-citations to Refs. [13] and [32] are contextual or methodological and are not the sole support of any central result. I find no step in which a claimed output reduces by construction to an input.

Axiom & Free-Parameter Ledger

1 free parameters · 4 axioms · 0 invented entities

The central spectrum and decay computation are self-contained given the numerically computed vortex background (eqs. 2.5-2.6) and carry no fitted constants: Γ(λ) follows from solving the inhomogeneous linear system (3.9)-(3.10) and is then compared with independent time-domain simulation. The two hand-chosen ingredients are the same-symmetry radiation truncation (3.1)-(3.3) and the phase assignment of the second zero mode in the extended CCM (Sec. 4); both are stated explicitly in the text. No new entities are introduced; the vector mode is a known fluctuation (Refs. [28-31]), and the 'waveguide' and 'Feshbach resonance' language is interpretive description of existing spectral states.

free parameters (1)
  • Relative phase δ of the y-zero-mode modulation in the 4-DoF CCM = π/2 (sin(k₀z) for Y vs cos(k₀z) for X)
    Introduced in eqs. (4.16)-(4.19) so that the cubic term X(t)Y(t)C(t) survives in (4.23); the sin modulation is adopted because field-theory snapshots show Y excited with that phase ('It is therefore natural to adopt this modulation', Sec. 4). The value is post-hoc-informed, but the model is subsequently tested against the full simulation (Fig. 12).
axioms (4)
  • standard math A localized attractive potential in 2D supports at least one bound state (Simon 1976; Yang & de Llano 1989)
    Invoked in Sec. 2.1 to conclude the vector mode remains bound for arbitrarily large λ; also the Sturm-comparison argument in Appendix B for the λ→0 count of bound states.
  • domain assumption The linearized A_z fluctuation decouples from the scalar and A_θ sectors (perturbations 2.18-2.20 keep ϕ and A_{x,y} at background values)
    Eqs. (2.21)-(2.23): at linear order f and a stay background, leaving the single-channel eigenvalue problem (2.25) for the vector mode.
  • domain assumption Radiation emitted in the decay has the same cylindrical symmetry as the vortex, and truncation at O(C²) captures the dominant channel
    Ansatz (3.1)-(3.3) drops non-symmetric radiation and higher-order sources; the resulting t^{-1/2} law (3.21) is cross-checked in Figs. 6-7.
  • ad hoc to paper The 4-DoF collective-coordinate ansatz (4.16)-(4.19) contains the modes needed to describe the parametric instability
    The 2-DoF ansatz fails (resonance at k₀=ω_z predicted, ω_z/2 observed, Fig. 9); the extended ansatz adds Y and S and fixes the phase; agreement is demonstrated only up to the first zero-mode oscillation (Fig. 12).

pith-pipeline@v1.3.0-alltime-deepseek · 28864 in / 19991 out tokens · 159865 ms · 2026-08-01T21:16:43.875602+00:00 · methodology

0 comments
read the original abstract

We investigate the dynamics of vector excitations localized on cosmic strings in the Abelian-Higgs model. These confined gauge field fluctuations behave as massive vector degrees of freedom propagating along the string. We show that they remain bounded for all values of the scalar self coupling $\lambda$, including parameter regimes in which the previously studied conventional shape mode is no longer present. After determining their spectra and spatial profiles, we analyze their decay through nonlinear coupling to bulk radiation. Their amplitudes exhibit the characteristic power law relaxation associated with non-linear radiation. For sufficiently large $\lambda$ the scalar radiation channel becomes kinematically inaccessible and the decay is instead mediated by quasinormal modes which can be interpreted as Feshbach resonances of the string. We then study the full $3+1$ dimensional dynamics of the vector-excited string. Field theory simulations reveal a parametric instability between the vector mode and the Goldstone sector. In contrast to the shape mode, the vector mode resonantly excites both transverse directions simultaneously. In addition, we build an effective model that captures the instability and identifies the nonlinear couplings responsible for it. Our results show that vector excitations constitute a long-lived massive degree of freedom with a non-trivial role in the dynamics of local strings.

Figures

Figures reproduced from arXiv: 2607.16120 by Jose J. Blanco-Pillado, Jose Queiruga, Sergio Alameda-Calvo.

Figure 1
Figure 1. Figure 1: FIG. 1: Profile functions [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Spectrum of the vector and shape modes as a function of [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Radial profiles of the normalized vector mode for different [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Radial profiles of the fields and current density for different values of [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Pictorial representation of the magnetic field lines of the excited string (blue tube). [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Decay constant of the vector mode amplitude as a function of the dimensionless [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Amplitude of the vector mode, [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Snapshots of the radial profiles of the scalar and gauge field perturbations extracted [PITH_FULL_IMAGE:figures/full_fig_p019_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Evolution of the extracted vector and zero mode amplitudes in the BPS case for [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Snapshots of the string configuration during the parametric resonance. The left [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: Time evolution of the amplitudes of the vector, shape and zero modes, extracted [PITH_FULL_IMAGE:figures/full_fig_p026_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Time evolution of the amplitudes of the vector, shape and zero modes. The blue [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Snapshots of the string configuration before and during the parametric resonance. [PITH_FULL_IMAGE:figures/full_fig_p028_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Amplitude of the vector mode, [PITH_FULL_IMAGE:figures/full_fig_p036_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Complementary snapshots to Fig. [PITH_FULL_IMAGE:figures/full_fig_p042_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16: Complementary snapshots to Fig. [PITH_FULL_IMAGE:figures/full_fig_p043_16.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 21 linked inside Pith

  1. [1]

    1+1 dimensions In section 3, the cylindrical symmetry of the problem at hand has allowed us to reduce system’s dynamics to the following set of radial equations: ∂2F ∂t2 − ∂2F ∂r 2 − 1 r ∂F ∂r + 1 r2 (1−B) 2F+F A 2 z − λ 2 (1−F 2)F= 0,(C.1) ∂2B ∂t2 − ∂2B ∂r 2 + 1 r ∂B ∂r −(1−B)F 2 = 0,(C.2) ∂2Az ∂t2 − ∂2Az ∂r 2 − 1 r ∂Az ∂r +A zF 2 = 0,(C.3) whereF(r, t) ...

  2. [2]

    CONVENTIONS AND SPECTRUM

    THE MODEL. CONVENTIONS AND SPECTRUM. We will consider the Abelian-Higgs model in 3+1 dimensions. In dimensionless variables, the model is defined by the following Lagrangian density L=− 1 4 FµνF µν + 1 2 (Dµϕ)∗Dµϕ− λ 8 (1−ϕ ∗ϕ)2 .(2.1) The covariant derivative for the complex scalar fieldϕis defined asD µ ≡∂ µ −iA µ and Fµν =∂ µAν −∂ νAµ is the field stre...

  3. [3]

    lump-like

    VECTOR MODE DECA Y In this section, we want to study both analytically and numerically the non-linear decay of the vector mode due to its higher order coupling to radiation modes. With that purpose in mind, let us extend the previous parametrizations of the fields in order to include radiation 12 fields ϕ(r, θ, t) =eiθf(r) +e iθRϕ(r, t),(3.1) Aθ(r, t) =a(...

  4. [4]

    cross coupling

    INTERACTION AND INST ABILITY: VECTOR VS. ZERO MODE Parametric instabilities are common phenomena, perhaps universal, in the dynamics of several types of (excited) strings [13, 19, 21], where a resonant energy transfer takes place between a massive internal mode of the string and the Goldstone sector (i.e. a transversal zero mode). This feature accelerates...

  5. [5]

    A way of making Europe

    CONCLUSIONS In this work, we have investigated the excitation spectrum of local strings in the Abelian–Higgs model, with particular emphasis on vector modes. This internal degree of freedom originates from gauge field fluctuations localized within the vortex core, where the suppression of the Higgs condensate reduces the effective gauge field mass relativ...

  6. [6]

    3+1 dimensions In this case, the setup is no longer reduced by exploiting any symmetry of the system, and the vortex string is no longer pinned down to the origin, allowing for a wider range of dynamics. 36 The set of equations that needs to be solved in the 3D cartesian grid is ∂2 t ϕ1 =∂ i∂iϕ1 +ϕ 2∂iAi + 2Ai∂iϕ2 −ϕ 1AiAi + λ 2 1−(ϕ 2 1 +ϕ 2 2) ϕ1 ,(C.8)...

  7. [7]

    For the zero mode IX2 =I ∥ X2 +I ⊥ X2 , I ⊥ X2 = 0,(D.3) since it is an eigenfunction of zero eigenvalue of the transverse fluctuation operator

    Effective model with 2 DoF The integrals appearing in the effective Lagrangian (4.8) are the following: IdX2 = 1 2 Z L/2 −L/2 cos2(k0z)dz Z 2π 0 Z ∞ 0 |δxϕ(r, θ)|2 +|δ xAr(r, θ)|2 +|δ xAθ(r, θ)|2 rdrdθ , (D.1) IdC2 = 1 2 Z L/2 −L/2 dz Z 2π 0 dθ Z ∞ 0 ψ2rdr .(D.2) For the integrals that are purely quadratic in their amplitudes, it is useful to separate the...

  8. [8]

    |δxϕ|2 α2 2 + λ 4 φ2 + φ2 2 (δxAr)2 + (δxAθ)2 + 2φαf ′a′ r cos2 θ+ λ 2 f ′2φ2 cos2 θ # ,(D.25) IY2S2 =− Z dzsin 2(k0z) cos2(ksz) Z 2π 0 dθ Z ∞ 0 rdr

    Effective model with 4 DoF We now include the second zero mode,Y(t), and the shape modeS(t). The coefficients al- ready defined in the 2 DoF model remain unchanged, so below we will just list the additional terms: IdY2 = 1 2 Z L/2 −L/2 sin2(k0z)dz Z 2π 0 Z ∞ 0 |δyϕ|2 +|δ yAr|2 +|δ yAθ|2 rdrdθ ,(D.10) IdS2 = 1 2 Z L/2 −L/2 dzcos 2(ksz) Z 2π 0 dθ Z ∞ 0 rdr(...

  9. [9]

    T. W. B. Kibble, J. Phys. A9, 1387 (1976)

  10. [10]

    Vilenkin and E

    A. Vilenkin and E. P. S. Shellard,Cosmic Strings and Other Topological Defects(Cambridge University Press, 2000), ISBN 978-0-521-65476-0

  11. [11]

    Abbott et al

    R. Abbott et al. (LIGO Scientific, Virgo, KAGRA), Phys. Rev. Lett.126, 241102 (2021), 2101.12248

  12. [12]

    J. J. Blanco-Pillado, Y. Cui, S. Kuroyanagi, M. Lewicki, G. Nardini, M. Pieroni, I. Y. Rybak, L. Sousa, and J. M. Wachter (LISA Cosmology Working Group), JCAP05, 006 (2025), 2405.03740

  13. [13]

    Afzal et al

    A. Afzal et al. (NANOGrav), Astrophys. J. Lett.951, L11 (2023), [Erratum: Astrophys.J.Lett. 971, L27 (2024), Erratum: Astrophys.J. 971, L27 (2024)], 2306.16219

  14. [14]

    Auclair et al., JCAP04, 034 (2020), 1909.00819

    P. Auclair et al., JCAP04, 034 (2020), 1909.00819

  15. [15]

    Abac et al

    A. Abac et al. (ET), JCAP03, 081 (2026), 2503.12263

  16. [16]

    Goto, Prog

    T. Goto, Prog. Theor. Phys.46, 1560 (1971)

  17. [17]

    Nambu, Symmetries and Quark Models, ed

    Y. Nambu, Symmetries and Quark Models, ed. R. Chand (Gordon and Breach, 1970)269 (1970)

  18. [18]

    Maeda and N

    K.-i. Maeda and N. Turok, Phys. Lett. B202, 376 (1988)

  19. [19]

    Gregory, Phys

    R. Gregory, Phys. Lett. B206, 199 (1988)

  20. [20]

    M. R. Anderson, F. Bonjour, R. Gregory, and J. Stewart, Phys. Rev. D56, 8014 (1997), hep-ph/9707324

  21. [21]

    J. C. Aurrekoetxea, J. J. Blanco-Pillado, A. Garc ´ ıa Mart ´ ın-Caro, and J. M. Queiruga (2026), 2603.05243

  22. [22]

    N. S. Manton and P. Sutcliffe,Topological solitons, Cambridge Monographs on Mathematical Physics (Cambridge University Press, 2004), ISBN 978-0-521-04096-9, 978-0-521-83836-8, 978- 0-511-20783-9

  23. [23]

    N. S. Manton and H. Merabet, Nonlinearity10, 3 (1997), hep-th/9605038

  24. [24]

    J. J. Blanco-Pillado, D. Jim´ enez-Aguilar, and J. Urrestilla, JCAP01, 027 (2021), 2006.13255

  25. [25]

    Sugiyama, Prog

    T. Sugiyama, Prog. Theor. Phys.61, 1550 (1979)

  26. [26]

    D. K. Campbell, J. F. Schonfeld, and C. A. Wingate, Physica D9, 1 (1983)

  27. [27]

    J. J. Blanco-Pillado, D. Jim´ enez-Aguilar, J. M. Queiruga, and J. Urrestilla, JCAP05, 011 (2023), 2209.12945

  28. [28]

    J. J. Blanco-Pillado, A. Garc ´ ıa Mart ´ ın-Caro, D. Jim´ enez-Aguilar, and J. M. Queiruga, Phys. Rev. D111, 056007 (2025), 2411.13521

  29. [29]

    J. J. Blanco-Pillado, D. Jim´ enez-Aguilar, J. M. Queiruga, and J. Urrestilla, JCAP04, 043 44 (2023), 2212.06194

  30. [30]

    Alonso-Izquierdo, J

    A. Alonso-Izquierdo, J. Mateos Guillarte, M. Rees, and A. Wereszczynski, Phys. Rev. D110, 065004 (2024), 2406.05725

  31. [31]

    Bachmaier and A

    M. Bachmaier and A. Wereszczynski, Phys. Lett. B875, 140324 (2026), 2510.17964

  32. [32]

    Alonso-Izquierdo, N

    A. Alonso-Izquierdo, N. S. Manton, J. Mateos Guilarte, M. Rees, and A. Wereszczynski, Phys. Rev. D111, 105021 (2025), 2502.15087

  33. [33]

    Alonso Izquierdo, M

    A. Alonso Izquierdo, M. Bachmaier, and A. Wereszczynski, Phys. Rev. D114, 015015 (2026), 2603.04495

  34. [34]

    Hindmarsh, J

    M. Hindmarsh, J. Lizarraga, A. Urio, and J. Urrestilla, Phys. Rev. D104, 043519 (2021), 2103.16248

  35. [35]

    J. J. Blanco-Pillado, K. D. Olum, and B. Shlaer, Phys. Rev. D89, 023512 (2014), 1309.6637

  36. [36]

    R. L. Davis, Nucl. Phys. B294, 867 (1987)

  37. [37]

    Arodz, Acta Phys

    H. Arodz, Acta Phys. Polon. B22, 511 (1991)

  38. [38]

    Arodz and L

    H. Arodz and L. Hadasz, Phys. Rev. D54, 4004 (1996), hep-th/9506021

  39. [39]

    Arodz and L

    H. Arodz and L. Hadasz, Phys. Rev. D55, 942 (1997), hep-th/9607116

  40. [40]

    Alonso-Izquierdo, J

    A. Alonso-Izquierdo, J. J. Blanco-Pillado, D. Migu´ elez-Caballero, S. Navarro-Obreg´ on, and J. Queiruga, Phys. Rev. D110, 065009 (2024), 2405.06030

  41. [41]

    A. A. Abrikosov, Sov. Phys. JETP5, 1174 (1957)

  42. [42]

    H. B. Nielsen and P. Olesen, Nucl. Phys. B61, 45 (1973)

  43. [43]

    Goodband and M

    M. Goodband and M. Hindmarsh, Phys. Rev. D52, 4621 (1995)

  44. [44]

    T. Kojo, H. Suganuma, and K. Tsumura, Phys. Rev. D75, 105015 (2007), hep-ph/0702014

  45. [45]

    Alonso-Izquierdo, W

    A. Alonso-Izquierdo, W. Garcia Fuertes, and J. Mateos Guilarte, Phys. Lett. B753, 29 (2016), 1509.06632

  46. [46]

    Alonso-Izquierdo, W

    A. Alonso-Izquierdo, W. Garcia Fuertes, and J. Mateos Guilarte, JHEP05, 074 (2016), 1602.09084

  47. [47]

    Simon, Annals Phys.97, 279 (1976)

    B. Simon, Annals Phys.97, 279 (1976)

  48. [48]

    Yang and M

    K. Yang and M. de Llano, Am. J. Phys.57, 85 (1989)

  49. [49]

    J. D. Jackson,Classical Electrodynamics(Wiley, 1998), ISBN 978-0-471-30932-1

  50. [50]

    Cheng and L.-F

    T.-P. Cheng and L.-F. Li,Gauge Theory of Elementary Particle Physics(Oxford University Press, Oxford, UK, 1984), ISBN 978-0-19-851961-4

  51. [51]

    Videos of the complete field-theory simulations are available as ancillary files in the arXiv submission

  52. [52]

    Mur, IEEE Transactions on Electromagnetic CompatibilityEMC-23, 377 (1981)

    G. Mur, IEEE Transactions on Electromagnetic CompatibilityEMC-23, 377 (1981). 45