Lecture-Notes Solid State Theory
Lecture-Notes Solid State Theory
Lecture-Notes Solid State Theory
Manfred Sigrist Institut f ur Theoretische Physik HIT K23.8 Tel.: 044-633-2584 Email: [email protected] Website: http://www.itp.phys.ethz.ch/research/condmat/strong/
Literature: N.W. Ashcroft and N.D. Mermin: Solid State Physics, HRW International Editions, 1976. C. Kittel: Einf uhrung in die Festk orperphysik, R. Oldenburg Verlag, 1983. C. Kittel: Quantentheorie der Festk orper, R. Oldenburg, 1970. O. Madelung: Introduction to solid-state theory, Springer 1981; auch in Deutsch in drei B anden: Festk operphysik I-III, Springer. J.M. Ziman: Principles of the Theory of Solids, Cambridge University Press, London, 1972. M.P. Marder: Condensed Matter Physics, John Wiley & Sons, 2000. G. Grosso & G.P. Parravicini: Solid State Physics, Academic Press, 2000. G. Czychol: Theoretische Festk orperphysik, Springer 2004. P.L. Taylor & O. Heinonen, A Quantum Approach to Condensed Matter Physics, Cambridge Press 2002. G.D. Mahan, Condensed Matter in a Nutshell, Princeton University Press 2011. numerous specialized books.
Contents
Introduction 1 Electrons in the periodic crystal - band structure 1.1 Symmetries of crystals . . . . . . . . . . . . . . . . . . . . . . . 1.1.1 Space groups of crystals . . . . . . . . . . . . . . . . . . 1.1.2 Reciprocal lattice . . . . . . . . . . . . . . . . . . . . . . 1.2 Blochs theorem and Bloch functions . . . . . . . . . . . . . . . 1.3 Nearly free electron approximation . . . . . . . . . . . . . . . . 1.4 Tight-binding approximation . . . . . . . . . . . . . . . . . . . 1.4.1 Linear combination of atomic orbitals - LCAO . . . . . 1.4.2 Band structure of s-orbitals . . . . . . . . . . . . . . . . 1.4.3 Band structure of p-orbitals . . . . . . . . . . . . . . . . 1.4.4 Wannier functions . . . . . . . . . . . . . . . . . . . . . 1.4.5 Tight binding model in second quantization formulation 1.5 Symmetry properties of the band structure . . . . . . . . . . . 1.6 Band-lling and materials properties . . . . . . . . . . . . . . . 1.6.1 Electron count and band lling . . . . . . . . . . . . . . 1.6.2 Metals, semiconductors and insulators . . . . . . . . . . 1.7 Semi-classical description . . . . . . . . . . . . . . . . . . . . . 1.7.1 Equations of motion . . . . . . . . . . . . . . . . . . . . 1.7.2 Bloch oscillations . . . . . . . . . . . . . . . . . . . . . . 1.7.3 Current densities . . . . . . . . . . . . . . . . . . . . . . 2 Semiconductors 2.1 The band structure in group IV . . 2.1.1 Crystal and band structure 2.2 Elementary excitations . . . . . . . 2.2.1 Electron-hole excitations . . 2.2.2 Excitons . . . . . . . . . . . 2.2.3 Optical properties . . . . . 2.3 Doping semiconductors . . . . . . . 2.3.1 Impurity state . . . . . . . 2.3.2 Carrier concentration . . . 2.4 Semiconductor devices . . . . . . . 2.4.1 pn-contacts . . . . . . . . . 2.4.2 Diodes . . . . . . . . . . . . 2.4.3 MOSFET . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 8 8 8 10 11 12 15 15 16 17 19 21 21 23 24 24 27 27 28 28 35 36 36 38 39 40 42 44 44 45 46 46 47 48 51 51 52 54
3 Metals 3.1 The Jellium model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.1.1 Theory of metals - Sommerfeld and Pauli . . . . . . . . . . . . . . . . . . 3.1.2 Stability of metals - a Hartree-Fock approach . . . . . . . . . . . . . . . . 2
3.2
3.3
Charge excitations . . . . . . . . . . . . . . . . . . . . . 3.2.1 Dielectric response and Lindhard function . . . . 3.2.2 Electron-hole excitation . . . . . . . . . . . . . . 3.2.3 Collective excitation - Plasmon . . . . . . . . . . 3.2.4 Screening . . . . . . . . . . . . . . . . . . . . . . Phonons . . . . . . . . . . . . . . . . . . . . . . . . . . . 3.3.1 Vibration of a isotropic continuous medium . . . 3.3.2 Phonons in metals . . . . . . . . . . . . . . . . . 3.3.3 Peierls instability in one dimension . . . . . . . . 3.3.4 Dynamics of phonons and the dielectric function . . . . . . . . . . . . . . . gas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
57 57 59 60 62 66 66 68 69 73 76 76 76 78 79 80 82 83 89 93 93 96 98 99 101 102 103 105 105 108 109 111 111 113 113 115 118 119 119 121 122 125 126 128 128 130 132 132 135 136
4 Itinerant electrons in a magnetic eld 4.1 The de Haas-van Alphen eect . . . . . . . . . . 4.1.1 Landau levels . . . . . . . . . . . . . . . . 4.1.2 Oscillatory behavior of the magnetization 4.1.3 Onsager equation . . . . . . . . . . . . . . 4.2 Quantum Hall Eect . . . . . . . . . . . . . . . . 4.2.1 Hall eect of the two-dimensional electron 4.2.2 Integer Quantum Hall Eect . . . . . . . 4.2.3 Fractional Quantum Hall Eect . . . . . . 5 Landaus Theory of Fermi Liquids 5.1 Lifetime of quasiparticles . . . . . . . . . . 5.2 Phenomenological Theory of Fermi Liquids 5.2.1 Specic heat . . . . . . . . . . . . . 5.2.2 Compressibility . . . . . . . . . . . . 5.2.3 Spin susceptibility . . . . . . . . . . 5.2.4 Galilei invariance . . . . . . . . . . . 5.2.5 Stability of the Fermi liquid . . . . . 5.3 Microscopic considerations . . . . . . . . . . 5.3.1 Landau parameters . . . . . . . . . . 5.3.2 Distribution function . . . . . . . . . 5.3.3 Fermi liquid in one dimension? . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
6 Transport properties of metals 6.1 Electrical conductivity . . . . . . . . . . . . . . 6.2 Transport equations and relaxation time . . . . 6.2.1 The Boltzmann equation . . . . . . . . 6.2.2 The Drude form . . . . . . . . . . . . . 6.2.3 The relaxation time . . . . . . . . . . . 6.3 Impurity scattering . . . . . . . . . . . . . . . . 6.3.1 Potential scattering . . . . . . . . . . . 6.3.2 Kondo eect . . . . . . . . . . . . . . . 6.4 Electron-phonon interaction . . . . . . . . . . . 6.5 Electron-electron scattering . . . . . . . . . . . 6.6 Matthiessens rule and the Ioe-Regel limit . . 6.7 General transport coecients . . . . . . . . . . 6.7.1 Generalized Boltzmann equation . . . . 6.7.2 Thermoelectric eect . . . . . . . . . . . 6.8 Anderson localization . . . . . . . . . . . . . . 6.8.1 Landauer Formula for a single impurity 6.8.2 Scattering at two impurities . . . . . . . 6.8.3 Anderson localization . . . . . . . . . . 3
7 Magnetism in metals 7.1 Stoner instability . . . . . . . . . . . . . . . . . . . . . . . 7.1.1 Stoner model within the mean eld approximation 7.1.2 Stoner criterion . . . . . . . . . . . . . . . . . . . . 7.1.3 Spin susceptibility for T > TC . . . . . . . . . . . . 7.2 General spin susceptibility and magnetic instabilities . . . 7.2.1 General dynamic spin susceptibility . . . . . . . . 7.2.2 Instability with nite wave vector Q . . . . . . . . 7.2.3 Inuence of the band structure . . . . . . . . . . . 7.3 Stoner excitations . . . . . . . . . . . . . . . . . . . . . . 8 Magnetism of localized moments 8.1 Mott transition . . . . . . . . . . . . . . . . . . . . . . . . 8.1.1 Hubbard model . . . . . . . . . . . . . . . . . . . . 8.1.2 Insulating state . . . . . . . . . . . . . . . . . . . . 8.1.3 The metallic state . . . . . . . . . . . . . . . . . . 8.1.4 Fermi liquid properties of the metallic state . . . . 8.2 The Mott insulator as a quantum spin system . . . . . . . 8.2.1 The eective Hamiltonian . . . . . . . . . . . . . . 8.2.2 Mean eld approximation of the anti-ferromagnet . 8.3 Collective modes spin wave excitations . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . .
137 138 138 139 142 143 143 146 147 149 152 153 153 154 155 157 159 159 160 162
Introduction
Solid state physics (or condensed matter physics) is one of the most active and versatile branches of modern physics that have developed in the wake of the discovery of quantum mechanics. It deals with problems concerning the properties of materials and, more generally, systems with many degrees of freedom, ranging from fundamental questions to technological applications. This richness of topics has turned solid state physics into the largest subeld of physics; furthermore, it has arguably contributed most to technological development in industrialized countries.
Figure 1: Atom cores and the surrounding electrons. Condensed matter (solid bodies) consists of atomic nuclei (ions), usually arranged in a regular (elastic) lattice, and of electrons (see Figure 1). As the macroscopic behavior of a solid is determined by the dynamics of these constituents, the description of the system requires the use of quantum mechanics. Thus, we introduce the Hamiltonian describing nuclei and electrons, H = He + Hn + Hne , with He =
i
(1)
p2 1 i + 2m 2
2 Pj
i= i
e2 , |r i r i | Zj Zj e2 , |R j R j | (2)
Hn =
j
2 Mj
1 2
j =j
Hne =
i,j
Zj e2 , |r i R j |
where He (Hn ) describes the dynamics of the electrons (nuclei) and their mutual interaction and Hne includes the interaction between ions and electrons. The parameters appearing are m e Mj Zj free electron mass elementary charge mass of j -th nucleus atomic (charge) number of j -th nucleus 9.1094 1031 kg 1.6022 1019 As 103 104 m
The characteristic scales known from atomic and molecular systems are 5
with the ne structure constant = e2 / c = 1/137. The energy scale of one Hartree is much less than the (relativistic) rest mass of an electron ( 0.5MeV), which in turn is considered small in particle physics. In fact, in high-energy physics even physics at the Planck scale is considered, at least theoretically. The Planck scale is an energy scale so large that even gravity is thought to be aected by quantum eects, as EPlanck = c2 c 1019 GeV, G lPlanck = G 1.6 1035 m, c3 (3)
where G = 6.673 1011 m3 kg1 s2 is the gravitational constant. This is the realm of the GUT (grand unied theory) and string theory. The goal is not to provide a better description of electrons or atomic cores, but to nd the most fundamental theory of physics. solid state physics 10 meV metals semiconductors magnets superconductors ferroelectrics ...... eective models known and established Figure 2: Energy scales in physics. In contrast, in solid state physics we are dealing with phenomena occurring at room temperature (T 300K) or below, i.e., at characteristic energies of about E kB T 0.03eV = 30meV, which is even much smaller than the energy scale of one Hartree. Correspondingly, the important length scales are given by the extension of the system or of the electronic wave functions. The focus is thus quite dierent from the one of high-energy physics. There, a highly successful phenomenological theory for low energies, the so-called standard model, exists, whereas the underlying theory for higher energies is unknown. In solid state physics, the situation is reversed. The Hamiltonian (1) describes the known high-energy physics (on the energy scale of Hartree), and one aims at describing the low-energy properties using reduced (eective, phenomenological) theories. Both tasks are far from trivial. Among the various states of condensed matter that solid state theory seeks to describe are metals, semiconductors, and insulators. Furthermore, there are phenomena such as magnetism, superconductivity, ferroelectricity, charge ordering, and the quantum Hall eect. All of these states share a common origin: Electrons interacting among themselves and with the ions through the Coulomb interaction. More often than not, the microscopic formulation in (1) is too complicated to allow an understanding of the low-energy behavior from rst principles. Consequently, the formulation of eective (reduced) theories is an important step in condensed matter theory. On the one hand, characterizing the ground state of a system is an important goal in itself. However, measurable quantities are determined by excited states, so that the concept of elementary excitations takes on a central role. Some celebrated examples are Landaus quasiparticles for 6 10 eV electrons, cores atom high-energy physics astrophysics and cosmology 1 MeV phenomenological particle physics standard model GUT string theory M-theory
Fermi liquids, the phonons connected to lattice vibrations, and magnons in ferromagnets. The idea is to treat the ground state as an eective vacuum in the sense of second quantization, with the elementary excitations as particles on that vacuum. Depending on the system, the vacuum may be the Fermi sea or some state with a broken symmetry, like a ferromagnet, a superconductor, or the crystal lattice itself. According to P. W. Anderson,1 the description of the properties of materials rests on two principles: The principle of adiabatic continuity and the principle of spontaneously broken symmetry. By adiabatic continuity we mean that complicated systems may be replaced by simpler systems that have the same essential properties in the sense that the two systems may be adiabatically deformed into each other without changing qualitative properties. Arguably the most impressive example is Landaus Fermi liquid theory mentioned above. The low-energy properties of strongly interacting electrons are the same as those of non-interacting fermions with renormalized parameters. On the other hand, phase transitions into states with qualitatively dierent properties can often be characterized by broken symmetries. In magnetically ordered states the rotational symmetry and the time-reversal invariance are broken, whereas in the superconducting state the global gauge symmetry is. In many cases the violation of a symmetry is a guiding principle which helps to simplify the theoretical description considerably. Moreover, in recent years some systems have been recognized as having topological order which may be considered as a further principle to characterize low-energy states of matter. A famous example for this is found in the context of the Quantum Hall eect. The goal of these lectures is to introduce these basic concepts on which virtually all more elaborate methods are building up. In the course of this, we will cover a wide range of frequently encountered ground states, starting with the theory of metals and semiconductors, proceeding with magnets, Mott insulators, and nally superconductors.
P.W. Anderson: Basic Notions of Condensed Matter Physics, Frontiers in Physics Lecture Notes Series, Addison-Wesley (1984).
Chapter 1
1.1
1.1.1
Symmetries of crystals
Space groups of crystals
Most solids consist of a regular lattice of atoms with perfectly repeating structures. The minimal repeating unit of such a lattice is the unit cell. The symmetries of a crystal are contained in the space group R, a group of symmetry operations (translations, rotations, the inversion or combinations) under which the crystal is left invariant. In three dimensions, there are 230 dierent space groups1 (cf. Table 1.1).
All symmetry transformations form together a set which has the properties of a group. A group G combined with a multipliation has the following properties: the product of two elements of G is also in G : a, b G multiplications are associative: a (b c) = (a b) c. a unit element e G exists with: e a = a e = a for all a G . for every element a G there is an inverse a1 G with a1 a = a a1 = e. A group with a b = b a for all pairs of element is called Abelian group, otherwise it is non-Abelian. A subset G G is called a subgroup of G , if it is a group as well. Guides to group theory in the context of condensed matter physics can be found in the textbooks - Mildred S. Dresselhaus, Gene Dresselhaus and Ado Jorio: Group Theory - Application to the Physics of Condensed Matter - Peter Y. Yu and Manuel Cardona: Fundamentals of Semiconductors, Springer. a b = c G.
1
We consider here a crystal displayed as a point lattice, each point either symbolizing an atom or a whole unit cell (Fig.1.1 for a square lattice). Translations in the space group are represented by linear combinations of a basic set of translation vectors {ai } (primitive lattice vectors) connecting lattice points. Any translation that maps the lattice onto itself is a linear combination of the {ai } with integer coecients, Rn = n1 a1 + n2 a2 + n3 a3 . (1.1)
a2 a1
Figure 1.1: Crystal point lattice in two dimensions: the vectors a1 and a2 form the basic set of translations and R = 2a1 + a2 in the gure. The shaded area is the Wigner-Seitz cell which is obtained by drawing perpendicular lines (planes in three dimensions) through the center of all lines connecting neighboring lattice points. The Wigner-Seitz cell also constitutes a unit cell of the lattice. General symmetry transformations including general elements of the space group may be written in the notation due to Wigner, r = g r + a = {g |a}r , (1.2)
where g represents a rotation, reection or inversion with respect to lattice points, axes or planes. The elements g form the generating point group P . In three dimensions there are 32 point groups. We distinguish the following basic symmetry operations: basic translations rotations, reections, inversions screw axes, glide planes {E |a}, {g |0}, {g |a},
where E is the unit element (identity) of P . A screw axis is a symmetry operation of a rotation followed by a translation along the rotation axis. A glide plane is a symmetry operation with reection at a plane followed by a translation along the same plane. The symmetry operations {g |a}, together with the associative multiplication {g |a}{g |a } = {gg |g a + a} (1.3)
form a group with unit element {E |0}. In general, these groups are non-Abelian, i.e., the group elements do not commute with each other. However, there is always an Abelian subgroup of R, the group of translations {E |a}. The elements g P do not necessarily form a subgroup, 9
because some of these elements (e.g., screw axes or glide planes) leave the lattice invariant only in combination with a translation. Nevertheless, the relation {g |a}{E |a }{g |a}1 = {E |g a } (1.4) (1.5)
holds generally. If P is a subgroup of R, then R is said to be symmorphic. In this case, the space group contains only primitive translations {E |a} and neither screw axes nor glide planes. The 14 Bravais lattices2 are symmorphic. Among the 230 space groups 73 are symmorphic and 157 are non-symmorphic. crystal system
(# point groups, # space groups)
point groups
Sch onies symbols
triclinic (2,2) monoclinic (3,13) orthorhombic (3,59) tetragonal (7,68) trigonal (5,25) hexagonal (7,27) cubic (5, 36)
C1 , C 1 C2 , Cs , C2h D2 , C2v , D2h C4 , S4 , C4h , D4 , C4v , D2d , D4h C3 , S6 , D3 , C3v , D3d C6 , C3h , C6h , D6 , C6v , D3h , D6h T, Th , O, Td , Oh
Table 1.1: List of the point and space groups for each crystal system in three dimensions.
1.1.2
Reciprocal lattice
We dene now the reciprocal lattice which is of importance for the electron band structure and x-ray diraction on a periodic lattice. The reciprocal lattice is also perfectly periodic with a translation symmetry with a basic set {bi } dening arbitrary reciprocal lattice vectors as Gm = m1 b1 + m2 b2 + m3 b3 , where mi are integers and ai bj = 2ij , such that bi = 2 aj ak ai (aj ak ) i, j = 1, 2, 3 , bj bk . bi (bj bk ) (1.7) (1.8) (1.6)
and ai = 2
The reciprocal lattice of a simple cubic lattice is simple cubic. However, a body centered cubic (bcc) lattice has a face centered cubic (fcc) reciprocal lattice and vice versa (see slides). It follows that any real space lattice vector Rn and any reciprocal lattice vector Gm have the property that Gm Rn = 2 (m1 n1 + m2 n2 + m3 n3 ) = 2N (1.9)
2
Crystal systems, crystals, Bravais lattices are discussed in more detail in - Czycholl, Theoretische Festk orperphysik, Springer
10
with N being an integer. This allows us to expand any function f (r ) periodically in the real lattice as fG eiGr (1.10) f (r ) =
G
(1.11)
where the integral runs over the unit cell of the periodic lattice with the volume UC . Finally, we dene the (rst) Brillouin zone as the Wigner-Seitz cell constructed in the reciprocal lattice (see Fig.1.1 and 1.2).
1.2
We consider a Hamiltonian H of electrons invariant under a discrete set of lattice translations {{E |a}}, a symmetry introduced by a periodic potential. This implies that the corresponding translation operator Ta on the Hilbert space commutes with the Hamiltonian H = He + Hie (purely electronic Hamiltonian He , interaction between electrons and ions Hie ), [Ta , H] = 0. (1.12)
This translation operator is dened through Ta |r = |r + a . Neglecting the interactions among electrons, which would be contained in He , we are left with a single particle problem H H0 = p2 + V (r ), 2m (1.13)
Vion (r Rj ),
(1.14)
describes the potential landscape of the single particle in the ionic background. With Rj being the position of the j -th ion, the potential V (r ) is by construction periodic, with V (r + a) = V (r ) for all lattice vectors a, and represents Hie . Therefore, H0 commutes with Ta . For a Hamiltonian H0 commuting with the translation operator Ta , the eigenstates of H0 are simultaneously eigenstates of Ta . Blochs theorem states that the eigenvalues of Ta lie on the unit circle of the complex plane, which ensures that these states are extended. This means Ta (r ) = (r a) = a (r ) ,
l Tla (r ) = (r la) = Ta (r ) = l a (r )
(1.15)
(1.16) (1.17)
requiring This condition is satised if we express the wave function as product of a plane wave eikr and a periodic Bloch function uk (r ) 1 n,k (r ) = eikr un,k (r ) . 11 (1.18) |a | = 1 a = eia .
The integer n is a quantum number called band index, k is the pseudo-momentum (wave vector) and represents the volume of the system. Note that the eigenvalue of n,k (r ) with respect to Ta , eika , implies periodicity in the reciprocal space, the k-space, because ei(k+G)a = eika for all reciprocal lattice vectors G. We may, therefore, restrict k to the rst Brillouin zone. Blochs theorem simplies the initial problem to the so-called Bloch equation for the periodic function uk , (p + k)2 + V (r ) uk (r ) = 2m
k uk (r ),
(1.22)
where we suppress the band index to simplify the notation. This equation follows from the relation p eikr uk (r ) = eikr (p + k)uk (r ), (1.23)
which can be used for more complex forms of the Hamiltonian as well. There are various numerical methods which allow to compute rather eciently the band energies k for a given Hamiltonian H.
1.3
We start here from the limit of free electrons assuming the periodic potential V (r ) is weak. Using Eqs.(1.10) and (1.11) we expand the periodic potential, V (r ) =
G
VG eiGr , 1 UC
UC
(1.24) (1.25)
VG =
d3 r V (r )eiGr .
The potential is real and we assume it also to be invariant under inversion (V (r ) = V (r )), =V leading to VG = V G . Note that the uniform component V0 corresponds to an irrelevant G energy shift and may be set to zero. Because of its periodicity, the Bloch function uk (r ) is expressed in the same way, uk (r ) =
G
cG eiGr ,
(1.26)
where the coecients cG = cG (k) are functions of k, in general. Inserting this Ansatz and the expansion (1.24) into the Bloch equation, (1.22), we obtain a linear eigenvalue problem for the band energies k ,
2
2m
(k G)2
cG +
G
VG G cG = 0.
(1.27)
This represents an eigenvalue problem in innite dimensions with eigenvectors cG (k) and eigenvalues k as band energies. These k include corrections to the bare parabolic dispersion, (0) 2 k2 /2m, due to the potential V (r ). Obviously, the dispersion (0) is naturally parabolic k = k 12
in absence of the potential V (r ) whereby the eigenstates would be simply plane waves. As a lowest order approach we obtain the approximative energy spectrum within the rst Brillouin zone, considering only all parabolic bands of the type k (G) = 2 (k G)2 /2m centered around the reciprocal wave vectors G (see dashed line in Fig.1.2)). We illustrate here the nearly free electron method using the case of a one-dimensional lattice. 2 G2 /2m (taking G Assuming that the periodic modulation of the potential is weak, say |VG | as a characteristic wavevector in the periodic system), the problem (1.27) can be simplied. Let us start with the lowest energy values around the center of the rst Brillouin zone, i.e. k 0 (|k | /a). For the lowest energy eigenvalue we solve Eq.(1.27) by 1 for G = 0 cG (1.28) 2mVG 1 for G = 0 2 {(k G)2 k 2 } leading to the energy eigenvalue
k 2 k2
2m
[(k G=0 2m E0 =
G)2 k 2 ] |G |2
2
|VG |2
2 k2
2m
+ E0
(1.29)
with
G=0
2m
G2
(1.30)
and
1 1 = 14 m m
G=0
|G |2
(1.31)
with G = VG /{ 2 G2 /2m} (|G | 1). We observe that k is parabolic with a slightly modied (eective) mass, m m > m. Note that this result resembles the lowest order corrections in the Rayleigh-Schr odinger perturbation theory for a non-degenerate state. This solution corresponds to the lowest branch of the band structure within this approach (see Fig.1.2). The parabolic approximation of the band structure at a symmetry point with an eective mass m , is a standard way to approximate band tops or bottoms. It is called k p-approximation (k-dot-p). We stay at the zone center and address the next eigenstates which are dominated by the parabola which cross for k = 0 at a value 2 G 2 /2m. Restricting originating from G = 2/a = G ourselves to these two components we obtain the two-dimensional eigenvalue equation system, 2 2 c V2G 2m (k G) k G = 0 . (1.32) 2 cG )2 k V (k + G
2G 2m
The eigenvalues are obtained through the secular equation, 2 )2 k (k G V2 G det 2m 2 2 V2G 2m (k + G) leading to
k
=0,
k
(1.33)
1 )2 + (k G )2 (k + G 2 2m
2
2m
)2 (k G )2 (k + G
2 + 4|V2G |
2m
2 |V | + G 2G
2m
k2 (1.34) 13
where m 1. We observe a energy band gap separating two bands | + > 0 and m < 0 as |G with opposite curvature (see Fig.1.2). Note that the curvature diverges for V2G 0 as m 0. The wavefunctions at k = 0 are given by ()
1 1 1 1 2|G | m m
(1.35)
cG =
V2 G |V2G |
c+G = c+G
()
()
(1.36)
where we have chosen V2G to be real and positive. Thus, uk=0 (x) = sin Gx for
k=0,+
(1.37)
k=0,
cos Gx for
one being even and the other odd under parity operation x x. A similar analysis can be done at the boundary of the rst Brillouin zone where two energy parabolas cross. For example at k = /a we nd the two dominant contributions originate from G = +2/a and 2/a, respectively, together with G = 0. Also here the energy eigenvalues show a band gap with parabolic bands centered at k = /a (boundary of the rst Brillouin zone in one-dimension) with positive and negative eective mass (see Fig.1.2). Analogous as for the band center we can distinguish the wavefunction with even and odd parity for the two bands at the Brillouin zone boundary. Indeed every crossing energy parabola centered around dierent reciprocal lattice points contributes to a band gap. By construction we can extend the band structure beyond the rst Brillouin zone and nd a periodic energy spectrum with
k +G
k k
(1.38) =
k
where G is a reciprocal lattice vector. Moreover, we nd in Fig.1.2 that as well as time reversal symmetry, like for free electrons.
due to parity
band gap