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.
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.
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 Tc. 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.
Hubbard model plays a critical role in the study of strongly correlated electron systems. The only available exact solutions of 1-D Hubbard Hamiltonian, obtained through nested Bethe ansatz, show a spin charge decoupling for infinite U when there is no possibility of double occupancy of an electron orbital state. It had been shown that these features can be directly obtained by employing orthofermion algebra. In the present communication, the thermodynamic properties, namely, specific heat and entropy are provided for a 1-D chain and 2-D square lattice using orthofermion approach. It is shown that in 1-D, our results are identical to the known exact results. We have also compared the thermodynamic properties of free fermions and orthofermions in both 1-D and 2-D. Our thermodynamic results, based on orthofermion approach, can be used to ascertain the accuracy of the approximate methods employed in the context of infinite U Hubbard model.
We present analytical results of fundamental properties of one-dimensional (1D) Hubbard model with a repulsive interaction, ranging from fractional excitations to universal thermodynamics, interaction-driven criticality, correlation functions, Contact susceptibilities and quantum cooling. Using the exact solutions of the Bethe Ansatz equations of the Hubbard model, we first rigorously calculate the gapless spin and charge excitations, exhibiting exotic features of fractionalized spinons and holons. Based on the analysis on the fractional charge and spin excitations, the spin-incoherent Luttinger liquid with only the charge propagation mode is elucidated by the asymptotic of the two-point correlation functions with the help of the conformal field theory. Near quadruple critical point, we then further analytically obtain the thermodynamical properties, dimensionless ratios and scaling functions near quantum phase transitions in terms of chemical potential, magnetic field and interaction. In particular, we determine additivity rules of spin and charge susceptibilities, and derive explicit forms of thermodynamics of spin-incoherent Luttinger liquid. Finally, in order to capture deeper insight into the Mott insulator and interaction driven criticality, we further study the double occupancy and its associated Contact and Contact susceptibilities through which an adiabatic cooling scheme upon the quantum criticality is introduced.
We investigate the momentum-resolved spin and charge susceptibilities, as well as the chemical potential and double occupancy in the two-dimensional Hubbard model as functions of doping, temperature, and interaction strength. Through these quantities, we identify a weak-coupling regime, a strong-coupling regime with short-range correlations and an intermediate-coupling regime with long magnetic correlation lengths. In the spin channel, we observe an additional crossover from commensurate to incommensurate correlations. In contrast, we find charge correlations to be only short ranged for all studied temperatures, which suggests that the spin and charge responses are decoupled. These findings were obtained by a connected determinant diagrammatic Monte Carlo algorithm for the computation of double expansions, which we introduce in this paper. This permits us to obtain numerically exact results at temperatures as low as T=0.067 and interactions up to U=7, while working on arbitrarily large lattices. Our method also allows us to gain physical insights from investigating the analytic structure of perturbative series. We connect to previous work by studying smaller lattice geometries and report substantial finite-size effects.
The one-dimensional (1D) Hubbard model, describing electrons on a lattice with an on-site repulsive interaction, provides a paradigm for the physics of quantum many-body phenomena. Here, by solving the thermodynamic Bethe ansatz equations, we study the universal thermodynamics, quantum criticality, and magnetism of the 1D attractive Hubbard model. We show that the compressibility and the susceptibility of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-like state obey simple additivity rules at low temperatures, indicating an existence of two free quantum fluids. The magnetic properties, such as magnetization and susceptibility, reveal three physical regions: quantum fluids at low temperatures, a non-Fermi liquid at high temperatures, and the quantum fluid to non-Fermi liquid crossover in between. The lattice interaction is seen to significantly influence the nature of the FFLO-like state in 1D. Furthermore, we show that the dimensionless Wilson ratio provides an ideal parameter to map out the various phase boundaries and to characterize the two free fluids of the FLLO-like state. The quantum scaling functions for the thermal and magnetic properties yield the same dynamic critical exponent z=2 and correlation critical exponent $\ensuremath{\nu}=1/2$ in the quantum critical region whenever a phase transition occurs. Our results provide a rigorous understanding of quantum criticality and free fluids of many-body systems on a 1D lattice.
We numerically solve the Hubbard model on the Bethe lattice with finite coordination number z = 3, and determine its zero-temperature phase diagram.For this purpose, we introduce and develop the "variational uniform tree state" (VUTS) algorithm, a tensor-network algorithm which generalizes the variational uniform matrix product state algorithm to tree tensor networks.Our results reveal an antiferromagnetic insulating phase and a paramagnetic metallic phase, separated by a first-order doping-driven metal-insulator transition.We show that the metallic state is a Fermi liquid with coherent quasiparticle excitations for all values of the interaction strength U , and we obtain the finite quasiparticle weight Z from the single-particle occupation function of a generalized "momentum" variable.We find that Z decreases with increasing U , ultimately saturating to a nonzero, doping-dependent value.Our work demonstrates that tensor-network calculations on tree lattices, and the VUTS algorithm in particular, are a platform for obtaining controlled results for phenomena absent in one dimension, such as Fermi liquids, while avoiding computational difficulties associated with tensor networks in two dimensions.We envision that future studies could observe non-Fermi liquids, interaction-driven metal-insulator transitions, and doped spin liquids using this platform.