Energy Bands: Kronig-Penney Model

[torehan]

In Kronig-Penney Model we seek solutions for time independent Scrödinger equation inside and outside of the periodic potential.

for 0 < x < a;​

V(x) =

V₀ for -b < x < 0​

My question is why solution of the KP Scrödinger equation in the region of [tex] -b < x < 0 [/tex] doesnot have complex exponentials as in the region of [tex] 0 < x < a [/tex] ?

 

[Niles]

In this region we have
[tex]
\frac{{d^2 \psi }}{{dx^2 }} = k^2 \psi,
[/tex]

where

[tex]
k^2 = \frac{{2m\left( {V_0 - E} \right)}}{{\hbar ^2 }}.
[/tex]
The solution to this 2nd order ODE is a linear combination of exponentials (with real arguments). If you think about it physically, the wavefunction will also decay in this area, which can only happen if the argument of the exponential is real.

 

[torehan]

In this region we have
[tex]
\frac{{d^2 \psi }}{{dx^2 }} = k^2 \psi,
[/tex]

where

[tex]
k^2 = \frac{{2m\left( {V_0 - E} \right)}}{{\hbar ^2 }}.
[/tex]
The solution to this 2nd order ODE is a linear combination of exponentials.

OK, so what about first region?

[tex]
\frac{{d^2 \psi }}{{dx^2 }} = q^2 \psi,
[/tex]

where
[tex]q^2 = \frac{{2m\left( {E} \right)}}{{\hbar ^2 }}.[/tex]

Isn't it also 2.nd order ODE ?

 

[Niles]

In the region with V=0 we have
[tex]
\frac{{d^2 \psi }}{{dx^2 }} = -\hat k^2 \psi,
[/tex]

where

[tex]
\hat k^2 = \frac{{2m\left( {E} \right)}}{{\hbar ^2 }}.
[/tex]

Solutions are traveling waves (going left and right).

 

[torehan]

Oh i see, I had some conceptual errors on writing the wave eq.
Thanks!
Torehan

 

[torehan]

I know that this topic doesn't take much attention of most of you as there are more interesting topics about paradoxes of physics but I need a little bit guidance about reproducing the band structure diagrams shown in Kittel , 170p in 8th Edition :)

So we have different wave solution for V₀ region and 0 region,

[tex]\psi_{1}(x)=Ae^{ikx}+Be^{-ikx}\ \ \ \ \ \ \ for \ \ \ \ \ \ \ \ 0<x<a[/tex]

[tex]\psi_{2}(x)=Ce^{qx}+De^{-qx}\ \ \ \ \ \ \ for \ \ \ \ \ \ \ \ \-b<x<0[/tex]

After we apply the boundary conditions (usual QM boundary condiions in in square well potentials ),

[tex] \psi_{1}(0) = \psi_{2}(0) [/tex] and [tex] \psi_{1}'(0) = \psi_{2}'(0) [/tex]

[tex] \psi_{1}(a) = \psi_{2}(-b) [/tex] and [tex] \psi_{1}'(a) = \psi_{2}'(-b) [/tex]

Then, we have four equations

[tex]A+B=C+D[/tex]

[tex]ik(A-B)=q(C-D)[/tex]

[tex]Ae^{ika}+Be^{-ika}=(Ce^{-qb}+De^{qb})e^{i\alpha(a+b)}[/tex]

[tex]ik(Ae^{ika}-Be^{-ika})=q(Ce^{-qb}-De^{qb})e^{i\alpha(a+b)}[/tex]

Determinant of the coefficients of this equation system mush vanish to have solutions. Than determinant yields,


[tex]([q^{2}-k^{2}]/2qk) sinh(qb)sin(ka) + cosh(qb)cosh(ka) = cos (\alpha(a+b)) [/tex]


Finally my issues

• How can I define the allowed and forbidden energy values?
• to be continued...

[tex]e^{i\alpha(a+b)}[/tex] ---->The term [tex]\alpha[/tex] came from Bloch's Theorem and (a+b) came from application of transitional symmetry operation to Bloch form wavefunction!