Construction of phase diagrams in anisotropic nanoscaled Ising thin films with diluted surfaces using cellular automata approach

The phase diagrams of the anisotropic three-layer Ising model on the square lattice with diluted surfaces have been constructed with high precision, using the probabilistic cellular automata with the Glauber algorithm. The thermal variation of magnetization is calculated for different values of nearest neighbor couplings: namely J ≠  Js  ≠  J1, where J and Js  are the nearest neighbor  couplings within inner-layer and surface-layer, respectively, and each magnetic site in the surface-layer is coupled with the nearest site in the inner-layer via the exchange interaction J1. In the case of the antiferromagnetic coupling between layers, simulation results show existence of the compensation point in the phase diagram.