Constructing the Critical Curve for the Two-Layer Potts Model Using   Cellular Automata
  The critical points of the 3-states two-layer Potts model on square lattice for different interlayer couplings (Kx, Ky,and Kz) are calculated with high precision using probabilistic cellular automata with Glauber algorithm, where Kx and Ky are the nearest-neighbor interactions within each layer in the x and y directions, respectively and Kz is the interlayer coupling. The obtained results are 0.726, 0.807, 0.928, 0.987 and 1.00 for the cases Kz per Kx = 1, 0.5, 0.1, 0.01 and 0.001, respectively. Then, the critical curve has been constructed for this model. 
Yazdan Asgari, Mehrdad Ghaemi