Mixed Ising models with two compensation temperatures

We calculate the free energy of a mixed Ising ferrimagnetic system from a mean field estimate of the total energy and the entropy. By minimizing the free energy we obtain the equilibrium magnetizations and the compensation temperatures. We found that compensation temperatures only appear in metastable systems, and that depending on the crystal field parameter, the system can present none, one or two compensation points.