Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1905.12478 · Physical review. B./Physical review. B

Thermal tensor network simulations of the Heisenberg model on the Bethe lattice

Dai-Wei Qu, Wei Li, Tao Xiang

Published Sep 10, 2019• 8 citations
arXivHTMLPDF

Abstract

We have extended the canonical tree tensor network (TTN) method, which was initially introduced to simulate the zero-temperature properties of quantum lattice models on the Bethe lattice, to finite temperature simulations. By representing the thermal density matrix with a canonicalized tree tensor product operator, we optimize the TTN and compute accurately the thermodynamic quantities, including the internal energy, specific heat, spontaneous magnetization, etc., at various temperatures. By varying the anisotropic coupling constant $\mathrm{\ensuremath{\Delta}}$, we obtain the phase diagram of the spin-1/2 Heisenberg XXZ model on the Bethe lattice, where three magnetic ordered phases, namely the ferromagnetic, XY, and antiferromagnetic phases, are found in low temperatures and separated from the high-$T$ paramagnetic phase by a continuous thermal phase transition at ${T}_{c}$. The XY phase is separated from the other two phases by two first-order phase transition lines at the symmetric coupling points $\mathrm{\ensuremath{\Delta}}=\ifmmode\pm\else\textpm\fi{}1$. We have also carried out a linear spin wave calculation on the Bethe lattice, showing that the low-energy magnetic excitations are always gapped, and find the obtained magnon gaps in very good agreement with those estimated from the TTN simulations. Despite the gapped excitation spectrum, Goldstone-like transverse fluctuation modes, as a manifestation of spontaneous continuous symmetry breaking, are observed in the ordered magnetic phases with $|\mathrm{\ensuremath{\Delta}}|\ensuremath{\le}1$. One remarkable feature there is that the transverse correlation length reaches ${\ensuremath{\xi}}_{c}=1/ln(z\ensuremath{-}1)$ for $T\ensuremath{\le}{T}_{c}$, the maximal value allowed on a $z$-coordinated Bethe lattice.

DOI
10.1103/physrevb.100.125121

Related papers

Thermal tensor renormalization group simulations of square-lattice quantum spin models

Han Li, Bin-Bin Chen, Ziyu Chen +3 more

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 ${T}_{c}$.

arXiv:1904.06273•Physical review. B./Physical review. B•Jul 9, 2019•39 citations

Efficient simulation of infinite tree tensor network states on the Bethe lattice

Wei Li, Jan von Delft, Tao Xiang

