Pith. sign in

REVIEW 3 major objections 5 minor 57 references

Networks of Z_N-charged cosmic strings with baryon-like junctions settle into a scaling regime rather than jamming, and the string density per horizon volume grows in proportion to N²−1.

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 06:56 UTC pith:HC343EDE

load-bearing objection New simulation evidence that global Z_N string networks with junctions scale with density ~(N^2−1), but the winding reconstruction that carries the result needs to be specified and validated before I would trust the normalization. the 3 major comments →

arxiv 2607.21701 v1 pith:HC343EDE submitted 2026-07-23 hep-ph astro-ph.COhep-th

Formation and scaling of mathbb{Z}_N strings for global SU(N)/mathbb{Z}_N symmetry

classification hep-ph astro-ph.COhep-th
keywords cosmic stringsZ_N symmetryglobal stringsscaling regimejunctionslattice simulationgravitational wavesPSU(N) vacuum
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.

This paper asks whether cosmic-string networks whose strings carry Z_N charges and meet at baryon-vertex-like junctions can settle into a scaling regime, or whether the junctions keep the network frustrated and unable to dilute. Using classical lattice simulations of a global SU(N)/Z_N scalar model in a radiation-dominated universe, it finds for N=2,3,4,5,8 that the networks do scale: the total string length per horizon volume approaches a plateau. The plateau value grows in proportion to the adjoint dimension N²−1, meaning the number of long strings per Hubble volume is of order N. This matters because a scaling network with this normalization gives a gravitational-wave background amplitude Ω_GW ∝ μ²(N²−1)², a concrete prediction for pulsar timing arrays.

Core claim

The central claim is that Z_N-charged global string networks with charge-conserving junctions are not frustrated: for N=2,3,4,5,8, the winding-based string density ζ approaches a plateau, and the plateau value divided by N²−1 is roughly constant (about 0.5 for fixed-comoving-core runs and 0.6 for physical-core runs). Non-minimal charge strings are subdominant, contributing 12–14% of the total length for N=4,5 and only a few percent for N=8. Because vertices can annihilate rather than permanently tie the network together, junction structure does not prevent scaling; instead the number of vertices per Hubble volume grows as O(N).

What carries the argument

The construction rests on a scalar-only model whose vacuum manifold is PSU(N)=SU(N)/Z_N. Three adjoint scalars with an F-term-like potential realize this vacuum, and π₁(PSU(N))=Z_N classifies strings by a charge q∈Z_N, allowing N unit-charge strings to meet at a center-neutral junction. The numerical diagnostic is a center-valued plaquette winding: at each lattice site a local SU(N) representative is fixed, link variables are projected onto the nearest integer modulo N, and a plaquette with nonzero winding is counted as a pierced string. The scaling parameter ζ counts such pierced plaquettes, converting Manhattan length to Euclidean length, and the empirical law ζ≈const×(N²−1) is the paper's

Load-bearing premise

The scaling observable ζ assumes that at every lattice site away from the string core one can reconstruct a single SU(N) representative U(x) from the three adjoint fields, but near cores the radial norm is not small (up to ~0.8 v²), so the reconstruction is ill-defined there and the paper does not specify how such sites are excluded or how the simultaneous U-reconstruction is solved numerically.

What would settle it

Compare the winding-based string length ζ with a core-based length obtained by counting sites where φ_r² falls below a chosen threshold, across N and over time; if the two measures disagree increasingly with N, or if ζ/(N²−1) drifts when the reconstruction cutoff is varied, the scaling law would be a numerical artifact.

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

If this is right

  • For N≤8, global Z_N string networks enter a scaling regime rather than a frustrated, string-dominated state.
  • The long-string density per horizon volume grows as N²−1, so the number of baryon-like junctions per Hubble volume is O(N), not O(1).
  • Non-minimal charge strings (q≥2) remain subdominant, so unit-charge strings dominate the network.
  • The gravitational-wave amplitude from such networks scales as Ω_GW ∝ μ²(N²−1)², with an additional (N−1)² factor if unit-string tension follows μ∝q(N−q).
  • Scaling behavior supports using standard scaling-network formulas for Yang-Mills-motivated confining string networks with junctions.

