Accurate determination of tensor network state of quantum lattice models in two dimensions.

We have proposed a novel numerical method to calculate accurately physical quantities of the ground state using the tensor network wave function in two dimensions. The tensor network wave function is determined by an iterative projection approach which uses the Trotter-Suzuki decomposition formula of quantum operators and the singular value decomposition of matrix. The norm of the wave function and the expectation value of a physical observable are evaluated by a coarse-grain tensor renormalization group approach. Our method allows a tensor network wave function with a high bond degree of freedom (such as D=8) to be handled accurately and efficiently in the thermodynamic limit. For the Heisenberg model on a honeycomb lattice, our results for the ground state energy and the staggered magnetization agree well with those obtained by the quantum Monte Carlo and other approaches.