In this article we discuss the convergence of the variational cumulant expansion (VCE) in the Ising model and lattice model. We find that in critical region VCE results do not converge so quickly as elsewhere and the analyticity of thermo-dynamical functions and the strong correlation in the critical region are responsible for this fact.