The NN-state chiral Potts model is exactly solvable — its Boltzmann weights satisfy the star-triangle (Yang–Baxter) relation — yet it lacks the rapidity-difference property that underlies the solution of the Ising and eight-vertex models, and its correlation functions have never been uniformized. The broken-rapidity-line method produces a two-rapidity function Gpq(r)G_{pq}(r) whose specialization gives the order parameter Mr=kr(N−r)/N2M_r = k^{r(N-r)/N^2}, conjectured in 1989 and proved by Baxter (2005) — a proof possible because MrM_r is rapidity-independent, so it suffices to know Gpq(r)G_{pq}(r) in the superintegrable case q=Vpq = Vp, the only case where it has been obtained. He then (2006) expressed (xp−1/3S(tp))3(x_p^{-1/3}S(t_p))^3 — the ingredient of Gpq(r)G_{pq}(r) in that case — as a product of functions of the hyperelliptic variables zp,wpz_p, w_p: the first thermodynamic quantity of the N>2N>2 chiral Potts model ever so expressed. The open question he posed: does such a representation exist for arbitrary independent rapidities p,qp, q? His own results prove an obstruction — the automorphism R6R^6 fixes zp,wpz_p, w_p but not Gpq(r)G_{pq}(r), so Gpq(r)G_{pq}(r) is not single-valued in these variables — while the combination Lpq(r)=Gpq(r) GRq,Rp(r)L_{pq}(r) = G_{pq}(r)\,G_{Rq,Rp}(r) is invariant and 'may be so expressible'; the automorphism MM also fixes zp,wpz_p, w_p but shifts rr, so residual rapidity monodromy must enter any answer. The problem: construct such a hyperelliptic representation of Lpq(r)L_{pq}(r) for arbitrary independent p,qp, q — a single-valued meromorphic function of zp,wp,zq,wqz_p, w_p, z_q, w_q (and the nome xx) times explicit residual monodromy factors of the fixed algebraic kind carried by the proven prefactor xp−1/3x_p^{-1/3} — or prove that Lpq(r)L_{pq}(r) admits none.

Background

