Variational method and duality in the 2D square Potts model

The ferromagnetic q-states Potts model on a square lattice is analysed, for q > 4, through an elaborate version of the operatorial variational method. In the variational approach proposed in this paper, the duality relations are exactly satisfied, involving at a more fundamental level, a duality relationship between variational parameters. Besides some exact predictions, the approach is very effective in the numerical estimates over a wide range of temperature and can be systematically improved.