Asymptotic scaling of multi-tension cosmic string networks with bound states: unequal tensions, bound-state abundance, and semi-analytic calibration

Cosmic superstrings are metastable macroscopic fundamental (F) and D-strings, formed in brane-inflation scenarios, which can bind into $(p,q)$ states joined at Y-junctions and must reach a scaling regime to be cosmologically viable and observationally predictable. Urrestilla and Vilenkin (2008) simulated a two-sector field-theory proxy for such networks (the Saffin model) and found that while network-level correlation lengths scale, the comoving lengths of the constituent-string segments showed a prolonged plateau, approaching linear growth only under an extreme extended-damping protocol of questionable physical reliability; unequal string tensions and reconnection probabilities below unity were not simulated. The subsequent literature shows: Nambu-Goto and semi-analytic studies found robust multi-tension scaling, and $4096^3$ field-theory simulations (Correia and Martins 2022) confirmed linear growth of the average bound-segment length from generic initial conditions, resolving the original scaling-status question for equal tensions. The surviving open problem is the quantitative completion of the scaling characterization: field-theory networks with unequal string tensions (the F/D hierarchy), a definitive asymptotic bound-state abundance and its dependence on cosmological epoch, and the calibration of the multi-tension velocity-dependent one-scale model against simulation, as required for quantitative CMB anisotropy and stochastic gravitational-wave-background predictions.

Background

Cosmic strings are line-like topological defects that can form at symmetry-breaking phase transitions in the early universe. In string-theoretic brane-inflation scenarios (inflationary models in which inflation ends when two higher-dimensional membranes collide), macroscopic fundamental strings (F-strings, the elementary strings of string theory) and D-strings (one-dimensional D-branes, extended objects on which F-strings can end) can survive as relic one-dimensional objects; because both species carry conserved charges, they can also form bound states known as $(p,q)$ strings, containing $p$ F-strings and $q$ D-strings, which meet at Y-shaped junctions where three string segments of different type join (Copeland, Myers, and Polchinski 2004). Their tensions form a hierarchy, $\mu_{pq}=\mu_F\sqrt{p^2+q^2/g_s^2}$, so for string coupling $g_s<1$ the D-string is heavier than the F-string, unlike in most field-theory analogues.

A string network only avoids dominating the universe's energy density if it reaches a scaling regime: the typical inter-string distance and other characteristic lengths grow in proportion to the horizon (linearly in conformal time, the time coordinate that rescales the expansion away), so that the network's energy density stays a fixed fraction of the total. Whether networks with junctions and bound states scale, or instead freeze into a frustrated static configuration, was historically in doubt, because junctions can lock segments together. The velocity-dependent one-scale (VOS) model — a semi-analytic framework evolving one correlation length and one root-mean-square velocity per string species — was extended to multi-tension networks with binding and unbinding: analytic solutions (Tye, Wasserman, and Wyman 2005; Pourtsidou et al. 2011) predict that such networks approach scaling solutions in which only a few lowest-tension states are substantially populated.

Direct field-theory simulations resolve the string microphysics but have limited dynamic range. The standard proxy is the Saffin model: two Abelian-Higgs sectors — each a complex scalar field coupled to a $U(1)$ gauge field, whose symmetry breaking produces gauge strings — with combined gauge group $U(1)_p\times U(1)_q$, coupled by an attractive interaction term, which supports stable bound strings carrying flux of both sectors — a field-theory stand-in for $(p,q)$ superstrings. Lattice simulations of it must also choose a width-evolution prescription, most commonly the Press–Ryden–Spergel (PRS) algorithm, which artificially keeps the comoving string width constant so that a single run spans more physical time; whether this distorts bound-state formation is itself a question.