The chiral Potts model assigns spins σ∈{0,…,N−1}\sigma \in \{0,\ldots,N-1\} to the sites of a lattice with ZNZ_N-symmetric interactions (ω=e2πi/N\omega = e^{2\pi i/N}). Its integrable version has Boltzmann weights satisfying the star-triangle (Yang–Baxter) relation, so transfer matrices T(p)T(p), T(q)T(q) commute; each matrix is labeled by a rapidity pp, a point on a high-genus algebraic curve with coordinates (xp,yp,μp)(x_p, y_p, \mu_p) (together with the dependent symmetric variable tp=xpypt_p = x_p y_p). Crucially, the weights are not functions of any difference of rapidities (no 'rapidity-difference property'), which is why the Bethe-ansatz and corner-transfer-matrix machinery that solves the Ising and eight-vertex models fails here. The order parameters are Mr=⟨ωrσ⟩M_r = \langle \omega^{r\sigma}\rangle; Albertini–McCoy–Perk–Tang (1989) conjectured Mr=kr(N−r)/N2M_r = k^{r(N-r)/N^2}, where the modulus kk and complementary modulus k′k' are fixed by the nome xx through −k′2=27x∏n≥1((1−x3n)/(1−xn))12-k'^2 = 27x\prod_{n\ge1}\big((1-x^{3n})/(1-x^n)\big)^{12}; for N=2N=2 this reduces to the Onsager–Yang magnetization k1/4k^{1/4}. The broken rapidity line method (adapted from Jimbo–Miwa–Nakayashiki) cuts one rapidity line of the lattice so the two sides carry different rapidities p,qp, q, producing a two-rapidity correlation function F~pq(r)\tilde F_{pq}(r); its ratios Gpq(r)=F~pq(r)/F~pq(r−1)G_{pq}(r) = \tilde F_{pq}(r)/\tilde F_{pq}(r-1) satisfy Gpq(1)⋯Gpq(N)=1G_{pq}(1)\cdots G_{pq}(N) = 1 and give Mr=Gpp(1)⋯Gpp(r)M_r = G_{pp}(1)\cdots G_{pp}(r) when the line is restored (q→pq \to p). In the superintegrable case q=Vpq = Vp — VV an automorphism of the rapidity curve with xq=xpx_q = x_p, yq=ωypy_q = \omega y_p, μq=μp\mu_q = \mu_p — Baxter (2005) proved the conjecture (under an analyticity assumption of the same type as in the free-energy derivation), obtaining Gpq(r)=Cr S(tp)ϵ(r)G_{pq}(r) = C_r\,S(t_p)^{\epsilon(r)} with Cr=k(N+1−2r)/N2C_r = k^{(N+1-2r)/N^2}, ϵ(0)=ϵ(N)=1−N\epsilon(0) = \epsilon(N) = 1-N and ϵ(r)=1\epsilon(r) = 1 otherwise; here SS is a function of the single variable tpt_p, living on the branched tpt_p-plane, normalized by S(0)=1S(0) = 1 and defined by a contour integral in Baxter (2005). For N=3N=3 the rapidity admits a hyperelliptic parametrization by variables (zp,wp)(z_p, w_p), tied to (xp,yp,μp)(x_p, y_p, \mu_p) through elliptic-like functions h,φh, \varphi; Baxter (2006) showed one must use his 'alternative' (section-5) parametrization, in which the curve automorphisms R,S,V,MR, S, V, M act by zRp=−xwpz_{Rp} = -xw_p, wRp=wp/zpw_{Rp} = w_p/z_p, zSp=−1/(xwp)z_{Sp} = -1/(xw_p), wSp=−1/(xzp)w_{Sp} = -1/(xz_p), zVp=−1/wpz_{Vp} = -1/w_p, wVp=zp/wpw_{Vp} = z_p/w_p, while the automorphism MM (which multiplies μp\mu_p by ω\omega, leaving xp,yp,tpx_p, y_p, t_p fixed) also fixes zp,wpz_p, w_p. Analytic continuation of zp,wpz_p, w_p at fixed tpt_p takes values xmzp0±1x^m z_{p_0}^{\pm1}, xnwp0±1x^n w_{p_0}^{\pm1} with m+n≡0(mod3)m + n \equiv 0 \pmod 3 on even sheets and ≡1(mod3)\equiv 1 \pmod 3 on odd sheets; the Riemann sheet graph is a honeycomb lattice. The main result of Baxter (2006) is (xp−1/3S(tp))3=k−1/3G(−1/wp)/G(zp)(x_p^{-1/3}S(t_p))^3 = k^{-1/3}G(-1/w_p)/G(z_p) with G(z)=F(ωz)/F(ω2z)G(z) = F(\omega z)/F(\omega^2 z) and F(z)=∏j≥1(1−xjz)j/(1−xjz−1)jF(z) = \prod_{j\ge1}(1-x^j z)^j/(1-x^j z^{-1})^j — infinite products generalizing elliptic functions in that factors carry the power jj. The R6R^6 obstruction: the automorphism rules give zR6p=zpz_{R^6p} = z_p and wR6p=wpw_{R^6p} = w_p, yet GR6p,q(r)≠Gpq(r)G_{R^6p,q}(r) \neq G_{pq}(r) and Gp,R6q(r)≠Gpq(r)G_{p,R^6q}(r) \neq G_{pq}(r); hence Gpq(r)G_{pq}(r) cannot be a single-valued function of zp,wp,zq,wqz_p, w_p, z_q, w_q. However the combination Lpq(r)=Gpq(r)GRq,Rp(r)L_{pq}(r) = G_{pq}(r)G_{Rq,Rp}(r), introduced in Baxter (1998), satisfies closed functional relations (LR2Np,q(r)=(−1)N−1Lpq(r)L_{R^{2N}p,q}(r) = (-1)^{N-1}L_{pq}(r), Lqp(r)=1/Lpq(N−r+1)L_{qp}(r) = 1/L_{pq}(N-r+1), LMp,q(r)=Lp,M−1q(r)=Lpq(r+1)L_{Mp,q}(r) = L_{p,M^{-1}q}(r) = L_{pq}(r+1), Lpq(1)⋯Lpq(N)=1L_{pq}(1)\cdots L_{pq}(N) = 1) and is unchanged by p→R6pp \to R^6p and q→R6qq \to R^6q, so it may admit such a representation. One further necessary condition follows from these same ingredients: MM also fixes zp,wpz_p, w_p for each rapidity, while the 1998 relations give Lp,M−1q(r)=Lpq(r+1)L_{p,M^{-1}q}(r) = L_{pq}(r+1), which differs from Lpq(r)L_{pq}(r) in general — so at fixed rr even Lpq(r)L_{pq}(r) is not a single-valued function of zp,wp,zq,wqz_p, w_p, z_q, w_q alone. The four hyperelliptic variables leave the μ\mu-branch of each rapidity undetermined (the choice among p,Mp,M2pp, Mp, M^2p), and an explicit answer must carry that residual monodromy by explicit factors — fixed algebraic functions of the rapidity coordinates, of the same kind as the prefactor xp−1/3x_p^{-1/3} carried by the proven single-rapidity identity — with all essential rapidity dependence residing in the single-valued meromorphic part.

