Hao, Nepomechie and Sommese (Phys. Rev. E 88, 052113 (2013)) used numerical homotopy continuation to enumerate every solution of the spin-1/2 periodic XXX Bethe equations for chains up to N=14N=14 sites, and conjectured that, once unphysical singular solutions are discarded, the number of remaining distinct-root solutions in the MM-down-spin sector equals the number of highest-weight states, (NM)−(NM−1)\binom{N}{M}-\binom{N}{M-1}. Later literature has essentially proven this combination identity through the polynomiality/Q-system framework of Mukhin-Tarasov-Varchenko, Granet-Jacobsen, Bajnok-Granet-Jacobsen-Nepomechie, and Hou-Jiang-Zhu, but only under an unproven assumption that Bethe roots are pairwise distinct. The finer questions the source explicitly posed remain open: closed-form expressions with proofs for the individual counting functions N(N,M)\mathcal{N}(N,M), Ns(N,M)\mathcal{N}_s(N,M), Nsp(N,M)\mathcal{N}_{sp}(N,M) for M≥3M\ge 3 (with the related Kirillov-Sakamoto rigged-configuration conjectures, formulated later in 2014, still unproven, and the odd-length periodicity patterns of Nsp\mathcal{N}_{sp} unexplained); a rigorous proof or refutation of the Pauli principle, i.e. the absence of physical repeated-root ("strange") solutions for spin 1/2; and the classification and thermodynamic-limit abundance of unphysical singular solutions, together with a precise account of why string-hypothesis-based thermodynamic computations nevertheless produce correct highest-weight counts.

Background

The isotropic Heisenberg (XXX) spin-1/21/2 chain on a ring of NN sites is the archetypal exactly solved quantum model. Its Hamiltonian H=∑n=1N(σ⃗n⋅σ⃗n+1−1)/4H=\sum_{n=1}^{N}\left(\vec\sigma_n\cdot\vec\sigma_{n+1}-1\right)/4 (periodic boundary conditions; σ⃗\vec\sigma are the Pauli matrices) conserves the total spin, so the Hilbert space splits into sectors with MM down spins. In the coordinate Bethe ansatz, eigenstates in the MM-down sector are parametrized by MM complex rapidities λ1,…,λM\lambda_1,\dots,\lambda_M solving the Bethe equations

(λj+i/2λj−i/2)N=∏k≠jMλj−λk+iλj−λk−i,j=1,…,M.\left(\frac{\lambda_j+i/2}{\lambda_j-i/2}\right)^{N}=\prod_{k\neq j}^{M}\frac{\lambda_j-\lambda_k+i}{\lambda_j-\lambda_k-i},\qquad j=1,\dots,M.

Completeness of the Bethe ansatz — the question whether these equations, used correctly, account for all 2N2^N eigenstates — is a long-standing foundational problem. Since the Hamiltonian commutes with the global sl(2)sl(2) spin algebra, every state descends from a highest-weight state (a state annihilated by the total spin-raising operator), and the number of highest-weight states in the MM-down sector is (NM)−(NM−1)\binom{N}{M}-\binom{N}{M-1}; completeness therefore reduces to counting highest-weight states. Two subtleties have made the problem resistant. First, the Bethe equations admit singular solutions containing the pair λ=±i/2\lambda=\pm i/2 (an "exact two-string"; a string is a group of roots spaced regularly by ii in the imaginary direction), for which the standard Bethe-vector and energy expressions are singular. Nepomechie and Wang (2013) derived a criterion, via a regularized limit of the equations, that decides which singular solutions are physical, i.e. correspond to genuine eigenvectors, and which are unphysical artifacts. Second, solutions may a priori have repeated roots; for higher-spin chains such "strange" solutions exist and can be physical (Hao-Nepomechie-Sommese 2014), but for spin 1/2 none has ever been found, an assertion called the Pauli principle for Bethe roots.

Hao, Nepomechie and Sommese (2013) settled the numerics: using homotopy continuation (a certified numerical method that tracks all isolated solutions of a polynomial system from a start system with known roots, guaranteeing that none is missed), they enumerated the full solution set for N≤14N\le 14. They found that the equations have "too many" distinct-root solutions, but that after discarding unphysical singular solutions the right number remains, leading to their completeness/Pauli-principle conjecture

N(N,M)−Ns(N,M)+Nsp(N,M)=(NM)−(NM−1).\mathcal{N}(N,M)-\mathcal{N}_s(N,M)+\mathcal{N}_{sp}(N,M)=\binom{N}{M}-\binom{N}{M-1}.

