Chapter 6 Free-Electrons in Solids: Fermi Gas The free electron model is the simplest possible model for a metal. We shall make two fundamental assumptions that are much more reasonable than they look: (i) the interactions between the electrons are irrelevant and can be ignored (the independent electron approximation); (ii) the electrons move in a constant potential and we can ignore everything about the structure of the material. These are heroic assumptions but they do give a fair representation of a simple metal. For the exam you should • be familiar with the free electron model where the potential energy of the electrons is zero and the electron-electron interactions are ignored. • be familiar with the solutions to the Schrοinger equation for the free electron model and know that the dispersion relation in this case. • be able to calculate the density of states for free electrons D(k) and D(E) in 1, 2, and 3 dimensions. • be familiar with the fermi function. • be able to start from the density of states and calculate the thermodynamic properties such as the heat capacity of electrons. 1 Contents 6.1 Fermi gas and free electron model....................................................................................................3 6.2 Free electrons in one Dimension.......................................................................................................5 6.2.1 Energy Levels in One Dimension ..........................................................................................5 6.2.2 Dispersion relations of one dimension...................................................................................6 6.2.3 Fermi Energy in One Dimension ...........................................................................................6 6.2.4 Density of states in one dimension ........................................................................................7 6.3 The Fermi-Dirac distribution ..........................................................................................................10 6.3.1 Fermi-Dirac distribution at T=0 Kelvin ...............................................................................10 6.3.2 Fermi-Dirac distribution at T>0 Kelvin ...............................................................................11 6.3.3 Temperature-dependent chemical potential .........................................................................12 6.4 Electrons in three dimensions .........................................................................................................14 6.4.1 Energy levels in 3D..............................................................................................................14 6.4.2 Dispersion relations in 3D....................................................................................................14 6.4.3 Density of states in 3D case .................................................................................................15 6.4.4 Fermi energy in 3D case ......................................................................................................17 6.4.5 Example ...............................................................................................................................19 6.5 Electron heat capacity .....................................................................................................................20 6.5.1 Qualitative Representation of Electron Specific Heat..........................................................20 6.5.2 Experimental Result of Electron Specific Heat....................................................................21 6.5.3 Quantitative Representation of Electron Specific Heat........................................................22 References.............................................................................................................................................24 A: 2-D Density of States ...............................................................................................................25 B: 0-D Density of States ...............................................................................................................29 Questions ..............................................................................................................................................30 2 6.1 Fermi gas and free electron model What is a Fermi gas? Fermi gas is a gas that consists of particles with half-integral spin and obeys Fermi-Dirac statistics. Fermi gas that consists of noninteracting particles is called ideal Fermi gas. Fermi gases include electrons in metals and semiconductors, electrons in atoms with large atomic numbers, nucleons in heavy atomic nuclei, and gases consisting of quasiparticles with half-integral spin. What is the Free electron (Fermi) gas model? Free electron (Fermi) gas is the simplest model of a metal that was proposed by Fermi. For this model, we make the following assumptions: (i) The crystal comprises a fixed background of N identical positively charge nuclei and N electrons, which can move freely inside the crystal without seeing any of the nuclei (monovalent case); and (ii) Coulomb interactions are negligible because the system is neutral overall. A free electron model is the simplest way to represent the electronic structure of metals. Although the free electron model is a great oversimplification of the reality, surprisingly in many cases it works pretty well, so that it is able to describe many important properties of metals. According to this model, the valence electrons of the constituent atoms of the crystal become conduction electrons and travel freely throughout the crystal. Therefore, within this model we neglect the interaction of conduction electrons with ions of the lattice and the interaction between the conduction electrons. In this sense we are talking about a free electron gas. However, there is a principle difference between the free electron gas and ordinary gas of molecules. First, electrons are charged particles. Therefore, in order to maintain the charge neutrality of the whole crystal, we need to include positive ions. This is done within the jelly model, according to which the positive charge of ions is smeared out uniformly throughout the crystal. This positive background maintains the charge neutrality but does not exert any field on the electrons. Ions form a uniform jelly into which electrons move. Second important property of the free electron gas is that it should meet the Pauli exclusion principle, which leads to important consequences. Figure 1. Schematic illustration of the free electron model. 3 This model works relatively well for alkali metals (Group 1 elements), such as Na, K, Rb and Cs. We would like to understand why electrons are only weakly scattered as they migrate through a metal. This model transfers the ideas of electron orbitals in atoms into a macroscopic object. Thus, the fundamental behaviour of a metal comes from the Pauli exclusion principle. For now, we will ignore the crystal lattice. By a free electron Fermi gas, we mean a gas of free electrons subject to the Pauli principle. One of the successes of free electron model is the verification of the Heat capacity functional form correct. 4 6.2 Free electrons in one Dimension 6.2.1 Energy Levels in One Dimension Consider a free electron gas in one dimension. An electron of mass m is confined to a length L by infinite barriers. We consider first a free electron gas in one dimension. The wavefunction ψn(x) of the electron is a solution of the Schrödinger equation Hψn(x) = Enψn(x), where En is the energy of electron orbital. Since we can assume that the potential lies at zero, the Hamiltonian H includes only the kinetic energy so that Hψ n ( x) = − h2 d 2 ψ n ( x) = E nψ n ( x) . 2m dx 2 Note that this is a one-electron equation, which means that we neglect the electron-electron interactions. We use the term orbital to describe the solution of this equation. Since the ψn(x) is a continuous function and is equal to zero beyond the length L, the boundary conditions for the wave function are ψn(0) =ψn(L). The solution of Eq.(6.1) is therefore nπ x L ψ n ( x) = A sin where A is a constant and n is an integer. Substituting (6.2) into (6.1) we obtain for the eigenvalues 2 h 2 nπ E n ( x) = . 2m L These solutions correspond to standing waves with a different number of nodes within the potential well as is shown in Fig.1. Figure 2. First three energy levels and wave-functions of a free electron of mass m confined to a 5 line of length L. The energy levels are labeled according to the quantum number n which gives the number of half-wavelengths in the wavefunction. The wavelengths are indicated on the wavefunctions. 6.2.2 Dispersion relations of one dimension The energy of the confined electron is related to its wavevector k. E n (k ) = h2k 2 2m So the dispersion relation of the confined electron is displayed beloe. Figure 3. Dispersion relation for one-dimensional electron system. 6.2.3 Fermi Energy in One Dimension Now we need to accommodate N valence electrons in these quantum states. According to the Pauli exclusion principle no two electrons can have their quantum number identical. That is, each electronic quantum state can be occupied by at most one electron. The electronic state in a 1D solid is characterized by two quantum numbers that are n and ms, where n describes the orbital ψn(x), and ms describes the projection of the spin momentum on a quantization axis. Electron spin is equal to S=1/2, so that there (2S+1)=2 possible spin states with ms = ±½. Therefore, each orbital labeled by the quantum number n can accommodate two electrons, one with spin up and one with spin down orientation. Let nF denote the highest filled energy level, where we start filling the levels from the bottom (n = 1) and continue filling higher levels with electrons until all N electrons are accommodated. It is convenient to suppose that N is an even number. The condition 2nF = N determines nF, the value of n for the uppermost 6 filled level. The energy of the highest occupied level is called the Fermi energy EF. For the one dimensional system of N electrons we find, using Eq. (6.3), 2 h 2 Nπ EF = . 2m 2 L In metals the value of the Fermi energy is of the order of 5 eV. The ground state of the N electron system is illustrated in Fig.x: All the electronic levels are filled up to the Fermi energy. All the levels above are empty. Figure 6. Occupation of energy levels according to the Pauli exclusion principle. 6.2.4 Density of states in one dimension To calculate various optical properties such as the rate of absorption or emission and how electrons and holes distribute themselves within a solid, we need to know the number of available states per unit volume per unit energy. We first calculate the available states in k-space and then use the energy-momentum relation in parabolic bands to give the density of states in terms of energy. In one dimensional case, we have derived its discrete wavevector kn kn = πn L This equation suggests that there is one quantum state within the momentum interval π/L. Consequently the number of quantum states dN in the momentum interval dk is dN = dk L = dk π /L π Multiplying by 2 for spin degeneracy, the number of states within the interval dk in 1-d, dN = 2 × dk 2L = dk π /L π For one dimensional quantum wires, we have h2k 2 E= 2m 7 dk = m dE 2h 2 E Therefore the density of states per unit length in 1-d and multiplying by 2 for spin degeneracy, dN = 2L dk = π DOS = 2L π m dE 2h 2 E dN L 2m 1 = dE π h 2 E For one dimension, the density of states per unit volume at energy E is given by DOS = dN L = dE π 2m 1 h2 E Figure 5. DOS of one dimensional system. Using more than the first energy level, the density of states function becomes DOS = L π 2 m ni H ( E − E i ) ∑ E−E h2 i i where once again, H(E-Ei) is the Heaviside function and ni is the degeneracy factor. For quantum structures with dimensions lower than 2, it is possible for the same energy level to occur for more than one arrangement of confined states. To account for this, a second factor n i (E) is introduced. 8 Figure 6. Density of states for quantum wire (1D green) which has multiple states. 9 6.3 The Fermi-Dirac distribution 6.3.1 Fermi-Dirac distribution at T=0 Kelvin Electrons are Fermions. Two electrons can occupy a state, one with spin up and one with spin down. Fermi function f(E) is is the probability that a state at energy E is occupied. Electrons obey the Fermi-Dirac distribution f (E) = 1 E−µ + 1 exp k T B , where µ is the chemical potential. At a temperature T = 0°K, an ideal Fermi gas is in the ground state, and the particles of the gas occupy all quantum states with energy up to a certain maximum value, which depends on the gas density and is called the Fermi level (EF). Quantum states with energy E > EF are empty; such a distribution of occupied and empty states corresponds to complete degeneracy of the Fermi gas. This is the ground state of the N electron system at absolute zero. Figure 7. Fermi distribution at zero temperature. It follows from the preceding discussion that the distribution function for electrons at T = 0°K has the form 10 1, E < µ . f (E) = 0, E > µ That is, all levels below EF are completely filled, and all those above EF are completely empty. This function is plotted in Fig. 2 (b), which shows the discontinuity at the Fermi energy. 6.3.2 Fermi-Dirac distribution at T>0 Kelvin What happens if the temperature is increased? The kinetic energy of the electron gas increases with temperature. Therefore, some energy levels become occupied which were vacant at zero temperature, and some levels become vacant which were occupied at absolute zero. The distribution of electrons among the levels is usually described by the distribution function, f(E), which is defined as the probability that the level E is occupied by an electron. Thus if the level is certainly empty, then, f(E) = 0, while if it is certainly full, then f(E) = 1. In general, f(E) has a value between zero and unity. Figure 8. Fermi-Dirac distribution function at various temperatures, for EF/kB = 50000 K. The results apply to a gas in three dimensions. The total number of particles is constant, independent of temperature. The chemical potential at each temperature may be read from the graph as the energy at which f = 0:5. At T ≠ 0°K, the mean occupation number for a quantum state of an ideal Fermi gas is described by the Fermi-Dirac distribution function. When the system is heated (T>0°K), thermal energy excites the electrons. However, all the electrons do not share this energy equally, as would be the case in the classical treatment, because the electrons lying well below the Fermi level EF cannot absorb energy. If they did so, they would move to a higher level, which would be already occupied, and hence the exclusion principle 11 would be violated. Recall in this context that the energy which an electron may absorb thermally is of the order kBT (=0.025 eV at room temperature), which is much smaller than EF, this being of the order of 5 eV. Therefore only those electrons close to the Fermi level can be excited, because the levels above EF are empty, and hence when those electrons move to a higher level there is no violation of the exclusion principle. Thus only these electrons - which are a small fraction of the total number - are capable of being thermally excited. The distribution function at non-zero temperature is given by the Fermi distribution function. The Fermi distribution function determines the probability that an orbital of energy E is occupied at thermal equilibrium. This function is also plotted in Fig. 2(b), which shows that it is substantially the same as the distribution at T = 0°K, except very close to the Fermi level, where some of the electrons are excited from below EF to above it. The quantity m is called the chemical potential. The chemical potential can be determined in a way that the total number of electrons in the system is equal to N. At absolute zero µ = EF . 6.3.3 Temperature-dependent chemical potential Figure 9. Variation of the chemical potential with temperature for free electron Fermi gases in one and three dimensions. A Fermi level also exists for a nonideal Fermi gas, although the particles of such a gas do not occupy specific quantum states. In a nonideal Fermi gas consisting of electrons in a metal, the formation of pairs of correlated electrons (the Cooper effect) and the transition of the metal to the superconducting state may 12 occur at very low temperatures because of the attraction of electrons with equal but oppositely directed momenta and spins. A Fermi gas consisting of electrons in heavy atoms is described by the Thomas-Fermi model. 13 6.4 Electrons in three dimensions 6.4.1 Energy levels in 3D The Schrödinger equation in the three dimensions takes the form h2 ∂2 ∂2 ∂2 Hψ n ( x, y, z ) = − + + 2m ∂x 2 ∂y 2 ∂z 2 ψ n ( x, y, z ) = E nψ n ( x, y, z ) If the electrons are confined to a cube of edge L, the solution is the standing wave nπ nπ nπ x sin y sin z L L L ψ n ( x, y, z ) = A sin where nx, ny, and nz are positive integers. In many cases, however, it convenient to introduce periodic boundary conditions, as we did for lattice vibrations. The advantage of this description is that we assume that our crystal is infinite and disregard the influence of the outer boundaries of the crystal on the solution. We require then that our wavefunction is periodic in x, y, and z directions with period L, so that ψ ( x + L, y , z ) = ψ ( x , y , z ) and similarly for the y and z coordinates. The solution of the Schrödinger equation (6.7) which satisfies these boundary conditions has the form of the traveling plane wave: r r r ψ k (r ) = A exp(ik • r ) provided that the component of the wavevector k are determined from kx = πn x L ,ky = πn y L , kz = πn z L where nx, ny, and nz are positive or negative integers. 6.4.2 Dispersion relations in 3D If we now substitute this solution to Eq.(6.7) we obtain for the energies of the orbital with the wavevector k h2 2 Ek = ( k x + k y2 + k z2 ) 2m The dispersion relation is Ek = h2 2 k 2m 14 Figure10. Dispersion relation for 3D. 6.4.3 Density of states in 3D case In three dimensional case, wavevector k has discrete values. Figure 11. Visualisation of k-space showing values of k as points. The number of allowed states is the number of these points contained in the shell of radius k and thickness dk. 15 The diagram shows the spherical shells. The problem of finding the number of allowed states amounts to finding the number of these allowed states between spheres of radius k and k +d k . In 3-dimensions, the volume between the two shells is given by: r dk = dk x dk y dk z = 4πk 2 dk Therefore, the number of states is given simply by dividing this volume by the volume of a single energy state. At this point is convenient to introduce an additional factor of two to account for the intrinsic angular momentum of the electrons or spin states. Figure 12. Discrete k values for 3D kx = 2πn y 2πn x 2πn z ,ky = ,kz = L L L This equation suggests that there is one quantum state within the momentum interval (2π/L)3. Consequently the number of quantum states dN in the momentum interval dk is r dk L 3 r L dN = = ( ) dk = ( ) 3 4πk 2 dk 3 2π 2π (2π / L) Multiplying by 2 for spin degeneracy, the number of states within the interval dk, dN = 2 × ( L 3 ) 4πk 2 dk 2π For 3 dimensional quantum case, we have 16 E= dk = h2k 2 2m m dE 2h 2 E Therefore the density of states per unit length in 1-d and multiplying by 2 for spin degeneracy, 3 L 2m 2 1 dN = 4π × ( ) 3 2 E 2 dE 2π h For 3 dimension, the density of states per unit volume at energy E is given by 3 dN L 2m 2 1 DOS = = 4π × ( ) 3 2 E 2 dE 2π h Figure 13 Density of states of three dimensional case. Figure 14. Density of states for bulk (3D blue). 6.4.4 Fermi energy in 3D case From the equation of DOS, we have 17 3 L 2m 2 1 dN = 4π × ( ) 3 2 E 2 dE 2π h 3 ∞ L 2m 2 1 L3 N = ∫ dN = ∫ 4π × ( ) 3 2 E 2 dE = 2π h 2π 2 0 [( 2 h2 2 N 3 h2 3nπ 2 EF = 3π 3 = 2m 2m L ) 1 3 2m 2 h 3 2 ], 2 where n is the electron density. Comparison to the equation. E= h2 2 k 2m leads to the Fermi wavevector k F = (3nπ 2 ) 3 1 Because the linear momentum of an electron is r r r p = mυ = hk , We can derive the Fermi velocity for an electron as υF = ( hk F h = 3nπ 2 m m ) 1/ 3 An alternative way to understand the Fermi energy is list as below. 18 EF L3 E dE = ∫0 3π 2 1 2 3 3 2m 2 E F2 h 2 Figure 15. Fermi surface. N = 2× 4πk F2 / 3 (2π / L) 3 We can get the same expressions of the Fermi wavevector and Fermi energy level. 6.4.5 Example Calculate the Fermi energy, Fermi wavevector, Fermi velocity and Fermi temperature for metal Li, whose electron density is 4.70×1022 cm-3. Solution: E F = [ 1 h2 ( 3nπ 2 ) 3 2m ] = 4.32eV , 2 where n is the electron density. The Fermi wavevector is given by k F = (3nπ 2 ) 3 = (3 × 4.7 × 10 28 × 3.14 2 )1 / 3 = 1.12 × 1010 m −1 = 1.12 × 10 8 cm −1 1 We can derive the Fermi velocity for an electron as υF = hk F 1.05 × 10 −34 × 1.12 × 1010 = 1.29 × 10 6 m / s = −30 m 0.91 × 10 The Fermi temperature is given by TF = E F / k = 4.32 × 1.6 × 10 −19 / 1.38 × 10 −23 = 5.0086 × 10 4 K Example 2. Calculate the number of states per unit energy in a 100 by 100 by 10 nm piece of silicon (m* = 1.08 m0) 100 meV above the conduction band edge. Write the result in units of eV-1. Solution: The density of states equals: 3 dN 1 2m 2 1 DOS = = ( 2 ) 2 E 2 dE 2π h DOS = ( 3 8π 2 1 )m * 2 ( E − E c ) 2 = 1.51 × 10 56 m −3 J −1 3 h So that the total number of states per unit energy equals: DOS = 1.51 × 10 56 × 10 22 J −1 = 2.41 × 10 5 eV −1 19 6.5 Electron heat capacity 6.5.1 Qualitative Representation of Electron Specific Heat At T=0 , the Fermi-Dirac distribution function drops precisely at E=EF. At low temperatures, the function remains roughly symmetric around E=EF. The number of electrons promoted is proportional to the area between the T=0 and finite temperature curves, which is proportional to T. The average increase of their energy is also proportional to T. So the finite-T system has an excess energy proportional to T2 , giving CV ∝T. Figure 16. Fermi-Dirac distribution for various temperatures, where T F ≡ F /kB = 50,000K . This result applies at a gas in three dimensions. The chemical potential at each temperature is the value at f=0.5. Only the electrons near the Fermi surface are excited by thermal energy kT. The number of excited electrons is roughly of the order of N′ ≈ N kT EF 20 The energy absorbed by the electrons is U ≈ N ′kT ≈ N (kT ) 2 EF The specific heat capacity of electrons is given by the following equation C el = ∂U T ≈ Nk ∂T TF In general, the Fermi temperature is about 50,000 K. For classic thermodynamics, the specific heat of a particle is C classic = ∂U 3 = Nk . ∂T 2 So the electronic heat of electrons is two orders less than expected in classic thermodynamics. Because T/TF ~0.01, therefore usually electron specific is much smaller than phonon specific heat. In general, C = C el + C ph = γT + AT 3 Cel becomes more important at very low T. 6.5.2 Experimental Result of Electron Specific Heat We expect the total specific heat to be of the form: C = C el + C ph = γT + AT 3 Dividing both sides by T gives: where the latter (phonon) term is from the Debye model. So if we plot CV/T as a function of T2, the intercept gives γ, whilst the gradient gives A. 21 Figure 17. The specific heat of metal potassium at low temperatures. 6.5.3 Quantitative Representation of Electron Specific Heat A major success of Fermi's consideration of electrons in metals is that it could explain the “excess” heat capacity in metals (beyond that given by the phonons). Combine DOS and thermal distribution, the total electron energy is given by the sum over all energies: ∞ U = ∫ D( E ) f ( E , T ) EdE 0 For T = 0, the Fermi function becomes a step, so we can define EF implicitly as: ∞ EF 0 0 N = ∫ D( E ) f ( E , T )dE = ∫ D( E )dE We need a few approximations. In general, T≪TF, where TF is the Fermi temperature. In other words, kBT ≪EF. In this case,T<<TF we can approximate N as in equation. This also allows us to approximate the occupation number as: f (E) ≈ 1 exp[(E − µ ) / k B T ] + 1 ( ) Let x = ( E − E F ) / k B T , then f = 1 / e + 1 . x Multiplying equation (6.17) by EF and subtract from equation (6.16) gives: ∞ ∆U = U − NE F = ∫ (E − E F )D( E ) F ( E )dE 0 Taking the partial differential over equation (6.19) with respect to temperature, and noticing that N and EF are independent of temperature gives an expression for the heat capacity: 22 ∞ CV = ∂U ∂F = ∫ (E − E F )D( E ) dE ∂T 0 ∂T where (E − E F ) ) . ∂f ex = x 2 ∂T (1 + e ) k BT 2 This only becomes larger near E = EF, corresponding to −1<x<1 , so we can approximate D(E)≈D(EF) = const. Thus ∞ (E − EF )2 ex dE x 2 k BT 2 0 (1 + e ) CV = D( E F ) ∫ And using the definition of x: dx=dE/kBT , gives: ∞ ∫ CV = k B2TD ( E F ) − EF k BT ex x2 dE (1 + e x ) 2 But, E F >> k B T , leads to − EF → −∞ , so the integral becomes k BT ∞ π2 x 2e x = dx ∫ x 2 3 − ∞ (e + 1) Thus the specific heat becomes: CV = π2 3 k B TD( E F ) = π2 2E F Nk 2B T To derive above equation we should have DOS expressed in N and E. As proven in previous problems, D( E ) = AE 1 / 2 , 3 NE F−3 / 2 . 2 3 D ( E ) = NE F−3 / 2 E 1 / 2 2 where A = D( E F ) = 3N 2E F 23 References 1 http://lamp.tu-graz.ac.at/~hadley/ss1/fermigas/fermigas.php 2 For DOS, see http://ecee.colorado.edu/~bart/book/book/chapter2/ch2_4.htm#fig2_4_4 3 For DOS, see http://www.tf.uni-kiel.de/matwis/amat/semi_en/kap_2/basics/b2_1_8.html 4 For semiconductors, see http://www.tf.uni-kiel.de/matwis/amat/semi_en/index.html 5 For DOS, see http://ecee.colorado.edu/~bart/book/book/chapter2/ch2_4.htm#fig2_4_4 6 For DOS, see http://en.wikipedia.org/wiki/Density_of_states (End of this page) 24 Appendix In solid-state and condensed matter physics, the density of states (DOS) of a system describes the number of states per interval of energy at each energy level that are available to be occupied. Unlike isolated systems, like atoms or molecules in gas phase, the density distributions are not discrete like a spectral density but continuous. A high DOS at a specific energy level means that there are many states available for occupation. A DOS of zero means that no states can be occupied at that energy level. In general a DOS is an average over the space and time domains occupied by the system. Local variations, most often due to distortions of the original system, are often called local density of states (LDOS). If the DOS of a undisturbed system is zero, the LDOS can locally be non-zero due to the presence of a local potential. A: 2-D Density of States Energy levels in 2D The Schrödinger equation in the two dimensions takes the form Hψ n ( x, y ) = − h2 ∂2 ∂2 2 + 2 2m ∂x ∂y ψ n ( x, y ) = E nψ n ( x, y ) If the electrons are confined to a rectangle of edge L, the solution is the standing wave nπ nπ x sin y L L ψ n ( x, y ) = A sin where nx and ny are positive integers. In many cases, however, it convenient to introduce periodic boundary conditions, as we did for lattice vibrations. The advantage of this description is that we assume that our crystal is infinite and disregard the influence of the outer boundaries of the crystal on the solution. We require then that our wavefunction is periodic in x, y, and z directions with period L, so that ψ ( x + L, y ) = ψ ( x , y ) and similarly for the y and y coordinates. The solution of the Schrödinger equation (6.7) which satisfies these boundary conditions has the form of the traveling plane wave: r r r ψ k (r ) = A exp(ik • r ) provided that the component of the wavevector k are determined from kx = 2πn y 2πn x ,ky = L L 25 where nx and ny are positive or negative integers. Dispersion relations in 2D If we now substitute this solution to Eq.(6.7) we obtain for the energies of the orbital with the wavevector k Ek = h2 2 (k x + k y2 ) 2m The dispersion relation is Ek = h2 2 k 2m Figure 18. Dispersion relation for 2D Density of states in 2D case In two dimensional case, wavevector k has discrete values. 26 Figure 19. k-space in 2D. The density of states at an energy E is the number of k-states per unit volume contained with the annulus of radius k and thickness dk. The diagram shows the spherical shells. The problem of finding the number of allowed states amounts to finding the number of these allowed states between areas of radius k and k +d k . In 2-dimensions, the area between the two shells is given by: r dk = dk x dk y = 2πkdk Therefore, the number of states is given simply by dividing this volume by the volume of a single energy state. At this point is convenient to introduce an additional factor of two to account for the intrinsic angular momentum of the electrons or spin states. kx = 2πn y 2πn x ,ky = L L This equation suggests that there is one quantum state occupies the area (π/L)2 in the momentum space. Consequently the number of quantum states dN in the momentum interval dk is r dk L 2 r L dN = = ( ) dk = ( ) 2 2πkdk 2 2π 2π (2π / L) Multiplying by 2 for spin degeneracy, the number of states within the interval dk, dN = 2 × ( L 2 ) 2πkdk 2π For 2 dimensional quantum case, we have E= h2k 2 2m 27 kdk = m dE h2 Therefore the density of states per unit area in 2-d and multiplying by 2 for spin degeneracy, dN = 4π × ( 1 2 m ) dE 2π h 2 For 2 dimension, the density of states per unit area at energy E is given by DOS = dN 1 m = dE π h 2 It is significant that the 2-d density of states does not depend on energy. Immediately, as the top of the energy-gap is reached, there is a significant number of available states. Taking into account the other energy levels in the quantum well, the density of states takes on a staircase like function given by: DOS = dN 1 m = ∑ H ( E − Ei )dE dE π h 2 i where H(E-Ei) is the Heaviside function. It takes the value of zero when E is less than Ei and 1, when E is equal to or greater than Ei . Ei is the i-th energy level within the quantum well. Figure 20. Density of states for bulk (3D blue), quantum well (2D red), quantum wire (1D green) and quantum Dot (0D black) 28 B: 0-D Density of States Finally we consider the density of states in a zero-dimensional system, the quantum box. No free motion is possible in such a quantum box, since the electron is confined in all three spatial dimensions. Consequently there is no k-space available which could be filled up with electrons. Each quantum state of a 0D system can therefore be occupied by two electrons. The density of states is therefore described by a delta function. DOS ( E ) = 2δ ( E − E 0 ) For more than one quantum state the density of states is given by DOS ( E ) = ∑ 2δ ( E − E n ) n Figure 21. Electronic density of states of semiconductors with 3, 2, 1, and 0 degrees of freedom for electron propagation. Systems with 2, 1, and 0 degrees of freedom are referred to as quantum wells, quantum wires, and quantum boxes, respectively. The density of states for one quantized level for 3D, 2D, 1D, and 0D electron system are schematically illustrated in Fig. 21. 29 Questions 3 NE F . 5 1 Calculate the kinetic energy of N electrons in 3D system zero temperature U 0 = 2 Calculate the Fermi energy of a metal. What is the relationship between the electron density and the Fermi energy in the free electron model? Explain how this relationship would be derived. Hint: The volume in k-space corresponding to one allowed k-state is ((2π)/L)³, where L is the length of the crystal. The total volume of the crystal is L³. There are twice as many electrons states as k-states due to spin. 3 Derive the relationship connecting pressure and volume of an electron gas at 0 K. The result may be written as p = 2(U0/V)/3. Show that the bulk modulus B = -V(dp/dV) of an electron gas at T = 0 is B = 5p/3 = 10U0/9V. 4 Average electron energy in 2-D. For a metal at a temperature of T = 0 K, the conduction electrons at the bottom of the band have an energy E = 0 and the conduction electrons with the highest energy have an energy E = EF. For free electrons in two-dimensions, what is the average energy of the conduction electrons? 5 Find the density of states as a function of energy for a non-interacting free electron gas in two dimensions. For this system it is possible to find an analytic expression for the temperature dependence of the chemical potential. Show that πnh 2 mkT µ (T ) = kT ln exp E − 1 = kT ln exp F − 1 , kT here n is the number of electrons per unit area. 6 Thermodynamic properties of a metal. A monovalent metal has a simple cubic Bravais lattice and a lattice constant of a = 0.15 nm. Calculate the chemical potential, the specific heat, the entropy, and the Helmholtz free energy of the electrons at temperatures of 10 K and 300 K assuming that the free electron model can be used. 7 Derive the density of states D(k) in two dimensions. The density of states for a free electron gas in two dimensions is, D ( E ) = m /(h 2 π ) . What is the heat capacity for a free electron gas in two dimensions? 30

© Copyright 2021 Paperzz