Kronig-Penneyev model kristala

Iz Fizika trdne snovi 2007 - 2008

(Primerjava redakcij)
Skoči na: navigacija, iskanje
Redakcija: 21:20, 21 januar 2008 (spremeni)
195.210.234.108 (Pogovor)

← Pojdi na prejšnje urejanje
Redakcija: 23:07, 21 januar 2008 (spremeni) (undo)
195.210.234.108 (Pogovor)

Novejše urejanje →
Vrstica 1: Vrstica 1:
Imamo delec v 1D Kronig-Penneyevem potencialu: Imamo delec v 1D Kronig-Penneyevem potencialu:
 +
<math>V\left(x\right)=\sum_n\lambda\delta\left(x-na\right)</math> <math>V\left(x\right)=\sum_n\lambda\delta\left(x-na\right)</math>
Vrstica 8: Vrstica 9:
Potencial je več kot očitno periodičen: <math>V\left(x+a\right)=V\left(x\right)</math> Potencial je več kot očitno periodičen: <math>V\left(x+a\right)=V\left(x\right)</math>
-Torej tudi za Hamiltonov operator: <math>\hat{H}=\hat{T}+\hat{V}=\frac{\hat{p_x}^2}{2m}+\hat{V}\left(x\right) =+Torej tudi za Hamiltonov operator velja, da je periodičen: <math>\hat{H}=\hat{T}+\hat{V}=\frac{\hat{p_x}^2}{2m}+\hat{V}\left(x\right) =
--\frac{\hbar^2}{2m} \frac{\partial^2}{\partial{x}^2} + \hat{V}(x)</math> velja, da je periodičen: <math>\hat{H}\left(x+a\right)=\hat{H}\left(x\right)</math>+-\frac{\hbar^2}{2m} \frac{\partial^2}{\partial{x}^2} + \hat{V}(x)</math>
 + 
 +<math>\hat{H}\left(x+a\right)=\hat{H}\left(x\right)</math>
Uvedimo operator translacije: <math>\hat{A}\psi(x)=\psi(x-a)</math> Uvedimo operator translacije: <math>\hat{A}\psi(x)=\psi(x-a)</math>
Vrstica 41: Vrstica 44:
<math>u_n(x)=u_{n+1}(x+a)</math> <math>u_n(x)=u_{n+1}(x+a)</math>
-<math>A_{n}e^{i(q-k)x}+B_{n}e^{-i(q+k)x}=A_{n+1}e^{i(q-k)(x+a)}+B_{n+1}e^{-i(q+k)(x+a)}</math>,+<math>A_{n}e^{i(q-k)x}+B_{n}e^{-i(q+k)x}=A_{n+1}e^{i(q-k)(x+a)}+B_{n+1}e^{-i(q+k)(x+a)}</math>
 + 
 +Ker to velja za vsak ''x'', sledi:
 + 
 +<math>A_{n+1}=A_{n}e^{-i(q-k)a}</math>
 + 
 +<math>B_{n+1}=B_{n}e^{i(q+k)a}</math>, torej
 + 
 +<math>A_{n}=Ae^{-i(q-k)na}</math>
 + 
 +<math>B_{n}=Be^{i(q+k)na}</math>, kjer sta: <math>A=A_0</math> in <math>B=B_0</math>
 + 
 +Torej lahko zapišemo rešitev Schrodingerjeve enačbe za ''n''-to področje kot:
 + 
 +<math>\psi_n(x)=e^{ikx}(Ae^{i(q-k)(x-na)}+Be^{-i(q+k)(x-na)})=</math>
 + 
 +<math>=Ae^{iqx}e^{-i(q-k)na}+Be^{-iqx}e^{i(q+k)na}</math>
 + 
 +Neznane količine: ''q'' (in torej tudi lastne energije ''E''), ''A'' in ''B'', dobimo iz začetnih in robnih pogojev.
 + 
 +----
 + 
 +Robni pogoji:
 + 
 +1.) <math>\psi_{n-1}(x=na)=\psi_{n}(x=na)</math>
 + 
 +2.) <math>(\frac{\partial\psi_{n}(x=na)}{\partial{x}}-\frac{\partial\psi_{n-1}(x=na)}{\partial{x}})=\frac{2m\lambda}{\hbar^2}\psi_{n}(x=na)=2Q\psi_{n}(x=na)</math>, kjer je <math>Q=\frac{m\lambda}{\hbar^2}</math>
 + 
 +Ko vstavimo valovno funkcijo v prvi robni pogoj, dobimo:
 + 
 +<math>\psi_{n-1}(na)=\psi{n}(na)</math>
 + 
 +<math>A_{n-1}e^{iqna}+B_{n-1}e^{-iqna}=A_{n-1}e^{-i(q-k)a}e^{iqna}+B_{n-1}e^{i(q+k)a}e{-iqna}</math>
 + 
 +<math>A_{n-1}e^{iqna}(1-e^{-iqa}e^{ika})=B_{n-1}e^{-iqna}(e^{iqa}e^{ika}-1)</math>

