“…Experimentally, several works have been published on the properties of the magnetic bilayer such as FePt/Fe [11], FM/TbFe (FM=Fe, Py-permalloy, FeCo) [12], SiGe, SiC and GeC [13], Ni/Au [14], Fe/Ni and Fe/Co [15]. On the other hand, theoretically, the equilibrium behaviors of bilayer Ising systems consisting of single ions, namely spin-1/2, spin-1 etc., have been widely studied by utilizing the wellknown theoretical methods used in equilibrium statistical physics, such as the mean field theory (MFT) [16,17], the effective field theory (EFT) [18,19], the exact recursion equation on the Bethe lattice [6,7], the renormalization-group (RG) method [20,21], and Monte Carlo (MC) simulations [22,23].…”