The spin-1/2 XX chain in a staggered longitudinal magnetic field (the Kontorovich–Tsukernik two-sublattice model) has been solvable by Jordan–Wigner free-fermion methods since 1968, but no R-matrix for it was known within the quantum inverse scattering method. In 2024, Bibikov conjectured an explicit 16×1616\times16 elliptic R-matrix R(λ,hst)R(\lambda,h_{\rm st}) whose entries are polynomials in three Jacobi elliptic functions α(λ)\alpha(\lambda), γ(λ)\gamma(\lambda), θ(λ)\theta(\lambda) with moduli determined by the staggered field hsth_{\rm st}; the matrix was found by a Taylor-expansion matching heuristic. The Yang–Baxter equation (YBE) for this matrix is rigorously proven only in the two trigonometric degenerations hst=0h_{\rm st}=0 and hst=ih_{\rm st}=i; for generic hsth_{\rm st} it is verified only by computer-algebra Taylor expansion up to total order 8, and the author states that the analytic machinery for a rigorous proof, and for the auxiliary identities among α\alpha, γ\gamma, θ\theta on which the construction rests, is unclear to him. The problem posed here is to settle this rigorously: prove, or refute, that R(λ,hst)R(\lambda,h_{\rm st}) satisfies the braid-form YBE for generic hsth_{\rm st}, and prove (or refute) the four auxiliary elliptic-function identities conjectured alongside it. A proof would give the first R-matrix construction of the staggered XX chain and add a new elliptic solution of the Yang–Baxter equation to the short list known since Baxter's eight-vertex model.

Background

The spin-1/2 XX chain, introduced by Lieb, Schultz and Mattis (1961), is the nearest-neighbor spin chain with Hamiltonian H^=∑n12(Sn+Sn+1−+Sn−Sn+1+)\hat H=\sum_n\frac12(S_n^+S_{n+1}^-+S_n^-S_{n+1}^+), where Sn±S_n^\pm and SnzS_n^z are spin-1/2 raising, lowering and zz operators at site nn. The Jordan–Wigner transformation, a rewrite of the spins as fermions that maps the XX chain onto free hopping fermions, makes its spectrum and thermodynamics exactly computable, and the model remains a canonical testing ground of statistical mechanics (Lieb, Schultz and Mattis 1961; Göhmann, Kozlowski, Sirker and Suzuki 2019). A stronger and more structural form of solvability is encoded in the quantum inverse scattering method: an R-matrix, i.e. a matrix R(λ)R(\lambda) depending on a spectral parameter λ\lambda, such that the local Hamiltonian density is proportional to dR/dλdR/d\lambda at λ=0\lambda=0 and such that RR satisfies the Yang–Baxter equation (YBE). In the braid form used throughout this problem the YBE reads

R12(λ−μ) R23(λ) R12(μ)=R23(μ) R12(λ) R23(λ−μ),R_{12}(\lambda-\mu)\,R_{23}(\lambda)\,R_{12}(\mu)=R_{23}(\mu)\,R_{12}(\lambda)\,R_{23}(\lambda-\mu),

where R12=R⊗IR_{12}=R\otimes I and R23=I⊗RR_{23}=I\otimes R act on three consecutive tensor factors. This is the difference form of the YBE: each factor depends on a single spectral parameter, obtained from the standard three-parameter form by setting the third parameter to zero. The YBE guarantees commuting transfer matrices and hence an extensive family of conserved charges. For the XX chain the R-matrix is a 4×44\times4 trigonometric matrix. Adding a staggered longitudinal magnetic field hst∑n(−1)nSnzh_{\rm st}\sum_n(-1)^nS_n^z breaks translation invariance to a two-site unit cell. Kontorovich and Tsukernik (1968) showed that this spin array with two sublattices is still solvable by the Jordan–Wigner free-fermion method, but no R-matrix for it was known within the quantum inverse scattering method. Because of the two-site cell, the Hamiltonian density becomes a 16×1616\times16 matrix HH acting on a block of four spins, viewed as C4⊗C4\mathbb C^4\otimes\mathbb C^4 where each C4\mathbb C^4 is a pair of spin sites. Bibikov (arXiv:2407.10395) verified that this HH satisfies the Reshetikhin condition, [H⊗I4+I4⊗H, [H⊗I4,I4⊗H]]=I4⊗K−K⊗I4[H\otimes I_4+I_4\otimes H,\,[H\otimes I_4,I_4\otimes H]]=I_4\otimes K-K\otimes I_4 for an explicitly computable KK, which is a necessary condition for the existence of an R-matrix with R(λ)=I+CHλ+o(λ)R(\lambda)=I+C H\lambda+o(\lambda), together with two higher-order necessary conditions, and constructed the R-matrix explicitly at the degenerate imaginary field hst=ih_{\rm st}=i. For generic hsth_{\rm st}, Bibikov (arXiv:2409.19775) conjectured a 16×1616\times16 R-matrix obtained by a Taylor-expansion matching heuristic (Bibikov 2003): expand the ratios of the unknown entries in Taylor series, guess closed forms, and verify. The resulting matrix is symmetric, has a sparse entry pattern, and all its entries are polynomials in three Jacobi elliptic functions

