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.
Classification of crystals and their excitations by symmetry is a general approach applicable to electronic states, vibrational states, and other properties. The first part of this chapter deals with translational symmetry which has the same universal form in all crystals, and which leads to the Bloch theorem that rigorously classifies excitations by their crystal momentum. (The discussion here follows Ashcroft and Mermin, [84], Chs. 4–8.) The other relevant symmetries are time reversal and point symmetries. The latter depend upon a specific crystal structure and are treated only briefly. Detailed classification can be found in many texts, and computer programs that deal with the symmetries can be found on-line at sites listed in Ch. 24.
Structures of crystals: lattice + basis
A crystal is an ordered state of matter in which the positions of the nuclei (and consequently all properties) are repeated periodically in space. It is completely specified by the types and positions of the nuclei in one repeat unit (primitive unit cell), and the rules that describe the repetition (translations).
• The positions and types of atoms in the primitive cell are called the basis. The set of translations, which generates the entire periodic crystal by repeating the basis, is a lattice of points in space called the Bravais lattice. Specification of the crystal can be summarized as:
The subject of this chapter is the role of plane waves and grids in modern electronic structure calculations, which builds upon the basic formulation of Ch. 12. Plane waves have played an important role from the early OPW calculations to widely used methods involving norm-conserving pseudopotentials. Plane waves continue to be the basis of choice for many new developments, such as quantum molecular dynamics simulations (Ch. 18), owing to the simplicity of operations. Efficient iterative methods (App. M) have made it feasible to apply plane waves to large systems, and recently developed approaches such as “ultrasoft” pseudopotentials and projector augmented waves (PAWs Ch. 11) have made it feasible to apply plane waves to difficult cases such as materials containing transition metals. Real-space grids are an intrinsic part of efficient planewave calculations and there is a growing development of real-space methods, including multigrids, finite elements, wavelets, etc.
Basic Schrödinger-like equations for eigenstates expanded in a plane wave basis can be found in Sec. 12.1 and related equations for real-space grids in Sec. 12.8. These methods are appropriate in cases where the potentials and wavefunctions are smooth. Thus application of these methods to real materials means that they must be combined with a transformation to remove the core states, such as OPWs, pseudopotentials, or PAWs (Ch. 11). Many aspects of pseudopotential calculations have been given in Sec. 12.6.
The fundamental tenet of density functional theory is that any property of a system of many interacting particles can be viewed as a functional of the ground state density n0(r); that is, one scalar function of position n0(r), in principle, determines all the information in the many-body wavefunctions for the ground state and all excited states. The existence proofs for such functionals, given in the original works of Hohenberg and Kohn and of Mermin, are disarmingly simple. However, they provide no guidance whatsoever for constructing the functionals, and no exact functionals are known for any system of more than one electron. Density functional theory (DFT)would remain a minor curiosity today if it were not for the ansatz made by Kohn and Sham, which has provided a way to make useful, approximate ground state functionals for real systems of many electrons. The subject of this chapter is density functional theory as a methodology for many-body systems; Ch. 7 describes the Kohn–Sham ansatz that replaces the interacting problem with an auxiliary independent-particle problem with all many-body effects included in an exchange–correlation functional; Ch. 8 deals with widely used approximations for the exchange–correlation functional; and Ch. 9 is devoted to solution of the Kohn–Sham independent-particle equations in a general form useful for all Kohn–Sham calculations. Following chapters in this volume are devoted to algorithms for actual calculations, and applications to problems in atomic, molecular, and condensed matter physics.
Many properties of materials – mechanical, electrostatic, magnetic, thermal, etc. – are determined by the variations of the total energy around the equilibrium configuration, defined by formulas such as (2.2)–(2.7). Experimentally, vast amounts of information about materials are garnered from studies of vibration spectra, magnetic excitations, and other responses to experimental probes. This chapter is devoted to the role of electronic structure in providing predictions and understanding of such properties, through the total energy and force methods described in previous chapters, as well as recent advances in efficient methods for calculation of response functions themselves. Through these developments, calculation of full phonon dispersion curves, dielectric functions, infrared activity, Raman scattering intensities, magnons, anharmonic energies to all orders, phase transitions, and many other properties have been brought into the fold of practical electronic structure theory.
The primary properties considered in previous chapters are the total energy and (generalized) forces. These are sufficient to treat a vast array of problems including stability of structures, phase transitions, surfaces and interfaces, spin polarization, “ab initio” molecular dynamics, etc. One can also use such direct methods to calculate all the derivatives of the energy with respect to perturbations, by carrying out full self-consistent calculations for various values of the perturbation, and extracting derivatives from finite difference formulas. This has been used very successfully, for example, in the “frozen phonon” method illustrated in Fig. 2.8 and described further in Sec. 19.2.
The subject of this appendix is formulations and explicit equations for the total energy that properly take into account the long-range effects of Coulomb interactions. We emphasize the Kohn–Sham independent-particle equations and expressions for total energy; however, the ideas and many of the equations also apply to many-body calculations. There are three main issues:
Identifying various convenient expressions that each yield properly the intrinsic total energy per formula unit for an extended bulk system.
Understanding and calculating the effect upon the average potential in a bulk material due to dipole terms at surfaces and interfaces.
Treating finite systems, where there is no essential difficulty, but where it is convenient to carry out the calculations in a periodic “supercell” geometry.
Basic issues
There is a simple set of guiding principles that must be followed to properly treat long-range Coulomb interactions in extended systems. If the calculations are carried out in a cell that represents an infinite system, i.e. the unit cell of a crystal, or a “supercell” constructed so that its limiting behavior represents a macroscopic system, then:
The cell must be chosen to be neutral;
The neutral cell can be used to define a proper thermodynamic “reference state” if in addition we require that there is no average (macroscopic) electric field;
Theoretical analysis of the electronic structure of matter provides understanding and quantitative methods that describe the great variety of phenomena observed. A list of these phenomena reads like the contents of a textbook on condensed matter physics, which naturally divides into ground state and excited state electronic properties. The aim of this chapter is to provide an introduction to electronic structure without recourse to mathematical formulas; the purpose is to lay out the role of electrons in determining the properties of matter and to present an overview of the challenges for electronic structure theory.
The properties of matter naturally fall into two categories determined, respectively, by the electronic ground state and by electronic excited states. This distinction is evident in the physical properties of materials and also determines the framework for theoretical understanding and development of the entire field of electronic structure. In essence, the list of ground state and excited state electronic properties is the same in most textbooks [84, 86, 88] on condensed matter physics:
Ground state: cohesive energy, equilibrium crystal structure, phase transitions between structures, elastic constants, charge density, magnetic order, static dielectric and magnetic susceptibilities, nuclear vibrations and motion (in the adiabatic approximation), and many other properties.
Excited states: low-energy excitations in metals involved in specific heat, Pauli spin susceptibility, transport, etc; higher energy excitations that determine insulating gaps in insulators, optical properties, spectra for adding or removing electrons, and many other properties.
A density is a field defined at each position r, for example the particle number density n(r), which is a well-defined, experimentally measurable function. It would be desirable to have expressions for other densities, in particular, energy and stress densities. However, energy and stress densities are not unique on a microscopic quantum scale, even though they are the basis of the theory of elasticity on a macroscopic scale. This appendix brings out three points: (1) certain integrals of energy and stress densities are unique and very useful; (2) there are important contributions to the energy or stress density that are completely unique – these include all terms that arise from the fact that electrons are a many-body system of fermions; (3) all other terms that are non-unique can be shown to involve only the single scalar number density – there are different possible choices for these terms, each involving only derivatives of the density n(r) or the classical Coulomb potential VCC(r) which is directly related to n(r). It follows that all the issues of non-uniqueness are exactly the same as in a one-particle problem.
Only one density is widely used in electronic structure – the particle density n(r). It is the fundamental measurable quantity in quantum mechanics and the fundamental density in density functional theory. Theoretical expressions for n(r) are well defined and lead to unique results. Here we emphasize that other densities have the potential to play a useful role in electronic structure theory.
Localized functions afford a satisfying description of electronic structure and bonding in an intuitive localized picture. They are widely used in chemistry and have been revived in recent years in physics for efficiency in large simulations, especially “order-N” methods (Ch. 23). The semi-empirical tight-binding method is particularly simple and instructive since the basis is not explicitly specified and one needs only the matrix elements of the overlap and the hamiltonian. This chapter starts with a definition of the problem of electronic structure in terms of localized orbitals, and considers various illustrative examples in the tight-binding approach. Many of the concepts and forms carry over to full calculations with localized functions that are the subject of the following chapter, Ch. 15.
The hallmark of the approaches considered in this chapter and the next is that the wavefunction is expanded in a linear combination of fixed energy-independent orbitals, each associated with a specific atom in the molecule or crystal. For example, the linear combination of atomic orbitals (LCAO) formulation denotes a basis of atomic or modified atomic-like functions. Such a basis provides a natural, physically motivated description of electronic states in materials; in fact, possibly the first theory of electrons in a crystal was the tight-binding method developed by Bloch [36] in 1928. The history of this approach is summarized nicely by Slater and Koster [589], who point out that the seminal work of Bloch considered only the simplest s-symmetry function and the first to consider a basis of different atomic orbitals were Jones, Mott, and Skinner [594] in 1934.
The fundamental idea of a “pseudopotential” is the replacement of one problem with another. The primary application in electronic structure is to replace the strong Coulomb potential of the nucleus and the effects of the tightly bound core electrons by an effective ionic potential acting on the valence electrons. A pseudopotential can be generated in an atomic calculation and then used to compute properties of valence electrons in molecules or solids, since the core states remain almost unchanged. Furthermore, the fact that pseudopotentials are not unique allows the freedom to choose forms that simplify the calculations and the interpretation of the resulting electronic structure. The advent of “ab initio norm-conserving” and “ultrasoft” pseudopotentials has led to accurate calculations that are the basis for much of the current research and development of new methods in electronic structure, as described in the following chapters.
Many of the ideas originated in the orthogonalized plane wave (OPW) approach that casts the eigenvalue problem in terms of a smooth part of the valence functions plus core (or core-like) functions. The OPW method has been brought into the modern framework of total energy functionals by the projector augmented wave (PAW) approach that uses pseudopotential operators but keeps the full core wavefunctions.
Scattering amplitudes and pseudopotentials
The scattering properties of a localized spherical potential at any energy ε can be formulated concisely in terms of the phase shift ηl(ε), which determines the cross-section and all properties of the wavefunction outside the localized region.
The subject of this appendix is the macroscopic stress that enters mechanical properties of matter in the form of stress–strain relations. The stress tensor is the generalization of pressure to all the independent components of dilation and shear, and the “stress theorem” is the generalization of the virial theorem for scalar pressure to all components of the stress tensor. In condensed matter, the state of the system is specified by the forces on each atom and the stress, which is an independent variable. The conditions for equilibrium are: (1) that the total force vanishes on each atom, and (2) that the macroscopic stress equals the externally applied stress.
Macroscopic stress and strain
Stress and strain are important concepts in characterizing the states of condensed matter [177, 721, 722, 890]. A body is in a state of stress if it is acted upon by external forces or if one part of the body exerts forces upon another part. If we consider two types of forces as illustrated in Fig. G.1: those acting interior to a volume element and those that act upon (or through) the surface of the element due to the surrounding material, which are shown as arrows in the figure. The latter forces (per unit area) are the stresses transmitted throughout the interior of the volume. Since these forces balance on any surface in equilibrium, the stress can be determined in terms of only the intrinsic internal forces; i.e. stress is an intrinsic property of a material in a given state.