Multirange Ising model on the square lattice

We study the Ising model on ${\mathbb{Z}}^2$ and show, via numerical simulation, that allowing interactions between spins separated by distances $1$ and $m$ (two ranges), the critical temperature, $T_c(m)$, converges monotonically to the critical temperature of the Ising model on ${\mathbb{Z}}^4$ as $m\to \mathrm{\infty }$. Only interactions between spins located in directions parallel to each coordinate axis are considered. We also simulated the model with interactions between spins at distances of $1$, $m$ and $u$ (three ranges), with $u$ a multiple of $m$; in this case our results indicate that $T_c(m,u)$ converges to the critical temperature of the model on ${\mathbb{Z}}^6$. For percolation, analogous results were proven for the critical probability $p_c$ [B. N. B. de Lima, R. P. Sanchis and R. W. C. Silva, Stochastic Process. Appl. {\bf 121}, 2043 (2011)].