A formulation is proposed to treat the ground-state properties of an assembly of electrons, interacting mutually through the Coulomb potential and confined by an external potential V whose maximum is finite so that both bound and extended continuum states appear as eigen states for the one-body Hamiltonian composed of the electron kinetic energy and V . These extended states are found to play a crucial role in incorporating the effect of electron correlation into the localized ground-state many-electron wave function. As an illustration of our formalism, the ground-state energy of a He atom is calculated in its simplest version and compared rather successfully with experiment as well as the results in other serious numerical methods. We also calculate the density and pair distribution functions in the same version. Special attention is paid to the latter function in this system in order to elucidate a qualitative difference from that in the homogeneous electron gas. This function shows clearly that the two electrons with opposite spins in He correlate strongly over quite a long range in contrast to the common expectation that the ground-state wave function for the He atom can be described accurately in the Hartree-Fock approximation.
A formulation is proposed to treat the ground-state properties of an assembly of electrons, interacting mutually through the Coulomb potential and confined by an external potential V whose maximum is finite so that both bound and extended continuum states appear as eigen states for the one-body Hamiltonian composed of the electron kinetic energy and V . These extended states are found to play a crucial role in incorporating the effect of electron correlation into the localized ground-state many-electron wave function. As an illustration of our formalism, the ground-state energy of a He atom is calculated in its simplest version and compared rather successfully with experiment as well as the results in other serious numerical methods. We also calculate the density and pair distribution functions in the same version. Special attention is paid to the latter function in this system in order to elucidate a qualitative difference from that in the homogeneous electron gas. This function shows clearly that the two electrons with opposite spins in He correlate strongly over quite a long range in contrast to the common expectation that the ground-state wave function for the He atom can be described accurately in the Hartree-Fock approximation.
Strong electron correlation in materials is explored within a class of model Hamiltonians that treat only pair interactions between electrons. The model is unique among typical spin Hamiltonians in that it does not have an effective mean-field reference wave function. The ground-state wave functions from all Hamiltonians in the model have the same one-electron reduced density matrix (1-RDM); consequently, one-electron theories such as the Hartree–Fock and density functional theories are inapplicable. In contrast, the ground-state two-electron reduced density matrix (2-RDM) has a one-to-one mapping to the ground-state wave function. For a range of lattices including linear, ladder, and square topologies, we variationally compute the 2-RDM subject to constraints, known as N -representability conditions, that are necessary for the 2-RDM to represent an N -electron ensemble density matrix. We find that for all model Hamiltonians the 2-RDM is accurately computed as long as the D, Q, and G N -representability conditions are supplemented with the T 2 condition. Energies, orbital and geminal occupation numbers, and correlation functions are computed. The model has applications to superconductivity as well as more general pairing phenomena in electronic systems. Effective methods are needed to treat strong electron correlation arising in the study of materials including the study of high-temperature superconductivity and phase transitions.
In parallel with standard density-functional theory, we study the energy of the ground state of a finite many-body system as a functional of the one-particle density matrix. We show that the formulation of a variational principle that is valid beyond the Hartree-Fock limit requires that two-body correlations be included not only in the ground-state energy but also in the constraints. As an illustration, we apply a linear-response argument to derive formulas for first-order corrections to the Hartree-Fock density matrix. Further analysis suggests an approach in terms of the density matrix of an independent-particle system, which can be introduced by the application of an alternative variational principle. This approach is reminiscent of Kohn-Sham theory, but the effective external potential is not required to be local. This variational method can be implemented in a systematic fashion by means of the linked-cluster expansion. In an appendix we study a variant of the Hohenberg-Kohn theorem for nonlocal potentials.
Real knowledge is to know the extent of one's ignorance. Confucius, 500 BCE Summary The topic of this chapter is a small selection of the vast array of experimentally observed phenomena chosen to exemplify crucial roles played by the electron– electron interaction. Examples in the present chapter bring out the effects of correlation in ground and excited states as well as in thermal equilibrium. These raise challenges for theory and quantitative many-body methods in treating interacting electrons, the topics of the following chapters. The title of this book is Interacting Electrons . Of course, there are no non-interacting electrons: in any system with more than one electron, the electron–electron interaction affects the energy and leads to correlation between the electrons. All first-principles theories deal with the electron–electron interaction in some way, but often they treat the electrons as independent fermions in a static mean-field potential that contains interaction effects approximately. As described in Ch. 4, the Hartree–Fock method is a variational approximation with a wavefunction for fermions that are uncorrelated, except for the requirement of antisymmetry. The Kohn–Sham approach to DFT defines an auxiliary system of independent fermions that is chosen to reproduce the ground-state density. It is exact in principle and remarkably successful in practice. However, many properties such as excitation energies are not supposed to be taken directly from the Kohn–Sham equations, even in principle. Various other methods attempt to incorporate some effect of correlation in the choice of the potential. This chapter is designed to highlight a few examples of experimentally observed phenomena that demonstrate qualitative consequences of electron–electron interactions beyond independent-particle approximations. Some examples illustrate effects that cannot be accounted for in any theory where electrons are considered as independent particles. Others are direct experimental measurements of correlation functions that would vanish if the electrons were independent. In yet other cases, a phenomenon can be explained in terms of independent particles in some effective potential, but it is deeply unsatisfying if one has to invent a different potential for every case, even for different properties in the same material. A satisfactory theory ultimately requires us to confront the problem of interacting, correlated electrons.
Abstract Density functional theory (DFT) is usually formulated in terms of the electron density as a function of position n(r). Here we discuss an alternative formulation of DFT in terms of the orbital occupation numbers {n α } associated with a local-orbital orthonormal basis set {ϕ α }. First, we discuss how the building blocks of DFT, namely the Hohenberg–Kohn theorems, the Levy–Lieb approach and the Kohn–Sham method, can be adapted for a description in terms of {n α }. In particular, the total energy is now a function of {n α }, E[{n α }], and a Kohn–Sham-like Hamiltonian is derived introducing the effects of the electron–electron interactions via effective potentials, { V α eff = ∂ E e e [ { n β } ] / ∂ n α } . In a second step we consider the Hartree and exchange energies and discuss how to describe them, in the spirit of a DFT approach, in terms of the orbital occupation numbers. In this contribution special attention is paid to the description of the (intra-atomic) correlation energy and corresponding correlation potentials {V corr,α }. For this purpose, a model system is analyzed in detail, whereby an atomic Hamiltonian interacts with the environment via a simplified model; the use of this model allows us to obtain the correlation energy and potentials (in terms of {n α }) for different cases corresponding to low, intermediate and high electron correlations.
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.
After a brief summary of some basic properties of ideal gases of bosons and of fermions, two many-body Hamiltonians are cited for which ground-state wave functions allow the generation of excited states. Because of the complexity of ground-state many-body wave functions, we then consider properties of reduced density matrices, and in particular, the diagonal element of the second-order density matrix. For both the homogeneous correlated electron liquid and for an assembly of charged bosons, the ground-state pair correlation function g(r) has fingerprints of the zero-point energy of the plasmon modes. These affect crucially the static structure factor, S(k), in the long wavelength limit. This is best understood by means of the Ornstein–Zernike direct correlation function c(r), which plays an important role throughout this article. Turning from such charged liquids, both boson and fermion, to superfluid 4He, the elevated temperature (T) structure factor S(k, T) is related, albeit approximately, to its zero-temperature counterpart, via the velocity of sound, reflecting the collective phonon excitations and the superfluid density. Finally, some future directions are pointed.