Consider a 2D Ising model on a square lattice Λ Λ with N sites and periodic boundary conditions in both the horizontal and vertical directions, which effectively reduces the topology of the model to a torus. Generally, the horizontal coupling J and the vertical coupling J ∗ ∗ are not equal. With β β = 1 k T and absolute temperature T and the Boltzmann constant k , the partition function