No Arabic abstract
Just as matrix product states represent ground states of one-dimensional quantum spin systems faithfully, continuous matrix product states (cMPS) provide faithful representations of the vacuum of interacting field theories in one spatial dimension. Unlike the quantum spin case however, for which the density matrix renormalization group and related matrix product state algorithms provide robust algorithms for optimizing the variational states, the optimization of cMPS for systems with inhomogeneous external potentials has been problematic. We resolve this problem by constructing a piecewise linear parameterization of the underlying matrix-valued functions, which enables the calculation of the exact reduced density matrices everywhere in the system by high-order Taylor expansions. This turns the variational cMPS problem into a variational algorithm from which both the energy and its backwards derivative can be calculated exactly and at a cost that scales as the cube of the bond dimension. We illustrate this by finding ground states of interacting bosons in external potentials, and by calculating boundary or Casimir energy corrections of continuous many-body systems with open boundary conditions.
We combine the Density Matrix Renormalization Group (DMRG) with Matrix Product State tangent space concepts to construct a variational algorithm for finding ground states of one dimensional quantum lattices in the thermodynamic limit. A careful comparison of this variational uniform Matrix Product State algorithm (VUMPS) with infinite Density Matrix Renormalization Group (IDMRG) and with infinite Time Evolving Block Decimation (ITEBD) reveals substantial gains in convergence speed and precision. We also demonstrate that VUMPS works very efficiently for Hamiltonians with long range interactions and also for the simulation of two dimensional models on infinite cylinders. The new algorithm can be conveniently implemented as an extension of an already existing DMRG implementation.
For the past twenty years, Matrix Product States (MPS) have been widely used in solid state physics to approximate the ground state of one-dimensional spin chains. In this paper, we study homogeneous MPS (hMPS), or MPS constructed via site-independent tensors and a boundary condition. Exploiting a connection with the theory of matrix algebras, we derive two structural properties shared by all hMPS, namely: a) there exist local operators which annihilate all hMPS of a given bond dimension; and b) there exist local operators which, when applied over any hMPS of a given bond dimension, decouple (cut) the particles where they act from the spin chain while at the same time join (glue) the two loose ends back again into a hMPS. Armed with these tools, we show how to systematically derive `bond dimension witnesses, or 2-local operators whose expectation value allows us to lower bound the bond dimension of the underlying hMPS. We extend some of these results to the ansatz of Projected Entangled Pairs States (PEPS). As a bonus, we use our insight on the structure of hMPS to: a) derive some theoretical limitations on the use of hMPS and hPEPS for ground state energy computations; b) show how to decrease the complexity and boost the speed of convergence of the semidefinite programming hierarchies described in [Phys. Rev. Lett. 115, 020501 (2015)] for the characterization of finite-dimensional quantum correlations.
We define matrix product states in the continuum limit, without any reference to an underlying lattice parameter. This allows to extend the density matrix renormalization group and variational matrix product state formalism to quantum field theories and continuum models in 1 spatial dimension. We illustrate our procedure with the Lieb-Liniger model.
A generic method to investigate many-body continuous-variable systems is pedagogically presented. It is based on the notion of matrix product states (so-called MPS) and the algorithms thereof. The method is quite versatile and can be applied to a wide variety of situations. As a first test, we show how it provides reliable results in the computation of fundamental properties of a chain of quantum harmonic oscillators achieving off-critical and critical relative errors of the order of 10^(-8) and 10^(-4) respectively. Next, we use it to study the ground state properties of the quantum rotor model in one spatial dimension, a model that can be mapped to the Mott insulator limit of the 1-dimensional Bose-Hubbard model. At the quantum critical point, the central charge associated to the underlying conformal field theory can be computed with good accuracy by measuring the finite-size corrections of the ground state energy. Examples of MPS-computations both in the finite-size regime and in the thermodynamic limit are given. The precision of our results are found to be comparable to those previously encountered in the MPS studies of, for instance, quantum spin chains. Finally, we present a spin-off application: an iterative technique to efficiently get numerical solutions of partial differential equations of many variables. We illustrate this technique by solving Poisson-like equations with precisions of the order of 10^(-7).
Any quantum process is represented by a sequence of quantum channels. We consider ergodic processes, obtained by sampling channel valued random variables along the trajectories of an ergodic dynamical system. Examples of such processes include the effect of repeated application of a fixed quantum channel perturbed by arbitrary correlated noise, or a sequence of channels drawn independently and identically from an ensemble. Under natural irreducibility conditions, we obtain a theorem showing that the state of a system evolving by such a process converges exponentially fast to an ergodic sequence of states depending on the process, but independent of the initial state of the system. As an application, we describe the thermodynamic limit of ergodic matrix product states and prove that the 2-point correlations of local observables in such states decay exponentially with their distance in the bulk. Further applications and physical implications of our results are discussed in the companion paper [11].