Application of the Monte Carlo method for evaluating the Niemeijer-van Leeuwen cumulant expansion

We show that the large cell size problem of the real space renormalization can be handled effectively by Monte Carlo methods. As a demonstration, the second-order cumulant expansion is calculated for the three-dimensional simple cubic Ising model, using a 3 × 3 × 3 cell.