Tensor networks, and in particular Projected Entangled Pair States (PEPS), are a powerful tool for the study of quantum many body physics, thanks to both their built-in ability of classifying and studying symmetries, and the efficient numerical calculations they allow. In this work, we introduce a way to extend the set of symmetric PEPS in order to include local gauge invariance and investigate lattice gauge theories with fermionic matter. To this purpose, we provide as a case study and first example, the construction of a fermionic PEPS, based on Gaussian schemes, invariant under both global and local U (1) gauge transformations. The obtained states correspond to a truncated U (1) lattice gauge theory in 2 + 1 dimensions, involving both the gauge field and fermionic matter. For the global symmetry (pure fermionic) case, these PEPS can be studied in terms of spinless fermions subject to a p-wave superconducting pairing. For the local symmetry (fermions and gauge fields) case, we find confined and deconfined phases in the pure gauge limit, and we discuss the screening properties of the phases arising in the presence of dynamical matter.arXiv:1507.08837v2 [quant-ph]

We propose a tensor network encoding the set of all eigenstates of a fully many-body localized system in one dimension. Our construction, conceptually based on the ansatz introduced in Phys. Rev. B 94, 041116(R) (2016), is built from two layers of unitary matrices which act on blocks of l contiguous sites. We argue that this yields an exponential reduction in computational time and memory requirement as compared to all previous approaches for finding a representation of the complete eigenspectrum of large many-body localized systems with a given accuracy. Concretely, we optimize the unitaries by minimizing the magnitude of the commutator of the approximate integrals of motion and the Hamiltonian, which can be done in a local fashion. This further reduces the computational complexity of the tensor networks arising in the minimization process compared to previous work. We test the accuracy of our method by comparing the approximate energy spectrum to exact diagonalization results for the random-field Heisenberg model on 16 sites. We find that the technique is highly accurate deep in the localized regime and maintains a surprising degree of accuracy in predicting certain local quantities even in the vicinity of the predicted dynamical phase transition. To demonstrate the power of our technique, we study a system of 72 sites, and we are able to see clear signatures of the phase transition. Our work opens a new avenue to study properties of the many-body localization transition in large systems.

We show that projected entangled-pair states (PEPS) in two spatial dimensions can describe chiral topological states by explicitly constructing a family of such states with a nontrivial Chern number. They are ground states of two different kinds of free-fermion Hamiltonians: (i) local and gapless; (ii) gapped, but with hopping amplitudes that decay according to a power law. We derive general conditions on topological free-fermionic projected entangled-pair states that show that they cannot correspond to exact ground states of gapped, local parent Hamiltonians and provide numerical evidence demonstrating that they can nevertheless approximate well the physical properties of topological insulators with local Hamiltonians at arbitrary temperatures.

Lessons from Anderson localization highlight the importance of dimensionality of real space for localization due to disorder. More recently, studies of many-body localization have focussed on the phenomenon in one dimension using techniques of exact diagonalization and tensor networks. On the other hand, experiments in two dimensions have provided concrete results going beyond the previously numerically accessible limits while posing several challenging questions. We present the first large-scale numerical examination of a disordered Bose-Hubbard model in two dimensions realized in cold atoms, which shows entanglement based signatures of many-body localization. By generalizing a low-depth quantum circuit to two dimensions we approximate eigenstates in the experimental parameter regimes for large systems, which is beyond the scope of exact diagonalization. A careful analysis of the eigenstate entanglement structure provides an indication of the putative phase transition marked by a peak in the fluctuations of entanglement entropy in a parameter range consistent with experiments.Many-body localization (MBL) is a paradigm shift in out-of-equilibrium quantum matter. This novel phase of matter is characterized by the absence of thermalization 1-6 . An MBL system retains a memory of its initial state and displays only a logarithmic growth of entanglement following quantum quenches 7 . By localizing the excitations, MBL can also protect certain forms of topological and symmetry-breaking orders in excited states and provides an opportunity to process quantum information in a system driven far from equilibrium 8-11 . The quantum phase transition separating the MBL and thermal phases poses a major challenge for developing a theory of dynamical critical phenomena described by manybody entanglement in highly excited states 12-18 .It is well-known that dimensionality of real space affects single particle Anderson localization where in one and two dimensions (without spin-orbit coupling and broken time-reversal symmetry), the entire spectrum of single particle eigenstates is localized for arbitrarily weak disorder 19 . Although the defining properties of MBL in one dimension are firmly established both theoretically and experimentally 20,21 , the existence of the phenomenon in two and higher dimensions is much debated 22-28 . Experiments in cold atoms measuring local and global equilibration have shown the persistence of quantum memory for long times providing indications of an MBL-like phase in two dimensions 28-30 . On the other hand, theoretical criteria suggest that the lifetime of local memory is finite, albeit extremely long 22,23 .In this article, we evaluate the eigenstates of bosons hopping in a disordered lattice in two dimensions with on-site interactions. We generalize the tensor network method developed earlier for one dimensional systems to approximate the eigenstates in the localized regime 31 .Because the system sizes accessible by our method are much larger than those that can be currently achieved with prior ...

Over the last years, Projected Entangled Pair States have demonstrated great power for the study of many body systems, as they naturally describe ground states of gapped many body Hamiltonians, and suggest a constructive way to encode and classify their symmetries. The PEPS study is not only limited to global symmetries, but has also been extended and applied for local symmetries, allowing to use them for the description of states in lattice gauge theories. In this paper we discuss PEPS with a local, SU(2) gauge symmetry, and demonstrate the use of PEPS features and techniques for the study of a simple family of many body states with a non-Abelian gauge symmetry. We present, in particular, the construction of fermionic PEPS able to describe both two-color fermionic matter and the degrees of freedom of an SU(2) gauge field with a suitable truncation. CONTENTS