Urrestilla and Vilenkin (2008) performed the first large simulations of this model (flat, radiation-era, and matter-era expansion; $512^3$ lattices). They found that an interconnected network does form and does not freeze: the network-level correlation length scales, and the average bound-segment length and junction spacing grow linearly. However, the average comoving lengths of the constituent A and B string segments stayed nearly constant — a 'transient' regime — and only an extreme protocol with a second, strongly damped evolution period (pushed to the edge of code stability) made these lengths approximately linear, a result the authors flagged as tentative because of likely spurious damping. Bound strings carried only about 2% of the total string length, they assumed equal constituent tensions and unit reconnection probability (crossing strings always reconnect rather than pass through), and they listed as open problems the partition of energy loss among loop production, direct particle radiation, and massive-field emission, the microphysical conditions for bound-segment formation in string collisions, and the matching to Nambu-Goto models (the idealization of strings as infinitesimally thin objects whose motion minimizes worldsheet area, ignoring their field-theoretic internal structure) and to analytic models, needed for observational predictions.

Subsequent work transformed this picture. Nambu-Goto-type simulations of F/D networks with controlled initial bound-state abundance found robust scaling of all three components (Sakellariadou and Stoica 2008). Kinematic constraints on junction formation were shown to extend the parameter range admitting scaling in the analytic treatment (Avgoustidis and Copeland 2010), and the full semi-analytic scaling configurations of F-D networks were worked out with implications for the cosmic microwave background (CMB) and for pulsar timing (Pourtsidou et al. 2011). In field theory, Lizarraga and Urrestilla (2016) showed, from two opposite initial conditions (with and without an initial bound-string population), that both cases scale with a low bound-string fraction and efficient 'unzipping' of bound segments (bound pairs separating into their constituents), with hints of a common asymptote that their dynamic range could not confirm. The averaged junction dynamics was incorporated into the multi-tension VOS model (Rybak, Avgoustidis, and Martins 2019), which awaits simulation calibration. On the energy-loss side, field-theory simulations of single-species Abelian-Higgs loops found that particle radiation — emitted mainly at collisions between kinks, the sharp corners a string acquires when segments reconnect — makes the loop decay time grow as $L^2$ in loop length $L$, dominating over gravitational radiation for small loops (Matsunami et al. 2019). Finally, Correia and Martins (2022) simulated the Saffin model on $4096^3$ lattices with $\Delta x=0.5$, under both PRS and physical evolution: they obtained robust scaling for the lightest strings, confirmed the linear growth of the average bound-segment length from generic (random-phase) initial conditions — resolving the tentative finding of the 2008 paper — and found tentative evidence for an asymptotically constant bound-state fraction of order 1–2%, possibly different in the radiation and matter eras, about which 'no firm conclusions can be drawn at present'. They simulated only equal constituent tensions, and explicitly identified the unequal-tension case (for example, one tension set to half the other, mimicking the F/D hierarchy), the definitive determination of the bound-state abundance, and the calibration of the multi-tension semi-analytic model as the outstanding tasks. Adjacent defect classes show a consistent picture: global $\mathbb{Z}_N$ string networks with multi-string junctions approach scaling rather than frustration for $N$ up to 8 (Yamada 2026).

Problem Statement

Complete the quantitative characterization of the asymptotic scaling regime of multi-tension string networks containing bound states and Y-junctions, in the regime relevant to cosmic superstrings. The problem has three linked components serving this single objective:

(1) Unequal tensions. In high-resolution field-theory simulations of two-sector (Saffin-type) models in which the two constituent string tensions are unequal — an F/D-type hierarchy, for instance with one tension set to twice the other — establish whether every string species, including the bound states, reaches a scaling regime, and measure the asymptotic slopes and values of all characteristic lengths (the correlation lengths of each species, the average bound-segment length, and the junction spacing) together with the root-mean-square velocities, in both radiation-dominated and matter-dominated expansion, with resolution convergence and width-evolution (PRS versus physical evolution) systematics under control.

(2) Bound-state abundance. Determine definitively the asymptotic fraction of total string length carried by bound states and its dependence on cosmological epoch; at equal tensions this fraction is known only tentatively (of order 1–2%, with a possible radiation-versus-matter difference that existing data cannot confirm), and its robustness against the threshold choices used to identify sparse bound-state segments on the lattice must be established.

(3) Model calibration. Use these measurements to calibrate the multi-tension velocity-dependent one-scale model — including its junction-formation kinematic constraints and effective binding/unbinding parameters — and test whether its predicted scaling densities and velocities agree with the field-theory results, thereby fixing the inputs required for quantitative predictions of CMB anisotropies and stochastic gravitational-wave backgrounds from superstring-like networks.

