An important method in many-particle physics is related to the Hubbard model which can be used to investigate the behavior of fermions, particularly to study exchange and correlation of electron- electron interaction. For a special cluster, the molecule consisting of two atoms, the electrons interact via one hopping integral, and for a small number of N = 1; 2; 3 electrons, the model has the advantage that it can be solved exactly. Therefore, some important physical quantities like the Green's function, the spectral function and the self-energy are calculated and discussed. Exchange and correlation are fully treated. The self-energies obtained by two approximate methods, namely the Hartree-Fock and the GW approximation, are compared with the exact self-energy. The comparison allows a discussion of the validity of the approximations versus the ratio of Hubbard parameter and hopping integral.
An important method in many-particle physics is related to the Hubbard model which can be used to investigate the behavior of fermions, particularly to study exchange and correlation of electron- electron interaction. For a special cluster, the molecule consisting of two atoms, the electrons interact via one hopping integral, and for a small number of N = 1; 2; 3 electrons, the model has the advantage that it can be solved exactly. Therefore, some important physical quantities like the Green's function, the spectral function and the self-energy are calculated and discussed. Exchange and correlation are fully treated. The self-energies obtained by two approximate methods, namely the Hartree-Fock and the GW approximation, are compared with the exact self-energy. The comparison allows a discussion of the validity of the approximations versus the ratio of Hubbard parameter and hopping integral.
We evaluate several approximations for the self-energy operator and dielectric function of systems of interacting electrons using a two-dimensional Hubbard cluster for which the self-energy, dielectric function, and one-particle Green's function may be calculated exactly. The results show the GW approximation (in the form in which it is commonly used in ab initio calculations for real materials) to be relatively successful in establishing the main features of the spectrum, even when the electron-electron interaction is not weak. It is also clear that improving the GandW used in this approximation without including vertex corrections in the self-energy does not lead to major improvements.
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].
Various methods and approximation schemes are used to study many-electron interacting systems. Two important many-particle models, the Anderson model and the Hubbard model, and their electromagnetic properties have been investigated in many parameter regimes, and applied to physical systems. An Anderson single-impurity model Hamiltonian based calculation of the magnetic susceptibility is performed for YbN in the presence of crystal fields using an alteration of the Non-Crossing Approximation proposed by Zwicknagl et.al., incorporating parameters obtained from ab initio band structure calculations. It yields good agreement with experimental data. For the Anderson lattice model, a variational scheme which uses specific many-electron wavefunctions as basis is applied to both one- and two-dimensional systems represented by symmetric Anderson lattice Hamiltonians. Without much computational effort, the ground state energy is well approximated, especially in strong-coupling limit. Some electronic properties are examined using the variational ground state wavefunction. The one-dimensional Hubbard model has been solved exactly for small-size clusters by diagonalizing the Hamiltonian in the basis of many-electron Bloch states. The results for the energy spectrum and eigenfunctions of the ground state and low-lying excited states are presented. Also, mean field calculations of the two-dimensional single-band Hubbard model and Cu-O lattice model (three-band Hubbard model) are carried out for various physical quantities including the energy, occupation probability, staggered magnetization, momentum distribution Fermi surface and density of states, by using a projection operator formalism. To develop a systematic approach to solving many-electron problems, the many-particle partition function for the free electron gas system is explored using a cumulant expansion scheme. Starting from the ground state, the partition function can be approximated to any order in terms of excitation energy. Its application to interacting systems such as the Anderson model and the Hubbard model is briefly discussed.
A density-functional theory for many-body lattice models is considered in which the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ is the basic variable. Eigenvalue equations are derived for solving Levy's constrained search of the interaction energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is expressed as the sum of Hartree-Fock energy ${E}_{\mathrm{HF}}[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ and the correlation energy ${E}_{C}[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ Exact results are obtained for ${E}_{C}({\ensuremath{\gamma}}_{12})$ of the Hubbard model on various periodic lattices, where ${\ensuremath{\gamma}}_{\mathrm{ij}}={\ensuremath{\gamma}}_{12}$ for all nearest neighbors i and j. The functional dependence of ${E}_{C}({\ensuremath{\gamma}}_{12})$ is analyzed by varying the number of sites Na, band filling Ne, and lattice structure. The infinite one-dimensional chain and one-, two-, or three-dimensional finite clusters with periodic boundary conditions are considered. The properties of ${E}_{C}({\ensuremath{\gamma}}_{12})$ are discussed in the limits of weak $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{0})$ and strong $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$ electronic correlations, and in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}<~{\ensuremath{\gamma}}_{12}<~{\ensuremath{\gamma}}_{12}^{0}).$ Using an appropriate scaling we observe that ${\ensuremath{\varepsilon}}_{C}{(g}_{12}{)=E}_{C}{/E}_{\mathrm{HF}}$ has a pseudo-universal behavior as a function of ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}).$ The fact that ${\ensuremath{\varepsilon}}_{C}{(g}_{12})$ depends weakly on Na,Ne, and lattice structure suggests that the correlation energy of extended systems could be obtained quite accurately from finite-cluster calculations. Finally, the behaviors of ${E}_{C}({\ensuremath{\gamma}}_{12})$ for repulsive (U>0) and attractive (U<0) interactions are contrasted.
Abstract. A simple and commonly employed approximate technique with which one can examine spatially disordered systems when strong electronic correlations are present is based on the use of real-space unrestricted self-consistent Hartree-Fock wave functions. In such an approach the disorder is treated exactly while the correlations are treated approximately. In this report we critique the success of this approximation by making comparisons between such solutions and the exact wave functions for the Anderson-Hubbard model. Due to the sizes of the complete Hilbert spaces for these problems, the comparisons are restricted to small one-dimensional chains, up to ten sites, and a 4x4 two-dimensional cluster, and at 1/2 filling these Hilbert spaces contain about 63,500 and 166 million states, respectively. We have completed these calculations both at and away from 1/2 filling. This approximation is based on a variational approach which minimizes the Hartree-Fock energy, and we have completed comparisons of the exact and Hartree-Fock energies. However, in order to assess the success of this approximation in reproducing ground-state correlations we have completed comparisons of the local charge and spin correlations, including the calculation of the overlap of the Hartree-Fock wave functions with those of the exact solutions. We find that this approximation reproduces the local charge densities to quite a high accuracy, but that the local spin correlations, as represented by 〈Si · Sj〉, are not as well represented. In addition to these comparisons, we discuss the properties of the spin degrees of freedom in the HF approximation, and where in the disorderinteraction phase diagram such physics may be important. PACS numbers: 71.27.+a,71.23.-k,71.10.-wCritiquing HF for the Anderson-Hubbard Model 2 1.
A simple and commonly employed approximate technique with which one can examine spatially disordered systems when strong electronic correlations are present is based on the use of real-space unrestricted self-consistent Hartree-Fock wave functions. In such an approach the disorder is treated exactly while the correlations are treated approximately. In this report we critique the success of this approximation by making comparisons between such solutions and the exact wave functions for the Anderson-Hubbard model. Due to the sizes of the complete Hilbert spaces for these problems, the comparisons are restricted to small one-dimensional chains, up to ten sites, and a 4x4 two-dimensional cluster, and at 1/2 filling these Hilbert spaces contain about 63,500 and 166 million states, respectively. We have completed these calculations both at and away from 1/2 filling. This approximation is based on a variational approach which minimizes the Hartree-Fock energy, and we have completed comparisons of the exact and Hartree-Fock energies. However, in order to assess the success of this approximation in reproducing ground-state correlations we have completed comparisons of the local charge and spin correlations, including the calculation of the overlap of the Hartree-Fock wave functions with those of the exact solutions. We find that this approximation reproduces the local charge densities to quite a high accuracy, but that the local spin correlations, as represented by < S_i . S_j >, are not as well represented. In addition to these comparisons, we discuss the properties of the spin degrees of freedom in the HF approximation, and where in the disorder-interaction phase diagram such physics may be important.