2D Potts Model Correlation Lengths: Numerical Evidence for $ξ_o = ξ_d$ at $β_t$

We have studied spin-spin correlation functions in the ordered phase of the two-dimensional $q$-state Potts model with $q=10$, 15, and 20 at the first-order transition point $\beta_t$. Through extensive Monte Carlo simulations we obtain strong numerical evidence that the correlation length in the ordered phase agrees with the exactly known and recently numerically confirmed correlation length in the disordered phase: $\xi_o(\beta_t) = \xi_d(\beta_t)$. As a byproduct we find the energy moments in the ordered phase at $\beta_t$ in very good agreement with a recent large $q$-expansion.