The ZpZ_p chiral Potts model — the most general known solution of the star-triangle (Yang–Baxter) relation, reducing to the Ising model at p=2p=2 — has order parameters MrM_r whose closed form Mr=(1−k′2)r(p−r)/(2p2)M_r=(1-k'^2)^{r(p-r)/(2p^2)} was conjectured in 1983–1989, derived by Baxter (2005) for the general model under unproven analyticity assumptions, and proven from the finite lattice only on the superintegrable subvariety. Grosjean–Maillet–Niccoli (2013, published 2015) expressed finite-size, inhomogeneous chiral Potts matrix elements — including the order-parameter element Mp−1M_{p-1} — as determinants via Sklyanin's separation of variables, and announced a future analysis of the homogeneous and thermodynamic limits aiming to prove the order-parameter formulae for general homogeneous ZpZ_p representations from the finite lattice. That announced program remains open. The open problem is the limit analysis of these determinant representations and the resulting unconditional, general-representation proof of the order-parameter formulae, together with the auxiliary goal of single-determinant formulae for the vnv_n generators.

Background

The ZpZ_p chiral Potts model is a two-dimensional lattice model in which each site carries a spin taking pp values (a clock variable ZZ with Zp=1Z^p=1), with Boltzmann weights invariant under the cyclic group ZpZ_p. Unlike the Ising model, its weights depend on two points ('rapidities') lying on a high-genus algebraic curve, and it is the most general known solution of the star-triangle relation — the form of the Yang–Baxter equation for edge-interacting spin models. Its transfer matrix also generates a one-dimensional quantum chain of pp-state clock variables. In the ordered (ferromagnetic) phase the ZpZ_p symmetry is spontaneously broken, and the order parameters (the spontaneous-magnetization components) are Mr=lim⁡L→∞⟨Zr⟩,r=1,…,p−1,M_r=\lim_{L\to\infty}\langle Z^r\rangle,\qquad r=1,\dots,p-1, the thermodynamic-limit ground-state expectations of the rr-th power of the local clock variable.

Howes–Kadanoff–den Nijs (1983) evaluated M1M_1 and M2M_2 for p=3p=3 by series expansion in the modulus and conjectured the closed form; Albertini–McCoy–Perk–Tang (1989) extended the series analysis to general pp and conjectured, for all rr, Mr=(1−k′2)r(p−r)/(2p2),M_r=(1-k'^2)^{r(p-r)/(2p^2)}, where k′k' is the modulus of the model's rapidity curve (0<k′<10<k'<1 in the ordered phase, k′→1k'\to 1 at criticality). In complementary-modulus notation k2+k′2=1k^2+k'^2=1 this is Mr=kr(p−r)/p2M_r=k^{r(p-r)/p^2}, the form used by Baxter; Grosjean–Maillet–Niccoli write it as Mr=(1−k2)r(p−r)/(2p2)M_r=(1-k^2)^{r(p-r)/(2p^2)}, their kk being Baxter's k′k'. For p=2p=2 the formula reduces to the Onsager–Yang spontaneous magnetization of the Ising model.

Baxter (2005) derived the conjecture for the general chiral Potts model (generic rapidities), generalizing the broken-rapidity-line functional-relation method introduced by Jimbo–Miwa–Nakayashiki for the eight-vertex model. His derivation is, by his own statement, 'subject to a plausible analyticity assumption' — two numbered assumptions in §7 of the long paper — concerning analytic properties of correlation functions with broken rapidity lines; these are supported by small-modulus evidence but unproven. Unconditional proofs starting from finite lattices exist only on the superintegrable subvariety, the special self-dual class of representations where the chiral Potts quantum chain satisfies quartic commutation relations implying an Onsager-algebra symmetry (the Lie-algebraic structure underlying Onsager's 1944 Ising solution). There, Baxter's generalized Onsager algebra led to a determinantal formula for the spontaneous magnetization, and the Au-Yang–Perk eigenvector constructions, together with the Iorgov–Pakuliak–Shadura–Tykhyy–von Gehlen proof of the conjectured factorized form of the spin-operator matrix elements, gave a second finite-lattice route. Au-Yang–Perk (2011) used ZZ-invariance (insensitivity of the model to deformations of the rapidity lattice) to reduce pair correlations of the general model to superintegrable eigenvectors, obtaining the order parameter 'as a sum almost identical to the one given by Baxter' — still importing Baxter's final evaluation, and resting on eigenvector constructions whose algebraic underpinnings are proven only for some states. Thus the formula's value is settled knowledge in the community, while a derivation that is both unconditional and valid for general representations, starting from the finite lattice, is missing.

Grosjean–Maillet–Niccoli (2013, published 2015) attacked the model with Sklyanin's quantum separation of variables (SOV), a framework in which the commuting transfer matrix acts as multiplication by the eigenvalues of an auxiliary separate-variable operator, so that scalar products of eigenstates become determinants. Working from the Bazhanov–Stroganov model (the six-vertex Yang–Baxter algebra with general cyclic representations, from which the chiral Potts transfer matrices are built) on finite inhomogeneous lattices (lattices whose sites carry independent rapidity parameters), they proved single-determinant formulae for scalar products of separate states, solved the quantum inverse problem — every local operator is reconstructed as a polynomial in the separate variables — and showed that the form factors (normalized matrix elements between transfer-matrix eigenstates) of any local operator are finite sums of determinants, with single-determinant formulae for a basis of elementary operators and for the generators un−1u_n^{-1} of the local Weyl algebra. The Weyl generators satisfy unp=vnp=1u_n^p=v_n^p=1 and unvn=ω vnunu_nv_n=\omega\,v_nu_n with ω\omega a primitive pp-th root of unity, and u1−1u_1^{-1} reduces to the clock-spin variable whose ground-state expectation gives the order parameter: Mp−1=⟨g.s.∣u1−1∣g.s.⟩⟨g.s.∣g.s.⟩,M_{p-1}=\frac{\langle g.s.|u_1^{-1}|g.s.\rangle}{\langle g.s.|g.s.\rangle}, for which their approach yields a simple determinant formula. They observe that these determinant formulae are exactly what 'will reproduce the chiral Potts order parameters in the thermodynamic limit', state that the known formulae were 'proven with techniques working only in the super-integrable case but they are expected to be true for the general homogeneous ZpZ_p chiral Potts models', and announce: 'In a future paper, we will analyze the homogeneous and thermodynamic limits focusing the attention on the derivation of the order parameter formulae for the general homogeneous ZpZ_p chiral Potts models.' They also pose the auxiliary problem of simplifying the form factors of the vnv_n generators — currently finite sums of determinants — to single determinants.

Problem Statement

Analyze the homogeneous and thermodynamic limits of the finite-lattice SOV determinant representations of ZpZ_p chiral Potts matrix elements, and prove from the finite lattice the order-parameter formulae Mr=(1−k′2)r(p−r)/(2p2),r=1,…,p−1,M_r=(1-k'^2)^{r(p-r)/(2p^2)},\qquad r=1,\dots,p-1, for general homogeneous ZpZ_p chiral Potts representations — that is, without restricting to the superintegrable subvariety and without postulating the unproven analyticity properties on which Baxter's broken-rapidity-line derivation rests. Concretely, this program entails constructing the homogeneous limit of the separate-variable basis and of the determinant formulae, controlling the coalescence of the inhomogeneity parameters that defines the separate-variable spectrum; isolating the ground-state sector for general cyclic representations and taking the thermodynamic limit of the relevant matrix elements — in particular the order-parameter element Mp−1=⟨g.s.∣u1−1∣g.s.⟩/⟨g.s.∣g.s.⟩M_{p-1}=\langle g.s.|u_1^{-1}|g.s.\rangle/\langle g.s.|g.s.\rangle (where g.s. denotes the ground state), which admits a single-determinant formula — with controlled interchange of the homogeneous and infinite-volume limits; and evaluating the limit in closed form, recovering the conjectured expression, with the p=2p=2 (Onsager–Yang Ising) and superintegrable specializations reproduced as consistency checks. As an auxiliary goal of the same determinant program, obtain single-determinant representations for the form factors of the Weyl generators vnv_n (currently available only as finite sums of determinants), a simplification that would also lighten the subsequent limit analysis. The common objective is to make the SOV determinant representation of chiral Potts matrix elements yield the physics of the homogeneous model — its order parameters first, and by the same apparatus its correlation functions, which are beyond current techniques for the general model. Proving the order-parameter formulae for all r=1,…,p−1r=1,\dots,p-1 in general homogeneous ZpZ_p representations resolves the task; the single-determinant form factors of the vnv_n generators and the extension to correlation functions are optional extensions and are not required for completion.

References

  1. N. Grosjean, J.-M. Maillet, G. Niccoli, On the Form Factors of Local Operators in the Bazhanov–Stroganov and Chiral Potts Models, Annales Henri Poincaré 16 (2015) 1103–1153, DOI 10.1007/s00023-014-0358-9, arXiv:1309.4701, https://arxiv.org/abs/1309.4701
  2. S. Howes, L. P. Kadanoff, M. den Nijs, Quantum model for commensurate-incommensurate transitions, Nuclear Physics B 215 (1983) 169–208, DOI 10.1016/0550-3213(83)90212-2, https://doi.org/10.1016/0550-3213(83)90212-2
  3. G. Albertini, B. M. McCoy, J. H. H. Perk, 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, https://doi.org/10.1016/0550-3213(89)90415-X
  4. L. Onsager, Crystal Statistics. I. A Two-Dimensional Model with an Order-Disorder Transition, Physical Review 65 (1944) 117–149, DOI 10.1103/PhysRev.65.117, https://doi.org/10.1103/PhysRev.65.117
  5. C. N. Yang, The Spontaneous Magnetization of a Two-Dimensional Ising Model, Physical Review 85 (1952) 808–816, DOI 10.1103/PhysRev.85.808, https://doi.org/10.1103/PhysRev.85.808
  6. M. Jimbo, T. Miwa, 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, https://doi.org/10.1088/0305-4470/26/9/015
  7. R. J. Baxter, Derivation of the Order Parameter of the Chiral Potts Model, Physical Review Letters 94 (2005) 130602, DOI 10.1103/PhysRevLett.94.130602, arXiv:cond-mat/0501227, https://arxiv.org/abs/cond-mat/0501227
  8. R. J. Baxter, The Order Parameter of the Chiral Potts Model, Journal of Statistical Physics 120 (2005) 1–36, DOI 10.1007/s10955-005-5534-3, arXiv:cond-mat/0501226, https://arxiv.org/abs/cond-mat/0501226
  9. R. J. Baxter, Algebraic Reduction of the Ising Model, Journal of Statistical Physics 132 (2008) 959–982, DOI 10.1007/s10955-008-9587-y, https://doi.org/10.1007/s10955-008-9587-y
  10. R. J. Baxter, Some Remarks on a Generalization of the Superintegrable Chiral Potts Model, Journal of Statistical Physics 137 (2009) 798–813, DOI 10.1007/s10955-009-9778-1, https://doi.org/10.1007/s10955-009-9778-1
  11. R. J. Baxter, Proof of the determinantal form of the spontaneous magnetization of the superintegrable chiral Potts model, ANZIAM Journal 51 (2010) 309–316, DOI 10.21914/anziamj.v51i0.2846, https://doi.org/10.21914/anziamj.v51i0.2846
  12. H. Au-Yang, J. H. H. Perk, Quantum loop subalgebra and eigenvectors of the superintegrable chiral Potts transfer matrices, Journal of Physics A: Mathematical and Theoretical 44 (2011) 025205, DOI 10.1088/1751-8113/44/2/025205, arXiv:0907.0362, https://arxiv.org/abs/0907.0362
  13. H. Au-Yang, J. H. H. Perk, Spontaneous magnetization of the integrable chiral Potts model, Journal of Physics A: Mathematical and Theoretical 44 (2011) 445005, DOI 10.1088/1751-8113/44/44/445005, arXiv:1003.4805, https://arxiv.org/abs/1003.4805
  14. N. Iorgov, S. Pakuliak, V. Shadura, Yu. Tykhyy, G. von Gehlen, Spin Operator Matrix Elements in the Superintegrable Chiral Potts Quantum Chain, Journal of Statistical Physics 139 (2010) 743–768, DOI 10.1007/s10955-010-9972-1, https://doi.org/10.1007/s10955-010-9972-1
  15. N. Jones, A. Prakash, P. Fendley, Pivoting through the chiral-clock family, SciPost Physics 18 (2025) 094, DOI 10.21468/SciPostPhys.18.3.094, arXiv:2406.01680, https://arxiv.org/abs/2406.01680
  16. I. Camp, N. G. Jones, Matrix-Product State Skeletons in Onsager-Integrable Quantum Chains, Journal of Statistical Physics 193 (2026), article 73, DOI 10.1007/s10955-026-03637-8, arXiv:2511.07212, https://arxiv.org/abs/2511.07212
  17. H. Zhu, A symmetry formula for correlation functions in the superintegrable chiral Potts spin chain, arXiv:2603.27930 (2026), https://arxiv.org/abs/2603.27930
  18. Z. Zhang, The Yang-Baxter Equation for the Chiral Potts Model and Integrable Parafermions, arXiv:2606.11146 (2026), https://arxiv.org/abs/2606.11146