International Journal of Quantum Chemistry

Theory of strongly correlated electron systems. I. Intersite Coulomb interaction and the approximation of renormalized fermions in total energy calculations

I. Sandalov, Urban Lundin, Olle Eriksson

Published Jan 1, 2005• 2 citations

Abstract

Abstract The diagrammatic strong‐coupling perturbation theory (SCPT) for correlated electron systems is developed for intersite Coulomb interaction and for a nonorthogonal basis set. The construction is based on iterations of exact closed equations for many‐electron Green functions (GFs) for Hubbard operators in terms of functional derivatives with respect to external sources. The graphs, which do not contain the contributions from the fluctuations of the local population numbers of the ion states, play a special role: a one‐to‐one correspondence is found between the subset of such graphs for the many‐electron GFs and the complete set of Feynman graphs of weak‐coupling perturbation theory (WCPT) for single‐electron GFs. This fact is used for formulation of the approximation of renormalized Fermions (ARF) in which the many‐electron quasi‐particles behave analogously to normal Fermions. Then, by analyzing: (a) Sham's equation, which connects the self‐energy and the exchange–correlation potential in density functional theory (DFT); and (b) the Galitskii and Migdal expressions for the total energy, written within WCPT and within ARF SCPT, a way we suggest a method to improve the description of the systems with correlated electrons within the local density approximation (LDA) to DFT. The formulation, in terms of renormalized Fermions LDA (RF LDA), is obtained by introducing the spectral weights of the many‐electron GFs into the definitions of the charge density, the overlap matrices, effective mixing and hopping matrix elements, into existing electronic structure codes, whereas the weights themselves have to be found from an additional set of equations. Compared with LDA+U and self‐interaction correction (SIC) methods, RF LDA has the advantage of taking into account the transfer of spectral weights, and, when formulated in terms of GFs, also allows for consideration of excitations and nonzero temperature. Going beyond the ARF SCPT, as well as RF LDA, and taking into account the fluctuations of ion population numbers would require writing completely new codes for ab initio calculations. The application of RF LDA for ab initio band structure calculations for rare earth metals is presented in part II of this study (this issue). © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005

DOI: 10.1002/qua.20210

Related papers

Theory of strongly correlated electron systems: Hubbard–Anderson models from an exact Hamiltonian, and perturbation theory near the atomic limit within a nonorthogonal basis set

I. Sandalov, Börje Johansson, Olle Eriksson

Abstract The theory of correlated electron systems is formulated in a form that allows the use as a reference point a density functional theory in the local density approximation (LDA DFT) for solids or molecules. The theory is constructed in two steps. As a first step, the total Hamiltonian is transformed into a correlated form. To elucidate the microscopic origin of the parameters of the periodic Hubbard–Anderson model (PHAM) the terms of the full Hamiltonian that have the operator structure of PHAM are separated. It is found that the matrix element of mixing interaction includes ion configuration‐ and number of particles‐dependent contributions from the Coulomb interaction. In a second step the diagram technique (DT) is developed by means of generalization of the Baym–Kadanoff method for correlated systems. The advantages of the method are that: (1) A nonorthogonal basis can be used, in particular the one generated by LDA DFT; (2) the equations for Green's functions (GFs) for the Fermi and Bose types of quasiparticles can be formulated in the form of a closed system of functional equations. The latter allows us to avoid the question of the nonunique decoupling procedure existing in previous versions of the DT and perform the expansion in terms of dressed GFs. Although the expressions for all interactions depend on the overlap matrix, it is shown that the theory is formally equivalent to one with orthogonal states with redefined interactions. When the PHAM is treated from the atomic‐limit side the vertexes are generated by kinematic interactions. The latter arise due to nontrivial commutation relations between X‐operators and come from the mixing, hopping, and overlap of states. The equations for GFs are derived within the nonorthogonal basis set in Hubbard‐I, one‐loop and random‐phase approximations with respect to kinematic interactions. The self‐consistent equation for “Hubbard Us” is derived. The technique developed is general, in particular its “bosonic” part can be used for description of spin systems with arbitrary anisotropy, systems with orbital ordering, or ordering of multipoles. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 113–143, 2003

arXiv:cond-mat/0011259International Journal of Quantum ChemistryJan 1, 200328 citations

Perturbation theory of the electronic properties in strongly correlated solids

M. M. Steiner, Robert C. Albers, Douglas J. Scalapino +1 more

In an attempt to obtain a better description than the local-density approximation (LDA) for the electronic structure of realistic narrow-band systems, we have studied a perturbation-theory method for the self-energy and have tested it on the electronic properties of the periodic one-dimensional Anderson model. For a number of models we have calculated a wide range of properties such as quasiparticle energies, spectral densities, reduced hybridizations, and various magnetic correlation functions. For two different models we were able to compare with Monte Carlo results and found good agreement. In general, we have obtained a description of narrow-band systems that is more reasonable than the one given by the LDA.

Physical review. B, Condensed matterJan 15, 199143 citations

A closed local-orbital unified description of DFT and many-body effects

Fernando Flores, Diego Soler‐Polo, José Ortega

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.

Journal of Physics Condensed MatterMay 11, 20223 citations

Many-body Theory for Excited-states Process

Meichun Huang