We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point ($q$ is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.

arXiv:1209.2387•Physical Review B•Nov 26, 2012•84 citations

Thermodynamics of the Hubbard model on the Bethe lattice

Chen Jia-lin, Zhen Fan, Bo Zhan +6 more

We investigate the thermodynamic properties of the Hubbard model on the Bethe lattice with a coordination number of 3 using the thermal canonical tree tensor network method. Our findings reveal two distinct thermodynamic phases: a low-temperature antiferromagnetic phase, where spin SU(2) symmetry is broken, and a high-temperature paramagnetic phase. A key feature of the system is the separation of energy scales for charge and spin excitations, which is reflected in the temperature dependence of thermodynamic quantities and the disparity between spin and charge gaps extracted from their respective susceptibilities. At the critical point, both spin and charge susceptibilities exhibit singularities, suggesting that charge excitations are not fully decoupled from their spin counterparts. Additionally, the double occupancy number exhibits a nonmonotonic temperature dependence, indicative of an entropy-driven Pomeranchuk effect. These results demonstrate that the loopless Bethe lattice effectively captures the essential physics of the Hubbard model while providing a computationally efficient framework for studying strongly correlated electronic systems.

arXiv:2503.11598•Physical review. B./Physical review. B•Aug 29, 2025•4 citations

Thermal Phase Diagram of the Square Lattice Ferro-antiferromagnetic $J_1-J_2$ Heisenberg Model

Olivier Gauthé, Frédéric Mila

Using state of the art tensor network computations combined with spin wave theory, we compute the finite temperature phase diagram of the spin 1/2 $J_1-J_2$ Heisenberg model on the square lattice with ferromagnetic $J_1 < 0$ and antiferromagnetic $J_2 > 0$. We localize the first order line originating from the first order zero temperature point as well as the second order thermal transition at large $J_2$ associated with spontaneous lattice rotation symmetry breaking. This second order line terminates at a critical end point where it meets the first order line. We found no evidence for an intermediate phase.

arXiv:2310.09344•Oct 13, 2023•0 citations

Tree Tensor Network States

Tao Xiang

This chapter discusses the properties of tree tensor network states and the methods for evaluating the ground state and thermodynamic properties of quantum lattice models on a Bethe lattice or, more generally, a Husimi lattice. It starts with a brief discussion of the canonical form of a tree tensor network state. Then, a canonicalization scheme is proposed. To calculate the ground state through the imaginary time evolution, the full and simple update methods are introduced to renormalize the local tensors. Finally, as the correlation length of a quantum system is finite even at a critical point, an accurate and efficient method is described to compute the thermodynamic quantities of quantum lattice models on the Bethe lattice.

Cambridge University Press eBooks•Aug 31, 2023•0 citations

Exponential Thermal Tensor Network Approach for Quantum Lattice Models

Bin-Bin Chen, Lei Chen, Ziyu Chen +2 more

We speed up thermal simulations of quantum many-body systems in both one-(1D) and twodimensional (2D) models in an exponential way by iteratively projecting the thermal density matrix ρ ¼ e -β Ĥ 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 (approximately ln β) are shown in both spin chains and squarelattice models with gapless towers of states.We also reveal that XTRG can be employed to accurately simulate the Heisenberg XXZ model 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.

arXiv:1801.00142•Physical Review X•Sep 26, 2018•95 citations
We have extended the canonical tree tensor network (TTN) method, which was initially introduced to simulate the zero-temperature properties of quantum lattice models on the Bethe lattice, to finite temperature simulations. By representing the thermal density matrix with a canonicalized tree tensor product operator, we optimize the TTN and compute accurately the thermodynamic quantities, including the internal energy, specific heat, spontaneous magnetization, etc., at various temperatures. By varying the anisotropic coupling constant \ensuremathΔ\mathrm{\ensuremath{\Delta}}\ensuremathΔ, we obtain the phase diagram of the spin-1/2 Heisenberg XXZ model on the Bethe lattice, where three magnetic ordered phases, namely the ferromagnetic, XY, and antiferromagnetic phases, are found in low temperatures and separated from the high-TTT paramagnetic phase by a continuous thermal phase transition at Tc{T}_{c}Tc​. The XY phase is separated from the other two phases by two first-order phase transition lines at the symmetric coupling points \ensuremathΔ=\ifmmode±\else\textpm\fi1\mathrm{\ensuremath{\Delta}}=\ifmmode\pm\else\textpm\fi{}1\ensuremathΔ=\ifmmode±\else\textpm\fi1. We have also carried out a linear spin wave calculation on the Bethe lattice, showing that the low-energy magnetic excitations are always gapped, and find the obtained magnon gaps in very good agreement with those estimated from the TTN simulations. Despite the gapped excitation spectrum, Goldstone-like transverse fluctuation modes, as a manifestation of spontaneous continuous symmetry breaking, are observed in the ordered magnetic phases with ∣\ensuremathΔ∣\ensuremath≤1|\mathrm{\ensuremath{\Delta}}|\ensuremath{\le}1∣\ensuremathΔ∣\ensuremath≤1. One remarkable feature there is that the transverse correlation length reaches \ensuremathξc=1/ln(z\ensuremath−1){\ensuremath{\xi}}_{c}=1/ln(z\ensuremath{-}1)\ensuremathξc​=1/ln(z\ensuremath−1) for T\ensuremath≤TcT\ensuremath{\le}{T}_{c}T\ensuremath≤Tc​, the maximal value allowed on a zzz-coordinated Bethe lattice.
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=10W=10W=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-TTT thermal data on $W\ensuremath{\le}8$ cylinders and $W\ensuremath{\le}10$ open strips, respectively. We extract an exponential divergence versus TTT of the structure factor S(M)S(M)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{T}_{c}Tc​.
We show that the simple update approach proposed by Jiang et al. [H. C. Jiang, Z. Y. Weng, and T. Xiang, Phys. Rev. Lett. 101, 090603 (2008)] is an efficient and accurate method for determining the infinite tree tensor network states on the Bethe lattice. Ground-state properties of the quantum transverse Ising model and the Heisenberg XXZ model on the Bethe lattice are studied. The transverse Ising model is found to undergo a second-order quantum phase transition with a diverging magnetic susceptibility but a finite correlation length which is upper bounded by $1/\mathrm{ln}(q\ensuremath{-}1)$ even at the transition point (qqq is the coordinate number of the Bethe lattice). An intuitive explanation on this peculiar ``critical'' phenomenon is given. The XXZ model on the Bethe lattice undergoes a first-order quantum phase transition at the isotropic point. Furthermore, the simple update scheme is found to be related with the Bethe approximation. Finally, by applying the simple update to various tree tensor clusters, we can obtain rather nice and scalable approximations for two-dimensional lattices.
We investigate the thermodynamic properties of the Hubbard model on the Bethe lattice with a coordination number of 3 using the thermal canonical tree tensor network method. Our findings reveal two distinct thermodynamic phases: a low-temperature antiferromagnetic phase, where spin SU(2) symmetry is broken, and a high-temperature paramagnetic phase. A key feature of the system is the separation of energy scales for charge and spin excitations, which is reflected in the temperature dependence of thermodynamic quantities and the disparity between spin and charge gaps extracted from their respective susceptibilities. At the critical point, both spin and charge susceptibilities exhibit singularities, suggesting that charge excitations are not fully decoupled from their spin counterparts. Additionally, the double occupancy number exhibits a nonmonotonic temperature dependence, indicative of an entropy-driven Pomeranchuk effect. These results demonstrate that the loopless Bethe lattice effectively captures the essential physics of the Hubbard model while providing a computationally efficient framework for studying strongly correlated electronic systems.
Thermal Phase Diagram of the Square Lattice Ferro-antiferromagnetic J1−J2J_1-J_2J1​−J2​ Heisenberg Model
Using state of the art tensor network computations combined with spin wave theory, we compute the finite temperature phase diagram of the spin 1/2 J1−J2J_1-J_2J1​−J2​ Heisenberg model on the square lattice with ferromagnetic J1<0J_1 < 0J1​<0 and antiferromagnetic J2>0J_2 > 0J2​>0. We localize the first order line originating from the first order zero temperature point as well as the second order thermal transition at large J2J_2J2​ associated with spontaneous lattice rotation symmetry breaking. This second order line terminates at a critical end point where it meets the first order line. We found no evidence for an intermediate phase.
This chapter discusses the properties of tree tensor network states and the methods for evaluating the ground state and thermodynamic properties of quantum lattice models on a Bethe lattice or, more generally, a Husimi lattice. It starts with a brief discussion of the canonical form of a tree tensor network state. Then, a canonicalization scheme is proposed. To calculate the ground state through the imaginary time evolution, the full and simple update methods are introduced to renormalize the local tensors. Finally, as the correlation length of a quantum system is finite even at a critical point, an accurate and efficient method is described to compute the thermodynamic quantities of quantum lattice models on the Bethe lattice.
We speed up thermal simulations of quantum many-body systems in both one-(1D) and twodimensional (2D) models in an exponential way by iteratively projecting the thermal density matrix ρ ¼ e -β Ĥ 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 (approximately ln β) are shown in both spin chains and squarelattice models with gapless towers of states.We also reveal that XTRG can be employed to accurately simulate the Heisenberg XXZ model 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.