Kronig-Penneyev model kristala

Iz Fizika trdne snovi 2007 - 2008

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

← Pojdi na prejšnje urejanje
Redakcija: 00:01, 22 januar 2008 (spremeni) (undo)
195.210.234.108 (Pogovor)

Novejše urejanje →
Vrstica 94: Vrstica 94:
<math>(iqe^{iqa}e^{ika}-iq+2Qe^{iqa}e^{ika}-iqe^{2ika}+iqe^{-iqa}e^{ika}-2Qe^{2ika})=(iqe^{2ika}-iqe^{iqa}e^{ika}-2Qe^{2ika}-iqe^{-iqa}e^{ika}+iq+2Qe^{-iqa}e^{ika})</math> <math>(iqe^{iqa}e^{ika}-iq+2Qe^{iqa}e^{ika}-iqe^{2ika}+iqe^{-iqa}e^{ika}-2Qe^{2ika})=(iqe^{2ika}-iqe^{iqa}e^{ika}-2Qe^{2ika}-iqe^{-iqa}e^{ika}+iq+2Qe^{-iqa}e^{ika})</math>
-Pomnožimo celotno enačbo s členom <math>e^{-ika}</math>+Pomnožimo celotno enačbo s členom <math>e^{-ika}</math>:
<math>(iqe^{iqa}-iqe^{-ika}+2Qe^{iqa}-iqe^{ika}+iqe^{-iqa}-2Qe^{ika})=(iqe^{ika}-iqe^{iqa}-2Qe^{ika}-iqe^{-iqa}+iqe^{-ika}+2Qe^{-iqa})</math> <math>(iqe^{iqa}-iqe^{-ika}+2Qe^{iqa}-iqe^{ika}+iqe^{-iqa}-2Qe^{ika})=(iqe^{ika}-iqe^{iqa}-2Qe^{ika}-iqe^{-iqa}+iqe^{-ika}+2Qe^{-iqa})</math>
 +
 +<math>2iqe^{iqa}+2iqe^{-iqa}-2iqe^{ika}-2iqe^{-ika}+2Qe^{iqa}-2Qe^{-iqa}=0</math>
 +
 +<math>4iq\cos(qa)-4iq\cos(ka)+4iQ\sin(qa)=0</math>
 +
 +<math>\cos(ka)=\cos(qa)+\frac{Q}{q}\sin(qa)</math>
 +
 +Prišli smo do enačbe, ki povezuje "valovni vektor" ''k'' s konstanto <math>q= \sqrt{\frac{2mE}{\hbar^2}}</math>
 +
 +----

Redakcija: 00:01, 22 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)

Iz drugega robnega pogoja pa dobimo:

ψn'(na) − ψn − 1'(na) = 2Qψn(na)

iq(An − 1ei(qk)aeiqnaBn − 1ei(q + k)aeiqnaAn − 1eiqna + Bn − 1eiqna) = 2Q(An − 1ei(qk)aeiqna + Bn − 1ei(q + k)aeiqna)

An − 1eiqna(iqeiqaeikaiq − 2Qeiqaeika) = Bn − 1eiqna(iqeiqaeikaiq + 2Qeiqaeika)

Sedaj lahko med seboj zdelimo obe enačbi, ki sledita iz robnih pogojev, ter se tako znebimo konstant A in B

(1 − eiqaeika)(iqeiqaeikaiq + 2Qeiqaeika) = (eiqaeika − 1)(iqeiqaeikaiq − 2Qeiqaeika)

(iqeiqaeikaiq + 2Qeiqaeikaiqe2ika + iqeiqaeika − 2Qe2ika) = (iqe2ikaiqeiqaeika − 2Qe2ikaiqeiqaeika + iq + 2Qeiqaeika)

Pomnožimo celotno enačbo s členom eika:

(iqeiqaiqeika + 2Qeiqaiqeika + iqeiqa − 2Qeika) = (iqeikaiqeiqa − 2Qeikaiqeiqa + iqeika + 2Qeiqa)

2iqeiqa + 2iqeiqa − 2iqeika − 2iqeika + 2Qeiqa − 2Qeiqa = 0

4iqcos(qa) − 4iqcos(ka) + 4iQsin(qa) = 0

\cos(ka)=\cos(qa)+\frac{Q}{q}\sin(qa)

Prišli smo do enačbe, ki povezuje "valovni vektor" k s konstanto q= \sqrt{\frac{2mE}{\hbar^2}}