Building on a previously introduced block Lanczos method, we demonstrate how to approximate any operator function of the form $\mathsf{Tr}f\left(A\right)$ when the argument $A$ is given as a Hermitian matrix product operator (MPO). This gives access to quantities that, depending on the full spectrum, are difficult to access for standard tensor network techniques, such as the von Neumann entropy and the trace norm of an MPO. We present a modified, more efficient strategy for computing thermal properties of short- or long-range Hamiltonians, and we illustrate the performance of the method with numerical results for the thermal equilibrium states of the Lipkin-Meshkov-Glick and Ising Hamiltonians.

• Received 9 May 2018
• Revised 26 July 2018

DOI:https://doi.org/10.1103/PhysRevB.98.075128