Redakcija: 23:07, 21 januar 2008

Imamo delec v 1D Kronig-Penneyevem potencialu:

V\left(x\right)=\sum_n\lambda\delta\left(x-na\right)

Iščemo lastne funkcije energije za ta sistem.


Potencial je več kot očitno periodičen: V\left(x+a\right)=V\left(x\right)

Torej tudi za Hamiltonov operator velja, da je periodičen: \hat{H}=\hat{T}+\hat{V}=\frac{\hat{p_x}^2}{2m}+\hat{V}\left(x\right) = -\frac{\hbar^2}{2m} \frac{\partial^2}{\partial{x}^2} + \hat{V}(x)

\hat{H}\left(x+a\right)=\hat{H}\left(x\right)

Uvedimo operator translacije: \hat{A}\psi(x)=\psi(x-a)

Za lastne funkcije tega operatorja torej sledi: \hat{A}\psi_a(x)=\alpha\psi_a(x)=\psi_a(x-a)

Ker je Hamiltonov operator periodičen s periodo a, komutira z operatorjem translacije, torej velja naslednja zveza: \hat{A}(\hat{H}\psi_a(x))=\hat{H}(\hat{A}\psi_a(x))=\alpha(\hat{H}\psi_a(x)).

Od tod sledi, da so lastne funkcije operatorja translacije hkrati tudi lastne fukcije Hamiltonovega operatorja: \hat{H}\psi_a(x)=\beta\psi_a(x)

Sedaj lahko izvedemo n-kratno translacijo tako, da n-krat delujemo na funkcijo z operatorjem translacije. Temu ustreza premik valovne fukcije za na. \hat{A}^n\psi_a(x)=\alpha^n\psi_a(x)=\psi_a(x-na)

Ker vemo, da je verjetnost, da tam najdemo delec končna in od nič različna in ker imamo neskončen kristal (torej lahko n izberemo poljubno velik), dobimo od tod pogoj za lastno vrednost operatorja translacije kot: | α | = 1, torej jo lahko zapišemo kot α = eiφ oziroma, če to izrazimo s periodo a, kot α = eika

Če sedaj povzamemo vse skupaj: ψa(xa) = eikaψa(x), oziroma ψa(x)eika = ψa(x + a)

Sedaj uporabimo sledeči nastavek:ψa(x) = eikxu(x), katerega nesemo v zgornjo enačbo: eik(x + a)u(x) = eik(x + a)u(x + a) in od tod dobimo pogoj za funkcijo u(x): u(x) = u(x + a). Toliko nam lahko da sama simetrija problema - torej periodičnost potenciala.


Lastne funkcije poiščemo s pomočjo Schrodingerjeve enačbe: i\hbar\frac{\partial\psi}{\partial{t}}=-\frac{\hbar^2}{2m} \frac{\partial^2\psi}{\partial{x}^2} + \hat{V}(x)\psi=\hat{H}\psi=E\psi

Rešitve v vmesnih področjih, kjer je V=0, so oblike

ψn(x) = Aneiqx + Bneiqx = ψn(x) = eikx(Anei(qk)x + Bnei(q + k)x) = eikxun(x), kjer je ψn(x) lastna funkcija energije v n-tem področju in q= \sqrt{\frac{2mE}{\hbar^2}}

Sedaj nesemo dobljeni izraz za u(x) v pogoj, dobljen iz periodičnosti potenciala: u(x)=u(x+a):

un(x) = un + 1(x + a)

Anei(qk)x + Bnei(q + k)x = An + 1ei(qk)(x + a) + Bn + 1ei(q + k)(x + a)

Ker to velja za vsak x, sledi:

An + 1 = Anei(qk)a

Bn + 1 = Bnei(q + k)a, torej

An = Aei(qk)na

Bn = Bei(q + k)na, kjer sta: A = A0 in B = B0

Torej lahko zapišemo rešitev Schrodingerjeve enačbe za n-to področje kot:

ψn(x) = eikx(Aei(qk)(xna) + Bei(q + k)(xna)) =

= Aeiqxei(qk)na + Beiqxei(q + k)na

Neznane količine: q (in torej tudi lastne energije E), A in B, dobimo iz začetnih in robnih pogojev.


Robni pogoji:

1.) ψn − 1(x = na) = ψn(x = na)

2.) (\frac{\partial\psi_{n}(x=na)}{\partial{x}}-\frac{\partial\psi_{n-1}(x=na)}{\partial{x}})=\frac{2m\lambda}{\hbar^2}\psi_{n}(x=na)=2Q\psi_{n}(x=na), kjer je Q=\frac{m\lambda}{\hbar^2}

Ko vstavimo valovno funkcijo v prvi robni pogoj, dobimo:

ψn − 1(na) = ψn(na)

An − 1eiqna + Bn − 1eiqna = An − 1ei(qk)aeiqna + Bn − 1ei(q + k)aeiqna

An − 1eiqna(1 − eiqaeika) = Bn − 1eiqna(eiqaeika − 1)