Where Pith is reading between the lines

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

  • If the N²−1 normalization persists for larger N, simulations with N=16 or N=32 could test whether the network remains in scaling or whether core-overlap and volume constraints change the behavior.
  • The tension's logarithmic growth in a global model may enhance gravitational-wave radiation relative to local strings, so a gauge-version simulation would separate the junction effect from the global-string effect.
  • The winding reconstruction assumes a well-defined SU(N) representative U(x) even near cores, where radial norms are large (up to ~0.8 v² for N=5); cross-checking with an energy-density-based core counter would test whether the measured scaling law is robust.
  • If scaling also holds for local PSU(N) flux tubes, pure Yang-Mills explanations of pulsar-timing-array backgrounds would need to include a (N²−1)² enhancement in the predicted signal amplitude.

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

3 major / 5 minor

Summary. The paper constructs a scalar field theory with three adjoint SU(N) fields whose Higgs vacuum manifold is PSU(N)=SU(N)/Z_N, motivated by the string-junction structure of mass-deformed N=4 SYM / pure Yang-Mills. It reports classical lattice simulations of the formation and evolution of global Z_N string networks in a radiation-dominated universe, reconstructing strings via a center-valued plaquette-winding diagnostic. For N=2,3,4,5,8 the manuscript finds that the dimensionless string-length scaling parameter ζ approaches a plateau rather than showing frustration, that non-minimal charge classes remain subdominant, and that the plateau normalization is approximately proportional to N^2−1. From this it infers an empirical scaling law for the long-string energy fraction and hence Ω_GW ∝ μ^2 (N^2−1)^2.

Significance. If the numerical result holds, it addresses a genuinely open dynamical question: Z_N string networks with charge-conserving junctions might frustrate rather than scale, and the simulations indicate scaling with a nontrivial N-dependence. The paper is careful to present the model as an effective global-scalar description rather than as a direct simulation of Yang-Mills flux tubes, and the s=0/s=1 comparison and the isolated-string profile calculation are useful elements. The N=3 ensemble and the fact that the N^2−1 normalization is an empirical fit rather than an input are positive features. However, the central observable ζ is obtained from the topological reconstruction in Appendix B, and that reconstruction is not specified or validated at the level required to support the headline claim. Because the scaling plateau and the N^2−1 law rest entirely on that observable, the paper needs additional numerical validation before the main conclusion can be considered established.

major comments (3)
  1. [Appendix B (Eqs. B1–B5), Table I, Sec. IV.D] The scaling observable ζ is computed from plaquette windings Q_p in Eq. (B5). This construction requires choosing an SU(N) representative U(x) satisfying Φ_i(x) ≃ U(x) Φ_i^(vac) U(x)^† for all three adjoint fields simultaneously 'away from the string core'; no fitting criterion, tolerance, or exclusion threshold is given. Eq. (B2) is overdetermined for generic field configurations, and Table I shows that the radial norm at the core is not small (φ_r^2(0)/v^2 = 0.40–0.83), so a radial-norm cut cannot separate core sites from vacuum-like sites. If q_ab in Eq. (B4) is evaluated on links touching or crossing cores, Q_p is not a protected topological quantity. Since ζ, R_q, Fig. 6, Table II, and the N^2−1 normalization in Eq. (58) all derive from Q_p, the headline result could be an artifact of the reconstruction. Please specify the numerical projection used to find U(x), define 'away from th
  2. [Table II, Fig. 6, Eq. (52)] The N-dependence claim is supported statistically only for N=3, which has a 10-run ensemble; N=2,4,5,8 are single production runs. The time series in Fig. 6 are shown only for N=2,3,4, and for the physical s=1 runs the curves are still slowly rising at η=70, so the quoted 'plateau' values include an extrapolation. Comparing single-run values at a fixed η can bias the N^2−1 pattern if the approach to scaling has different speeds for different N. Please provide multiple seeds for at least N=4 and N=5, report the N=5 and N=8 time series, and give a convergence or extrapolation estimate for the asymptotic ζ̃ rather than the η=70 snapshot.
  3. [Sec. IV.D, Table II (N=8)] The N=8 run uses 512^3 lattice points while the other main runs use 1024^3, so the core resolution is coarser by a factor of two in each direction. The N=8 non-minimal charge fractions are at the few-percent level (R_{q=2}=5%, R_{q=3}=2%, R_{q=4}=2%), precisely the regime where reconstruction noise or finite-resolution effects could dominate. A single coarse run is insufficient to establish that non-minimal charges remain subdominant for N=8. Please quantify the core-width resolution for N=8 and, if feasible, verify with a 1024^3 run or with a different fat-string parameter s.
