Efficient approximation for global functions of matrix product operators

Building on a previously introduced block Lanczos method, we demonstrate how to approximate any operator function of the form Tr f (A ) 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.