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 accurately evaluate the thermodynamic quantities, including the internal energy, specific heat, and the spontaneous magnetization, etc, at various temperatures. By varying the anisotropic coupling constant $\Delta$, we obtain the phase diagram of the spin-1/2 Heisenberg XXZ model on the Bethe lattice, where three kinds of magnetic ordered phases, namely the ferromagnetic, XY and antiferromagnetic ordered 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 $ \Delta=\pm 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 $|\Delta|\le 1$. One remarkable feature there is that the prominent transverse correlation length reaches $\xi_c=1/\ln{(z-1)}$ for $T\leq T_c$, the maximal value allowed on a $z$-coordinated Bethe lattice.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to use the tensor network method to simulate the thermodynamic properties of the Heisenberg model at a finite temperature on the Bethe lattice. Specifically, the author extends the classical tree - tensor network (TTN) method so that it can be applied to the simulation of quantum lattice models at a finite temperature. By optimizing TTN and accurately evaluating thermodynamic quantities (such as internal energy, specific heat, and spontaneous magnetization, etc.), the phase diagrams and low - energy excitation characteristics of the system at different temperatures are studied. The main focus is on the phase - transition behavior of the spin - 1/2 Heisenberg XXZ model on the Bethe lattice, especially the magnetic - order phases and their transition characteristics under different anisotropic coupling constants \(\Delta\).
### Main research contents:
1. **Method extension**: Extend the classical tree - tensor network (TTN) method from zero - temperature to finite - temperature for simulating quantum lattice models on the Bethe lattice.
2. **Thermodynamic property calculation**: Accurately calculate thermodynamic quantities such as internal energy, specific heat, and spontaneous magnetization by optimizing TTN.
3. **Phase - diagram plotting**: By changing the anisotropic coupling constant \(\Delta\), the phase diagram of the spin - 1/2 Heisenberg XXZ model on the Bethe lattice is plotted. Three magnetic - order phases, namely ferromagnetic, XY, and antiferromagnetic, are discovered, and the continuous thermal - phase - transition critical temperatures \(T_c\) between these phases and the high - temperature paramagnetic phase are determined.
4. **Low - energy excitation analysis**: Calculate the energy gap of low - energy magnetic excitations by the linear spin - wave theory (LSWT) and compare it with the TTN simulation results.
5. **Transverse correlation length**: Study the transverse correlation length in the magnetic - order phases and find that even in the presence of an energy gap, the transverse fluctuation mode still exhibits quasi - critical behavior similar to the Goldstone mode.
### Key findings:
- **Phase diagram**: At low temperatures, the system can be in the ferromagnetic, XY, or antiferromagnetic phase, and these phases undergo continuous phase transitions with the high - temperature paramagnetic phase at \(T_c\).
- **First - order phase transition**: There is a first - order phase transition between the XY phase and the ferromagnetic and antiferromagnetic phases at the symmetric coupling point \(\Delta=\pm1\).
- **Low - energy excitation**: Low - energy magnetic excitations always have an energy gap, but the transverse fluctuation mode exhibits quasi - critical characteristics in the magnetic - order phase, and the transverse correlation length reaches a maximum value of \(\xi_c = 1/\ln(z - 1)\).
- **Quasi - Goldstone mode**: Although the system has an energy gap, a transverse fluctuation mode similar to the Goldstone mode is observed in the magnetic - order phase, and these modes are activated at certain finite energy scales.
Through these studies, the author not only verifies the effectiveness and accuracy of the TTN method at a finite temperature but also deeply explores the thermodynamic properties and phase - transition behaviors of quantum many - body systems on the Bethe lattice.