For the two-dimensional 4-state Potts model — the marginal member of the Potts family, whose continuous transition is dressed with logarithmic corrections — the published determinations of the universal critical-amplitude ratios are mutually inconsistent: the two-kink form-factor predictions give Γ+/ΓL=4.013\Gamma_+/\Gamma_L = 4.013, ΓT/ΓL=0.129\Gamma_T/\Gamma_L = 0.129, RC+=0.0204R_C^+ = 0.0204 and RC−=0.00508R_C^- = 0.00508, while Monte Carlo and series-expansion estimates for the Potts model itself give Γ+/ΓL=3.14\Gamma_+/\Gamma_L = 3.14–6.936.93 and ΓT/ΓL=0.148\Gamma_T/\Gamma_L = 0.148–0.1670.167, with a Baxter–Wu Monte Carlo simulation in the same universality class yielding Γ+/ΓL=3.9(1)\Gamma_+/\Gamma_L = 3.9(1). The problem is to determine these universal ratios definitively, at the few-per-cent level, either by an analytical computation with controlled accuracy beyond the two-kink approximation or by numerical analyses that resolve the background and logarithmic-correction systematics, and thereby to settle whether the two-kink form-factor approximation fails at the marginal point q=4q=4. As of 2026 the conflict remains unadjudicated in the literature — in particular the ΓT/ΓL\Gamma_T/\Gamma_L and RC+R_C^+ discrepancies have never been explained — so the problem is open.

Background

The two-dimensional qq-state Potts model assigns to each lattice site a spin sis_i taking qq values, with nearest-neighbour Hamiltonian H=−J∑⟨ij⟩δsi,sjH=-J\sum_{\langle ij\rangle}\delta_{s_i,s_j}; on the square lattice it is exactly solvable at criticality and has a continuous transition for q≤4q\le 4 and a first-order transition for q>4q>4. Near the critical temperature TcT_c the singular parts of thermodynamic quantities behave as M≃B(−τ)βM\simeq B(-\tau)^\beta, χ±≃Γ±∣τ∣−γ\chi_\pm\simeq\Gamma_\pm|\tau|^{-\gamma}, C±≃(A±/α)∣τ∣−αC_\pm\simeq (A^\pm/\alpha)|\tau|^{-\alpha}, with reduced temperature τ=(T−Tc)/Tc\tau=(T-T_c)/T_c; below TcT_c the susceptibility splits into longitudinal (ΓL≡Γ−\Gamma_L\equiv\Gamma_-) and transverse (ΓT\Gamma_T) parts. The amplitudes B,Γ±,A±B,\Gamma_\pm,A^\pm are non-universal, but combinations such as Γ+/ΓL\Gamma_+/\Gamma_L, ΓT/ΓL\Gamma_T/\Gamma_L and RC±=A±Γ±/B2R_C^\pm=A^\pm\Gamma_\pm/B^2 (with the normalization conventions used in the cited works, e.g. RC−=αA−Γ−/B2R_C^-=\alpha A^-\Gamma^-/B^2 in Shchur–Janke) are universal numbers characterizing the universality class, on the same footing as critical exponents. The q=4q=4 case is marginal: a marginally irrelevant operator (an operator whose coupling dies out only logarithmically along the renormalization-group flow) generates multiplicative and additive logarithmic corrections (e.g. χ≃Γ∣τ∣−7/6∣ln⁡∣τ∣∣3/4\chi\simeq\Gamma|\tau|^{-7/6}|\ln|\tau||^{3/4}), together with confluent corrections such as ∣τ∣2/3|\tau|^{2/3}, which cancel in the ratios in principle but complicate their numerical extraction. On the analytical side, Delfino and Cardy (1998) computed amplitude ratios in the scaling (continuum, massive) field theory using form factors — matrix elements of local operators between multiparticle states — of the exact factorized scattering theory of Chim and Zamolodchikov for the qq-state Potts model (two-dimensional integrability makes scattering elastic and factorized into two-particle amplitudes), retaining only one- and two-kink intermediate states (the kinks interpolate the qq degenerate vacua below TcT_c; for 3<q≤43<q\le4 a kink bound state must also be included). This two-kink (two-particle) approximation is essentially exact for q=2q=2 (at Ising it gives Γ+/Γ−=37.699\Gamma_+/\Gamma_-=37.699 against the exact 37.6936…37.6936\ldots) and works to about 1% for q=3q=3, but at q=4q=4 it has no controlled small parameter: it already reproduces the central charge — the number characterizing the conformal field theory of the critical point — as c=0.985c=0.985 instead of the exact c=1c=1. Delfino, Barkema and Cardy (2000) extended the computation to ΓT/ΓL\Gamma_T/\Gamma_L (claimed accuracy ∼1%\sim1\%) and noted that the results for Γ+/ΓL\Gamma_+/\Gamma_L and RCR_C are sensitive to the relative normalization of the order- and disorder-operator form factors — a specifically delicate point at q=3,4q=3,4. Numerically, Monte Carlo (MC) simulations and series expansions (SE) must disentangle the leading amplitudes from logarithmic corrections, confluent singularities and large regular background terms (notably in the high-temperature susceptibility), and the published estimates disagree well beyond their error bars. The 4-state Potts, Baxter–Wu (three-spin interaction on the triangular lattice) and a particular coupling of the Ashkin–Teller model (two coupled Ising spins per site, with a line of critical points) are believed to belong to the same universality class, so their universal ratios must coincide — the Baxter–Wu model having the practical advantage of no logarithmic corrections — which provides an independent cross-check and an alternative numerical route.

