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.
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.
The Hubbard model is investigated in the framework of lattice density functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable of the many-body problem. A new approximation to the interaction-energy functional $W[\ensuremath{\gamma}]$ is proposed which is based on its scaling properties and which recovers exactly the limit of strong electron correlations at half-band filling. In this way, a more accurate description of W is obtained throughout the domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}},$ including the crossover from weak to strong correlations. As examples of applications results are given for the ground-state energy, charge-excitation gap, and charge susceptibility of the Hubbard model in one-, two-, and three-dimensional lattices. The performance of the method is demonstrated by comparison with available exact solutions, with numerical calculations, and with LDFT using a simpler dimer ansatz for W. Goals and limitations of the different approximations are discussed.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
A methodology is introduced for the systematic treatment of electron-electron correlations in solids and other interacting quantum N-particle systems. The method is developed within the framework of electronic structure theory (band theory) but, in contrast to conventional approaches, which are based on the single-particle picture, it is formulated within a many-particle picture in which n particles in d-dimensional phase space are treated as a single particle in a phase space of nd dimensions. In this phase space, interparticle interactions appear as external potentials allowing the treatment of the system of particles through the use of single-particle methods, while at the same time allowing a systematic, direct, and nonperturbative treatment of interparticle interactions. The method makes use of the invariance of the Hamiltonian describing an interacting-particle system under partitioning into subsystems of n particles. This treatment leads to exact results in the limit $n\ensuremath{\rightarrow}N$. Based on such partitioning, we propose a generalization of density functional theory and an appropriately defined local density approximation to treat the interactions between the n-particle units in a system of N>n particles. This approach yields n-particle correlated densities and n-particle states which can be used in an analysis of the electronic properties of materials, such as total energy, phase stability, electronic transport, and others. We use the formal construct of multiple-scattering theory to develop the method for the calculation of the two-particle electronic structure of a solid and the corresponding total energy of the ground state. We also illustrate some of the properties of the method in terms of a Hubbard model Hamiltonian on a linear ring. Various features of the method and further possible applications are presented in a discussion section.
In the general frame of a modified many-body perturbation theory, the single-band Hubbard model is treated in the limit of weak electron-electron interaction. In this modification the decomposition of the Hamiltonian into an effective free system and an effective perturbation is rearranged iteratively in each order of the formalism, resulting in a successive reduction of the interaction part. In comparison with conventional perturbation techniques, this means an extended range of validity. Imposing thermodynamic self-consistency, we calculate the wave-vector-dependent electronic self-energy up to second order and analyze the influence of the Coulomb correlations. Several well-known physical properties of the Hubbard model are recovered in characteristic single-particle quantities. Therefore, the effect of the band occupation n and the electron-electron repulsion U on spectral densities, band structures, densities of states, magnetization, etc. is investigated in detail. Some of the most remarkable results are: first, a splitting of the free Bloch band into two quasiparticle subbands when the coupling strength reaches moderate values. These subbands are linked by a small damping-induced background in the density of states and are separated roughly by an energy of the order of U. In addition, the upper subband is shifted to higher energies in the spectrum proportionally to U while the lower one is nearly fixed. Second, we found a paramagnetic-ferromagnetic phase diagram that has a critical coupling strength Uc and exhibits strong indications of a lower critical occupation number nc. There is no spontaneous magnetization in the system for parameters smaller than these values. Compared with the Stoner model, one thus obtains striking improvements. In many cases our conclusions for the weak-coupling limit are qualitatively in good agreement with well-known results in the strong-correlation regime.
In this work we investigate some aspects of the physics of strongly correlated systems by taking into account both electron-electron and electron-phonon interactions as basic mechanisms for reproducing electronic correlations in real materials. The relevance of the electron-electron interactions is discussed in the first part of this thesis in the framework of a self-consistent theoretical approach, named Composite Operator Method (COM), which accounts for the relevant quasi-particle excitations in terms of a set of composite operators that appear as a result of the modification imposed by the interactions on the canonical electronic fields. We show that the COM allows the calculation of all the relevant Green’s and correlation functions in terms of a number of unknown internal parameters to be determined self-consistently. Therefore, depending on the balance between unknown parameters and self-consistent equations, exact and approximate solutions can be obtained. By way of example, we discuss the application of the COM to the extended t-U- J-h model in the atomic limit, and to the two-dimensional single-band Hubbard model. In the former case, we show that the COM provides the exact solution of the model in one dimension. We study the effects of electronic correlations as responsible for the formation of a plethora of different charge and/or spin orderings. We report the phase diagram of the model, as well as a detailed analysis of both zero and finite temperature single-particle and thermodynamic properties. As far as the single-band Hubbard model is concerned, we illustrate an approximated selfconsistent scheme based on the choice of a two-field basis. We report a detailed analysis of many unconventional features that arise in single-particle properties, thermodynamics and system’s response functions. We emphasize that the accuracy of the COM in describing the effects of electronic correlations strongly relies on the choice of the basis, paving the way for possible multi-pole extensions to the twofield theory. To this purpose, we also study a three-field approach to the single-band Hubbard model, showing a significant step forward in the agreements with numerical data with respect to the two-pole results. The role of the electron-phonon interaction in the physics of strongly correlated systems is discussed in the second part of this thesis. We show that in highly polarizable lattices the competition between unscreened Coulomb and Fröhlich interactions results in a short-range polaronic exchange term Jp that favours the formation of local and light pairs of bosonic nature, named bipolarons, which condense with a critical temperature well in excess of hundred kelvins. These findings, discussed in the framework of the so-called polaronic t-Jp model, are further investigated in the presence of a finite on-site potential ~U , coming from the competition between on-site Coulomb and Fröhlich interactions. We discuss the role of ~U as the driving parameter for a small-to-large bipolaron transition, providing a possible explanation of the BEC-BCS crossover in terms of the properties of the bipolaronic ground state. Finally, we show that a hard-core bipolarons gas, studied as a charged Bose-Fermi mixture, allows for the description of many non Fermi liquid behaviours, allowing also for a microscopic explanation of pseudogap features in terms of a thermal-induced recombination of polarons and bipolarons, without any assumption on preexisting order or broken symmetries. [edited by Author]
The problem of the ground state of the electronic system in the Hubbard model for U=∞ is discussed. The author investigates the normal (singlet or nonmagnetic) N state of the electronic system over the entire range of electron densities n⩽1. It is shown that the energy of the N state ɛ 0 (1) (n) in a one-particle approximation, such as (e.g.) the extended Hartree-Fock approximation, is lower than the energy of the saturated ferromagnetic FM state ɛ FM(n) for all n. The dynamic magnetic susceptibility is calculated in the random phase approximation, and it is shown that the N state is stable over the entire range of electron densities: The static susceptibility (ω=0) does not have a band singularity in the zero-wave vector limit q→0. A formally exact representation is obtained for the mass operator of the one-particle Green’s function, and an approximation of this operator is proposed: M k(E)⋍λF(E), where λ=n(1−n)/(1−n/2)z is the kinematic interaction parameter, z is the number of nearest neighbors, and F(E) is the total single-site Green’s function. For an elliptical density of states the integral equation for F(E) is solved exactly, ad it is shown that the spectral intensity rigorously satisfies the sum rule. The calculated energy of the strongly correlated N state ɛ 0(n)<ɛ FM(n) for all n, and in light of this relationship the author discusses the hypothesis that the ground state of the system is the normal (singlet) state in the thermodynamic limit. The electron distribution function at T=0 differs significantly from the Fermi step; it is “smeared” along the entire energy spectrum, and discontinuities do not occur in the region of the chemical potential m.