“…This flip will then happen for 2J z = (q − q )J xy , where again q and q are the numbers of xy neighbors bonded to the flipping spin that are, respectively, +1 and −1. The possible values are (q,q ) = (4,0), (3,1), giving the threshold values of J z /J xy = 1,2, calculationally seen in the bottom panels of Figs. 2 and 3.…”