In this paper, the Wang-Landau method, one of the modern Monte Carlo methods was used. A detailed description of the method can be found both in the original work by Wang [1]}, and in our previous work [2]}.