Problem Statement

Determine the universal critical-amplitude ratios of the two-dimensional 4-state Potts universality class — at minimum the susceptibility ratio Γ+/ΓL\Gamma_+/\Gamma_L (with ΓL≡Γ−\Gamma_L\equiv\Gamma_-), the ordered-phase transverse-to-longitudinal ratio ΓT/ΓL\Gamma_T/\Gamma_L, and the combinations RC±=A±Γ±/B2R_C^\pm=A^\pm\Gamma_\pm/B^2 in explicitly stated normalization conventions (the cited works differ, e.g. in whether RC±R_C^\pm carries the extra factor α\alpha) — with errors at the few-per-cent level, i.e. small enough to discriminate between the conflicting published values, and with an account that settles the standing discrepancy: either (a) an analytical determination with controlled accuracy, going beyond or rigorously bounding the error of the two-kink form-factor approximation, which establishes whether that approximation is insufficient at the marginal point q=4q=4 (where the two-kink central charge is 0.985 instead of the exact 1 and the order/disorder form-factor normalization is uncontrolled), or (b) high-precision numerical determinations (Monte Carlo, series expansions, transfer matrices, or other) that resolve the systematic issues — logarithmic corrections, confluent singularities and regular background terms, especially in the high-temperature susceptibility — and reconcile the mutually inconsistent published estimates, including specifically the unexplained ΓT/ΓL\Gamma_T/\Gamma_L tension (0.129 vs 0.148–0.167) and the RC+R_C^+ tension (0.0204 vs ≈0.033), which the background-term explanation of the Γ+\Gamma_+ discrepancy proposed by Shchur and Janke [7] cannot address. The single overarching objective is a mutually consistent set of values for this universality class, determined to the few-per-cent accuracy stated above; a single determination by either route (a) or (b) suffices to complete the task, and agreement across its realizations (4-state Potts on the square lattice, Baxter–Wu model, Ashkin–Teller model at the Potts point) and between independent methods is corroboration to be reported where available rather than a separately required deliverable.

