We present a compendium of numerical simulation techniques, based on tensor network methods, aiming to address problems of many-body quantum mechanics on a classical computer. The core setting of this anthology are lattice problems in low spatial dimension at finite size, a physical scenario where tensor network methods, both Density Matrix Renormalization Group and beyond, have long proven to be winning strategies. Here we explore in detail the numerical frameworks and methods employed to deal with low-dimension physical setups, from a computational physics perspective. We focus on symmetries and closed-system simulations in arbitrary boundary conditions, while discussing the numerical data structures and linear algebra manipulation routines involved, which form the core libraries of any tensor network code. At a higher level, we put the spotlight on loop-free network geometries, discussing their advantages, and presenting in detail algorithms to simulate low-energy equilibrium states. Accompanied by discussions of data structures, numerical techniques and performance, this anthology serves as a programmer's companion, as well as a self-contained introduction and review of the basic and selected advanced concepts in tensor networks, including examples of their applications.

We investigate the finite-density phase diagram of a non-abelian SU(2) lattice gauge theory in (1+1)-dimensions using tensor network methods. We numerically characterise the phase diagram as a function of the matter filling and of the matterfield coupling, identifying different phases, some of them appearing only at finite densities. For weak matter-field coupling we find a meson BCS liquid phase, which is confirmed by second-order analytical perturbation theory. At unit filling and for strong coupling, the system undergoes a phase transition to a charge density wave of single-site (spin-0) mesons via spontaneous chiral symmetry breaking. At finite densities, the chiral symmetry is restored almost everywhere, and the meson BCS liquid becomes a simple liquid at strong couplings, with the exception of filling two-thirds, where a charge density wave of mesons spreading over neighbouring sites appears. Finally, we identify two tri-critical points between the chiral and the two liquid phases which are compatible with a SU (2) 2 Wess-Zumino-NovikovWitten model. Here we do not perform the continuum limit but we explicitly address the global U (1) charge conservation symmetry.

We study the equilibrium properties of the one-dimensional disordered Bose-Hubbard model by means of a gauge-adaptive tree tensor network variational method suitable for systems with periodic boundary conditions. We compute the superfluid stiffness and superfluid correlations close to the superfluid to glass transition line, obtaining accurate locations of the critical points. By studying the statistics of the exponent of the power-law decay of the correlation, we determine the boundary between the superfluid region and the Bose glass phase in the regime of strong disorder and in the weakly interacting region, not explored numerically before. In the former case our simulations are in agreement with previous Monte Carlo calculations.

We investigate the ground state phase diagram of square ice -a U(1) lattice gauge theory in two spatial dimensions -using gauge invariant tensor network techniques. By correlation function, Wilson loop, and entanglement diagnostics, we characterize its phases and the transitions between them, finding good agreement with previous studies. We study the entanglement properties of string excitations on top of the ground state, and provide direct evidence of the fact that the latter are described by a conformal field theory. Our results pave the way to the application of tensor network methods to confining, two-dimensional lattice gauge theories, to investigate their phase diagrams and low-lying excitations.Contents arXiv:1807.00826v2 [cond-mat.stat-mech]