Search

We present a numerical scheme for efficiently extracting the higher-order moments and cumulants of various operators on spin systems represented as tensor product states, for both finite and infinite systems, and present several applications for such quantities. For example, the second cumulant of the energy of a state, (ΔH²), gives a straightforward method to check the convergence of numerical ground-state approximation algorithms. Additionally, we discuss the use of moments and cumulants in the study of phase transitions. Of particular interest is the application of our method to calculate the so-called Binder cumulant, which we use to detect critical points and study the critical exponent of the correlation length with only small finite numerical calculations. We apply these methods to study the behavior of a family of one-dimensional models (the transverse Ising model, the spin-1 Ising model, and the spin-1 Ising model in a crystal field), as well as the two-dimensional Ising model on a square lattice. Our results show that in one dimension, cumulant-based methods can produce precise estimates of the critical points at a low computational cost. The method shows promise for two-dimensional systems as well. [ABSTRACT FROM AUTHOR]

Using exact diagonalization and tensor network techniques, we compute the gap for the Affleck-Kennedy-Lieb-Tasaki (AKLT) Hamiltonian in one and two spatial dimensions. Tensor network methods are used to extract physical properties directly in the thermodynamic limit, and we support these results using finite-size scalings from exact diagonalization. Studying the AKLT Hamiltonian perturbed by an external field, we show how to obtain an accurate value of the gap of the original AKLT Hamiltonian from the field value at which the ground state verifies e0 < 0, which is a quantum critical point. With the tensor network renormalization group methods we provide direct evidence of a finite gap in the thermodynamic limit for the AKLT models in the one-dimensional chain and two-dimensional hexagonal and square lattices. This method can be applied generally to Hamiltonians with rotational symmetry, and we also show results beyond the AKLT model. [ABSTRACT FROM AUTHOR]