Problem Statement

Let p,qp, q be arbitrary independent rapidities of the N=3N=3 integrable chiral Potts model, and let Gpq(r)G_{pq}(r) (r=1,…,Nr = 1, \ldots, N with N=3N = 3; Gpq(1)⋯Gpq(N)=1G_{pq}(1)\cdots G_{pq}(N) = 1) be the generalized order-parameter ratios of the broken rapidity line method, Gpq(r)=F~pq(r)/F~pq(r−1)G_{pq}(r) = \tilde F_{pq}(r)/\tilde F_{pq}(r-1), where F~pq(r)\tilde F_{pq}(r) is the spin correlation at the broken rapidity line. Known: (i) in the superintegrable case q=Vpq = Vp, Gpq(r)=k(N+1−2r)/N2S(tp)ϵ(r)G_{pq}(r) = k^{(N+1-2r)/N^2}S(t_p)^{\epsilon(r)} with ϵ(0)=ϵ(N)=1−N\epsilon(0) = \epsilon(N) = 1-N, ϵ(r)=1\epsilon(r) = 1 otherwise, and (xp−1/3S(tp))3=k−1/3G(−1/wp)/G(zp)(x_p^{-1/3}S(t_p))^3 = k^{-1/3}G(-1/w_p)/G(z_p) with G(z)=F(ωz)/F(ω2z)G(z) = F(\omega z)/F(\omega^2 z), F(z)=∏j≥1(1−xjz)j/(1−xjz−1)jF(z) = \prod_{j\ge1}(1-x^j z)^j/(1-x^j z^{-1})^j; (ii) Gpq(r)G_{pq}(r) cannot be a single-valued meromorphic function of the hyperelliptic variables (zp,wp,zq,wq)(z_p, w_p, z_q, w_q) of the section-5 parametrization, because R6R^6 fixes all four variables while changing GG; (iii) the automorphism MM, which likewise fixes zp,wpz_p, w_p (of each rapidity) while sending μ→ωμ\mu \to \omega\mu, shifts rr through the 1998 relation Lp,M−1q(r)=Lpq(r+1)L_{p,M^{-1}q}(r) = L_{pq}(r+1), so at fixed rr even Lpq(r)L_{pq}(r) is not a single-valued function of these four variables alone — the variables leave the μ\mu-branch of each rapidity undetermined, and residual factors of the kind specified below must carry it, just as the proven identity (i) carries the prefactor xp−1/3x_p^{-1/3}. Open question: can Lpq(r)=Gpq(r) GRq,Rp(r)L_{pq}(r) = G_{pq}(r)\,G_{Rq,Rp}(r) — the combination proposed by Baxter (2006), unchanged by p→R6pp \to R^6p and q→R6qq \to R^6q — be simply expressed for arbitrary independent p,qp, q in terms of the section-5 hyperelliptic variables, that is, as a single-valued meromorphic function of zp,wp,zq,wqz_p, w_p, z_q, w_q (and the nome xx) times explicit residual monodromy factors, where these factors are fixed algebraic functions of the rapidity coordinates (xp,yp,μp)(x_p, y_p, \mu_p), (xq,yq,μq)(x_q, y_q, \mu_q) — of the same kind as the prefactor xp−1/3x_p^{-1/3} in (i) — which carry only the monodromy under the honeycomb sheet structure (zp→xmzp±1z_p \to x^m z_p^{\pm1}, wp→xnwp±1w_p \to x^n w_p^{\pm1}, m+n≡0,1(mod3)m+n \equiv 0,1 \pmod 3) and under the automorphism fibers (R6R^6, MM) that fix the hyperelliptic variables, while all essential rapidity dependence resides in the single-valued meromorphic function (so a factorization with constant single-valued part, absorbing Lpq(r)L_{pq}(r) itself into the residual factor, is not admissible)? Construct such a representation of Lpq(r)L_{pq}(r), or prove that Lpq(r)L_{pq}(r) admits none (for instance by showing that its monodromy under the honeycomb sheet structure and the automorphism fibers is incompatible with any such factorization). Any such representation must reduce to (i) at q=Vpq = Vp, give Gpp(r)=k(N+1−2r)/N2G_{pp}(r) = k^{(N+1-2r)/N^2} and Mr=kr(N−r)/N2M_r = k^{r(N-r)/N^2} when the line is restored, and satisfy the 1998 functional relations for LpqL_{pq} (including LR2Np,q(r)=(−1)N−1Lpq(r)L_{R^{2N}p,q}(r) = (-1)^{N-1}L_{pq}(r), Lqp(r)=1/Lpq(N−r+1)L_{qp}(r) = 1/L_{pq}(N-r+1), LMp,q(r)=Lp,M−1q(r)=Lpq(r+1)L_{Mp,q}(r) = L_{p,M^{-1}q}(r) = L_{pq}(r+1)).

