The theory of entanglement has provided a new language to describe interacting quantum many body systems, opening up new possibilities in determining static and dynamical properties of such systems. This proposal aims at using the theory of entanglement and more specifically the theory of tensor networks to simulate gauge theories. The big advantage of this method is the fact that tensor networks do not suffer from the sign problem, thereby giving access to properties and quantitites not accessible by quantum Monte Carlo methods. It also provides a new ways of constructing order parameters for topologically ordered phases and of characterizing elementary excitations. We will focus on studying gauge theories in 2+1 dimensions, which is a territory as yet almost untouched upon by tensor network methods, and open up the way of simulating gauge theories in 3+1 dimensions.