When one modeling the magnetic systems, one is surely faced with the question of optimizing the code and the used model. Since the dipole-dipole interaction is long-range, it is computationally expensive to simulate the magnetic systems without a cutoff radius, and many researchers cut off the dipole-dipole interaction [1]}, [2]}, [3]}, [4]}, [5]}. The question naturally becomes, how much can the dipole interaction be cutoff, so that it does not affect the basic thermodynamic characteristics of the model? This is the question that was asked in [6]}. Rougemaille in his work write, that macrospins in artificial spin ice are coupled via the magnetostatic interaction, which extends beyond nearest neighbors and that including this dipole interaction to the model remove the ground state degeneracy [7]}. His work shows the differences between short-range and long-range interaction models. This problem was also considered in [8]}, [9]}. Also our scientific group in previous works describes this problem in [6]}, [11]}. It was shown [6]} that the dipole-dipole interaction cutoff to the nearest neighbors (\(R_1\)  in Fig. REF b, includes 4 nearest neighbors) leads to a loss of the low-temperature phase transition (see Fig. REF ).
<FIGURE>