The one-dimensional t 1 − t 2 half-filled Hubbard model is considered at finite temperatures T within a dynamical cluster approximation (DCA) with N c = 24 in combination with a semiclassical approximation (SCA) impurity solver. The SCA approach accounts for long-range spatial fluctuations, where exact numerical impurity solvers can not capture due to computational expense, even though dynamical fluctuations are freezing. Therefore, it can consider both frequency- and momentum-resolved physical properties beyond the DCA with small cluster in combination with exact impurity solvers. By the computation of the static spin-spin correlation, the density of states, and the double occupancy, we examine the description of the frustrated one-dimensional systems at finite T within given approximations. We confirm not only the interaction-driven metal-insulator transition in the regions of t 2/t 1 > 0.5, but also the commensurate-incommensurate transition by tunning t 2/t 1 in the strong interaction region. We also observe finite T-driven metal-insulator transition.
We formulate a finite-temperature scheme of the variational cluster approximation (VCA) particularly suitable for an exact-diagonalization cluster solver. Based on the analytical properties of the single-particle Green's function matrices, we explicitly show the branch-cut structure of logarithm of the complex determinant functions appearing in the self-energy-functional theory (SFT) and whereby construct an efficient scheme for the finite-temperature VCA. We also derive the explicit formulas for entropy and specific heat within the framework of the SFT. We first apply the method to explore the antiferromagnetic order in a half-filled Hubbard model by calculating the entropy, specific heat, and single-particle excitation spectrum for different values of on-site Coulomb repulsion
U and temperature
T. We also calculate the
T dependence of the single-particle excitation spectrum in the strong coupling region, and discuss the overall similarities to and the fine differences from the spectrum obtained by the spin-density-wave mean-field theory at low temperatures and the Hubbard-I approximation at high temperatures. Moreover, we show a necessary and sufficient condition for the third law of thermodynamics in the SFT. On the basis of the thermodynamic properties, such as the entropy and the double occupancy, calculated via the
T and/or
U derivative of the grand potential, we obtain a crossover diagram in the
(U,T) plane, which separates a Slater-type insulator and a Mott-type insulator. Next, we demonstrate the finite-temperature scheme in the cluster-dynamical-impurity approximation (CDIA), i.e., the VCA with noninteracting bath orbitals attached to each cluster, and study the paramagnetic Mott metal-insulator transition in the half-filled Hubbard model. Formulating the finite-temperature CDIA, we first address a subtle issue regarding the treatment of the artificially introduced bath degrees of freedom which are absent in the originally considered Hubbard model. We then apply the finite-temperature CDIA to calculate the finite-temperature phase diagram in the
(U,T) plane. Metallic, insulating, coexistence, and crossover regions are distinguished from the bath-cluster hybridization-variational-parameter dependence of the grand-potential functional. We find that the Mott transition at low temperatures is discontinuous, and the coexistence region of the metallic and insulating states persists down to zero temperature. The result obtained here by the finite-temperature CDIA is complementary to the previously reported zero-temperature CDIA phase diagram.
We investigate the Hubbard model on the triangular lattice at half filling using the dynamical cluster approximation (DCA) and dual fermion (DF) methods in combination with continuous-time quantum Monte Carlo (CT QMC) and semiclassical approximation (SCA) methods. We study the one-particle properties and nearest-neighbor spin correlations using the DCA method. We calculate the spectral functions using the CT QMC and SCA methods. The spectral function in the SCA and obtained by analytic continuation of the Pad\'e approximation in CT QMC are in good agreement. We determine the metal-insulator transition and the hysteresis associated with a first-order transition in the double-occupancy and nearest-neighbor spin-correlation functions as functions of temperature. As a further check, we employ the DF method and discuss the advantages and limitation of the dynamical mean-field theory, DCA, and recently developed DF methods by comparing the Green's functions. We find an enhancement of antiferromagnetic correlations and provide evidence for magnetically ordered phases by calculating the spin susceptibility.
Abstract We investigate the finite temperature properties of the half filled Hubbard model in two dimensions, with onsite interaction ( U ), in presence of (frustrating) next nearest neighbor hopping ( t ′) using a semiclassical approximation scheme. We show that introduction of t ′ results in a finite temperature pseudogapped (PG) phase that separates the small U Fermi liquid and large U Mott insulator. We map out the PG to normal metal crossover temperature scale ( T *) as a function of U and t ′. We demonstrate that in the PG phase, the quadratic dependence of resistivity on temperature is violated due to thermally induced spin fluctuations. We conclude with exact diagonalization calculations, that complement our finite temperature results, and indicate the presence of a frustration driven PG state between the Fermi liquid and the Mott insulator at zero temperature as well.
Via a dynamical cluster approximation with N c = 4 in combination with a semiclassical approximation (DCA+SCA), we study the doped two-dimensional Hubbard model. We obtain a plaquette antiferromagnetic (AF) Mott insulator, a plaquette AF ordered metal, a pseudogap (or d-wave superconductor) and a paramagnetic metal by tuning the doping concentration. These features are similar to the behaviors observed in copper-oxide superconductors and are in qualitative agreement with the results calculated by the cluster dynamical mean field theory with the continuous-time quantum Monte Carlo (CDMFT+CTQMC) approach. The results of our DCA+SCA differ from those of the CDMFT+CTQMC approach in that the d-wave superconducting order parameters are shown even in the high doped region, unlike the results of the CDMFT+CTQMC approach. We think that the strong plaquette AF orderings in the dynamical cluster approximation (DCA) with N c = 4 suppress superconducting states with increasing doping up to strongly doped region, because frozen dynamical fluctuations in a semiclassical approximation (SCA) approach are unable to destroy those orderings. Our calculation with short-range spatial fluctuations is initial research, because the SCA can manage long-range spatial fluctuations in feasible computational times beyond the CDMFT+CTQMC tool. We believe that our future DCA+SCA calculations should supply information on the fully momentum-resolved physical properties, which could be compared with the results measured by angle-resolved photoemission spectroscopy experiments.
Abstract We study metal-insulator phase diagram in the half-filled disordered Hubbard model at finite temperature. We calculate the averaged local density of states at the Fermi level and the site occupation as a function of the on-site energy at finite temperature by using typical medium theory with an impurity solver of equation of motion method. Although the metallicity in the correlated metallic region changes with temperature change, the phase diagram is almost temperature independent. We find fairly good agreement between our site occupation and those obtained by the typical medium method with other impurity solver.
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.