where N\mathcal{N} counts solutions with pairwise distinct roots, Ns\mathcal{N}_s the singular ones among them, and Nsp\mathcal{N}_{sp} the singular ones that are physical. They could conjecture only this combination, not the three functions separately; closed forms are known only for M=1M=1 (N=N−1\mathcal{N}=N-1) and M=2M=2 (the formula consistent with their Table 2 for every N≤14N\le 14 is N(N,2)=(N2−3N+1−(−1)N)/2\mathcal{N}(N,2)=(N^2-3N+1-(-1)^N)/2, whereas the printed equation (5.2) of the paper is a typo that contradicts the table). The same enumeration revealed arithmetic structure in the singular sector: for odd NN almost all singular solutions are unphysical, with striking exceptions repeating with period 6 (Nsp(9,3)=Nsp(15,3)=2\mathcal{N}_{sp}(9,3)=\mathcal{N}_{sp}(15,3)=2), and for M∼N/2M\sim N/2 the number of unphysical singular solutions is comparable to the number of highest-weight states — suggesting that the string hypothesis (the assumption, traditional in thermodynamic analyses, that roots organize into strings), which underlies the classical completeness arguments for N→∞N\to\infty, must implicitly project out unphysical solutions for reasons that were never made precise.

Since 2013, the completeness question proper has been largely settled by a chain of rigorous work in a different, algebraic-geometric formulation. The eigenvalue T(λ)T(\lambda) of the transfer matrix (the generating function of commuting conserved quantities, whose logarithmic derivative is the Hamiltonian) obeys the TQTQ relation T(λ)Q(λ)=(λ−i/2)NQ(λ+i)+(λ+i/2)NQ(λ−i)T(\lambda)Q(\lambda)=(\lambda-i/2)^N Q(\lambda+i)+(\lambda+i/2)^N Q(\lambda-i) with Q(λ)=∏k(λ−λk)Q(\lambda)=\prod_k(\lambda-\lambda_k). Mukhin, Tarasov and Varchenko (2009) proved that the Bethe algebra of the homogeneous XXX model has simple spectrum, with the byproduct that there are exactly (nl)−(nl−1)\binom{n}{l}-\binom{n}{l-1} two-dimensional spaces of polynomial solutions of the associated Wronskian difference equation. Granet and Jacobsen (2020) proved that a Bethe solution is physical if and only if two zero-remainder conditions hold (two nested TQTQ-type relations have polynomial right-hand sides), reconciling this framework with the Nepomechie-Wang regularization of singular solutions. Bajnok-Granet-Jacobsen-Nepomechie (2020) recast the physicality conditions as Q-systems — triangular systems of polynomial QQQQ-relations whose solutions directly yield the admissible, i.e. physical, solutions — for the closed XXZ, open XXX, and open quantum-group-invariant XXZ chains, proving the equivalence of admissibility and physicality while leaving the solution counts as conjectures; Hou-Jiang-Zhu (2024) formulated the spin-ss rational Q-system and proved its completeness for the periodic chain. However, for the spin-1/21/2 chain these arguments assume the Bethe roots are pairwise distinct, citing precisely the unproven spin-1/21/2 Pauli principle; Hou-Jiang-Zhu state explicitly that for the spin-1/21/2 XXX chain this has good numerical support but no rigorous proof. Kirillov and Sakamoto (2014) proposed conjectural formulas, in terms of rigged configurations (Young-diagram data with integer labels counting fermionic excitations, from the Kyoto-school theory of solvable lattice models), for the number of solutions with nonzero Bethe vectors and for the number of physical singular solutions at even NN; these conjectures remain unproven, and no closed forms exist at odd NN.

Problem Statement

