A new nonperturbative method calculating the one-particle Green's function is introduced and applied to the Hubbard model. We show the metal–insulator transition in the paramagnetic ground state of the half-filled Hubbard model on a Bethe lattice with infinite connectivity. This method describes the transition regime appropriately. The analytic result using the continued fraction formalism which has a close relationship with present formalism is provided for comparison. As for the Hubbard model with doping, we obtain the optical conductivities and show the interesting features, such as midinfrared band and isobestic point. We also show that the sharp frequency-dependent behavior at large doping is the dispersed Drude type in this model.
We report on a nonperturbative approach to the one-dimensional (1D) and two-dimensional (2D) Hubbard models that is capable of recovering both strong- $(U\ensuremath{\gg}t)$ and weak-coupling $(U\ensuremath{\ll}t)$ limits, with U the on-site Coulomb repulsion and t the kinetic energy. Dynamical corrections to the electron self-energy in the single particle Green function are explicitly included by expanding in terms of the 16 eigenstates that characterize two nearest neighbor sites. We first show that even when U is much smaller than the bandwidth, the Mott-Hubbard gap never closes at half filling in both 1D and 2D. Consequently, the Hubbard model at half filling is always in the strong-coupling nonperturbative regime. For both large and small U, we find that the population of nearest-neighbor singlet states approaches a value of order unity as $\stackrel{\ensuremath{\rightarrow}}{T}0$ as would be expected for antiferromagnetic order. We also find that the double occupancy is a smooth monotonic function of U and approaches the anticipated noninteracting limit of
1/4 as $\stackrel{\ensuremath{\rightarrow}}{U}0$ and vanishes as $\stackrel{\ensuremath{\rightarrow}}{U}\ensuremath{\infty}.$ Finally, we compute the heat capacity
[C(T,U)] for both 1D and 2D. Our results for 1D at moderate to high temperatures are in quantitative agreement with those of the exact Bethe ansatz solution, differing by no more than
1 In addition, we find that in 2D, the
C(T,U) curves vs T for different values of U exhibit a universal crossing point at two characteristic temperatures $T\ensuremath{\approx}1.7t\ifmmode\pm\else\textpm\fi{}0.1t$ and $T\ensuremath{\approx}0.4\ifmmode\pm\else\textpm\fi{}0.1t$ as is seen universally in Hubbard models and experimentally in a wide range of strongly correlated systems such as
3He, UBe3, and ${\mathrm{CeCu}}_{6\ensuremath{-}x}{\mathrm{Al}}_{x}.$ The success of this method in recovering well-established results that stem fundamentally from the Coulomb interaction suggests that local dynamics are at the heart of the physics of strongly correlated systems.
Abstract We use the cumulant Green’s functions method (CGFM) to study the single-band Hubbard model. The starting point of the method is to diagonalize a cluster (‘seed’) containing N correlated sites and employ the cumulants calculated from the cluster solution to obtain the full Green’s functions for the lattice. All calculations are done directly; no variational or self-consistent process is needed. We benchmark the one-dimensional results for the gap, the double occupancy, and the ground-state energy as functions of the electronic correlation at half-filling and the occupation numbers as functions of the chemical potential obtained from the CGFM against the corresponding results of the thermodynamic Bethe ansatz and the quantum transfer matrix methods. The particle-hole symmetry of the density of states is fulfilled, and the gap, occupation numbers, and ground-state energy tend systematically to the known results as the cluster size increases. We include a straightforward application of the CGFM to simulate the singles occupation of an optical lattice experiment with lithium-6 atoms in an eight-site Fermi-Hubbard chain near half-filling. The method can be applied to any parameter space for one, two, or three-dimensional Hubbard Hamiltonians and extended to other strongly correlated models, like the Anderson Hamiltonian, the t − J, Kondo, and Coqblin-Schrieffer models.
The single-band Hubbard Hamiltonian [1] is the simplest model of interacting electrons in a lattice.It was developed by John Hubbard, who introduced the model to take into account the local electronic correlations in narrow energy bands.The model has a kinetic term that allows hopping of electrons between atomic sites and a term that considers the on-site interaction between electrons.This model is largely used in solid state physics to study magnetic properties of solids, insulator-metal transition (Mott transition) and high temperature superconductors.However, this is a problem that has exact analytical solution in a very few limiting cases and must be treated approximately or numerically.In this work, we propose a new methodology to analytically solve the Hubbard Hamiltonian, mapping it into a two-site model (the Hubbard dimer).To obtain the Green's function for the lattice, we employed the cumulant expansion technique, using as a 'seed' the exact two-site Green's function (the atomic method) [2].From the Green's function, we obtained the density of states [DOS] and the occupation numbers as functions of the external parameters of the model, and compared the results with the Hubbard I approximation [1].
Results for the single-particle Green's function for the strong-coupling limit of the single-band Hubbard model ($t\ensuremath{-}J$ model) obtained by numerical and analytical methods are discussed. The numerical calculation of the spectral function is performed for small systems of up to 18 sites using the Lanczos method. A detailed analysis of the spectral functions for the $t\ensuremath{-}J$ model in one and two dimensions is given. We also investigate the effect of the next-nearest-neighbor hopping terms, which appear in order
Ut2 and are omitted in the standard $t\ensuremath{-}J$ model. The role of quantum fluctuations on the shape of the spectra and the fundamental difference between one and two dimensions is elucidated. In 2D we find evidence for a coherent quasiparticle band of width
2J with dispersion $^{\ensuremath{\omega}}\mathrm{k}\ensuremath{\sim}(\frac{J}{2}){(cos{k}_{x}+cos{k}_{y})}^{2}$ and spin
S=21, which emerges at the low-energy edge of the incoherent spectrum. This form of dispersion is due to effective next-nearest-neighbor processes related to spin flips and agrees with our extension of the approach of Brinkman and Rice to finite
J.
To explore correlated electrons in the presence of local and nonlocal disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single particle spectral functions on the real-frequency axis. In the absence of off-diagonal hopping, we establish exact bounds of the spectral function of the noninteracting Bethe lattice with coordination number
Z. In the presence of interaction, the Mott insulating paramagnetic phase of the one-band Hubbard model is computed at zero temperature in alloys with site- and off-diagonal disorder. When the Hubbard
U parameter is increased, transitions from an alloy band insulator through a correlated metal into a Mott insulating phase are found to take place.