The remaining open tasks, which jointly characterize the scaling regime, are: (1) the scaling characterization for at least one unequal-tension case alongside the equal-tension benchmark, with controlled systematics and in both radiation-dominated and matter-dominated expansion; (2) the definitive asymptotic fraction of total string length carried by bound states and its dependence on cosmological epoch; and (3) validation or falsification of the multi-tension semi-analytic model against those results. The original equal-tension scaling-status question posed by the source paper has been affirmatively resolved by later work and is not part of this problem.

References

1. J. Urrestilla, A. Vilenkin, 'Evolution of cosmic superstring networks: a numerical simulation', JHEP 02 (2008) 037, arXiv:0712.1146, doi:10.1088/1126-6708/2008/02/037, https://arxiv.org/abs/0712.1146
2. E. J. Copeland, R. C. Myers, J. Polchinski, 'Cosmic F- and D-strings', JHEP 06 (2004) 013, arXiv:hep-th/0312067, doi:10.1088/1126-6708/2004/06/013, https://arxiv.org/abs/hep-th/0312067
3. P. M. Saffin, 'A practical model for cosmic (p,q) superstrings', JHEP 09 (2005) 011, arXiv:hep-th/0506138, doi:10.1088/1126-6708/2005/09/011, https://arxiv.org/abs/hep-th/0506138
4. S.-H. H. Tye, I. Wasserman, M. Wyman, 'Scaling of Multi-Tension Cosmic Superstring Networks', Phys. Rev. D 71 (2005) 103508, arXiv:astro-ph/0503506, doi:10.1103/PhysRevD.71.103508, https://arxiv.org/abs/astro-ph/0503506
5. M. Sakellariadou, H. Stoica, 'Dynamics of F/D networks: the role of bound states', JCAP 08 (2008) 038, arXiv:0806.3219, doi:10.1088/1475-7516/2008/08/038, https://arxiv.org/abs/0806.3219
6. A. Avgoustidis, E. J. Copeland, 'The effect of kinematic constraints on multi-tension string network evolution', Phys. Rev. D 81 (2010) 063517, arXiv:0912.4004, doi:10.1103/PhysRevD.81.063517, https://arxiv.org/abs/0912.4004
7. A. Pourtsidou, A. Avgoustidis, E. J. Copeland, L. Pogosian, D. A. Steer, 'Scaling configurations of cosmic superstring networks and their cosmological implications', Phys. Rev. D 83 (2011) 063525, arXiv:1012.5014, doi:10.1103/PhysRevD.83.063525, https://arxiv.org/abs/1012.5014
8. J. Lizarraga, J. Urrestilla, 'Survival of pq-superstrings in field theory simulations', JCAP 04 (2016) 053, arXiv:1602.08014, doi:10.1088/1475-7516/2016/04/053, https://arxiv.org/abs/1602.08014
9. I. Yu. Rybak, A. Avgoustidis, C. J. A. P. Martins, 'Dynamics of junctions and the multi-tension velocity-dependent one-scale model', Phys. Rev. D 99 (2019) 063516, arXiv:1812.04584, doi:10.1103/PhysRevD.99.063516, https://arxiv.org/abs/1812.04584
10. D. Matsunami, L. Pogosian, A. Saurabh, T. Vachaspati, 'Decay of Cosmic String Loops Due to Particle Radiation', Phys. Rev. Lett. 122 (2019) 201301, arXiv:1903.05102, doi:10.1103/PhysRevLett.122.201301, https://arxiv.org/abs/1903.05102
11. J. R. C. C. C. Correia, C. J. A. P. Martins, 'Multitension strings in high-resolution $U(1)\times U(1)$ simulations', Phys. Rev. D 106 (2022) 043521, arXiv:2208.01525, doi:10.1103/PhysRevD.106.043521, https://arxiv.org/abs/2208.01525
12. M. Yamada, 'Formation and scaling of $\mathbb{Z}_N$ strings for global $\mathrm{SU}(N)/\mathbb{Z}_N$ symmetry', arXiv:2607.21701 (2026), https://arxiv.org/abs/2607.21701