Using the above resource, show allowed energy states of
a particle in this potential form bands and gaps.
LUL V(x) -b 0 a FIGURE 1. THE KRONIG-PENNEY POTENTIAL. it is usually developed in a way that makes it fairly abstruse to most chemistry under- graduate students. The Kronig-Penney model is represented by the one-dimensional periodic potential shown in Figure 1. Even though the model is one-dimensional, it is the periodicity of the potential that is the crucial property that yields electronic band structure. The mathematical form of the repeating unit of the potential is V(x) = V. -b <x< 0 (1) = 0 0<x<a As shown in Figure 1, the potential has a period of c = a +b. The Schrödinger equation for this model is dºt , 2mE + V = 0 0<x<a (2a) dx2 72 + dy 2m(E - V) dx2 y = 0 - b<x< 0 (2b) P2 The coefficients in both of these equations are constants, so the solutions are similar to that of a particle in a box. The solutions to Equations 2 are V = Alax + Be-iax 0<x<a (3a) and V2 = CeBx + De-Bx -b<x< 0 (3b)
where Q= (2mE)/2 ħ (4) and [2m(Vo - E)]1/2 B = (5) ħ are real quantities if we assume for the time being that E < V. The fact that the potential is periodic imposes a restriction on the form of the wave function y(x). Define an operator Î that shifts y(x) by a distance c. The resulting wave function y(x + c) must be a multiple of y(x), and so we write Ty(x) = ty(x) = y(x +c) (6) where t is a constant that is of magnitude unity. We can apply this translation operator repeatedly to obtain Î"y(x) = +"y(x) = y(x + nc) (7) At this point we shall apply a periodic boundary condition that pictures the one-dimen- sional lattice to be a ring of N atoms, so that y(x) = y(x + Nc). Clearly if N is large, this condition has no physical consequence. This (standard) boundary condition is called the Born-von Karmen boundary condition [3]. If we apply this boundary condition to Equation 7, then we have ÎNy(x) = +"y(x) = y(x + Nc) = y(x) (8) or that tº = 1. Thus we find that t must be one of the N roots of unity, or that = t = exil/N 1 = 0, 1, 2, ..., N - 1 (9) Defining k by k = 2^l/Nc, t becomes t = elke and so we can write = = V(x + c) = elke y(x) (10)
Now, any function y(x) satisfying this equation can be written in the form V(x) = u(x)eike (11) where u(x) is a periodic function with period c (u(x + c) = u(x)). To see that this must be so, we write = (x + c) = u(x + c)eľk(x+C) = eikeu(x + c)eikx = eikeu(x)eikx =e = eike y(x) (12) as in Equation 10. Equation 11 is a fundamental result of solid-state physics and is known as Bloch's theorem [4]. Any wave function for a periodic potential of period c must necessarily have that form. The solutions to Equations 2 yield four integration constants, and so we require four boundary conditions to determine them. Two of these boundary conditions can be taken to be that y(x) and its first derivative must be continuous across any discontinuity in the potential, and we shall impose this condition at the point x = 0. If V1(x) denotes y (x) for –b < x < 0 and 42(x) denotes 42(x) for 0 < x < a, then these two boundary conditions read V1(0) = 4200 and V10 = V20) (13) where the prime denotes the first derivative. Furthermore, from Bloch's theorem, u(x) = u(x + a) and so not only is u(x) periodic, but its first derivative must also be periodic. Therefore, we have the two conditions u(-5) = u(a) and u'(-b) = u'(a) (14) These four boundary conditions are sufficient to determine y(x) completely. Equations 13 give the two equations A+B=C+D BC - BD = ia A-iaB (15)
Ael(-k) + Be-ila+k)a and Equations 14 give Ce-(-ik)b + De Btik)b (16) ila – k) Ael(a-k)a – i (a + k) Be-1(a+k)a = (B - ik)Ce-(-ik)b – (B + ik) De(B+ik)b In order for there to be a nontrivial solution to the above four equations for A, B, C, and D, the determinant of their coefficients must vanish. This gives us the determinantal equation = 1 1 1 1 B -B ia -ia elik-B). elik+B) pla(-k) piw(a+k) (B – ik)elik-Bb -(B+ ik)elik+B)b i la – k)ela(a–k) -i(a + k)eia(@+k) = 0 (17) It is straightforward, albeit lengthy, to show that Equation 17 is equivalent to B2-a? -sinh(Bb) sin(aa) - cosh(Bb) cos(aa) = cos(a + b)k (18) 2αβ (The step-by-step procedure is outlined in the Appendix for completeness, although it is certainly not really necessary to carry it out in class, or anywhere else for that matter.) Equation 18 may appear to be hopeless, but it can be cast into a more transparent form. From Equations 4 and 5 we have 2m V. E 2m V a? (19) 12 V -E R2 and B2 2m(Vo - E) h2 2m V. -(1 - E) h2 (20) where ε = E/V. Therefore B2 - 2 2αβ 11 = 1-28 2€1/2(1 - €)1/2 (21)
Furthermore, 1/2 Bb = (20%)" act2 =rA(1 – €)12 (22) and Vo 1/2 aa = a81/2 Ag1/2 (23) p2 = where r = b/a and A = (2m Vo/ti2)1/2a are parameters whose values are fixed by the potential. Using Equations 21 through 23, Equation 18 becomes (1 - 28) sinh[r A(1 - €)1/2] sin(A&1/2) + (24) 2€1/2(1 - €)1/2 cosh[r A(1 – €)1/2] cos(Aɛ1/2) = cos[ka(1+r)] << 1 = Equation 24 is valid for the case in which € < 1. For the case in which € > 1, we replace the factor (1 – €)1/2 by i(€ – 1)1/2 and use the identities sinh iz = i sin z and cosh iz = cos z for € > 1. For € > 1, Equation 18 becomes (1 – 2€) (25) 281/2(€ – 1)1/2 sin[r A(E – 1)1/2]sin(A&1/2) + cos[r A(€ – 1)1/2] cos(Aɛ1/2) = cos[ka(1+r)] € > 1 At the point ε = 1 both Equations 24 and 25 give TA COS A - sin A = cos[ka(1+r)] E = 1 (26) 2 For notational convenience, we shall write Equations 24 through 26 collectively as FE) = cos[ka(1+r)] (27) where F(8) is given by the left-hand sides of Equations 24 for € < 1 and by the left-hand side of Equation 25 for ε > 1. Equation 27 gives ε as a function of ka for all values of ε for fixed values of r and A, or, in other words, it gives e(ka), at least in principle, for a fixed potential. Although it is not possible to solve Equation 27 analytically, it can
2.0 in 0.0 F(8) -2.0 -4.0 0.0 1.0 3.0 4.0 2.0 E FIGURE 2. A PLOT OF F(C), THE LEFT-HAND SIDE OF EQUATION 27, AGAINST . THE ONLY ALLOWED VALUES OF E ARE THOSE FOR WHICH F(8) LIES BETWEEN +1 AND -1. be easily solved graphically. The function F(E), which is just a function of ε for a fixed potential, is plotted against e in Figure 2 for A = 10 and r A = 2. The value of cos[ka(1+r)], the right-hand side of Equation 27, must lie between +1 and -1; consequently, the only allowed values of ε are those for which F(E), the left- hand side of Equation 27 lies between +1 and -1. The two horizontal lines in Figure 2 represent the two extreme values of cos[ka(1+r)]. The only allowed values of e are those for which F(8) lies between the two horizontal lines. These values are given in Table 1. As either Figure 2 or Table 1 indicates, the allowed values of ε fall into continuous regions separated by gaps. This distribution of allowed energies illustrates the band structure of crystalline solids. A solid will act as a conductor, an insulator, or a semiconductor, depending upon the electronic population of the bands and their separations. It is customary to present the results in Table 1 as a plot of ε against the wave vector k, as shown in Figure 3. To see how this plot is obtained, let's invert Equation 27 and write k a cosF(E). OC
E 1 Зл ka(1+r) 0 AT 2π 40 FIGURE 3. A PLOT OF ENERGY E AGAINST WAVE NUMBER K SHOWING THE BAND GAPS THAT APPEAR WHEN KA(1 + R)= N1. THE DASHED LINE IS THE PARABOLIC RESULT FOR FREE ELECTRONS. According to Figure 2 and Table 1, for values of ε from 0 to 0.0616, F(E) is greater than one and so there is no real value of k ( or really ka(1+r)) that satisfies Equation 27. At ε = 0.0616, k = 0. As e varies from 0.0616 to 0.0744, cos[ka(1+r)] varies from +1 to -1, and so ka(1+r) = p. As e varies from 0.0744 to 0.2425, there is no real = TABLE 1. The allowed values of s = E/V, for the potential shown in Figure 1 with r = b/a = 0.20 and A = (2mV)/2/h = 10. band range of 1 2 3 4 0.0616-0.0744 0.2425-0.2965 0.5325-0.6652 0.9192-1.1776 1.3980-1.8305 1.9791-2.6190 2.6850– 5 6 7
= value of k that satisfies Equation 27. At ε = 0.2425, ka(1+r) once again equals 1. At ε = 0.2965, ka(1+r) now equals 21, and another band gap appears in Figure 3. Continuing this argument, we get the series of discontinuities in the plot of e against k (or really ka(1+r)) as shown in Figure 3. It is instructive to compare the result shown in Figure 3 with that of a free electron. We can obtain this result by letting Vo = 0 in Equation 18. In that case, B = ia and so (B2 – 2%)/2aB = i. Using the fact that sinh iab = i sin ab and that cosh iab = cos ab, Equation 18 becomes - sin ab sinaa + cos ab cosaa = cos a (a+b) = cosk(a+b) = (28) or that a = k. Therefore, from Equation 4 we find that ti?k2 E= (29) 2m for a free particle. Equation 29 is shown as the parabola in Figure 3. The other extreme limit is to let V, be large, in which case we should obtain the energies of a particle in a box. If V, is large, then ε is small, so that the factor multiplying the first term in Equation 24 becomes 1/281/2, which is a large quantity. Therefore, the first term in Equation 24 dominates the second term and we have erA 41/2 sin A£1/2 = cos[ka(1+r)] (30) The factor in front of sin A£1/2 approaches infinity for large values of V, and so sin Aɛ1/2 must approach zero in order that the value of left-hand side of Equation 30 lie between +1 and -1. In the limit, then, sin Aɛ1/2, and so Aɛ1/2 = Inn, or, using the definition of A, n'h2 E= 8ma? (31) which represent the allowed energies of a particle in a box.
V2 - Vi a 0
LUL V(x) -b 0 a FIGURE 1. THE KRONIG-PENNEY POTENTIAL. it is usually developed in a way that makes it fairly abstruse to
-
answerhappygod
- Site Admin
- Posts: 899604
- Joined: Mon Aug 02, 2021 8:13 am
LUL V(x) -b 0 a FIGURE 1. THE KRONIG-PENNEY POTENTIAL. it is usually developed in a way that makes it fairly abstruse to
Join a community of subject matter experts. Register for FREE to view solutions, replies, and use search function. Request answer by replying!