Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Energy Bands: Kronig-Penney Model

  1. Feb 14, 2010 #1
    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) =
    V0 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] ?
  2. jcsd
  3. Feb 14, 2010 #2
    In this region we have
    \frac{{d^2 \psi }}{{dx^2 }} = k^2 \psi,


    k^2 = \frac{{2m\left( {V_0 - E} \right)}}{{\hbar ^2 }}.
    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.
  4. Feb 14, 2010 #3

    OK, so what about first region?

    \frac{{d^2 \psi }}{{dx^2 }} = q^2 \psi,

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

    Isn't it also 2.nd order ODE ?
  5. Feb 14, 2010 #4
    In the region with V=0 we have
    \frac{{d^2 \psi }}{{dx^2 }} = -\hat k^2 \psi,


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

    Solutions are travelling waves (going left and right).
  6. Feb 14, 2010 #5
    Oh i see, I had some conceptual errors on writing the wave eq.
  7. Feb 19, 2010 #6
    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 V0 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





    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!
    Last edited: Feb 19, 2010
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook