In this study, we explore the hysteresis behavior of a bilayer graphene system consisting of mixed spins \((S = 7/2, \sigma = 1/2)\) in ferrimagnetic interaction. The model considered is based on a modified Blume-Capel Hamiltonian, incorporating intra- and inter-layer exchange interactions, a uniaxial crystal field, and an external magnetic field. The analysis is performed using Monte Carlo simulations based on the Metropolis algorithm. We examined the influence of physical parameters, in particular the exchange interactions \(J_1\) , \(J_2\) , \(J_3\) , the crystal field D and the temperature T, on the hysteretic behavior of the system. The results obtained reveal a strong dependence of the topology of the hysteresis loops on these parameters. For example, increasing the temperature causes a gradual decrease in the surface area of the loops, until they disappear above a critical temperature. These observations are consistent with the results reported in other multilayer magnetic systems and open up prospects for the development of multistate magnetic memory devices.