α(λ)=iκ sn(λ/(iκ), κ2),γ(λ)=ξ sn(ζλ, ξ2),θ(λ)=cn(iηλ, 1/η),\alpha(\lambda)=i\kappa\,{\rm sn}\bigl(\lambda/(i\kappa),\,\kappa^2\bigr),\qquad \gamma(\lambda)=\xi\,{\rm sn}\bigl(\zeta\lambda,\,\xi^2\bigr),\qquad \theta(\lambda)={\rm cn}\bigl(i\eta\lambda,\,1/\eta\bigr),

where sn and cn are the Jacobi elliptic functions, meromorphic doubly periodic generalizations of sine and cosine taken in the convention of the paper's Maple appendix (Maple's JacobiSN/JacobiCN), in which the second argument is the modulus kk (the parameter is m=k2m=k^2): the moduli of α\alpha, γ\gamma, θ\theta are κ2\kappa^2, ξ2\xi^2, 1/η1/\eta respectively. The four quantities κ\kappa, ξ\xi, η\eta, ζ\zeta are fixed by the single physical parameter through hst=(1−κ2)/(2κ)h_{\rm st}=(1-\kappa^2)/(2\kappa), ξ=(1−κ)/(1+κ)\xi=(1-\kappa)/(1+\kappa), η=(1+κ2)/(2κ)\eta=(1+\kappa^2)/(2\kappa), ζ=(1+κ)2/(2κ)\zeta=(1+\kappa)^2/(2\kappa). Elliptic R-matrices are rare: since Baxter's 1972 solution of the eight-vertex lattice model, the known elliptic difference-form solutions of the YBE essentially form one small family, so a putative new one is noteworthy. The verification status is asymmetric. At the two degenerate values of the field the elliptic functions collapse to trigonometric or hyperbolic ones (hst=0h_{\rm st}=0: α=tan⁡λ\alpha=\tan\lambda, γ=0\gamma=0, θ=1/cos⁡λ\theta=1/\cos\lambda; hst=ih_{\rm st}=i: α=tanh⁡λ\alpha=\tanh\lambda, γ=itanh⁡λ\gamma=i\tanh\lambda, θ=1\theta=1) and the YBE is proven exactly by direct calculation. For generic hsth_{\rm st}, the paper verifies the YBE only by computer-algebra Taylor expansion of the matrix entries up to total order 8 (a Maple file whose truncation order is set by the author's computer, with an invitation to the reader to increase it), and verifies four auxiliary identities among α\alpha, γ\gamma, θ\theta only up to order 30. The author states that a complete rigorous verification 'seems to be far from the elementary theory of elliptic functions' and that the corresponding machinery, like the proofs of the auxiliary identities, is unclear to him. The question is whether this finite-order evidence can be replaced by rigorous analytic proof.

Problem Statement

Fix a complex staggered field hsth_{\rm st} away from the two degenerate values 00 and ii, and define the moduli hst=(1−κ2)/(2κ)h_{\rm st}=(1-\kappa^2)/(2\kappa), ξ=(1−κ)/(1+κ)\xi=(1-\kappa)/(1+\kappa), η=(1+κ2)/(2κ)\eta=(1+\kappa^2)/(2\kappa), ζ=(1+κ)2/(2κ)\zeta=(1+\kappa)^2/(2\kappa), with the explicit branch choice of Eq. (20) of Bibikov, arXiv:2409.19775: κ=1+hst2−hst\kappa=\sqrt{1+h_{\rm st}^2}-h_{\rm st}, ξ=1+1/hst2−1/hst\xi=\sqrt{1+1/h_{\rm st}^2}-1/h_{\rm st}, η=1+hst2\eta=\sqrt{1+h_{\rm st}^2}, ζ=hst/ξ\zeta=h_{\rm st}/\xi, together with the Jacobi elliptic functions

α(λ)=iκ sn(λ/(iκ), κ2),γ(λ)=ξ sn(ζλ, ξ2),θ(λ)=cn(iηλ, 1/η).\alpha(\lambda)=i\kappa\,{\rm sn}\bigl(\lambda/(i\kappa),\,\kappa^2\bigr),\qquad \gamma(\lambda)=\xi\,{\rm sn}\bigl(\zeta\lambda,\,\xi^2\bigr),\qquad \theta(\lambda)={\rm cn}\bigl(i\eta\lambda,\,1/\eta\bigr).

