“…From the experimental point of view, single layer, double layer and few (3 to 10) layer honeycomb structures are classified as three different types of 2D crystals, and thin film limit is reached for thicker systems [7]. In this regard, investigation of magnetic properties of graphenelike multilayers gained particular attention, and a wide variety of such systems have been successfully modeled within the framework of Ising model and its variants [19,20,21,22,23]. For instance, using the effective field theory (EFT) formalism, Jiang and coworkers [21] investigated the magnetic properties such as magnetization and the magnetic susceptibility of a nanographene bilayer.…”