Almost 50 years after the invention of the Hubbard model [1–3] and despite modern petaflop supercomputers, a precise analysis of the criticality of this most basic model for electronic correlations has not been achieved so far, at least not in three dimensions. Dynamical mean field theory (DMFT) [4–6] was a big step forward to quantitatively calculate the three dimensional Hubbard model since the major contribution of electronic correlations, i.e., the local ones, are well captured within this theory. Local correlations give rise to quasiparticle renormalization, the Mott-Hubbard transition, magnetism, and even more subtle issues such as kinks in a purely electronic models [7]. However, non-local spatial correlations are also naturally generated by a purely local Hubbard interaction, and, as it is well known, they become of essential importance in the vicinity of secondorder phase transitions. As these correlations are neglected in DMFT, this scheme can provide only for a conventional mean-field (MF) description of the critical properties. To overcome this shortcoming cluster extensions to DMFT such as the dynamical cluster approximation (DCA) and cluster-DMFT have been proposed, see Ref. 8 for a review. In these approaches non-local correlations beyond DMFT are taken into account, however only within the range of the cluster size; and due to computational limitations the actual size of the cluster is severely restricted to about 3-5 sites for each direction. Hence, short-range correlations are included by these approaches, whereas long-range ones are not. Nonetheless, Kent et al. [9] were able to extrapolate the cluster size of so-called Betts clusters to infinity, albeit assuming from the very beginning the critical exponents to be those of the Heisenberg model. This way they extrapolated the critical temperature of the paramagnetic-toantiferromagnetic phase transition which was found to be in agreement with earlier quantum Monte Carlo (QMC) simulations of Staudt et al. [10] on finite lattices.
Magnetic properties of the half-filled Hubbard model on three-dimensional bipartite lattices have been studied by using the linked-cluster series expansion method. The series for the free energy and the staggered susceptibility have been calculated as a function of \ensuremath{\mathrm{B}}t for three-dimensional bipartite lattices for an arbitrary value of the Coulomb interaction U. We obtain a transition from a paramagnetic to an antiferromagnetic U/t-
TN /t phase diagram from the divergence of the staggered susceptibility series. Comparisons with Monte Carlo simulation, Onsager reaction field result of an effective spin Hamiltonian, and large-U Heisenberg high-temperature solutions are made.
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.
We study the pseudogaps in the spectra of the 2D Hubbard model using both finite-size and dynamical cluster approximation (DCA) quantum Monte Carlo calculations. At half-filling, a charge pseudogap, accompanied by non-Fermi-liquid behavior in the self-energy, is shown to persist in the thermodynamic limit. The DCA (finite-size) method systematically underestimates (overestimates) the width of the pseudogap. A spin pseudogap is not seen at half-filling. At finite doping, a divergent d-wave pair susceptibility is observed.
We employ the dynamical cluster approximation (DCA) in conjunction with the fluctuation exchange approximation (FLEX) to study the Hubbard model. The DCA is a technique to systematically restore the momentum conservation at the internal vertices of Feynman diagrams relinquished in the dynamical mean field approximation. The FLEX is a perturbative diagrammatic approach in which classes of Feynman diagrams are summed over analytically using geometric series. The FLEX is used as a tool to investigate the complementarity of the DCA and the finite size lattice technique with periodic boundary conditions by comparing their results for the Hubbard model. We also study the microscopic theory underlying the DCA in terms of compact (skeletal) and noncompact diagrammatic contributions to the thermodynamic potential independent of a specific model. The significant advantages of the DCA implementation in momentum space suggests the development of the same formalism for the frequency space. However, we show that such a formalism for the Matsubara frequencies at finite temperatures leads to acausal results and is not viable. However, a real frequency approach is shown to be feasible.
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.