Optimizing compensation in Ising spin-$1/2$ site-diluted trilayered square ferrimagnets: Insights from Monte Carlo Simulations

This study employs Monte Carlo simulations using the Metropolis algorithm to investigate the thermomagnetic properties of a site-diluted trilayered ferrimagnetic system with spin-$1/2$ particles placed on square monolayers. The system is composed of two different types of theoretical atoms, labeled $A$ and $B$, resulting in two distinct configurations, $ABA$ and $AAB$. Atoms of the same type $(A-A$ and $B-B)$ exhibit ferromagnetic interactions, while unlike atoms $(A-B)$ display antiferromagnetic interactions. The $A$-layers are randomly site-diluted with non-magnetic atoms and the percentage of dilution is varied from 5\% to 45\%. Depending upon the concentration of these impurities, both the compensation and critical points shift, leading to different equilibrium ferrimagnetic behaviors. The phase area in the Hamiltonian parameter space, \textit{without} compensation, scales according to the relation: $\ln |A(\rho)/\tilde{A}|=ae^{b\rho}$ . Attention is paid to the mathematical relationships between the compensation temperature and the concentration of nonmagnetic impurities. Suggested mathematical formulas reveal how the threshold impurity concentration is linked to the Hamiltonian parameters.