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 →
Formation and scaling of mathbb{Z}_N strings for global SU(N)/mathbb{Z}_N symmetry
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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
- 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.
Referee Report
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)
- [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
- [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.
- [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)
- [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.
- [References] Reference [12] is a duplicate of reference [8]; one of the two EPTA/INPTA entries should be removed or merged.
- [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.
- [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.
- [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
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
free parameters (5)
- m (mass scale) =
1
- g (three-scalar coupling) =
1
- λ_r v² (radial stabilizing term) =
2 (m_core = m)
- η_0 (initial conformal time) =
10
- initial fluctuation amplitude =
0.1
axioms (5)
- standard math π_1(PSU(N)) = Z_N
- domain assumption The vacuum manifold is PSU(N) after SSB
- ad hoc to paper The scalar F-term potential of three adjoint fields captures the confining-string junction structure of N=1*/pure YM
- ad hoc to paper Global strings (no gauge fields) have the same scaling behavior as the local flux tubes
- standard math Clebsch-Gordan decomposition of the adjoint of SU(N) under principal SU(2)
invented entities (1)
-
Effective PSU(N) scalar model with three adjoint fields (Eqs. 6–15)
no independent evidence
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
Reference graph
Works this paper leans on
-
[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]
Vilenkin and E
A. Vilenkin and E. P . S. Shellard,Cosmic Strings and Other Topological Defects. Cambridge University Press, 7, 2000
2000
-
[3]
T . W . B. Kibble,Topology of Cosmic Domains and Strings,J. Phys. A9(1976) 1387
1976
-
[4]
Allen and E
B. Allen and E. P . S. Shellard,Cosmic string evolution: a numerical simulation,Phys. Rev. Lett.64(1990) 119
1990
-
[5]
D. P . Bennett and F . R. Bouchet,Evidence for a Scaling Solution in Cosmic String Evolution,Phys. Rev. Lett.60(1988) 257
1988
-
[6]
Vachaspati and A
T . Vachaspati and A. Vilenkin,Gravitational Radiation from Cosmic Strings,Phys. Rev. D31(1985) 3052
1985
-
[7]
Vilenkin,Gravitational radiation from cosmic strings,Phys
A. Vilenkin,Gravitational radiation from cosmic strings,Phys. Lett. B107(1981) 47
1981
-
[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]
Pith/arXiv arXiv 2023
-
[10]
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]
Pith/arXiv arXiv 2023
-
[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]
Pith/arXiv arXiv 2023
-
[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]
Pith/arXiv arXiv 2023
-
[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]
Pith/arXiv arXiv 2023
-
[14]
M. Yamada and K. Yonekura,Cosmic F- and D-strings from pure Yang–Mills theory,Phys. Lett. B838(2023) 137724 [2204.13125]
Pith/arXiv arXiv 2023
-
[15]
M. Yamada and K. Yonekura,Cosmic strings from pure Yang–Mills theory,Phys. Rev. D106(2022) 123515 [2204.13123]
Pith/arXiv arXiv 2022
-
[16]
Witten,Cosmic Superstrings,Phys
E. Witten,Cosmic Superstrings,Phys. Lett. B153(1985) 243
1985
-
[17]
M. Yamada and K. Yonekura,Dark baryon from pure Yang-Mills theory and its GW signature from cosmic strings, JHEP09(2023) 197 [2307.06586]
Pith/arXiv arXiv 2023
-
[18]
G. Dvali and A. Vilenkin,Formation and evolution of cosmic D strings,JCAP03(2004) 010 [hep-th/0312007]
Pith/arXiv arXiv 2004
-
[19]
Polchinski,Collision of Macroscopic Fundamental Strings, Phys
J. Polchinski,Collision of Macroscopic Fundamental Strings, Phys. Lett. B209(1988) 252
1988
-
[20]
M. G. Jackson, N. T . Jones and J. Polchinski,Collisions of cosmic F and D-strings,JHEP10(2005) 013 [hep-th/0405229]
Pith/arXiv arXiv 2005
-
[21]
E. J. Copeland, R. C. Myers and J. Polchinski,Cosmic F and D strings,JHEP06(2004) 013 [hep-th/0312067]
Pith/arXiv arXiv 2004
-
[22]
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]
Pith/arXiv arXiv 2023
-
[23]
A. Hanany and K. Hashimoto,Reconnection of colliding cosmic strings,JHEP06(2005) 021 [hep-th/0501031]
Pith/arXiv arXiv 2005
-
[24]
A. Hanany, M. J. Strassler and A. Zaffaroni,Confinement and strings in MQCD,Nucl. Phys. B513(1998) 87 [hep-th/9707244]
Pith/arXiv arXiv 1998
-
[25]
M. R. Douglas and S. H. Shenker,Dynamics of SU(N) supersymmetric gauge theory,Nucl. Phys. B447(1995) 271 [hep-th/9503163]
Pith/arXiv arXiv 1995
-
[26]
J. Polchinski and M. J. Strassler,The String dual of a confining four-dimensional gauge theory,hep-th/0003136
-
[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]
Pith/arXiv arXiv 1998
-
[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]
Pith/arXiv arXiv 2001
-
[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]
Pith/arXiv arXiv 2000
-
[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]
Pith/arXiv arXiv 1998
-
[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]
Pith/arXiv arXiv 2001
-
[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]
Pith/arXiv arXiv 1998
-
[33]
D. Spergel and U.-L. Pen,Cosmology in a string dominated universe,Astrophys. J. Lett.491(1997) L67 [astro-ph/9611198]
Pith/arXiv arXiv 1997
-
[34]
Vachaspati and A
T . Vachaspati and A. Vilenkin,Evolution of cosmic networks, Phys. Rev. D35(1987) 1131
1987
-
[35]
A. Avgoustidis and E. P . S. Shellard,Velocity-Dependent Models for Non-Abelian/Entangled String Networks,Phys. Rev. D78(2008) 103510 [0705.3395]
Pith/arXiv arXiv 2008
-
[36]
E. J. Copeland and P . M. Saffin,On the evolution of cosmic-superstring networks,JHEP11(2005) 023 [hep-th/0505110]
Pith/arXiv arXiv 2005
-
[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]
Pith/arXiv arXiv 2008
-
[38]
J. Urrestilla and A. Vilenkin,Evolution of cosmic superstring networks: A Numerical simulation,JHEP02(2008) 037 [0712.1146]
Pith/arXiv arXiv 2008
-
[39]
M. Hindmarsh and P . M. Saffin,Scaling in a SU(2)/Z3 model of cosmic superstring networks,JHEP08(2006) 066 [hep-th/0605014]
Pith/arXiv arXiv 2006
-
[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
1978
-
[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
1976
-
[42]
N. Seiberg and E. Witten,Monopoles, duality and chiral symmetry breaking in N=2 supersymmetric QCD,Nucl. Phys. B431(1994) 484 [hep-th/9408099]
Pith/arXiv arXiv 1994
-
[43]
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]
Pith/arXiv arXiv 1994
-
[44]
N. Dorey,An Elliptic superpotential for softly broken N=4 supersymmetric Yang-Mills theory,JHEP07(1999) 021 [hep-th/9906011]
Pith/arXiv arXiv 1999
-
[45]
R. Donagi and E. Witten,Supersymmetric Yang-Mills theory and integrable systems,Nucl. Phys. B460(1996) 299 [hep-th/9510101]
Pith/arXiv arXiv 1996
-
[46]
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]
Pith/arXiv arXiv 1998
-
[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]
Pith/arXiv arXiv 2001
-
[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]
Pith/arXiv arXiv 1999
-
[49]
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]
Pith/arXiv arXiv 1998
-
[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
1989
-
[51]
Aryal, L
M. Aryal, L. H. Ford and A. Vilenkin,Cosmic Strings and Black Holes,Phys. Rev. D34(1986) 2263
1986
-
[52]
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]
Pith/arXiv arXiv 2017
-
[53]
L. Fleury and G. D. Moore,Axion dark matter: strings and their cores,JCAP01(2016) 004 [1509.00026]
Pith/arXiv arXiv 2016
-
[54]
M. Gorghetto, E. Hardy and G. Villadoro,Axions from Strings: the Attractive Solution,JHEP07(2018) 151 [1806.04677]
Pith/arXiv arXiv 2018
-
[55]
M. Hindmarsh, J. Lizarraga, A. Lopez-Eiguren and J. Urrestilla,Scaling Density of Axion Strings,Phys. Rev. Lett. 124(2020) 021301 [1908.03522]
Pith/arXiv arXiv 2020
-
[56]
K. Strobl,Improvements for Vachaspati-Vilenkin type algorithms for cosmic string and disclination formation, hep-lat/9608085
-
[57]
Vachaspati and A
T . Vachaspati and A. Vilenkin,Formation and Evolution of Cosmic Strings,Phys. Rev. D30(1984) 2036
1984
-
[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...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.