With the eigenfunctional theory, we study the correlation of the local electrons in the Kondo lattice model, and give a general formulation to calculate the correlation functions of the local electrons where it ends in order to calculate the partition function of the conduction electrons with time-dependent external potentials. This formulation may be very useful to study the low-energy physical property of the Kondo lattice model by numerical simulations, especially by the quantum Monte Carlo method. In order to justify this formulation, we apply it for one and two magnetic impurities scattering problems in the continuum limit, and obtain reliable results that are in agreement with previous exact results.
By studying the Falicov--Kimball lattice model, we demonstrate that the correlation functions of local electrons can be written directly in the form of the partition function of conduction electrons with a time dependent external potential. These general expressions of the correlation functions are equivalent to the exact ones obtained by the dynamical mean field theory in infinite dimension limit. Moreover, we also prove that the asymptotic power-law behavior of the local electron correlation functions is universal in any dimensions, as long as the conduction electron excitation spectrum is gapless, regardless whether they are Landau Fermi liquid or not. This conclusion is also valid for other quantum many-particle systems, if the coupling between the conduction and local electrons is on-site or short-range Coulomb interaction.
Ground-state properties of the one-dimensional Kondo lattice model are investigated using the quantum Monte Carlo (QMC) simulation method. To reduce the difficulty of the negative sign problem, we employ the variationally optimized singlet wave function. Since the improvement of Monte Carlo sampling efficiency is achieved in the half-filling case, we can perform the simulations on large systems (up to 26-site system). The QMC data on the local spin correlations are closely compared with those of the variational Monte Carlo (VMC) calculations as well as with the exact diagonalization data of the Ruderman-Kittel-Kasuya-Yoshida (RKKY) lattice model. We find that the local spin correlation decays exponentially with the correlation length ${\ensuremath{\xi}}_{\mathit{f}\mathit{f}}$\ensuremath{\sim}2.54 in the QMC data, while ${\ensuremath{\xi}}_{\mathit{f}\mathit{f}}$\ensuremath{\sim}1.08 in the VMC results (J/t=1). This implies the absence of the long-range order (the existence of the spin gap) and the importance of intersite correlation effects which are not taken into account in the VMC results. On the other hand, the system size dependence is almost absent in the Kondo lattice model results, which contrasts the behavior in the RKKY data.
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.
We present, in this dissertation, a numerical simulation method to study interacting fermion systems. The general simulation procedures are discussed in connection with a description of a Quantum Monte Carlo simulation algorithm for interacting electrons on lattices. The algorithm presented here has been used to simulate interacting electrons on lattices, and it makes possible the study of substantially larger systems than can be studied by other numerical methods. As long as certain limits of applicability are respected, model Hamiltonians of interacting electrons can be studied without resort to uncontrolled approximations. The method then provides nearly exact solutions to model Hamiltonians of many-body systems, in the sense that the degree of error can be controlled. We also discuss some results obtained from simulations of the extended Hubbard model and and the periodic Anderson model. For the extended Hubbard model in two dimensions, different regions are identified where electron correlation produces antiferromagnetism, charge-density-waves, and singlet pairing superconducting behaviors. We also find regions where transitions from one type of correlation to another occur. For the symmetric periodic Anderson model, we observe the formation of local spin moments at high temperatures and their quenching at low temperatures, as in the single impurity Anderson model. In addition, we find antiferromagnetic interaction between f local moments at low temperatures which are not present in a single impurity problem.
This thesis discusses mainly two Fermionic lattice systems, first a Kondo lattice with additional Hubbard interaction and second a Hubbard Hamiltonian augmented with additional spin and charge interactions. We first introduce the Quantum Monte Carlo technique, which is then employed to study the two respective systems. We present an innovation that allows to calculate time displaced Greens functions more efficiently. Compared with previously used numerically stable algorithms the new method gains an order of magnitude in speed, but is just as precise, and very simple to implement. In the second chapter we consider the Kondo lattice model in two dimensions at half filling. In addition to the Fermionic hopping integral t and the superexchange coupling J the role of a Coulomb repulsion U in the conduction band is investigated. We find the model to display a magnetic order-disorder transition in the U-J plane with a critical value of Jc which is decreasing as a function of U. The single-particle spectral function A(k,ω) is computed across this transition. We conclude that (i) the local screening of impurity spins determines the low-energy behavior of the spectral function and (ii) one cannot deform continuously the spectral function of the half-filled Hubbard model at J=0 to that of the Kondo insulator at J>Jc. In the third chapter we investigate the phase diagram of a new model that exhibits a first order transition between s-wave superconducting and antiferromagnetic phases. The model, a generalized Hubbard model augmented with competing spin-spin and pair-pair interactions, was investigated using the projector quantum Monte Carlo method. Upon varying the Hubbard U from attractive to repulsive, we find a first order phase transition between superconducting and antiferromagnetic states.
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.