References

  1. B. Berche, P. Butera, W. Janke, L. Shchur, 'Universal ratios of critical amplitudes in the Potts model universality class', Comput. Phys. Commun. 180 (2009) 493-496. DOI: 10.1016/j.cpc.2008.12.022; arXiv:0809.5045; https://arxiv.org/abs/0809.5045
  2. G. Delfino, J. Cardy, 'Universal amplitude ratios in the two-dimensional q-state Potts model and percolation from quantum field theory', Nucl. Phys. B 519 (1998) 551-578. DOI: 10.1016/S0550-3213(98)00144-8; arXiv:hep-th/9712111; https://arxiv.org/abs/hep-th/9712111
  3. G. Delfino, G. T. Barkema, J. L. Cardy, 'Susceptibility amplitude ratios in the two-dimensional Potts model and percolation', Nucl. Phys. B 565 (2000) 521-534. DOI: 10.1016/S0550-3213(99)00629-X; arXiv:cond-mat/9908453; https://arxiv.org/abs/cond-mat/9908453
  4. M. Caselle, R. Tateo, S. Vinti, 'Universal amplitude ratios in the 2D four-state Potts model', Nucl. Phys. B 562 (1999) 549-566. DOI: 10.1016/S0550-3213(99)00521-0; arXiv:cond-mat/9902146; https://arxiv.org/abs/cond-mat/9902146
  5. I. G. Enting, A. J. Guttmann, 'Susceptibility amplitudes for the three- and four-state Potts models', Physica A 321 (2003) 90-107. DOI: 10.1016/S0378-4371(02)01790-9; https://doi.org/10.1016/S0378-4371(02)01790-9
  6. L. N. Shchur, B. Berche, P. Butera, 'Numerical revision of the universal amplitude ratios for the two-dimensional 4-state Potts model', Nucl. Phys. B 811 (2009) 491-518. DOI: 10.1016/j.nuclphysb.2008.10.024; arXiv:0809.4553; https://arxiv.org/abs/0809.4553
  7. L. N. Shchur, W. Janke, 'Critical amplitude ratios of the Baxter-Wu model', Nucl. Phys. B 840 (2010) 491-512. DOI: 10.1016/j.nuclphysb.2010.07.016; arXiv:1007.1838; https://arxiv.org/abs/1007.1838
  8. G. Delfino, P. Grinza, 'Universal ratios along a line of critical points. The Ashkin-Teller model', Nucl. Phys. B 682 (2004) 521-550. DOI: 10.1016/j.nuclphysb.2004.01.007; arXiv:hep-th/0309129; https://arxiv.org/abs/hep-th/0309129
  9. G. Delfino, 'Fields, particles and universality in two dimensions', Ann. Phys. 360 (2015) 477-519. DOI: 10.1016/j.aop.2015.05.020; arXiv:1502.05538; https://arxiv.org/abs/1502.05538
  10. B. Berche, P. Butera, L. N. Shchur, 'The two-dimensional 4-state Potts model in a magnetic field', J. Phys. A 46 (2013) 095001. DOI: 10.1088/1751-8113/46/9/095001; arXiv:1301.5219; https://arxiv.org/abs/1301.5219
  11. J. Salas, A. D. Sokal, 'Logarithmic Corrections and Finite-Size Scaling in the Two-Dimensional 4-State Potts Model', J. Stat. Phys. 88 (1997) 567-615. DOI: 10.1023/B:JOSS.0000015164.98296.85; arXiv:hep-lat/9607030; https://arxiv.org/abs/hep-lat/9607030
  12. L. Chim, A. B. Zamolodchikov, 'Integrable field theory of the q-state Potts model with 0<q<4', Int. J. Mod. Phys. A 7 (1992) 5317-5335. DOI: 10.1142/S0217751X9200243X; https://doi.org/10.1142/S0217751X9200243X
  13. A. Raju, C. B. Clement, L. X. Hayden, J. P. Kent-Dobias, D. B. Liarte, D. Z. Rocklin, J. P. Sethna, 'Normal Form for Renormalization Groups', Phys. Rev. X 9 (2019) 021014. DOI: 10.1103/PhysRevX.9.021014; arXiv:1706.00137; https://arxiv.org/abs/1706.00137