  In this work, Monte Carlo simulation was implemented to study the phase
diagrams and magnetic properties of a mixed spin (1-1/2-1) three-layer system
with a cubic structure. In particular, the effects of the exchange and coupling
parameters on the system compensation temperature were explored. More
importantly, the effect of the crystal field on the system compensation,
critical temperatures, and magnetic and thermodynamic properties of the system
have been explored in detail. The crystal field was found to have a profound
impact on the system compensation and critical temperatures. Our results
suggest that the presence of the crystal field will cause the compensation
temperature to exist for a short range of negative crystal field values.
Furthermore, the threshold values of the crystal field were found to determine
the change in both the compensation and critical temperatures with the
variation of the crystal field. The double-peak phenomena confirmed these
results in the susceptibility and specific heat curves. Finally, the system
magnetic hysteresis loops were studied for different temperatures and crystal
field values, revealing the presence of two-step hysteresis loops for several
values of these parameters.