References

  1. R. J. Baxter, 'Hyperelliptic parametrisation of the generalised order parameter of the N=3 chiral Potts model', The ANZIAM Journal 47 (2006) 289-307. arXiv:cond-mat/0601325. DOI: 10.1017/S1446181100009858. URL: https://doi.org/10.1017/S1446181100009858
  2. R. J. Baxter, 'The order parameter of the chiral Potts model', Journal of Statistical Physics 120 (2005) 1-36. arXiv:cond-mat/0501226. DOI: 10.1007/s10955-005-5534-3. URL: https://doi.org/10.1007/s10955-005-5534-3
  3. R. J. Baxter, 'Derivation of the order parameter of the chiral Potts model', Physical Review Letters 94 (2005) 130602. arXiv:cond-mat/0501227. DOI: 10.1103/PhysRevLett.94.130602. URL: https://doi.org/10.1103/PhysRevLett.94.130602
  4. R. J. Baxter, 'Functional relations for the order parameters of the chiral Potts model', Journal of Statistical Physics 91 (1998) 499-524. DOI: 10.1023/A:1023096408679. URL: https://doi.org/10.1023/A:1023096408679
  5. R. J. Baxter, 'Functional relations for the order parameters of the chiral Potts model: low-temperature expansions', Physica A 260 (1998) 117-130. arXiv:cond-mat/9808120. DOI: 10.1016/S0378-4371(98)00307-0. URL: https://doi.org/10.1016/S0378-4371(98)00307-0
  6. R. J. Baxter, 'Some hyperelliptic function identities that occur in the chiral Potts model', Journal of Physics A: Mathematical and General 31 (1998) 6807-6818. DOI: 10.1088/0305-4470/31/32/005. URL: https://doi.org/10.1088/0305-4470/31/32/005
  7. R. J. Baxter, 'Free energy of the three-state τ2(tq)\tau_2(t_q) model as a product of elliptic functions', Proceedings of the Royal Society A 462 (2006) 3323-3341. arXiv:cond-mat/0602367. DOI: 10.1098/rspa.2006.1734. URL: https://doi.org/10.1098/rspa.2006.1734
  8. R. J. Baxter, 'The challenge of the chiral Potts model', Journal of Physics: Conference Series 42 (2006) 11-24. arXiv:cond-mat/0510683. DOI: 10.1088/1742-6596/42/1/003. URL: https://doi.org/10.1088/1742-6596/42/1/003
  9. R. J. Baxter, 'Corner transfer matrices in statistical mechanics', Journal of Physics A: Mathematical and Theoretical 40 (2007) 12577-12588. arXiv:cond-mat/0611167. DOI: 10.1088/1751-8113/40/42/S05. URL: https://doi.org/10.1088/1751-8113/40/42/S05
  10. R. J. Baxter, 'Some comments on developments in exact solutions in statistical mechanics since 1944', arXiv:1010.0710 (2010). URL: https://arxiv.org/abs/1010.0710
  11. G. Albertini, B. M. McCoy, J. H. H. Perk and S. Tang, 'Excitation spectrum and order parameter for the integrable N-state chiral Potts model', Nuclear Physics B 314 (1989) 741-763. DOI: 10.1016/0550-3213(89)90415-X. URL: https://doi.org/10.1016/0550-3213(89)90415-X
  12. M. Jimbo, T. Miwa and A. Nakayashiki, 'Difference equations for the correlation functions of the eight-vertex model', Journal of Physics A: Mathematical and General 26 (1993) 2199-2209. DOI: 10.1088/0305-4470/26/9/015. URL: https://doi.org/10.1088/0305-4470/26/9/015
  13. H. Au-Yang and J. H. H. Perk, 'Spontaneous magnetization of the integrable chiral Potts model', Journal of Physics A: Mathematical and Theoretical 44 (2011) 445005. arXiv:1003.4805. DOI: 10.1088/1751-8113/44/44/445005. URL: https://doi.org/10.1088/1751-8113/44/44/445005