A numerical study for two-dimensional spin 1/2 antiferromagnets: A generalization of Entanglement Perturbation Theory to two-dimensional systems

Two-dimensional spin 1/2 antiferromagnetic Heisenberg models are numerically studied using Entanglement Perturbation Theory, where the ground state wave function is described by a product of local matrices defined at every site, and each matrix is optimized variationally to minimize the energy. We first apply this method for the spin 1/2 antiferromagnetic Heisenberg model on the square lattices to reproduce correctly the known ground state energy and the spin structure factors. Then, we study the spin 1/2 antiferromagnetic Heisenberg model on the triangular lattice with spatially anisotropic nearest neighbor couplings, J and J', to establish the ground state phase diagram as a function of J'/ J.