The matrix product state formalism is used to simulate Hamiltonian lattice gauge theories. To this end, we define matrix product state manifolds which are manifestly gauge invariant. As an application, we study 1+1 dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model, and are able to determine very accurately the ground state properties and elementary one-particle excitations in the continuum limit. In particular, a novel particle excitation in the form of a heavy vector boson is uncovered, compatible with the strong coupling expansion in the continuum. We also study full quantum non-equilibrium dynamics by simulating the real-time evolution of the system induced by a quench in the form of a uniform background electric field.
PACS numbers:Gauge theories hold a most prominent place in physics. They appear as effective low energy descriptions at different instances in condensed matter physics and nuclear physics. But far and foremost they lie at the root of our understanding of the four fundamental interactions that are each mediated by the gauge fields corresponding to a particular gauge symmetry. At the perturbative quantum level, this picture translates to the Feynman diagrammatic approach that has produced physical predictions with unlevelled precision, most famously in quantum electrodynamics (QED). However the perturbative approach miserably fails once the interactions become strong. This problem is most pressing for quantum chromodynamics (QCD), where all low energy features like quark confinement, chiral symmetry breaking and mass generation are essentially non-perturbative.Lattice QCD, which is based on Monte Carlo sampling of Wilson's Euclidean lattice version of gauge theories, has historically been by far the most successful method in tackling this strongly coupled regime. Using up a sizable fraction of the global supercomputer time, state of the art calculations have now reached impressive accuracy, for instance in the ab initio determination of the light hadron masses [1]. But in spite of its clear superiority, the lattice Monte Carlo sampling also suffers from a few drawbacks. There is the infamous sign problem that prevents application to systems with large fermionic densities. In addition, the use of Euclidean time, as opposed to real time, presents a serious barrier for the understanding of dynamical non-equilibrium phenomena. Over the last few years there has been a growing experimental and theoretical interest in precisely these elusive regimes, e.g. in the study of heavy ion collisions or early time cosmology.In this letter we study the application of tensor network states (TNS) as a possible complementary approach to the numerical simulation of gauge theories. This is highly relevant as this Hamiltonian method is free from the sign problem and allows for real-time dynamics. As a first application we focus on the massive Schwinger model. For this model the TNS approach has been studied before by Byrnes et al [2] and Bañuls et al [3]. By integrating out the ...