Establish a complete quantitative theory of the solution set of the spin-1/21/2 periodic XXX Bethe equations, continuing the HNS program beyond the now-essentially-proven combination identity. Three questions, all for the homogeneous periodic spin-1/21/2 chain, constitute the open core, and all three are mandatory: each of them is an independently required deliverable, so a result that settles only one or two of them counts as progress but does not complete the task. Together they support the single objective of an exact arithmetic characterization of the solution set — how many solutions exist, which are physical, and how the unphysical ones behave at large NN.

  1. Closed-form counting. This sub-question comprises three separately required deliverables: (a) find closed-form expressions, with proofs, for the individual counting functions N(N,M)\mathcal{N}(N,M), Ns(N,M)\mathcal{N}_s(N,M), and Nsp(N,M)\mathcal{N}_{sp}(N,M) for general NN and M≥3M\ge 3 (in the source's convention M≤⌊N/2⌋M\le\lfloor N/2\rfloor, the range in which these functions are defined and tabulated), covering both parities of NN; (b) prove or refute the Kirillov-Sakamoto rigged-configuration conjectures (formulated for even NN); (c) explain the parity and periodicity structure of Nsp\mathcal{N}_{sp} at odd NN (in the M=3M=3 sector the existence of singular solutions precisely for odd N=3(2k+1)N=3(2k+1) is proven, and the tables give Nsp=2\mathcal{N}_{sp}=2 on this sequence, but no general-MM statement exists). Each of (a), (b), and (c) is itself a required deliverable.

  2. Pauli principle. Prove, or refute by an explicit example, that no solution of the spin-1/21/2 Bethe equations with repeated roots is physical — i.e. that the number Nstrange(N,M)\mathcal{N}_{strange}(N,M) of repeated-root ("strange") solutions that nevertheless yield eigenvectors of the Hamiltonian vanishes for all N,MN,M — thereby closing the distinct-roots hypothesis on which all current completeness proofs rest.

  3. Thermodynamic limit of the singular sector. Classify the unphysical singular solutions and determine their asymptotic abundance as N→∞N\to\infty — in particular, whether Ns−Nsp\mathcal{N}_s-\mathcal{N}_{sp} is of the same order of magnitude as (NM)−(NM−1)\binom{N}{M}-\binom{N}{M-1} for M∼N/2M\sim N/2, i.e. whether there exist constants 0<c≤C0<c\le C with c≤(Ns−Nsp)/[(NM)−(NM−1)]≤Cc\le\bigl(\mathcal{N}_s-\mathcal{N}_{sp}\bigr)\big/\bigl[\binom{N}{M}-\binom{N}{M-1}\bigr]\le C as N→∞N\to\infty, as the finite-NN data suggest — and give a precise account of the mechanism by which string-hypothesis-based thermodynamic computations nevertheless produce the correct number of highest-weight states.

References

  1. W. Hao, R. I. Nepomechie, A. J. Sommese, "On the completeness of solutions of Bethe's equations", Phys. Rev. E 88, 052113 (2013), DOI 10.1103/PhysRevE.88.052113, arXiv:1308.4645, https://arxiv.org/abs/1308.4645
  2. E. Mukhin, V. Tarasov, A. Varchenko, "Bethe Algebra of Homogeneous XXX Heisenberg Model Has Simple Spectrum", Commun. Math. Phys. 288, 1 (2009), DOI 10.1007/s00220-009-0733-4, arXiv:0706.0688, https://arxiv.org/abs/0706.0688
  3. R. I. Nepomechie, C. Wang, "Algebraic Bethe ansatz for singular solutions", J. Phys. A 46, 325002 (2013), DOI 10.1088/1751-8113/46/32/325002, arXiv:1304.7978, https://arxiv.org/abs/1304.7978
  4. R. I. Nepomechie, C. Wang, "Twisting singular solutions of Bethe's equations", J. Phys. A 47, 505004 (2014), DOI 10.1088/1751-8113/47/50/505004, arXiv:1409.7382, https://arxiv.org/abs/1409.7382
  5. A. N. Kirillov, R. Sakamoto, "Singular solutions to the Bethe ansatz equations and rigged configurations", J. Phys. A 47, 205207 (2014), DOI 10.1088/1751-8113/47/20/205207, arXiv:1402.0651, https://arxiv.org/abs/1402.0651
  6. P. R. Giri, T. Deguchi, "Singular eigenstates in the even(odd) length Heisenberg spin chain", J. Phys. A 48, 175207 (2015), DOI 10.1088/1751-8113/48/17/175207, arXiv:1411.5839, https://arxiv.org/abs/1411.5839
  7. T. Deguchi, P. R. Giri, "Exact quantum numbers of collapsed and non-collapsed two-string solutions in the spin-1/2 Heisenberg spin chain", J. Phys. A 49, 174001 (2016), DOI 10.1088/1751-8113/49/17/174001, arXiv:1509.00108, https://arxiv.org/abs/1509.00108
  8. C. Marboe, D. Volin, "Fast analytic solver of rational Bethe equations", J. Phys. A 50, 204002 (2017), DOI 10.1088/1751-8121/aa6b88, arXiv:1608.06504, https://arxiv.org/abs/1608.06504
  9. Y. Jiang, Y. Zhang, "Algebraic geometry and Bethe ansatz. Part I. The quotient ring for BAE", JHEP 03 (2018) 087, DOI 10.1007/JHEP03(2018)087, arXiv:1710.04693, https://arxiv.org/abs/1710.04693
  10. E. Granet, J. L. Jacobsen, "On zero-remainder conditions in the Bethe ansatz", JHEP 03 (2020) 178, DOI 10.1007/JHEP03(2020)178, arXiv:1910.07797, https://arxiv.org/abs/1910.07797
  11. Z. Bajnok, E. Granet, J. L. Jacobsen, R. I. Nepomechie, "On Generalized Q-systems", JHEP 03 (2020) 177, DOI 10.1007/JHEP03(2020)177, arXiv:1910.07805, https://arxiv.org/abs/1910.07805
  12. J. Hou, Y. Jiang, R.-D. Zhu, "Spin-s Rational Q-system", SciPost Phys. 16, 113 (2024), DOI 10.21468/SciPostPhys.16.4.113, arXiv:2303.07640, https://arxiv.org/abs/2303.07640
  13. W. Hao, R. I. Nepomechie, A. J. Sommese, "Singular solutions, repeated roots and completeness for higher-spin chains", J. Stat. Mech. P03024 (2014), DOI 10.1088/1742-5468/2014/03/P03024, arXiv:1312.2982, https://arxiv.org/abs/1312.2982
  14. D. Chernyak, S. Leurent, D. Volin, "Completeness of Wronskian Bethe equations for rational gl(m|n) spin chains", Commun. Math. Phys. 391, 969 (2022), DOI 10.1007/s00220-021-04275-9, arXiv:2004.02865, https://arxiv.org/abs/2004.02865