minor comments (5)
  1. [Figs. 3–5] The figure labels use ω(x)>ω_cut and ω_r^2(x)<c_r v^2 where the text defines ρ(x)>ρ_cut and φ_r^2(x)<c_r v^2. The notation should be made consistent.
  2. [References] Reference [12] is a duplicate of reference [8]; one of the two EPTA/INPTA entries should be removed or merged.
  3. [Eq. (52)] The definition of ζ_q would be clearer if written as ζ_q = t_phys^2 (a ℓ_com,q)/(a^3 L^3), explicitly displaying the physical length a ℓ_com,q and physical volume a^3 L^3.
  4. [Sec. IV.B] The initial-condition prescription specifies a Gaussian variance for the field values but does not state the initial field derivatives. Please state whether the fields are initially at rest or given a thermal-like momentum distribution.
  5. [Table II] The table lists uncertainties only for the N=3, s=0 row. For single-run entries, please mark them explicitly as single-run and, if possible, give a rough run-to-run scatter from a small number of additional seeds for at least one value of N.

Circularity Check

0 steps flagged

No circularity: the scaling and GW conclusions are numerical measurements and an algebraic consequence, not definitions or fitted inputs.

full rationale

The paper's central claims are based on direct lattice measurements. The scaling parameter ζ is reconstructed from plaquette windings (Appendix B) and the N²−1 normalization is an empirical fit (Eq. (58) and Table II); the gravitational-wave scaling is then an algebraic consequence of Eq. (61) using the measured ζ, with no parameter fitted to the GW amplitude. Self-citations to Yamada & Yonekura (Refs. [13–15]) only motivate the effective model and do not enter the simulation, scaling analysis, or GW derivation. The Appendix B reconstruction, while relying on a local PSU(N) representative 'away from the string core' without an explicit exclusion threshold, is a numerical diagnostic; any possible bias near cores is a correctness/robustness concern, not a logical reduction of the result to its inputs. No step was found in which a prediction is equivalent by construction to an input.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 1 invented entities

The central scaling result rests on a specific choice of scalar potential parameters (m,g,λ_r), an initial conformal time η_0, and an instantaneous-quench initial condition. Only one point in this parameter space is simulated, so the N²−1 law is an empirical observation at that point, not a demonstrated universal property. The chain from the toy model to real Yang-Mills flux tubes rests on two explicit ad-hoc extrapolations: the vacuum-manifold structure of N=1* and the neglect of gauge fields. The winding reconstruction itself assumes a well-defined PSU(N) orientation away from cores, which is not automatically satisfied given the nonzero core values in Table I.

free parameters (5)
  • m (mass scale) = 1
    Sets the overall mass scale and core width δ~1/m_core; the dimensionless scaling parameter ζ and its N-dependence are measured at this point. Scanning over m (or g) is not performed, so the N²−1 result may be specific to this parameter slice.
  • g (three-scalar coupling) = 1
    Together with m sets the Higgs vacuum v² = (m/g)² N(N²−1)/2 and the mass spectrum. The ratio of string tension to core width depends on this choice.
  • λ_r v² (radial stabilizing term) = 2 (m_core = m)
    Chosen so the tachyonic curvature at the origin equals m². This sets the radial mass M_rad² = 5 m² and the lightest massive mode M_light = 2m. The core structure and junction dynamics may depend on this ratio; no variation is reported.
  • η_0 (initial conformal time) = 10
    Sets the initial horizon-to-core ratio. The scaling normalization depends logarithmically on this choice (Refs 48,51,54), and the paper quotes ζ at η=70 for this η_0; a different η_0 would shift ζ and hence the GW normalization.
  • initial fluctuation amplitude = 0.1
    Initial condition for the quench; paper says late-time scaling is insensitive, but no systematic scan is documented.