The most of the experimental measurable quantities in many-body electron system, such as the optical absorption, luminescent spectra and exciton effect, are related to the excitation states description. The density functional theory with local density approximation (DFT-LDA) as a ground state theory, i.e. the solutions based on the Kohn-Sham equation is a powerful tool for studying the ground state properties in many-particle systems. However, the first principles computation of excited states is more complexity than ground-state calculations. A key problem is that the exchange-correlation interaction in excited states is differ from the ground-states. In recent years, nevertheless, several electronic-excitation theories have been developed. The most important theoretical and computational method include the many-body perturbation theory which is based on the quasi-particle concept and the Green function equation, the time-dependent DFT and the Bethe-Salpeter equation for describing the electron-hole interaction. Among them the central ingredient is an electron's self-energy operator Σ that describes the exchange correlation effect beyond Hartree approximation. The implementations of the above theories are unavoidable to introduce some of approximations. A good approximation for Σ is the GW approach by Hedin. It is shown from computation simulations for many real condensed matter systems that the GW approach is a successful method for electronic-excitation problems. Another effective method for excited-states is so-called screen-exchange local-density approximation (sX-LDA). It is a combination of LDA and Hartree-Fock (HF)theory, in which the intrinsic local screen-exchange interaction has been replaced by a non-local one and a generalized Kohn-Sham equation (GKS) is given. Based on a PWscf package, the sX-LDA-PWscf method has been tested and indicates that an available approach for treating band-gap problem in semiconductors. In this paper, the many-body theories for the excitation process, its development and signification are reviewed. The relation and differences consist in different theoretical methods are discussed. On the based of above discussions, the applications of many-body theory to the band-band transition (band-gap underestimation problem) and exciton effect in semiconductors are presented.

Chinese Journal of LuminescenceJan 1, 20050 citations

Self-consistent density-functional approach to the correlated ground states and an unrestricted many-body perturbation theory

Behnam

We introduce a scheme through which the ground-state electronic density n and total energy Ey[n] of inhomogeneous interacting systems can be determined in a self-consistent (SC) way. It is based on a SC perturbation expansion for the single-particle Green function G. Our employed expansion is argued to be valid under all conditions provided that the ground-state charge density of the interacting system is pure-state non-interacting v representable. The formalism discussed in this paper has been applied to a quasi-one-dimensional periodic system. It has even been possible to carry out SC calculations with the dynamically screened exchange diagram (the so-called G W diagram, where W is the dynamically screened electron4ectron interaction function) for the selfenergy operator C. A coherent mode of interpretation of a number of our results entails that the overwhelming success of the G W framework in relatively accurately reproducing the experimental energy gaps of a number of semiconductors and insulators may be a direct consequence of non-selfconsistency effects; whereas a non-self-consistent calculation in our tests increases the value of the local-density approximation (LDA) gap in the range of 40-65% of the ‘experimental’ gaps, self-consistency effects reduce this range to 1 4 1 7 % . The same line of argument leads us to conclude that, in our model, 2635% of the LDA error in the gap energies is due to the absence of a well known discontinuity in the LDA exchange-correlation potential, and the remaining part is due to the Kohn-Sham eigenvalues not being identifiable with one-particle excitation energies. A natural conclusion that may be drawn from these observations is that a correct description of excitation energies (e.g. of the gap energy) cannot be achieved without including, at least approximately, the vertex part of C, which is entirely absent in the GW scheme. 5 1 . INTRODUCTION The ground-state (GS) density-functional theory (DFT) due to Hohenberg and Kohn (1964) provides a rigorous foundation for determining properties of systems in their ground states (for reviews see Lundqvist and March (1983), Jones and Gunnarsson (1989) and Dreizler and Gross (1990)). The central position in this theory is taken by the total-energy functional EY[n] which, owing to its unique absolute minimum a t the GS electronic density n, can be used to obtain not only the GS energy but also the GS charge density. tPresent address: Max-Planck-Institut fur Festkorperforschung, HeisenbergstraRe 1, 70569 Stuttgart, Federal Republic of Germany. 0141-8637/97 $12 .00

Jan 1, 20060 citations

Self-consistent density-functional approach to the correlated ground states and an unrestricted many-body perturbation theory

Behnam Farid

We introduce a scheme through which the ground-state electronic density n and total energy Ev [n] of inhomogeneous interacting systems can be determined in a self-consistent (SC) way. It is based on a SC perturbation expansion for the single-particle Green function G. Our employed expansion is argued to be valid under all conditions provided that the ground-state charge density of the interacting system is pure-state non-interacting v representable. The formalism discussed in this paper has been applied to a quasi-one-dimensional periodic system. It has even been possible to carry out SC calculations with the dynamically screened exchange diagram (the so-called GW diagram, where W is the dynamically screened electron-electron interaction function) for the self-energy operator Σ. A coherent mode of interpretation of a number of our results entails that the overwhelming success of the GW framework in relatively accurately reproducing the experimental energy gaps of a number of semiconductors and insulators may be a direct consequence of non-self-consistency effects; whereas a non-self-consistent calculation in our tests increases the value of the local-density approximation (LDA) gap in the range of 40–65% of the ‘experimental’ gaps, self-consistency effects reduce this range to 14–17%. The same line of argument leads us to conclude that, in our model, 26–35% of the LDA error in the gap energies is due to the absence of a well known discontinuity in the LDA exchange-correlation potential, and the remaining part is due to the Kohn-Sham eigenvalues not being identifiable with one-particle excitation energies. A natural conclusion that may be drawn from these observations is that a correct description of excitation energies (e.g. of the gap energy) cannot be achieved without including, at least approximately, the vertex part of Σ, which is entirely absent in the GW scheme.

Philosophical Magazine BAug 1, 199718 citations