We describe how to implement the time-dependent variational principle for matrix product states in the thermodynamic limit for nonuniform lattice systems. This is achieved by confining the ...nonuniformity to a (dynamically expandable) finite region with fixed boundary conditions. The suppression of nonphysical quasiparticle reflections from the boundary of the nonuniform region is also discussed. Using this algorithm we study the dynamics of localized excitations in infinite systems, which we illustrate in the case of the spin-1 antiferromagnetic Heisenberg model and the varphi super(4) model.
A natural way to generalize tensor network variational classes to quantum field systems is via a continuous tensor contraction. This approach is first illustrated for the class of quantum field ...states known as continuous matrix-product states (cMPS). As a simple example of the path-integral representation we show that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. A completeness argument is also provided that shows that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity. Beyond this, we obtain a well-behaved field limit of projected entangled-pair states (PEPS) in two dimensions that provide an abstract class of quantum field states with natural symmetries. We demonstrate how symmetries of the physical field state are encoded within the dynamics of an auxiliary field system of one dimension less. In particular, the imposition of Euclidean symmetries on the physical system requires that the auxiliary system involved in the class' definition must be Lorentz-invariant. The physical field states automatically inherit entropy area laws from the PEPS class, and are fully described by the dissipative dynamics of a lower dimensional virtual field system. Our results lie at the intersection many-body physics, quantum field theory and quantum information theory, and facilitate future exchanges of ideas and insights between these disciplines.
We extend the recently introduced continuous matrix product state variational class to the setting of (1+1)-dimensional relativistic quantum field theories. This allows one to overcome the ...difficulties highlighted by Feynman concerning the application of the variational procedure to relativistic theories, and provides a new way to regularize quantum field theories. A fermionic version of the continuous matrix product state is introduced which is manifestly free of fermion doubling and sign problems. We illustrate the power of the formalism by studying the momentum occupation for free massive Dirac fermions and the chiral symmetry breaking in the Gross-Neveu model.
We develop a method based on tensor networks to create localized single-particle excitations on top of strongly correlated quantum spin chains. In analogy to the problem of creating localized Wannier ...modes, this is achieved by optimizing the gauge freedom of momentum excitations on top of matrix product states. The corresponding wave packets propagate almost dispersionlessly. The time-dependent variational principle is used to scatter two such wave packets, and we extract the phase shift from the collision data. We also study reflection and transmission coefficients of a wave packet impinging on an impurity.
An essential primitive in quantum tensor network simulations is the problem of approximating a matrix product state with one of a smaller bond dimension. This problem forms the central bottleneck in ...algorithms for time evolution and for contracting projected entangled pair states. We formulate a tangent-space based variational algorithm to achieve this goal for uniform (infinite) matrix product states. The algorithm exhibits a favourable scaling of the computational cost, and we demonstrate its usefulness by several examples involving the multiplication of a matrix product state with a matrix product operator.
We show that the time-dependent variational principle provides a unifying framework for time-evolution methods and optimization methods in the context of matrix product states. In particular, we ...introduce a new integration scheme for studying time evolution, which can cope with arbitrary Hamiltonians, including those with long-range interactions. Rather than a Suzuki-Trotter splitting of the Hamiltonian, which is the idea behind the adaptive time-dependent density matrix renormalization group method or time-evolving block decimation, our method is based on splitting the projector onto the matrix product state tangent space as it appears in the Dirac-Frenkel time-dependent variational principle. We discuss how the resulting algorithm resembles the density matrix renormalization group (DMRG) algorithm for finding ground states so closely that it can be implemented by changing just a few lines of code and it inherits the same stability and efficiency. In particular, our method is compatible with any Hamiltonian for which ground-state DMRG can be implemented efficiently. In fact, DMRG is obtained as a special case of our scheme for imaginary time evolution with infinite time step.
The efficient evaluation of tensor expressions involving sums over multiple indices is of significant importance to many fields of research, including quantum many-body physics, loop quantum gravity, ...and quantum chemistry. The computational cost of evaluating an expression may depend strongly on the order in which the index sums are evaluated, and determination of the operation-minimizing contraction sequence for a single tensor network (single term, in quantum chemistry) is known to be NP-hard. The current preferred solution is an exhaustive search, using either an iterative depth-first approach with pruning or dynamic programming and memoization, but these approaches are impractical for many of the larger tensor network ansätze encountered in quantum many-body physics. We present a modified search algorithm with enhanced pruning which exhibits a performance increase of several orders of magnitude while still guaranteeing identification of an optimal operation-minimizing contraction sequence for a single tensor network. A reference implementation for matlab, compatible with the ncon() and multienv() network contractors of arXiv:1402.0939 and Evenbly and Pfeifer, Phys. Rev. B 89, 245118 (2014), respectively, is supplied.