axioms (5)
  • standard math π_1(PSU(N)) = Z_N
    Used to classify string charges q ∈ Z_N (Eq. 21, Sec. II.B).
  • domain assumption The vacuum manifold is PSU(N) after SSB
    Schur's lemma implies trivial stabilizer in PSU(N) for the Higgs vacuum (Eq. 12). This is a property of the model, not of the real YM theory.
  • ad hoc to paper The scalar F-term potential of three adjoint fields captures the confining-string junction structure of N=1*/pure YM
    The paper explicitly states the model 'is not intended to reproduce all microscopic properties' (Sec. I) but is a tractable effective description. The central claim about scaling of real YM strings depends on this.
  • ad hoc to paper Global strings (no gauge fields) have the same scaling behavior as the local flux tubes
    Sec. V: 'Gauge fields would screen the angular gradient energy...' and direct treatment is left to future work. The GW implication assumes this extrapolation.
  • standard math Clebsch-Gordan decomposition of the adjoint of SU(N) under principal SU(2)
    Used in Appendix A to compute the mass spectrum at the Higgs vacuum.
invented entities (1)
  • Effective PSU(N) scalar model with three adjoint fields (Eqs. 6–15) no independent evidence
    purpose: Tractable effective description of Z_N string networks with baryon-vertex-like junctions, motivated by N=1* vacua
    The model is constructed for this paper; it is not derived from pure YM and makes no new falsifiable prediction outside the paper. The paper states it is not intended to reproduce all microscopic properties of pure YM.

pith-pipeline@v1.3.0-alltime-deepseek · 23117 in / 18258 out tokens · 168321 ms · 2026-08-01T06:56:03.920649+00:00 · methodology

0 comments
read the original abstract

We numerically investigate networks of global $\mathbb{Z}_N$ strings with multi-string junctions in a scalar field model whose vacuum manifold is $\mathrm{PSU}(N)=\mathrm{SU}(N)/\mathbb{Z}_N$. The construction of the model is motivated by the Higgs vacua of mass-deformed $\mathcal{N}=4$ supersymmetric Yang-Mills theory, commonly known as $\mathcal{N}=1^*$ theory, and provides a tractable effective description of string networks with baryon-vertex-like junctions. We perform classical lattice simulations of the formation and evolution of these networks in a radiation-dominated universe. For $N=2,3,4,5,$ and $8$, we find that the networks approach a scaling regime rather than becoming frustrated. The normalization of the string density grows in proportion to the dimension of the adjoint representation, $N^2-1$, while non-minimal-charge components remain subdominant. These results imply that, at least for $N\lesssim 8$, the amplitude of the gravitational-wave energy density generated by the cosmic-string network scales as $\Omega_{\rm GW}\propto \mu^2 (N^2-1)^2$, where $\mu$ is the tension of a unit-charge string.

Figures

Figures reproduced from arXiv: 2607.21701 by Masaki Yamada.

