Exponential Thermal Tensor Network Approach for Quantum Lattice Models

We speed up thermal simulations of quantum many-body systems in both one- (1D) and two-dimensional (2D) models in an exponential way by iteratively projecting the thermal density matrixρ^=e−βH^onto itself. We refer to this scheme of doublingβin each step of the imaginary time evolution as the expone...

Full description

Saved in:
Bibliographic Details
Published in:Physical review. X 2018-09, Vol.8 (3), p.031082, Article 031082
Main Authors: Chen, Bin-Bin, Chen, Lei, Chen, Ziyu, Li, Wei, Weichselbaum, Andreas
Format: Article
Language:English
Subjects:
Accuracy
Chains
Computational efficiency
Computing time
CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY
Critical temperature
Density
Heisenberg theory
High temperature superconductors
Iterative methods
Low temperature
Magnetic properties
Mathematical analysis
Matrices (mathematics)
Phase transitions
Quantum entanglement
Quantum theory
Series expansion
Simulation
Statistical models
Strong interactions (field theory)
Temperature
Temperature scales
Tensors
Thermal simulation
Thermodynamic properties
Transition points
Two dimensional materials
Two dimensional models
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:We speed up thermal simulations of quantum many-body systems in both one- (1D) and two-dimensional (2D) models in an exponential way by iteratively projecting the thermal density matrixρ^=e−βH^onto itself. We refer to this scheme of doublingβin each step of the imaginary time evolution as the exponential tensor renormalization group (XTRG). This approach is in stark contrast to conventional Trotter-Suzuki-type methods which evolveρ^on a linear quasicontinuous grid in inverse temperatureβ≡1/T. As an aside, the large steps in XTRG allow one to swiftly jump across finite-temperature phase transitions, i.e., without the need to resolve each singularly expensive phase-transition point right away, e.g., when interested in low-energy behavior. A fine temperature resolution can be obtained, nevertheless, by using interleaved temperature grids. In general, XTRG can reach low temperatures exponentially fast 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 approach benefits in both efficiency and precision, from the logarithmic temperature scale setup. We work in an (effective) 1D setting exploiting matrix product operators (MPOs), which allows us to fully and uniquely implement non-Abelian and Abelian symmetries to greatly enhance numerical performance. We use our XTRG machinery to explore the thermal properties of Heisenberg models on 1D chains and 2D square and triangular lattices down to low temperatures approaching ground-state properties. The entanglement properties, as well as the renormalization-group flow of entanglement spectra in MPOs, are discussed, where logarithmic entropies (approximatelylnβ) are shown in both spin chains and square-lattice models with gapless towers of states. We also reveal that XTRG can be employed to accurately simulate the HeisenbergXXZmodel on the square lattice which undergoes a thermal phase transition. We determine its critical temperature based on thermal physical observables, as well as entanglement measures. Overall, we demonstrate that XTRG provides an elegant, versatile, and highly competitive approach to explore thermal properties, including finite-temperature thermal phase transitions as well as the different ordering tendencies at various temperature scales for frustrated systems.
ISSN:2160-3308
2160-3308
DOI:10.1103/PhysRevX.8.031082