“…We use a self-consistent local field algorithm, [19][20][21] and, at each temperature, the equilibrium state ͓͑m i ͒ , i =1, 2, ... ,N͔ is found by seeking a magnetic configuration in which the torque is zero in all the cells ͑m ជ i ϫ H ជ eff i = 0 for i =1,2, ... ,N͒, where the effective field is H ជ…”