Here sn and cn are taken in the convention of the paper's Maple appendix (Maple's JacobiSN/JacobiCN), where the second argument is the modulus kk (the parameter is m=k2m=k^2); the moduli of α\alpha, γ\gamma, θ\theta are κ2\kappa^2, ξ2\xi^2, 1/η1/\eta respectively. Let R(λ,hst)R(\lambda,h_{\rm st}) be the symmetric 16×1616\times16 matrix acting on C4⊗C4\mathbb C^4\otimes\mathbb C^4 whose entries are specified by Eqs. (14)–(18) of Bibikov, arXiv:2409.19775 (entries Rij=Aij+BijR_{ij}=A_{ij}+B_{ij}, with AijA_{ij} and BijB_{ij} the polynomials in α(λ)\alpha(\lambda), γ(λ)\gamma(\lambda), θ(λ)\theta(\lambda), hsth_{\rm st} given there, the overall normalization A11≡1A_{11}\equiv1 of the paper's Maple file 'dataR', and the reflection symmetries A17−i,17−j=AijA_{17-i,17-j}=A_{ij}, B17−i,17−j=−BijB_{17-i,17-j}=-B_{ij}). (1) Prove, or rigorously refute, the conjecture that for generic hsth_{\rm st} (as an identity in the indeterminate hsth_{\rm st}, or for every complex hsth_{\rm st} outside the degenerate values) this matrix satisfies the braid-form Yang–Baxter equation

R12(λ−μ) R23(λ) R12(μ)=R23(μ) R12(λ) R23(λ−μ),R12=R⊗I4,    R23=I4⊗R,R_{12}(\lambda-\mu)\,R_{23}(\lambda)\,R_{12}(\mu)=R_{23}(\mu)\,R_{12}(\lambda)\,R_{23}(\lambda-\mu),\qquad R_{12}=R\otimes I_4,\;\; R_{23}=I_4\otimes R,

as an identity of meromorphic functions of the spectral parameters λ\lambda and μ\mu on (C4)⊗3(\mathbb C^4)^{\otimes3}. A refutation may consist of a finite total Taylor order at which some entry of the difference of the two sides fails to vanish. (2) Prove, or rigorously refute, the auxiliary identities conjectured in the same paper for the above parametrization:

γ (1+α2)=hst α (1−γ2),γ θ2=hst α,\gamma\,(1+\alpha^2)=h_{\rm st}\,\alpha\,(1-\gamma^2),\qquad \gamma\,\theta^2=h_{\rm st}\,\alpha, dαdλ=1+α2+2hstαγ,dγdλ=hst(1−γ2)−2αγ.\frac{d\alpha}{d\lambda}=1+\alpha^2+2h_{\rm st}\alpha\gamma,\qquad \frac{d\gamma}{d\lambda}=h_{\rm st}(1-\gamma^2)-2\alpha\gamma.

The single scientific objective is to replace the finite-order computer-algebra evidence for this conjectural R-matrix by a complete rigorous analytic foundation (or by a definitive refutation); the auxiliary identities are the conjectured elliptic-function backbone of the construction and are part of the same objective, not an independent add-on. No method is prescribed. The problem is complete when both (1) and (2) are settled, each by proof or by refutation.

References

  1. P. N. Bibikov, R-matrix for the XX spin chain in a staggered magnetic field, arXiv:2409.19775 (2024), DOI: 10.48550/arxiv.2409.19775, https://arxiv.org/abs/2409.19775
  2. P. N. Bibikov, R-matrix for the XX spin chain at the special imaginary value of staggered magnetic field, arXiv:2407.10395 (2024), DOI: 10.48550/arxiv.2407.10395, https://arxiv.org/abs/2407.10395
  3. P. N. Bibikov, Thermodynamics of the XX spin chain in the QTM approach, arXiv:2503.21831 (2025), DOI: 10.48550/arxiv.2503.21831, https://arxiv.org/abs/2503.21831
  4. P. N. Bibikov, How to solve Yang–Baxter equation using the Taylor expansion of R-matrix, Physics Letters A 314, 209–213 (2003), DOI: 10.1016/S0375-9601(03)00818-1, arXiv:nlin/0112001, https://arxiv.org/abs/nlin/0112001
  5. Elliott Lieb, Theodore Schultz, Daniel Mattis, Two soluble models of an antiferromagnetic chain, Annals of Physics 16, 407–466 (1961), DOI: 10.1016/0003-4916(61)90115-4, https://doi.org/10.1016/0003-4916(61)90115-4
  6. V. M. Kontorovich, V. M. Tsukernik, Magnetic properties of a spin array with two sublattices, Soviet Physics JETP 26, 687–691 (1968).
  7. Rodney J. Baxter, Partition function of the Eight-Vertex lattice model, Annals of Physics 70, 193–228 (1972), DOI: 10.1016/0003-4916(72)90335-1, https://doi.org/10.1016/0003-4916(72)90335-1
  8. Frank Göhmann, Karol K. Kozlowski, Jesko Sirker, Junji Suzuki, Equilibrium dynamics of the XX chain, Physical Review B 100, 155428 (2019), DOI: 10.1103/PhysRevB.100.155428, https://doi.org/10.1103/PhysRevB.100.155428