LAUSR

Abstract Motivated by the recent investigations on the fullerene nanostructures, we study the critical behavior of the ferromagnetic classical Ising model on fullerene-like systems in the absence of any magnetic field by the Monte Carlo simulation. Since the traditional Metropolis algorithm of filliping is often time-consuming near the reduced critical temperature, a novel technique named Pika is introduced which solves this problem. Since the magnetic properties of fullerene-like materials are greatly affected by the particle size and its structure, the effect of various parameters, including system size, spin value, and vacancy defects on the thermal behavior of the systems, are studied. Based on the results, the reduced critical temperature of the system is proportional to its spin values. Furthermore, the reduced critical temperature strongly depends on the number of vacancies in the system. Increment in the number of vacancies in the system near the reduced critical temperature causes a reduction in the heat capacity value. Moreover, in the temperature T = 0, the total magnetization per site reduces when increasing the number of vacancies.