We propose a novel algorithm with a modified Tucker decomposition for tensor network that allows for efficiently and precisely calculating the ground state and thermodynamic properties of two-dimensional (2D) quantum spin lattice systems, and is coined as the fast linearized tensor renormalization group (fLTRG). Its amazing efficiency and precision are examined by studying the spin-1/2 anisotropic Heisenberg antiferromagnet on a honeycomb lattice, and the results are found to be fairly in agreement with the quantum Monte Carlo calculations. It is also successfully applied to tackle a quasi-2D spin-1/2 frustrated bilayer honeycomb Heisenberg model, where a quantum phase transition from an ordered antiferromagnetic state to a gapless quantum spin liquid phase is found. The thermodynamic behaviors of this frustrated spin system are also explored. The present fLTRG algorithm could be readily extended to other quantum lattice systems.
Thermal tensor networks constitute an efficient and versatile representation for quantum lattice models at finite temperatures. By Trotter-Suzuki decomposition, one obtains a
D+1 dimensional TTN for the
D-dimensional quantum system and then employs efficient renormalizaton group (RG) contractions to obtain the thermodynamic properties with high precision. The linearized tensor renormalization group (LTRG) method, which can be used to contract TTN efficiently and calculate the thermodynamics, is briefly reviewed and then generalized to a bilayer form. We dub this bilayer algorithm as LTRG++ and explore its performance in both finite- and infinite-size systems, finding the numerical accuracy significantly improved compared to single-layer algorithm. Moreover, we show that the LTRG++ algorithm in an infinite-size system is in essence equivalent to transfer-matrix renormalization group method, while reformulated in a tensor network language. As an application of LTRG++, we simulate an extended fermionic Hubbard model numerically, where the phase separation phenomenon, ground-state phase diagram, as well as quantum criticality-enhanced magnetocaloric effects, are investigated.
We review the development of numerical renormalization group methods,including Wilson's numerical renormalization group,White's density matrix renormalization group,and several recent rapidly developing tensor network state-based algorithms.Among these,particular emphasis is laid on the linearized tensor renormalization group method,which was recently proposed to accurately calculate the thermodynamic properties of quantum many-body lattice systems,and which may have broad applications in both one-and two-dimensional quantum systems.
In this work, we benchmark the well-controlled and numerically accurate exponential thermal tensor renormalization group (XTRG) in the simulation of interacting spin models in two dimensions. Finite temperature introduces a finite thermal correlation length $\ensuremath{\xi}$, such that for system sizes $L\ensuremath{\gg}\ensuremath{\xi}$ finite-size calculations actually simulate the thermodynamic limit. In this paper, we focus on the square lattice Heisenberg antiferromagnet (SLH) and quantum Ising models (QIM) on open and cylindrical geometries up to width
W=10. We explore various one-dimensional mapping paths in the matrix product operator (MPO) representation, whose performance is clearly shown to be geometry dependent. We benchmark against quantum Monte Carlo (QMC) data, yet also the series-expansion thermal tensor network results. Thermal properties including the internal energy, specific heat, and spin structure factors, etc. are computed with high precision, obtaining excellent agreement with QMC results. XTRG also allows us to reach remarkably low temperatures. For SLH, we obtain an energy per site ${u}_{g}^{*}\ensuremath{\simeq}\ensuremath{-}0.6694(4)$ and a spontaneous magnetization ${m}_{S}^{*}\ensuremath{\simeq}0.30(1)$ already consistent with the ground-state properties, which is obtained from extrapolated low-
T thermal data on $W\ensuremath{\le}8$ cylinders and $W\ensuremath{\le}10$ open strips, respectively. We extract an exponential divergence versus
T of the structure factor
S(M), as well as the correlation length $\ensuremath{\xi}$, at the ordering wave vector $M=(\ensuremath{\pi},\ensuremath{\pi})$, which represents the renormalized classical behavior and can be observed over a narrow but appreciable temperature window, by analyzing the finite-size data by XTRG simulations. For the QIM with a finite-temperature phase transition, we employ several thermal quantities, including the specific heat, Binder ratio, as well as the MPO entanglement to determine the critical temperature
Tc.
The recently developed tensor renormalization-group (TRG) method provides a highly precise technique for deriving thermodynamic and critical properties of lattice Hamiltonians. The TRG is a local coarse-graining transformation, with the elements of the tensor at each lattice site playing the part of the interactions that undergo the renormalization-group flows. These tensor flows are directly related to the phase diagram structure of the infinite system, with each phase flowing to a distinct surface of fixed points. Fixed-point analysis and summation along the flows give the critical exponents, as well as thermodynamic functions along the entire temperature range. Thus, for the ferromagnetic triangular lattice Ising model, the free energy is calculated to better than 10(-5) along the entire temperature range. Unlike previous position-space renormalization-group methods, the truncation (of the tensor index range D) in this general method converges under straightforward and systematic improvements. Our best results are easily obtained with D=24, corresponding to 4624-dimensional renormalization-group flows.