We propose a subvolume method to study the $\theta$ dependence of free energy density four-dimensional SU($N$) Yang-Mills theory on lattice. As an attempt, is first applied SU(2) at $T=1.2\,T_c$ understand systematics method. then proceed calculation vacuum and obtain qualitatively different from high temperature case. The numerical results combined with theoretical requirements provide evidenc...