Figure 1
Figure 1. Figure 1: for a schematic illustration). For N = 3, this is closely analogous to the familiar baryonic configuration in QCD, where three quarks are connected by color flux tubes meet￾ing at a Y-shaped junction. We are interested in the case without quarks, but the junction structure can still form. In the cosmological context, the existence of such junctions implies that the string network is not simply a collection… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Numerical isolated-string radial profiles for the field-space [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Snapshot of an [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Scaling parameters as functions of conformal time. Top: [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Schematic of the plaquette-winding diagnostic used in the [PITH_FULL_IMAGE:figures/full_fig_p011_7.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

57 extracted references · 41 linked inside Pith

  1. [1]

    SmYTpbSPuWYoDvdQYf1WF+HZOSE=

    Although the topological charge group is the finite group ZN , the amount of field-space structure that can relax and radiate grows rapidly withN. To factor out this empirical <latexit sha1_base64="SmYTpbSPuWYoDvdQYf1WF+HZOSE=">AAAKTXicfVbNjts2EFbSn8Rum27aYy9GjaQbIFjYRdrk0ARBd4M6i12vu8huAliGQdEjmTBFaUkqsUuwb9Cn6bV9hp77IL0VRUlbsiVKWR3s4Xwz35Azo6GClBIhe72/b...

  2. [2]

    Vilenkin and E

    A. Vilenkin and E. P . S. Shellard,Cosmic Strings and Other Topological Defects. Cambridge University Press, 7, 2000

  3. [3]

    T . W . B. Kibble,Topology of Cosmic Domains and Strings,J. Phys. A9(1976) 1387

  4. [4]

    Allen and E

    B. Allen and E. P . S. Shellard,Cosmic string evolution: a numerical simulation,Phys. Rev. Lett.64(1990) 119

  5. [5]

    D. P . Bennett and F . R. Bouchet,Evidence for a Scaling Solution in Cosmic String Evolution,Phys. Rev. Lett.60(1988) 257

  6. [6]

    Vachaspati and A

    T . Vachaspati and A. Vilenkin,Gravitational Radiation from Cosmic Strings,Phys. Rev. D31(1985) 3052

  7. [7]

    Vilenkin,Gravitational radiation from cosmic strings,Phys

    A. Vilenkin,Gravitational radiation from cosmic strings,Phys. Lett. B107(1981) 47

  8. [9]

    Agazie et al.,The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background,Astrophys

    NANOGRAVcollaboration, G. Agazie et al.,The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background,Astrophys. J. Lett.951(2023) L8 [2306.16213]

  9. [10]

    Xu et al.,Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I,Res

    H. Xu et al.,Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I,Res. Astron. Astrophys.23(2023) 075024 [2306.16216]

  10. [11]

    D. J. Reardon et al.,Search for an Isotropic Gravitational-wave Background with the Parkes Pulsar Timing Array,Astrophys. J. Lett.951(2023) L6 [2306.16215]

  11. [12]

    Antoniadis et al.,The second data release from the European Pulsar Timing Array - III

    EPTA, INPTA: collaboration, J. Antoniadis et al.,The second data release from the European Pulsar Timing Array - III. 12 Search for gravitational wave signals,Astron. Astrophys.678 (2023) A50 [2306.16214]

  12. [13]

    Afzal et al.,The NANOGrav 15 yr Data Set: Search for Signals from New Physics,Astrophys

    NANOGRAVcollaboration, A. Afzal et al.,The NANOGrav 15 yr Data Set: Search for Signals from New Physics,Astrophys. J. Lett.951(2023) L11 [2306.16219]

  13. [14]

    Yamada and K

    M. Yamada and K. Yonekura,Cosmic F- and D-strings from pure Yang–Mills theory,Phys. Lett. B838(2023) 137724 [2204.13125]

  14. [15]

    Yamada and K

    M. Yamada and K. Yonekura,Cosmic strings from pure Yang–Mills theory,Phys. Rev. D106(2022) 123515 [2204.13123]

  15. [16]

    Witten,Cosmic Superstrings,Phys

    E. Witten,Cosmic Superstrings,Phys. Lett. B153(1985) 243

  16. [17]

    Yamada and K

    M. Yamada and K. Yonekura,Dark baryon from pure Yang-Mills theory and its GW signature from cosmic strings, JHEP09(2023) 197 [2307.06586]

  17. [18]

    Dvali and A

    G. Dvali and A. Vilenkin,Formation and evolution of cosmic D strings,JCAP03(2004) 010 [hep-th/0312007]

  18. [19]

    Polchinski,Collision of Macroscopic Fundamental Strings, Phys

    J. Polchinski,Collision of Macroscopic Fundamental Strings, Phys. Lett. B209(1988) 252

  19. [20]

    M. G. Jackson, N. T . Jones and J. Polchinski,Collisions of cosmic F and D-strings,JHEP10(2005) 013 [hep-th/0405229]

  20. [21]

    E. J. Copeland, R. C. Myers and J. Polchinski,Cosmic F and D strings,JHEP06(2004) 013 [hep-th/0312067]

  21. [22]

    Ellis, M

    J. Ellis, M. Lewicki, C. Lin and V . Vaskonen,Cosmic superstrings revisited in light of NANOGrav 15-year data, Phys. Rev. D108(2023) 103511 [2306.17147]

  22. [23]

    Hanany and K

    A. Hanany and K. Hashimoto,Reconnection of colliding cosmic strings,JHEP06(2005) 021 [hep-th/0501031]

  23. [24]

    Hanany, M

    A. Hanany, M. J. Strassler and A. Zaffaroni,Confinement and strings in MQCD,Nucl. Phys. B513(1998) 87 [hep-th/9707244]

  24. [25]

    M. R. Douglas and S. H. Shenker,Dynamics of SU(N) supersymmetric gauge theory,Nucl. Phys. B447(1995) 271 [hep-th/9503163]

  25. [26]

    Polchinski and M

    J. Polchinski and M. J. Strassler,The String dual of a confining four-dimensional gauge theory,hep-th/0003136

  26. [27]

    Witten,Anti-de Sitter space, thermal phase transition, and confinement in gauge theories,Adv

    E. Witten,Anti-de Sitter space, thermal phase transition, and confinement in gauge theories,Adv. Theor . Math. Phys.2 (1998) 505 [hep-th/9803131]

  27. [28]

    J. M. Maldacena and C. Nunez,Towards the large N limit of pure N=1 superYang-Mills,Phys. Rev. Lett.86(2001) 588 [hep-th/0008001]

  28. [29]

    I. R. Klebanov and M. J. Strassler,Supergravity and a confining gauge theory: Duality cascades and chi SB resolution of naked singularities,JHEP08(2000) 052 [hep-th/0007191]

  29. [30]

    Witten,Baryons and branes in anti-de Sitter space,JHEP07 (1998) 006 [hep-th/9805112]

    E. Witten,Baryons and branes in anti-de Sitter space,JHEP07 (1998) 006 [hep-th/9805112]

  30. [31]

    Vafa,Superstrings and topological strings at large N,J

    C. Vafa,Superstrings and topological strings at large N,J. Math. Phys.42(2001) 2798 [hep-th/0008142]

  31. [32]

    McGraw,Evolution of a nonAbelian cosmic string network, Phys

    P . McGraw,Evolution of a nonAbelian cosmic string network, Phys. Rev. D57(1998) 3317 [astro-ph/9706182]

  32. [33]

    Spergel and U.-L

    D. Spergel and U.-L. Pen,Cosmology in a string dominated universe,Astrophys. J. Lett.491(1997) L67 [astro-ph/9611198]

  33. [34]

    Vachaspati and A

    T . Vachaspati and A. Vilenkin,Evolution of cosmic networks, Phys. Rev. D35(1987) 1131

  34. [35]

    Avgoustidis and E

    A. Avgoustidis and E. P . S. Shellard,Velocity-Dependent Models for Non-Abelian/Entangled String Networks,Phys. Rev. D78(2008) 103510 [0705.3395]

  35. [36]

    E. J. Copeland and P . M. Saffin,On the evolution of cosmic-superstring networks,JHEP11(2005) 023 [hep-th/0505110]

  36. [37]

    Y. Ng, T . W . B. Kibble and T . Vachaspati,Formation of Non-Abelian Monopoles Connected by Strings,Phys. Rev. D78 (2008) 046001 [0806.0155]

  37. [38]

    Urrestilla and A

    J. Urrestilla and A. Vilenkin,Evolution of cosmic superstring networks: A Numerical simulation,JHEP02(2008) 037 [0712.1146]

  38. [39]

    Hindmarsh and P

    M. Hindmarsh and P . M. Saffin,Scaling in a SU(2)/Z3 model of cosmic superstring networks,JHEP08(2006) 066 [hep-th/0605014]

  39. [40]

    ’t Hooft,On the Phase Transition Towards Permanent Quark Confinement,Nucl

    G. ’t Hooft,On the Phase Transition Towards Permanent Quark Confinement,Nucl. Phys. B138(1978) 1

  40. [41]

    Mandelstam,Vortices and Quark Confinement in Nonabelian Gauge Theories,Phys

    S. Mandelstam,Vortices and Quark Confinement in Nonabelian Gauge Theories,Phys. Rept.23(1976) 245

  41. [42]

    Seiberg and E

    N. Seiberg and E. Witten,Monopoles, duality and chiral symmetry breaking in N=2 supersymmetric QCD,Nucl. Phys. B431(1994) 484 [hep-th/9408099]

  42. [43]

    Seiberg and E

    N. Seiberg and E. Witten,Electric - magnetic duality, monopole condensation, and confinement in N=2 supersymmetric Yang-Mills theory,Nucl. Phys. B426(1994) 19 [hep-th/9407087]

  43. [44]

    Dorey,An Elliptic superpotential for softly broken N=4 supersymmetric Yang-Mills theory,JHEP07(1999) 021 [hep-th/9906011]

    N. Dorey,An Elliptic superpotential for softly broken N=4 supersymmetric Yang-Mills theory,JHEP07(1999) 021 [hep-th/9906011]

  44. [45]

    Donagi and E

    R. Donagi and E. Witten,Supersymmetric Yang-Mills theory and integrable systems,Nucl. Phys. B460(1996) 299 [hep-th/9510101]

  45. [46]

    Vincent, N

    G. Vincent, N. D. Antunes and M. Hindmarsh,Numerical simulations of string networks in the Abelian Higgs model, Phys. Rev. Lett.80(1998) 2277 [hep-ph/9708427]

  46. [47]

    S. G. Naculich, H. J. Schnitzer and N. Wyllard,Vacuum states of N=1* mass deformations of N=4 and N=2 conformal gauge theories and their brane interpretations,Nucl. Phys. B609 (2001) 283 [hep-th/0103047]

  47. [48]

    Yamaguchi,Scaling property of the global string in the radiation dominated universe,Phys

    M. Yamaguchi,Scaling property of the global string in the radiation dominated universe,Phys. Rev. D60(1999) 103511 [hep-ph/9907506]

  48. [49]

    Yamaguchi, J

    M. Yamaguchi, J. Yokoyama and M. Kawasaki,Numerical analysis of formation and evolution of global strings in ( 2+1)-dimensions,Prog. Theor . Phys.100(1998) 535 [hep-ph/9808326]

  49. [50]

    W . H. Press, B. S. Ryden and D. N. Spergel,Dynamical Evolution of Domain Walls in an Expanding Universe, Astrophys. J.347(1989) 590

  50. [51]

    Aryal, L

    M. Aryal, L. H. Ford and A. Vilenkin,Cosmic Strings and Black Holes,Phys. Rev. D34(1986) 2263

  51. [52]

    Hindmarsh, J

    M. Hindmarsh, J. Lizarraga, J. Urrestilla, D. Daverio and M. Kunz,Scaling from gauge and scalar radiation in Abelian Higgs string networks,Phys. Rev. D96(2017) 023525 [1703.06696]

  52. [53]

    Fleury and G

    L. Fleury and G. D. Moore,Axion dark matter: strings and their cores,JCAP01(2016) 004 [1509.00026]

  53. [54]

    Gorghetto, E

    M. Gorghetto, E. Hardy and G. Villadoro,Axions from Strings: the Attractive Solution,JHEP07(2018) 151 [1806.04677]

  54. [55]

    Hindmarsh, J

    M. Hindmarsh, J. Lizarraga, A. Lopez-Eiguren and J. Urrestilla,Scaling Density of Axion Strings,Phys. Rev. Lett. 124(2020) 021301 [1908.03522]

  55. [56]

    Strobl,Improvements for Vachaspati-Vilenkin type algorithms for cosmic string and disclination formation, hep-lat/9608085

    K. Strobl,Improvements for Vachaspati-Vilenkin type algorithms for cosmic string and disclination formation, hep-lat/9608085

  56. [57]

    Vachaspati and A

    T . Vachaspati and A. Vilenkin,Formation and Evolution of Cosmic Strings,Phys. Rev. D30(1984) 2036

  57. [70]

    Cf/9YpptQMwY3KV5d0BcVoPrs1s=

    ForN=3, the chargeq=2 is the conjugate ofq=1 and is not listed as an independent unoriented type. ForN=8, the higher-charge classes are individually at the few-percent 8 <latexit sha1_base64="Cf/9YpptQMwY3KV5d0BcVoPrs1s=">AAAKTXicfVZfb9s2EFe7P629rU23x70YM9qlQBHYQ7DtYS2KJcXcIHG8oEkLWIZB0SeZMEUpJNXaI7hvsE+z1+0z7HkfZG/DMFKWbIlyowf7eL+735F3p6OClBIhe72/b93+4MO...