Infinite lattice models by expansion with a non-Gaussian initial approximation

Recently, a convergent series employing a non-Gaussian initial approximation was constructed and shown to be an effective computational tool for the finite size lattice models with a polynomial interaction. Here we numerically examine the applicability of the convergent series method to models defined on infinite lattices. The comparison of the convergent series computations and the infinite lattice extrapolations of the Monte Carlo simulations reveals an agreement between two approaches.