*

The classical Ising model on the triangular lattice (we will call it I-3 model below), while simple in the nearest-neighbors (nn) only approximation, becomes increasingly richer and more complex when further interactions are incorporated. However, the studies so far have not been exhaustive, nor have any attempts been made to estimate how realistic are the parameter ranges that generate strong metamagnetism with a large number of magnetization steps. In this study, we identify one such candidate, ErGa₂ , a material known to have one strong magnetization step, albeit some narrow steps below and above cannot be confidently excluded. It has been established, and we can confirm the same computationally, to have an easy axis perpendicular to the triangular Er plane, with a strong anisotropy and with a large magnetic moment of 9.5 μB, making it a perfect implementation of the classical I-3 model. In the first part of the analysis, we present the I-3 model with up to the 3rd nearest neighbors in a range of parameters J2 and J3 ( J 1 in this part is set to 1), and in some cases adding a rather small J 4 in order to reveal new phase otherwise degenerate with some others. The richest phase diagram is, as expected, observed when all interactions are antiferromagnetic (AF). Subsequently, a more realistic case, inspired by RKKY and by our calculations for ErGa₂, where J1, J 2 >0 (antiferromagnetic) and J 3, J 4 <0 (ferromagnetic), is presented. Finally, we present our first-principles calculations of J 1-4 in ErGa₂ and compare the phase diagram in the regime corresponding to the calculated values.