Energy Scales and Exponential Speedup in Thermal Tensor Network Simulations
Bin-Bin Chen,Lei Chen,Ziyu Chen,Wei Li,Andreas Weichselbaum
2017-01-01
Abstract:In this work, we exponentially speed up thermal simulations of quantum many-body systems in both one- (1D) and two-dimensional (2D) models using matrix product operators (MPOs). Instead of evolving the density operator {$hat{rho}=e^{-beta hat{H}}$} linearly in inverse temperature $beta equiv 1/T$ as conventional Trotter-Suzuki methods do, we cool down the system {by iteratively projecting the MPO representation of $hat{rho}$ to itself, i.e., doubling $beta$} in the process of imaginary time evolution. This exponential tensor renormalization group (XTRG) scheme, compared to linear evolution schemes, reaches low temperatures much faster, and thus not only saves computational time but also merits better accuracy due to significantly fewer truncation steps. For similar reasons, we also find that the series expansion thermal tensor network (SETTN) approach benefits in both efficiency and precision, from the logarithmic temperature scale setup. For both thermal algorithms, XTRG as well as SETTN, we fully implement non-Abelian and Abelian symmetries to greatly enhance their numerical performance. We employ these cutting-edge techniques for finite temperature simulations to explore low-temperature thermal states of both 1D and 2D Heisenberg models. The entanglement properties, as well as the renormalization group flow of entanglement spectra in MPOs, are discussed, where logarithmic entropies ($sim ln{beta}$) are shown in both spin chains and square lattice models with gapless tower of states.