To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
There are nine and sixty ways of constructing tribal lays, And every single one of them is right!
Rudyard Kipling, In the Neolithic Age
Overview of Chapters 12–17
There are three basic approaches to the calculation of independent-particle electronic states in materials. There are no fundamental disagreements: all agree when applied carefully and taken to convergence. Indeed, each of the approaches leads to instructive, complementary ways to understand electronic structure and each can be developed into a general framework for accurate calculations.
Each method has its advantages: each is most appropriate for a range of problems and can provide particularly insightful information in its realm of application.
Each method has its pitfalls: the user beware. It is all too easy to make glaring errors or over-interpret results if the user does not understand the basics of the methods.
The three types of methods and their characteristic pedagogical values are:
1. Plane wave and grid methods provide general approaches for solution of differential equations, including the Schrödinger and Poisson equations. At first sight, plane waves and grids are very different, but in fact each is an effective way of representing smooth functions. Furthermore, grids are involved in modern efficient plane wave calculations that use fast Fourier transforms.
Chapter 12 is devoted to the basic concepts and methods of electronic structure. Plane waves are presented first because of their simplicity and because Fourier transforms provide a simple derivation of the Bloch theorem.
In this appendix are given representative forms for the exchange–correlation energy and potential in the LSDA and GGA approximations. The forms given here are chosen because they are widely used and are relatively simple. Actual programs that provide energies and potentials for these and other forms can be found on-line (see Ch. 24).
Local spin density approximation (LSDA)
The local density approximation is based upon the exact expressions for the exchange energy, Eq. (5.15), and various approximations and fitting to numerical correlation energies for the homogeneous gas. Comparison of the forms is shown in Fig. 5.4. The first functions were the Wigner interpolation formula, Eq. (5.22), and the Hedin–Lundqvist [220] form; the latter is derived from many-body perturbation theory and is given below. As described in Ch. 5, the quantum Monte Carlo (QMC) calculations of Ceperley and Alder [297], and more recent work [298, 299, 303] provide essentially exact results for unpolarized and fully polarized cases. These results have been fitted to analytic forms for εc(rs), where rs is given by Eq. (5.1), leading to two widely used functionals due to Perdew and Zunger (PZ) [300] and Vosko, Wilkes, and Nusiar (VWN) [301], which are very similar quantitatively. Both functionals assume an interpolation form for fractional spin polarization, and Ortiz and Balone [298] report that their QMC calculations at intermediate polarization are somewhat better described by the VWN form.
The “TBPW” code is a modular code for Slater–Koster orthogonal tight-binding (TB, Sec. 14.4) and plane wave empirical pseudopotential (EPM, Sec. 12.6) calculations of electron energies in finite systems or bands in crystals. This appendix is a schematic description. A much more complete description, sample files, and the full source codes are available on-line at the site given in Ch. 24.
TBPW is meant to be an informative, instructional code that can bring out much of the physics of electronic structure. It has many of the main features of full density functional codes, but is much simpler. A main characteristic is that the codes are modular and are organized to separate the features common to all electronic structure calculations from the aspects that are specialized to a given method.
A sample input file is given below in Sec. N.4. The input file uses keywords that are recognized by the input routines, so that the same file can be used for either TB or PW calculations.
Density functional theory is the most widely used method today for electronic structure calculations because of the success of practical, approximate functionals. The crucial quantity in the Kohn–Sham approach is the exchange–correlation energy which is expressed as a functional of the density Exc[n]. This chapter is devoted to relevant approximate functionals, in particular, the local density approximation (LDA) and examples of generalized-gradient approximations (GGAs). Explicit formulas for certain widely used functionals are given in App. B. Non-local formulations are an active area of research leading to new classes of functionals, in particular, orbital-dependent functionals including the “optimized effective potential” (OEP) method and “hybrid functionals.” Important features are illustrated by a few selected results on atoms and molecules.
As emphasized in the previous chapter, the genius of the Kohn–Sham approach is two-fold: first, the construction of an auxiliary system leads to tractable independent-particle equations that hold the hope of solving interacting many-body problems. The famous Kohn–Sham equations are given in (7.11)–(7.13). Second, and perhaps more important, by explicitly separating out the independent-particle kinetic energy and the long-range Hartree terms, the remaining exchange–correlation functional Exc[n] can be reasonably approximated as a local or nearly local functional of the density. Even though the exact functional Exc[n] must be very complex, great progress has been made with remarkably simple approximations. This chapter is devoted to those approximations.
Excitation spectra reveal the properties of matter in terms of the response to time- or frequency-dependent perturbations. Particularly important examples are the dielectric function ε(ω) and the inverse function ε−1(ω) defined in App. E. The basic formulas relating the response to the electronic structure are rooted in perturbation theory and response functions (Sec. 3.7 and App. D). This chapter is devoted to dynamic response functions for electrons in self-consistent field methods, such as the Kohn–Sham approach, and to the alternative approach of solving directly the time-dependent Kohn–Sham equation to find the solution to all orders. The formal structure is based upon time-dependent density functional theory (TDDFT) (Sec. 7.6), which provides an exact framework in principle. In practice, simple approximations are remarkably successful in many cases and there is active research to develop new functionals.
As emphasized in the overview, Sec. 2.10, two types of excitations are of primary importance for electronic structure: excitations in which an electron is added or subtracted from the system, and excitations in which the number of electrons remains fixed. The former are of great interest as the “one-electron excitations” in an interacting many-body system; however, in independent-particle theories, such as the Hartree–Fock or Kohn–Sham approaches, these excitations are just the eigenstates of the independent-particle hamiltonian. In a crystal, the eigenvalues form well-defined bands εik with none of the renormalization and broadening that can only be included in a many-body theory.
Electric polarization is one of the basic quantities in physics, essential to the theory of dielectrics, effective charges in lattice dynamics, piezoelectricity, ferroelectricity, and other phenomena. However, descriptions in widely used texts are often based upon oversimplified models that are misleading or incorrect. The basic problem is that the expression for a dipole moment is ill defined in an extended system, and there is no unique way to find the moment as a sum of dipoles by “cutting” the charge density into finite regions. For extended matter such as crystals, a theory of polarization formulated directly in terms of the quantum mechanical wavefunction of the electrons has only recently been derived, with an elegant formulation in terms of a Berry's phase and alternative expressions using Wannier functions. The other essential property of insulators is “localization” of the electrons. Although the concept of localization is well known, recent theoretical advances have provided new quantitative approaches and demonstrated that localization is directly measurable by optical experiments. This chapter is closely related to Ch. 21 on Wannier functions, in particular to the gauge-invariant center of mass and contribution to the spread of Wannier functions ΩI of Sec. 21.3.
The theory of electrodynamics of matter [448, 790] (see App. E) is cast in terms of electric fields E(r, t) and currents j(r′, t′). (Here we ignore response to magnetic fields.) In metals, there are real currents and, in the static limit, electrons flow to screen all macroscopic electric fields.
If you don't like the answer, change the question.
Summary
Density functional theory is the most widely used method today for electronic structure calculations because of the approach proposed by Kohn and Sham in 1965: to replace the original many-body problem by an auxiliary independent-particle problem. This is an ansatz that, in principle, leads to exact calculations of properties of many-body systems using independent-particle methods; in practice, it has made possible approximate formulations that have proved to be remarkably successful. As a self-consistent method, the Kohn–Sham approach involves independent particles but an interacting density, an appreciation of which clarifies the way the method is used. The present chapter is devoted to the basic formulation of the Kohn–Sham approach and the ideas behind the crucial ingredient, the exchange–correlation energy functional Exc[n]. Information on approximate functionals in widespread use is deferred to Ch. 8, and methods for solution of the Kohn–Sham equations using the functionals are the subjects of Ch. 9 and much of the remainder of this tome.
Replacing one problem with another
The Kohn–Sham approach is to replace the difficult interacting many-body system obeying the hamiltonian (3.1) with a different auxiliary system that can be solved more easily. Since there is no unique prescription for choosing the simpler auxiliary system, this is an ansatz that rephrases the issues. The ansatz of Kohn and Sham assumes that the ground state density of the original interacting system is equal to that of some chosen non-interacting system.
It is not appropriate to summarize or conclude this volume on the basic theory and methods of electronic structure. The field is evolving rapidly with new advances in basic theory, algorithms, and computational methods. New developments and applications are opening unforeseen vistas. Volumes of information are now available on-line at thousands of sites.
It is more appropriate to provide a resource for information that will be updated in the future, on-line information available at a site maintained at the University of Illinois:
http://ElectronicStructure.org
A link is maintained at the Cambridge University Press site:
http://books.cambridge.org/0521782856.htm
Resources for materials computation are maintained by the Materials Computation Center at the University of Illinois, supported by the National Science Foundation:
http://www.mcc.uiuc.edu/
Additional sites include the Department of Physics, the electronic structure group at the University of Illinois, and the author's home page:
Additional material coordinated with descriptions in this book, as well as future updates, corrections, additions, and convenient feedback forms.
Information related to many-body methods beyond the scope of this volume.
Links to courses, tutorials, and codes maintained at the University of Illinois. Specific codes are meant for pedagogical use, teaching, or individual study, and are coordinated with descriptions in this book, for example, the general purpose empirical pseudopotential and tight-binding code (TBPW) in App. N.
Links to codes for electronic structure calculations. This will include resources at the Materials Computation Center and many other sites.
The concept of localization can be imbedded directly into the methods of electronic structure to create new algorithms that take advantage of locality or “nearsightedness” as coined by W. Kohn. As opposed to the textbook starting point for describing crystals in terms of extended Bloch eigenstates, many physical properties can be calculated from the density matrix ρ(r, r′), which is exponentially localized in an insulator or a metal at finite T. For large systems, this fact can be used to make “order-N” or O(N) methods where the computational time scales linearly in the size of the system. There are two aspects of the problem: “building” the hamiltonian and “solving” the equations. Here we emphasize the second aspect, which is more fundamental, and describe representative O(N) approaches that either treat ρ(r, r′) directly or work in terms of Wannier-like localized orbitals.
The reader should be aware (beware) that O(N) methods are under development; there are problems and shortcomings in actual practice.
Every textbook on solid state physics begins with the symmetry of crystals and the entire subject of electronic structure is cast in the framework of the eigenstates of the hamiltonian classified in k space by the Bloch theorem. So far this volume is no exception. However, the real goal is to understand the properties of materials from the fundamental theory of the electrons and this is not always the best approach, either for understanding or for calculations.