Magnetic phase diagram of the Ising model with the long-range RKKY interaction

The standard Metropolis algorithm and the parallel tempering method are used to examine magnetization processes in the Ising model with the long-range RKKY interaction on the Shastry-Sutherland lattice. It is shown that the Ising model with RKKY interaction exhibits, depending on the value of the Fermi wave vector k(F), the reach spectrum of magnetic solutions, which is manifested in the appearance of new magnetization plateaus on the magnetization curve. In particular, we have found the following set of individual magnetization plateaus with fractional magnetization m / m(s) = 1 / 18, 1/9, 1/8, 1/5, 1/4, 1/3, 3/8, 5/12, 1/2, 3/5, 2/3, which for different values of k(F) form various sequences of plateaus, changing from very complex, appearing near the point k(F) = 2 pi / 1.2, to very simple appearing away this point. Since the change of k(F) can be induced by doping (the substitution of rare-earth ion by other magnetic ion that introduces the additional electrons to the conduction band) the model is able to predict the complete sequences of magnetization plateaus, which could appear in the tetraboride solid solutions.