Ground state approximation for strongly interacting systems in arbitrary dimension

We introduce a variational method for the approximation of ground states of strongly interacting spin systems in arbitrary geometries and spatial dimensions. The approach is based on weighted graph states and superpositions thereof. These states allow for the efficient computation of all local observables (e.g. energy) and include states with diverging correlation length and unbounded multi-particle entanglement. As a demonstration we apply our approach to the Ising model on 1D, 2D and 3D square-lattices. We also present generalizations to higher spins and continuous-variable systems, which allows for the investigation of lattice field theories.