Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-$N$ tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-$N$ tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.
Numerical simulations of strongly correlated fermions at finite temperature are essential for studying high-temperature superconductivity and other quantum many-body phenomena. The recently developed tangent-space tensor renormalization group (tanTRG) provides an efficient and accurate framework by representing thermal density operators as matrix product operators. However, the particle number generally varies during the cooling process. The conventional strategy of fine-tuning chemical potentials to reach a target filling is computationally demanding. Here we propose a fixed-N tanTRG algorithm that stabilizes the average particle number by adaptively tuning the chemical potential within the imaginary-time evolution. We benchmark its accuracy on exactly solvable free fermions, and further apply it to the square-lattice Hubbard model. For hole-doped cases, we study the temperature evolution of charge and spin correlations, identifying several characteristic temperature scales for stripe formation. Our results establish fixed-N tanTRG as an efficient and reliable tool for finite-temperature studies of correlated fermion systems.
Accurate simulations of the two-dimensional (2D) Hubbard model constitute one of the most challenging problems in condensed matter and quantum physics. Here we develop a tangent space tensor renormalization group (tanTRG) approach for the calculations of the 2D Hubbard model at finite temperature. An optimal evolution of the density operator is achieved in tanTRG with a mild O(D^{3}) complexity, where the bond dimension D controls the accuracy. With the tanTRG approach we boost the low-temperature calculations of large-scale 2D Hubbard systems on up to a width-8 cylinder and 10×10 square lattice. For the half-filled Hubbard model, the obtained results are in excellent agreement with those of determinant quantum Monte Carlo (DQMC). Moreover, tanTRG can be used to explore the low-temperature, finite-doping regime inaccessible for DQMC. The calculated charge compressibility and Matsubara Green's function are found to reflect the strange metal and pseudogap behaviors, respectively. The superconductive pairing susceptibility is computed down to a low temperature of approximately 1/24 of the hopping energy, where we find d-wave pairing responses are most significant near the optimal doping. Equipped with the tangent-space technique, tanTRG constitutes a well-controlled, highly efficient and accurate tensor network method for strongly correlated 2D lattice models at finite temperature.
We propose a series-expansion thermal tensor network (SETTN) approach for efficient simulations of quantum lattice models. This continuous-time SETTN method is based on the numerically exact Taylor series expansion of the equilibrium density operator ${e}^{\ensuremath{-}\ensuremath{\beta}H}$ (with H the total Hamiltonian and $\ensuremath{\beta}$ the imaginary time), and is thus Trotter-error free. We discover, through simulating XXZ spin chain and square-lattice quantum Ising models, that not only the Hamiltonian H, but also its powers Hn, can be efficiently expressed as matrix product operators, which enables us to calculate with high precision the equilibrium and dynamical properties of quantum lattice models at finite temperatures. Our SETTN method provides an alternative to conventional Trotter-Suzuki renormalization-group (RG) approaches, and achieves a very high standard of thermal RG simulations in terms of accuracy and flexibility.
In this Letter we propose a continuous-time thermal tensor network (TTN) approach for efficient simulations of quantum lattice models. The continuous-time TTN method is based on numerically exact Taylor series expansions of equilibrium density operators e−βH (with H the total Hamiltonian and β the imaginary time), as well as unitary time-evolution operators e−iHt, and is thus Trotter-error free. We discover, through simulating a typical XXZ spin chain, that not only the Hamiltonian H, but also its powers Hn, can be efficiently expressed as matrix product operators, which enables us to calculate with high precision the equilibrium and dynamical properties of quantum lattice models at finite temperatures. Our series-expansion TTN method provides an alternative to conventional Trotter-Suzuki renormalization group (RG) approaches, and achieves an unprecedented standard of thermal RG simulations in terms of accuracy and flexibility.
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.
Various successful approaches using numerical annealing and renormalization groups have been conceived to study the thermodynamics of strongly correlated systems where perturbation or expansion theories fail to work. As the process of lowering the temperatures is usually involved in different manners, these approaches in general become much less efficient or accurate at the low temperatures. In this work, we propose to access the finite-temperature properties from the tensor network (TN) representing the zero-temperature partition function. We propose to ``scissor'' a finite part from such an infinite-size TN and ``stitch'' it to possess the periodic boundary condition along the imaginary-time direction. We dub this approach ``TN tailoring.'' Exceptional accuracy is achieved with a fine-tuning process, surpassing the previous methods including the linearized tensor renormalization group [W. Li, S.-J. Ran, S.-S. Gong, Y. Zhao, B. Xi, F. Ye, and G. Su, Phys. Rev. Lett. 106, 127202 (2011)], the continuous matrix product operator [W. Tang, H.-H. Tu, and L. Wang, Phys. Rev. Lett. 125, 170604 (2020)], etc. High efficiency is demonstrated, where the time cost is nearly independent of the target temperature, including the extremely low temperatures. The proposed idea can be extended to higher-dimensional systems of bosons and fermions.
Understanding quantum many-body states of correlated electrons is one main theme in modern condensed-matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, was recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Chen et al., Phys. Rev. X 8, 031082 (2018)], complemented by independent determinant quantum Monte Carlo, offers a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to an unprecedented low temperature of a few percent of the tunneling energy. We observe excellent agreement with ultracold fermion measurements at both half filling and finite doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behaviors of the antimoments.