Thermodynamic description of the Ising antiferromagnet on a triangular lattice with selective dilution by a modified pair-approximation method

The pair-approximation method is modified in order to describe systems with geometrical frustration. The Ising antiferromagnet on a triangular lattice with selective dilution (Kaya-Berker model) is considered and a self-consistent thermodynamic description of this model is obtained. For this purpose, the Gibbs free energy as a function of temperature, concentration of magnetic atoms on the selected sublattice, and external magnetic field is derived. In particular, the phase diagram is constructed and a comparison of different methods is presented. The thermodynamic quantities are discussed in the context of their physical validity, and the improvement in the description introduced by the modified method is emphasized.