functii de distributie in fizica starii solide - unibuc.ro · intrinsec), functia de distributie...
TRANSCRIPT
Functii de distributie in fizica starii solide
In cursul de fizica solidului s-a folosit frecvent functia de distributie Fermi-Dirac (FD)
pentru a descrie probabilitatea de ocupare cu electroni (in general, fermioni) a starilor intr-un
metal, functia de distributie Maxwell-Boltzmann (MB) pentru a descrie probabilitatea de
ocupare cu electroni a starilor dintr-un semiconductor (in particular, un semiconductor
intrinsec), functia de distributie Bose-Einstein BE (Planck) in cazul fononilor si o functie de
distributie diferita de toate acestea pentru starile impuritatilor donoare/acceptoare. Toate
aceste functii de distributie sunt de echilibru (definite la echilibru termodinamic).
Statistica MB mai este cunoscuta si ca statistica clasica pentru ca se aplica doar
particulelor care interactioneaza foarte slab intre ele (densitatea de particule este mica) si in
situatii in care numarul particulelor permise pe o stare cuantica nu este restrictionat de
principiul Pauli (ex: molecule de gaz, electroni si goluri intr-un semiconductor nedegenerat).
BE si FD sunt statistici cuantice pentru ca sunt gasite din principii ale mecanicii cuantice.
Particulele pentru care se aplica au o mai mare densitate si interactii mai puternice intre ele
decat particulele clasice. Particulele care se supun statisticii BE se numesc bosoni (fotoni,
fononi), cele care se supun statisticii FD sunt fermioni (electroni in metal, electroni si goluri
in semiconductor degenerat). Diferenta intre BE si FD este ca ocuparea starii cuantice pentru
fermioni este restrictionata de principiul Pauli, acest principiu ne-actionand pentru bosoni.
Principiul de excludere Pauli: o stare cuantica nu poate fi ocupata decat de cel mult
doua particule cu spini opusi, astfel incat numarul total de particule cu acelasi spin este egal
sau mai mic decat numarul total al starilor cuantice pe care particulele le pot ocupa intr-un
solid.
Toate aceste functii de distributie se gasesc folosind principiile mecanicii statistice,
care este bazata pe presupunerea comportarii aleatorii a unui numar foarte mare de particule.
In mecanica statistica se defineste:
microstare: stare a sistemului in care toti parametrii particulelor constituente sunt specificati
(pozitia si impulsul sau energia fiecarei molecule intr-un gaz)
macrostare: stare a sistemului specificata doar de parametri care sunt medii pe constituentii
microscopici ai sistemului (distributia particulelor pe nivele de energie)
ansamblu microcanonic: sistem izolat cu enerrgie fixa E si numar particule fix, N
ansamblu canonic: sistem la temperatura constanta (baie de caldura), cu N fixat
ansamblu macrocanonic: nimic nu este fixat
Functii de distributie de echilibru 2
Postulatul fundamental al mecanicii statistice este ca toate microstarile cu o energie
data E sunt egal probabile. Aceasta presupunere este adevarata in sisteme perfect aleatorii, dar
trebuie aplicata cu grija cand exista forte intre particule.
Cea mai probabila distributie este cea cu cel mai mare numar de numarul de microstari
W corespunzatoare unei macrostari. Acest numar este foarte mare pentru macrostarea de
echilibru, in care entropia este maxima. WkS B ln=
Functia de distributie MB In statistica clasica se presupune ca toate particulele din sistem pot fi distinse (le pot identifica
si urmari evolutia), ceea ce implica ca schimbarea pozitiei a doua particule duce la o alta
configuratie a sistemului. In plus, nu exista restrictii privind plasarea a mai mult de o particula
pe orice stare a sistemului.
Consider un sistem izolat cu N particule pe care le pot distinge unele de altele, cu
energie totala E. Problema este: cum pot distribui cele N particule pe q nivele de energie E1,
E2,.., Eq, fiecare nivel fiind ocupat de n1, n2,..., nq particule. Exista doua restrictii:
1) conservarea energiei ∑== i iiq EnEnnnC ),...,,( 211
2) conservarea numarului de particule ∑== i iq nNnnnC ),...,,( 212 .
Nu exista restrictii legate de numarul de particule pe care le pot pune pe un nivel energetic,
dar starile de energie pot avea degenerari g1, g2,..., gq.
Exemplu: problema distribuirii a n1 si n2 particule pe doua nivele nedegenerate. Daca
este numarul total de moduri independente de a aranja n1 particule pe nivelul E1 si
n2 particule pe nivelul E2, este numarul de moduri independente in care nivelul E1
este liber si nivelul E2 este ocupat. In acest caz exista doar o singura posibilitate, adica
= 1. Daca vreau sa am o particula pe nivelul E1 si
),( 21 nnW
),0( 2nW
),0( 2nW 1−N pe nivelul E2, exista N
moduri diferite de a pune o particula pe nivelul E1, deci NnW =),1( 2 . Daca vreau doua
particule pe nivelul E1 si pe nivelul E2, avem N moduri de a aranja prima particula pe
nivelul E1 si moduri de a pune cea de-a doua particula pe nivelul E2, astfel incat
numarul de moduri de a aranja 2 particule pe nivelul E1 si
2−N
1−N
2−N pe nivelul E2 este
unde factorul 2! tine cont de permutarile intre doua particule
identice. In general, pentru a aranja n1 particule din N pe nivelul E1
!2/)1()2,2( −=− NNNW
!!!
)!(!!
!)1)...(2)(1(),(
21111
121 nn
NnNn
Nn
nNNNNnnW =−
=+−−−
=
Functii de distributie de echilibru 3
iar pentru a aranja ulterior n2 particule pe nivelul E2 avem
)!(!)!(),(
212
121 nnNn
nNnnW−−
−=
Similar, numarul total de moduri de a aranja n1, n2,...,nq particule pe q nivele nedegenerate
este
=),...,,( 21 qnnnW ...)!(!
)!()!(!
!212
1
11×
−−−
×− nnNn
nNnNn
N
∏==
=
q
ii
q n
Nnnn
N
1
21 !
!!!...!
!
cu numarul total de particule. qnnnN +++= ...21
Similar, intr-un solid cu stari cuantice degenerate g1, g2,..., gq si n1, n2,..., nq particule pe
fiecare nivel energetic E1, E2,..., Eq, singura diferenta este ca am gi stari cuantice in plus in
fiecare stare energetica Ei (am gi stari cuantice cu aceeasi valoare a energiei). In acest caz
numarul de moduri in care pot aranja n1 particule in g1 stari cuantice pe nivelul energetic E1
este , etc., astfel incat numarul total de a aranja ni particule in fiecare din cele gi stari
cuantice cu energi Ei este
1)( 1ng
∏ ⎟⎟⎠
⎞⎜⎜⎝
⎛==
=
q
i i
ni
q
nq
nn
q ngN
nnngggN
nnnWiq
121
2121 !
)(!!!...!
)...()()(!),...,,(
21
sau
∑ −+==
q
iiiiq ngnNnnnW
121 )]!ln()ln([)!ln(),...,,(ln
Deoarece ni >> 1 putem folosi aproximatia Stirling: 1,ln!ln >>−≅ xxxxx , astfel incat
∑ +−+−≅=
q
iiiiiiq nnngnNNNnnnW
121 ])ln()ln([)ln(),...,,(ln
Folosind metoda multiplicatorilor Lagrange,
ii
q
i i
ii
ii
q
dndC
dndC
ngnNN
dnd
dnnnnWd 21
1
21 lnln),...,,(ln
βα +=⎥⎦
⎤⎢⎣
⎡∑ ⎟⎟
⎠
⎞⎜⎜⎝
⎛⎟⎠
⎞⎜⎝
⎛+==
,
Functii de distributie de echilibru 4
sau
ii
i Eng βα +=−⎟⎠
⎞⎜⎝
⎛ 1ln ,
cu α si β constante ce trebuie determinate.
Functia de distributie a particulelor clasice este definita prin
)]1(exp[)( ii
ii E
gnEf βα ++−==
sau, renuntand la indicele i, obtinem functia de distributie MB:
)exp()( EAEf β−= ,
unde constantele A si B se determina, de exemplu, din distributia moleculelor intr-un gaz ideal
in care energia particulelor este pur cinetica: si energia totala a gazului continand
N molecule este . Pentru o distributie izotropa a vitezelor, in spatiul vitezelor,
numarul particulelor cu viteze v, aflate pe suprafata unei sfere de raza v, este
2/2mvE =
TNkU B)2/3(=
)2/exp(4)(4 222 mvAvvfvN βππ −==
iar numarul particulelor cu viteze intre v si dvv + , aflate pe o coaja sferica de latime dv este
dvmvAvdvvfvdN )2/exp(4)(4 222 βππ −==
Punand conditiile ca numarul total de particule sa fie N:
∫ −=∫=∞
0
22 )2/exp(4 dvmvAvdNN βπ
si ca
TNkdvmvAvmvEdNU B23)2/exp(4
20
222
=∫ −⎟⎟⎠
⎞⎜⎜⎝
⎛=∫=∞
βπ
Functii de distributie de echilibru 5
rezulta 2/3
2⎟⎠⎞
⎜⎝⎛=
TkmNA
Bπ, TkB/1=β .
Viteza medie a unei particule clasice care satisface statistica MB se defineste ca
∫
∫=⟩⟨= ∞
∞
0
0
)(
)(
dvvN
dvvvNvv .
In general ∫
∫=⟩⟨ ∞
∞
0
0
)(
)(
dvvN
dvvNvv
n
n (in particular energia cinetica medie este ). ⟩⟨=⟩⟨ 2)2/1( vmE
Functii de distributie cuantice Desi in functia de distributie MB se presupune ca particulele nu interactioneaza si pot fi
distinse, in realitate este practic imposibil sa distingi electronii dintr-un metal sau
semiconductor degenerat datorita densitatii mari de electroni (intre 1019 si 1023 cm-3). In
statistica cuantica particulele nu pot fi distinse unele de altele, astfel incat schimband intre ele
pozitiile a doua particule nu se modifica configuratia sistemului.
Distributia BE se aplica particulelor cu spin intreg, pentru care dupa interschimbarea a
doua particule functia de unda a sistemului este aceeasi si nu exista restrictii asupra numarului
de particule pe o stare. Pentru un sistem de N particule cu numere cuantice n1, n2, ..., nN,
starea cuantica total simetrica fata de schimbarea a oricarei doua particule este
Functii de distributie de echilibru 6
∑ ⟩⟩⟩∏
=⟩p
Npppj j
SN nnnN
nnnn )()2()1(21 |...||
!!
...|
unde suma este peste toate starile diferite obtinute in urma permutarilor p a celor N elemente
si nj este numarul de ori in care fiecare stare uni-particule apare in starea a N particule.
Functia de distributie FD se aplica particulelor cu spin semiintreg pentru care, dupa
interschimbarea a doua particule functia de unda totala schimba semnul. Ca si consecinta. o
stare poate fi ocupata de cel mult o particula. Starea cuantica antisimetrica la schimbarea a
doua particule se scrie
∑ ⟩⟩⟩=⟩p
NpppAN nnnpN
nnn )()2()1(21 |...||)sgn(!
1...|
unde sgn(p) = semnul fiecarei permutari (+1 daca p este un numar par de transpozitii, -1 daca
numar impar).
Avem
1...|... 2121 =⟩⟨ SNNS nnnnnn
1...|... 2121 =⟩⟨ ANNA nnnnnn
Functiile de unda avnd expresiile
SNNSN
Snnn nnnxxxxxx
N⟩⟨=Ψ ...|...)...( 212121...21
∑∏
=p
NNpppj j xxxN
n)()...()(
!!
)(2)2(1)1( ψψψ
Functii de distributie de echilibru 7
ANNANA
nnn nnnxxxxxxN
⟩⟨=Ψ ...|...)...( 212121...21
)(...)()(............
)(...)()()(...)()(
!1)()...()()sgn(
!1
21
21
21
)(2)2(1)1(222
111
Nnnn
Nnnn
Nnnn
pNNppp
xxx
xxxxxx
Nxxxp
N
NNN ψψψ
ψψψψψψ
ψψψ =∑=
Functia de distributie FD
In acest caz, in afara de cele doua restrictii C1 si C2 considerate si in cazul MB (energia totala
si numarul de particule fix), principiul Pauli impune si a treia restrictie: . ii gn ≤
Daca calculez intai distributia celor ni particule pe nivelul de energie Ei, in cele gi stari
cuantice exista gi moduri de a aranja prima particula, 1−ig moduri de a aranja a doua
particula, etc., astfel incat numarul total de moduri de a aranja ni particule in gi stari cuantice
in nivelul de energie Ei este
)!(!)1)...(1(
ii
iiiii ng
gnggg−
=+−−
Deoarece particulele nu pot fi distinse, orice permutare intre ele nu reprezinta un aranjament
independent, astfel incat
)!(!!)(
iii
ii ngn
gnW−
= .
•••−−••−•−•−••−−•••−••−−••−•−•−••−••−−••−•−••−−•••
10987654321
Exemplu 1: gi = 5, ni = 3: numarul modurilor de a aranja 3 electroni
in 5 stari la energia Ei este
10!3)!35(
!5=
−
(vezi figura din dreapta)
Functii de distributie de echilibru 8
Exemplu 2: Din numarul total al modurilor de a aranja 5 electroni pe 3 nivele energetice cu
diferite degenerari, astfel incat energia totala sa fie 12, unele sunt mai probabile decat altele
Analog,
∏ ⎟⎟⎠
⎞⎜⎜⎝
⎛−
=−−−
==
q
i iii
i
qqq
qq ngn
gngngngnnn
gggnnnW
1221121
2121 )!(!
!)!)!...(()!(!!...!
!!...!),...,,(
din care obtinem
∑ ⎟⎟⎠
⎞⎜⎜⎝
⎛−
==
q
i iii
iq ngn
gnnnW1
21 )!(!!ln),...,,(ln
si folosind aproximatia Stirling
)]ln()(lnln[),...,,(ln1
21 iiiiiiiq
iiq ngngnnggnnnW −−−−∑≅
=
Cea mai probabila functie de distributie in statistica FD se obtine diferentiind W in raport cu ni
si aplicand metoda multiplicatorilor Lagrange asupra constrangerilor. Rezulta
iiii
ii
i
q EdndC
dndC
nng
dnnnnWd
βηβη +=+=⎟⎠⎞
⎜⎝⎛ −
= 2121 ln),...,,(ln
,
de unde se obtine functia de distributie FD
Functii de distributie de echilibru 9
)exp(11)(
ii
ii Eg
nEfβη ++
==
unde TkE BF /−=η este energia Fermi redusa si TkB/1=β . Lasand deoparte indicele i:
]/)exp[(11)(
TkEEEf
BF−+=
)(Ef = probabilitatea ca o stare cu energia E sa fie ocupata cu electroni
)(1 Ef− = probabilitatea ca starea cu energie E sa fie libera.
La T = 0 K, , toate starile cuantice de sub nivelul Fermi fiind
ocupate complet, toate cele de deasupra fiind libere.
⎩⎨⎧
><
=F
F
EEEE
Ef,0,1
)(
Acest rezultat este in contradictie cu MB, care prezice ca energia cinetica a
electronilor (=(3/2)kBT) este zero la T = 0 K. De fapt, energia cinetica a electronilor in metal
la T = 0 K este , cu EF(0) nivelul Fermi la T = 0 K. )0()5/3( FE
La T > 0K, , starile sub EF fiind partial ocupate, cele deasupra
partial umplute. Aceste stari sunt responsabile de conductia electrica in semiconductori. Doar
starile aflate intr-un interval de cativa kBT in jurul EF contribuie la conductia electrica in
metale sau semiconductori.
⎪⎩
⎪⎨
⎧
><=<>
=
F
F
F
EEEEEE
Ef,2/1
,2/1,2/1
)(
Functii de distributie de echilibru 10
Nivelul de energie Fermi
In general, este definit ca energia pana la care, la T = 0 K starile electronice in metale sau
semiconductori sunt ocupate complet, deasupra nivelului Fermi starile electronice fiind goale.
Trebuie precizat ca notiunea de energie Fermi are un inteles mai larg : la T = 0, EF este
energia degajata la adaugarea unui electron sau energia care trebuie data sistemului pentru a
scoate un electron, in conditiile in care starile initiale si finale sunt in echilibru. Daca
numarul initial de electroni este Ne si energia la echilibru este , )( eNE
)1()()()1( −−=−+= eeeeF NENENENEE
Presupunerea ca )1()()()1( −−=−+ eeee NENENENE este valabila in cazul metalelor dar
nu si al semiconductorilor intrinseci, pentru care, daca maximul BV este la energie 0 si
minimul BC este la Eg, energia de scoatere a unui electron = 0, cea de adaugare = Eg. In acest
caz EF se poate defini la jumatatea Eg. In general, EF depinde de dopanti sau impuritati in
semiconductori extrinseci.
Functii de distributie de echilibru 11
In teoria de benzi EF este un nivel ipotetic de energie al unui electron a.i. in echilibtru
termodinamic acest nivel de energie are probabilitate 50% de a fi ocupat. EF nu trebuie sa
corespunda unui nivel real de energie si poate fi definit si daca nu exista o structura de benzi
(LUMO, HOMO).
Potential chimic
Concept foarte similar cu energia Fermi. Definit ca: energia castigata sau pierduta cand o
particula este adaugata sau scoasa din sistem, daca sistemul este in echilibru in starile initiale
si finale.
Spre deosebire de EF, potentialul chimic μ este definit si la temperaturi finite, energia
din definitie fiind de fapt energia libera TSEF −= , cu S entropia. La temperaturi finite
starea de echilibru se obtine minimizand energia libera si nu doar E.
Importanta potentialului chimic: daca doua sisteme continand acelasi tip de particule
sunt in contact si in echilibru unul cu celalat, potentialul chimic al particulelor (care pot fi
bozoni sau fermioni) pentru cele doua sisteme trebuie sa fie identic.
Orice doua sisteme in contact schimba constant unul cu altul particule datorita
fluctatiilor termice, de exemplu. Daca potentialele chimice, μ1 si μ2 ale sistemelor cu numar
de particule N10 si N20 ar fi diferite, s-ar crea un flux de particule de la potentialul chimic mai
mare la cel mai mic, pentru ca astfel energia libera totala ar scadea. Acest flux se opreste doar
cand potentialele chimice se egalizeaza. Deci, daca initial μ1 > μ2, in final voi avea un
potential μ, a.i. μ1 > μ > μ2. In starea finala N1 < N10, N2 > N20.
Daca corpul schimba electroni si energie cu un electrod (rezervor), este descris de
ansamblul macrocanonic (grand canonical ensemble); μ e definit de electrozi si numarul de
electroni poate fluctua. In acest caz, TNF )/( ∂∂=μ cu N = numar mediu de electroni, T = la
temperatura constanta.
Functii de distributie de echilibru 12
Daca numarul de electroni in corp este fix (dar conductorul este conectat termic, de
exemplu, la un rezervor de caldura), avem ansamblul canonic. In acest caz lucrul necesar
pentru a adauga un electron este
),(),1(' TNFTNF −+=μ ,
si
),1('),1(),(),('' TNTNFTNFTN −=−−= μμ .
''' μμμ == doar in limita termodinamica.
Rezumat:
μ este definit la orice temperatura si pentru orice tip de particule
EF definit la T = 0, doar pentru fermioni
EF = μ la T = 0.
Prin exetensie in solid EF = μ , indiferent de temperatura
Cand se masoara o cadere de tensiune intr-un circuit, masuram lucrul total per unitate
de sarcina care poate fi obtinut permitand sarcinii sa curga de la un punct la celalat (daca
scurtcircuitez). Mai exact, ceea ce masor este eVV /)( 2121 μμ −−=− unde –e = sarcina
electronului. Electronii curg de la un corp cu μ mare (cadere de potential mica) la unul cu μ
mic (tensiune mare). In acest proces μ mic creste (datorita incarcarii), si μ mai mare descreste
astfel incat la echilibru termodinamic nivelul Fermi/μ este acelasi in tot circuitul. Deci si
tensiunea masurata cu voltmetru = 0 la echilibru.
Curentii curg doar daca f1 si f2 sunt diferiti!!
Functii de distributie de echilibru 13
EF/μ nu pot fi definite in afara echilibrului termodinamic. Atunci se defineste un cvasi-
EF pentru o anumita pozitie, care ar corespunde ocuparii starilor din acel punct in termeni de
functie de distributie. Cvasi-EF nu este acelasi de-a lungul circuitului si se defineste cand
- am sursa de tensiune exterioara
- sistemul este expus la campuri electromagnetice variabile (ca in capacitori, inductori
sau transformatori)
- la iluminare de la o sursa de lumina cu temperatura diferita (soare)
- cand temperatura nu este constanta in circuit/dispozitiv (ca in termocuple)
- cand materialul a fost supus unui stres si nu s-a intors inca in starea de echilibru, ca in
materiale piezoelectrice
Daca sistemul este perturbat foarte tare distributia electronilor nu poate fi descrisa de nici o
distributie termica si se spune ca electronii „nu sunt termalizati”. In situatii mai putin
dramatice, se pot defini μ si T diferite pentru banda de conductie si valenta.
In termodinamica potential chimic = energia potentiala care poate fi absorbita sau
eliberata in cursul unei reactii chimice. μ se poate schimba si in cursul unei tranzitii de faza.
Potentialul chimic al unei specii intr-o mixtura = panta energiei libere in sistem fata de o
variatie in numarul de moli al speciei respective. μ este o generalizare a „potentialului” in
fizica, pentru ca particulele se misca de la μ mare la μ mic. (in camp gravitational, de
exemplu, particula se misca in jos de la inaltime mai mare la mai mica). In timpul acestei
miscari numarul de particule se modifica; numarul de particule este marimea conjugata
potentialului chimic.
Exemple: Pentru o temperatura data o molecula are un potential chimic mai mare intr-o
regiune in care este concentratia este mai mare decat in una in care concentratia este mai
mica. Miscarea moleculelor de la μ mare la mic este insotita de eliberare de energie libera si
deci este un proces spontan.
Deasupra 0 oC moleculele de H2O in stare lichida au μ mai mic decat in stare solida, a.i.
gheata se topeste si H2O migreaza din stare solida in lichida. Sub 0 oC, moleculele in gheata
au potential chimic mai mic. La 0 oC potentialele chimice ale apei si ghetii sunt aceleasi =
sistem in echilibru.
Functii de distributie de echilibru 14
Functia de distributie BE
Valabila pentru fotoni si fononi. Presupuneri: 1) particulele nu pot fi distinse unele de altele,
fiecare avand o energie si un numar de unda cuantificat, 2) numarul de ocupare al fiecarei
stari cuantice nu este restrans de principiul de excludere Pauli.
Restrictiile C1 si C2 sunt in continuare valide (restrictia referitoatre la numarul fix de particule
in total nu este valabila pentru fononi (particule de tip Planck)).
Consider ni particule si partitii necesare sa divid aceste particule in gi stari
cuantice. Numarul de moduri de aranjare a ni particule in gi stari cuantice = numarul de
permutari independente a particulelor si partitiilor. Deoarece sunt particule si
partitii, ele pot fi aranjate liniar in
1−ig
1−+ ii gn
)!1( −+ ii gn moduri. Dar, permutarea particulelor intre ele,
in numar de , si a partitiilor intre ele, in numar de !in )!1( −ig , nu sunt aranjamente
independente. Numarul total de moduri independente de a aranja ni particule in gi stari
cuantice pe nivelul de energie Ei este
)!1(!)!1(
−−+
=ii
iii gn
gnW
Spre deosebire de statistica FD, ni poate fi mai mare ca gi in BE.
∑ ⎟⎟⎠
⎞⎜⎜⎝
⎛−−+
==
q
i ii
iiq gn
gnnnnW1
21 )!1(!)!1(ln),...,,(ln
ii
i
i
q Eng
dnnnnWd
βα +=⎥⎦
⎤⎢⎣
⎡+⎟
⎠
⎞⎜⎝
⎛≅ 1ln),...,,(ln 21
Astfel incat obtinem functia de distributie BE
1)exp(1)(
−+==
ii
ii Eg
nEfβα
Functii de distributie de echilibru 15
cu TkB/1=β si α arbitrar (poate fi chiar = 0 pentru ca nu exista nici o restrictie asupra lui).
De fapt α = 0 corespunde cazului cand numarul particulelor nu este fixat; distributia fiind
numita in acest caz Planck. Acesta este si cazul fononilor pentru care functia de distributie are
forma
1)/exp(1)(
−=
TkEf
Bωh.
Ca si distributia FD, si BE se reduce la MB la temperaturi mari.
Statistica impuritatilor putin adanci in semiconductor Functia de distributie a electronilor pe nivele donoare putin adanci sau a golurilor pe nivele
acceptoare putin adanci este diferita de FD (care descrie distributia electronilor si golurilor in
BC si BV) deoarece starile de impuritate putin adanci pot fi ocupate doar de un singur
electron/gol cu spin sus sau jos. In momentul in care este ocupata, starea impuritatii donoare
devine neutra, si nu mai este ocupata de un al doilea electron.
Presupun stari donoare la cativa kBT sub Ec cu densitatea Nd si energia de ionizare Ed. Daca
Wd este numarul de moduri in care pot aranja nd electroni in cele Nd stari, exista 2Nd moduri
de a pune primul electron in starile Nd tinand cont de degenerarea de spin (care este 2),
2(Nd – 1) moduri de a pune al doilea electron in celelalte Nd – 1 stari, etc. Deoarece electronii
nu pot fi distinsi, permutarile intre ei nu sunt aranjamente independente.
Numarul total de a aranja nd electroni pe Nd stari donoare este
)!(!!2!/)1(2)...1(2)(2
ddd
dn
dddddd nNnNnnNNNW
d
−=+−−=
Factorul 2nd tine cont de natura Coulombiana a starilor de impuritate donoare putin adanci.
Similar cu cazurile de mai sus,
Functii de distributie de echilibru 16
]/)exp[(211
1)(TkEEN
nEfBFd
d
ddd
−+==
In general,
]/)exp[(11)(
1 TkEEgEf
BFdddd
−+=
−,
cu gd = factorul de degenerare a starilor donoare. (gd = 2 in semiconductori cu banda de
conductie sferica, si = 12 in Si, cu banda de conductie cu mai multe minime.)
Pentru functia de distributie a starilor acceptoare putin adanci se gaseste, similar,
]/)exp[(11)(
TkEEgNpEf
BaFaa
aaa
−+==
In cazul semicond cu benzi de valenta a golurilor usoare si grele ga = 4.
Observatie: )(1)( adaa EfEf −=
Exemplu 1:
Comparatie intre MB, FD si BE
Functii de distributie de echilibru 17
Densitate energie u pentru spectrul Planck
Exemplu 2:
Avem un sistem cu 6 electroni cu energia totala 9 eV.
Numarul total de configuratii cu aceste doua restrictii este 26 (vezi mai jos) pentru particule
care nu pot fi distinse (doar 5 pentru FD – verzui), dar starile diferite sunt 2002 daca
particulele pot fi distinse (numerele de deasupra indica microstarile compatibile cu fiecare
configuratie)
E Numar mediu particule: MB BE FD
0 2.143
(5x6+4x4x30+3x3x120+3x3x60+60+120+4x2x180+120+30+120+60+180+30)/2002
2.269 (5+4x4+7x3+6x2+5x1)
/26
1.8 (2x4+1)/5
Functii de distributie de echilibru 18
1
1.484 (30+120x2+2x60+60+3x60+2x2x180+2x180+4x30+3x120+3x60+2x180+30+5x6+4x3
0+3x20)/2002
1.538 (7x1+4x2+4x3+2x4+5)
/26
1.6 (3x2+2)/5
2 0.989
(1x30+2x120+2x60+4x180+3x60+1x120+2x180+4x30+30+3x20)/2002
0.885 (7x1+3x2+2x3+4)/26
1.2 (2x2+2)/5
3 0.629
(1x30+2x120+3x20+3x180+3x60+1x180+1x30)/2002
0.538 (7x1+2x2+3)/26
0.8 (2+2)/5
4 0.378 (1x30+1x120+2x60+2x180+1x120+6)/2002
0.307 (6x1+2)/26
0.4 2/5
5 0.210 (1x30+1x120+1x60+1x180+1x30)/2002
0.192 5/26
0.2 1/5
6 0.105 (1x30+1x120+1x60)/2002
0.115 3/26 0
7 0.045 (1x60+1x30)/2002
0.077 2/26 0
8 0.015 30/2002
0.038 1/26 0
9 0.003 6/2002
0.038 1/26
0
Total: 6 Total: 6 Total: 6
Exemplu 3:
8 din cele 24 configuratii posibile in care 20 Probabilitate versus enrgie mediata pe electroni pot fi plasati astfel incat energia cele 24 configuratii (cercuri), fitata cu o totala sa fie 106 eV. functie FD cu kBT =1.447 eV si
EF = 9.998 eV
Operatori de creare si anihilare. Cuantizarea a doua
In mecanica cuantica, pentru un Hamiltonian independent de timp exista trei formalisme:
1) Schrödinger, in care operatorii nu depind de timp, iar vectorii de stare depind de timp
⟩Ψ=∂
⟩Ψ∂ )(|ˆ)(| tHt
ti sSsh ,
⟩Ψ⎟⎠⎞
⎜⎝⎛−=⟩Ψ⎟⎟
⎠
⎞⎜⎜⎝
⎛−≡⟩Ψ⎟
⎠⎞
⎜⎝⎛
∑=
=∫ ⎟⎟⎠
⎞⎜⎜⎝
⎛⟩Ψ∫⎟
⎠⎞
⎜⎝⎛+⟩Ψ+⟩Ψ=
⟩Ψ∫+⟩Ψ=⟩Ψ
∞
=
<
)0(|exp)0(|ˆ
exp)0(|ˆ!
1
....'")"(|ˆˆ1)0(|ˆ)0(|
')'(|ˆ1)0(|)(|
,0
0
'
0
20
iSi
SS
SnS
n
n
tS
ttSSSSS
St
SSS
itEHitHit
n
dtdttHHi
Hit
dttHi
t
hhh
hh
h
daca este o functie proprie a Hamiltonianului (definitia exponentialei unui
operator este dezvoltarea in serie a operatorului).
iS ,Ψ
2) Heisenberg, in care operatorii depind de timp si vectorii de stare sunt independenti de
timp
⟩Ψ⎟⎟⎠
⎞⎜⎜⎝
⎛≡⟩Ψ )(|
ˆexp| tHit
sS
Hh
De fapt, SH HH ˆˆ =
0)0(|)0(|ˆ
expˆ
exp)(|ˆ
exp|=
∂⟩Ψ∂
=⎥⎥⎦
⎤
⎢⎢⎣
⎡⟩Ψ⎟⎟
⎠
⎞⎜⎜⎝
⎛−⎟⎟
⎠
⎞⎜⎜⎝
⎛
∂∂
=⎥⎥⎦
⎤
⎢⎢⎣
⎡⟩Ψ⎟⎟
⎠
⎞⎜⎜⎝
⎛
∂∂
≡∂
⟩Ψ∂t
HitHitt
tHittt
Sss
H
hhh
⎟⎟⎠
⎞⎜⎜⎝
⎛−⎟⎟
⎠
⎞⎜⎜⎝
⎛=
hh
HitOHittO SHˆ
expˆˆexp)(ˆ
]ˆ,ˆ[ˆ
exp]ˆ,ˆ[ˆ
expˆ
HOHitHOHitt
Oi HSH =⎟⎟
⎠
⎞⎜⎜⎝
⎛−⎟⎟
⎠
⎞⎜⎜⎝
⎛=
∂∂
hhh
Relatie: ⟩Ψ−Ψ⟨=⟩ΨΨ⟨ HSHSSS HitOHittOt |)/ˆexp(ˆ)/ˆexp(|')(|ˆ|)(' hh
3) de interactie: si operatorii si vectorii de stare depind de timp. In general Hamiltonianul
este impartit in doua parti: unde corespunde unei probleme cu 10 ˆˆˆ HHH += 0H
solutie stiuta si se doreste aflarea efectelor . 1H
Operatori de creare si anihilare 2
⟩Ψ⎟⎟⎠
⎞⎜⎜⎝
⎛≡⟩Ψ )(|
ˆexp| 0 tHit
sIh
⟩Ψ≡⟩Ψ⎟⎟⎠
⎞⎜⎜⎝
⎛−⎟⎟
⎠
⎞⎜⎜⎝
⎛=
∂⟩Ψ∂
⎟⎟⎠
⎞⎜⎜⎝
⎛+⟩Ψ⎟⎟
⎠
⎞⎜⎜⎝
⎛−=
∂⟩Ψ∂
)(|)(ˆ)(|ˆ
expˆˆexp
|ˆexp)(|
ˆexpˆ|
,10
10
000
ttHtHitHHit
tiHittHitH
ti
III
SS
I
hh
hhh
h
⎟⎟⎠
⎞⎜⎜⎝
⎛−⎟⎟
⎠
⎞⎜⎜⎝
⎛=
hh
00 ˆexpˆˆ
exp)(ˆ HitOHittO SI
In toate aceste trei formalisme, cunoscute si ca prima cuantizare, particulelor clasice li se
asociaza functii de unda.
In a doua cuantizare functiile de unda sunt cuantizate, astfel incat descriu cuante sau
particule cuantice iar procesele/interactiile se descriu ca si schimburi (exchange) de particule
reale sau virtuale.
Pentru fermioni
Operatori de creare si anihilare 3
Operatorii de creare si anihilare sunt folositi in diverse domenii cum ar fi in studiul
oscilatorului cuantic armonic si in sisteme cu mai multe particule. Operatorul de anihilare
scade numarul particulelor dintr-o stare data cu o unitate. Operatorul de creare creste numarul
particulelor intr-o stare data cu unu, adjunctul sau fiind operatorul de anihilare. Folosirea
acestor operatori in locul functiilor de unda se mai numeste si cuantizarea a doua.
Prima si a doua cuantizare trebuie sa duca la un acelasi rezultat, dar este mai
convenabil sa folosim a doua cuantizare pentru calculul proprietatilor sistemelor multi-
electronice deoarece proprietatile antisimetrice sunt inglobate in operatorii de creare si
anihilare si nu in determinanti Slater. Acesti operatori sunt diferiti pentru bozoni (spin intreg)
si fermioni (spin semi-intreg) deoarece functiile de unda au simetrii diferite.
Bozoni: ∑ ⟩⟩⟩∏
=⟩p
Npppj j
SN nnnN
nnnn )()2()1(21 |...||
!!
...|
SNNSNS
nnn nnnxxxxxxN
⟩⟨=Ψ ...|...)...( 212121...21
∑∏
=p
NNpppj j xxxN
n)()...()(
!!
)(2)2(1)1( ψψψ
Fermioni: ∑ ⟩⟩⟩=⟩p
NpppAN nnnpN
nnn )()2()1(21 |...||)sgn(!
1...|
ANNANA
nnn nnnxxxxxxN
⟩⟨=Ψ ...|...)...( 212121...21
)(...)()(............
)(...)()()(...)()(
!1)()...()()sgn(
!1
21
21
21
)(2)2(1)1(222
111
Nnnn
Nnnn
Nnnn
pNNppp
xxx
xxxxxx
Nxxxp
N
NNN ψψψ
ψψψψψψ
ψψψ =∑=
Operatori de creare si anihilare 4
In chimia cuantica sau in teoria many-body operatorii de creare si anihilare actioneaza
asupra starilor electronice, pe cand in cazul oscilatorului armonic cuantic acestia (numiti si
„ladder operators”) cresc sau scad numarul cuantelor de energie a sistemului; pot fi folositi sa
reprezinte si fononi. Matematica operatorilor de creare si anihilare pentru bozoni este aceeasi
ca si in cazul ladder operators a oscilatorului, comutatorul acestor operatori asociat cu aceeasi
stare bozonica fiind egal cu 1, si 0 in celelalte cazuri. Pentru fermioni matematica este
diferita, si comutatorii sunt inlocuiti de anticomutatori.
Oscilatorul armonic cuantic. Ladder operations
⟩Ψ=∂
⟩Ψ∂ ),(|ˆ),(| txHt
txih , ),(,ˆ xpHxx
iH cl=⎟⎠⎞
⎜⎝⎛
∂∂
− h
222
2222
2
21
221
2xm
xmxm
mp ωω +
∂
∂−→+h
In acest caz operatorii de creare sau anihilare aduna la sau scad din sistemul oscilatorului o
cuanta fixa de energie. Pentru un sistem bozonic de fononi ai unui oscilator armonic cuantic,
daca in ecuatia Schrödinger pentru un oscilator 1D independent de timp
)()(21
222
2
22xExxm
dxd
mψψω =⎟
⎟⎠
⎞⎜⎜⎝
⎛+−
h
se face schimbarea de variabila qm
xωh
= aceasta devine
Operatori de creare si anihilare 5
)()(2
22
2qEqq
dqd ψψω
=⎟⎟⎠
⎞⎜⎜⎝
⎛+−
h
Parametrul νω h=h este energia unei cuante de lumina. Termenii din paranteza de mai sus
pot fi aranjati ca
dqdqq
dqdq
dqdq
dqdq
dqd
−+⎟⎟⎠
⎞⎜⎜⎝
⎛+⎟⎟
⎠
⎞⎜⎜⎝
⎛+−=+− 2
2
2
ultimii doi termeni putand fi simplificati daca se considera efectul lor asupra unei functii
arbitrare : )(qf
)()()]([)( qfdq
qdfqqqfdqdqf
dqdqq
dqd
=−=⎟⎟⎠
⎞⎜⎜⎝
⎛−
adica
1=−dqdqq
dqd
sau
122
2+⎟⎟
⎠
⎞⎜⎜⎝
⎛+⎟⎟
⎠
⎞⎜⎜⎝
⎛+−=+− q
dqdq
dqdq
dqd
si ecuatia Schrödinger se poate scrie
)()(21
21
21 qEqq
dqdq
dqd ψψω =⎥
⎦
⎤⎢⎣
⎡+⎟⎟
⎠
⎞⎜⎜⎝
⎛+⎟⎟
⎠
⎞⎜⎜⎝
⎛+−h .
Daca definim
⎟⎟⎠
⎞⎜⎜⎝
⎛+−=+ q
dqda
21 ca si operator de creare si
⎟⎟⎠
⎞⎜⎜⎝
⎛+= q
dqda
21 ca si operator de anihilare,
ecuatia Schrödinger pentru oscilator devine
Operatori de creare si anihilare 6
)()(21 qEqaa ψψω =⎟⎠⎞
⎜⎝⎛ ++h
Aceasta ecuatie este semnificativ mai simpla decat cea originala, dar poate fi simplificata in
continuare.
Introducand operatorul impulsului dqdip −= avem
ipq =],[
si
)(2
12
1 ipqdqdqa +=⎟⎟⎠
⎞⎜⎜⎝
⎛+= , )(
21
21 ipq
dqdqa −=⎟⎟⎠
⎞⎜⎜⎝
⎛−=+
ceea ce implica
1]),[],([2
]),[],([21],[
21],[ =+
−=+−=−+=+ pqpqiqipipqipqipqaa
adica operatorii q si p nu au functii proprii sau distributii proprii comune.
In operatorii originali, cu dxdip h−=ˆ ,
⎟⎠⎞
⎜⎝⎛ −=+
ωω
mpixmaˆˆ
2h, ⎟
⎠⎞
⎜⎝⎛ +=
ωω
mpixmaˆˆ
2h
)(2
ˆ ++= aam
xωh , )(
2ˆ +−= aamip ωh
Operatorul Hamiltonian poate fi exprimat ca
⎟⎠⎞
⎜⎝⎛ −=⎟
⎠⎞
⎜⎝⎛ += ++
21
21ˆ aaaaH ωω hh
cu
aaH ωh−=],ˆ[ , ++ = aaH ωh],ˆ[
Operatori de creare si anihilare 7
Aceste relatii pot fi folosite pentru a gasi starile proprii de energie ale oscilatorului armonic
cuantic, nψ . Daca , din relatiile de comutare obtinem nnn EH ψψ =ˆ
nnn aEaH ψωψ )(ˆ h−= , nnn aEaH ψωψ ++ += )(ˆ h
ceea ce arata ca naψ si sunt si stari proprii ale hamiltonianului cu valori proprii na ψ+
ωh−nE si ωh+nE . Aceste relatii identifica operatorii a si ca operatori de anihilare si
creare a unor cuante de energie, diferenta in energie intre doua stari proprii fiind
+a
ωh=ΔE .
Fiecare excitatie (cuanta de energie) este reprezentata prin operatorul care are
propria sa ecuatie de miscare (in reprezentare Heisenberg):
+a
]ˆ,[ Hat
ai HH ++=
∂∂
h .
Prin introducerea operatorilor de creare si anihilare excitatiile unui sistem sunt considerate ca
particule individuale, cu propria ecuatie de miscare.
Starea fundamentala se defineste ca cea pe care operatorul de anihilare colapseaza,
adica 00 =ψa . Apoi, folosind Hamiltonianul scris in termenii acestor operatori,
00
00 21
21ˆ
0 ψω
ψω
ψ ⎟⎠⎞
⎜⎝⎛ −=⎟⎟
⎠
⎞⎜⎜⎝
⎛−==+
hh
EHaa
Obtinem ca termenul in paranteza trebuie sa fie zero pentru ca functia de unda nu poate fi,
adica energia starii fundamentale este 2/0 ωh=E . Aceasta relatie permite identificarea
valorii proprii a energiei pentru orice stare proprie nψ ca
ωh⎟⎠⎞
⎜⎝⎛ +=
21nEn
In plus, se poate introduce operatorul numar de particule , astfel incat ,
, si
aaN += 1+=+ Naa
0],[ =HN
)()(21 qEqN ψψω =⎟⎠⎞
⎜⎝⎛ +h
Operatori de creare si anihilare 8
Daca , starea fundamentala are proprietatea ⟩=⟩ nnnN || 00| =⟩a , respectiv
⟩=⟩⎟⎠⎞
⎜⎝⎛ +=⟩ + 0|
210|
210| ωω hhaaH ,
Toate celelalte stari putand fi obtinute din starea fundamentala prin aplicarea repetata a
operatorului de creare:
⟩⎟⎠⎞
⎜⎝⎛ +=⟩⎟
⎠⎞
⎜⎝⎛ +=⟩⎟
⎠⎞
⎜⎝⎛ +=⟩=⟩ ++++ nnanaaaaHnH nnn |
210|)(
210|)(
210|)(| ωωω hhh
cu ⟩=⟩+
0|!)(|
nan
n
⟩−=⟩ 1|| nnna
⟩++=⟩+ 1|1| nnna ⟩=⟩+ nnnaa ||
⟩⟨= 0|)(0 qqψ ⟩⟨= + 0|)(|!
1)( nn aq
nqψ
Operatorii de anihilare/anihilare extind conceptul operatorilor de creare/anihilare la un sistem
de mai multe particule in interactie.
Acesti operatori tin cont de simetria/statistica corecta pentru bozoni sau fermioni si permit
definirea starilor cu numar variabil de particule (Fock)
Functii de unda
Exprimand starea fundamentala )(0 qψ a oscilatorului armonic cuantic prin conditia
0)(0 =qaψ scrisa ca
000 =+
dqdq ψψ
obtinem solutia
)2/exp()( 2
0 qCq −=ψ
Operatori de creare si anihilare 9
constanta 4/1 π=C fiind obtinuta din conditia . 10*0 =∫
∞
∞−dqψψ
In general
)2/exp(),(!2
1)( 2qqnHn
qn
n −=π
ψ
),( qnH functii Hermite, care formeaza o baza ortonormata: nmmn dqqq δψψ =∫∞
∞−)()(
)(1 qψ )(2 qψ )(3 qψ
q q q
)()()( 322 qqdqdqq ψψψ ≈− )()()( 122 qq
dqdqq ψψψ ≈+
q q
Reprezentare matriciala
In formularea matriciala a mecanicii cuantice functiile proprii ale oscilatorului armonic sunt
vectori:
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
=⟩
0
01
| 0M
ψ , | , . . . |
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
=⟩
0
10
1M
ψM
Nψ
⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜
⎝
⎛
=⟩
1
00
Operatori de creare si anihilare 10
Acesti vectori formeaza o baza ortonormala. Expresiile matriciale ale operatorilor de creare si
anihilare obtinuti din modelul oscilatorului armonic cuantic sunt
⎟⎟⎟⎟⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜⎜⎜⎜⎜
⎝
⎛
=+
OMMMMM
K
KMOMMM
KK
KK
KK
KK
n
a
0000
0300002000010000
ˆ
⎟⎟⎟⎟⎟⎟⎟⎟⎟⎟
⎠
⎞
⎜⎜⎜⎜⎜⎜⎜⎜⎜⎜
⎝
⎛
=
OMMMMMM
OK
KOMMMM
KMO
KK
KK
KK
00000
0000030000020000010
ˆ
na
Din aceste expresii se pot obtine elementele de matrice ale laddering operators folosind
, ⟩⟨= ++jiij aa ψψ |ˆ| ⟩⟨= jiij aa ψψ |ˆ| , unde functiile de unda sunt cele ale oscilatorului
armonic cuantic, numite si „baza numarului de particule” („number basis”).
Detalii matematice
Operatorii derivati mai sus sunt un exemplu dintr-o clasa mai generala de operatori de creare
si anihilare, care satisfac proprietatile de mai jos.
Daca H este spatiul Hilbert uni-particula, pentru a gasi algebra bozonica ne uitam intai
la algebra generata de pentru orice f in H. Operatorul este numit operator de
anihilare si harta (map) a(.) este antilineara. Adjunctul operatorului este , care este
linear in H.
)( fa )( fa
)( fa+
Pentru bozoni,
0)](),([)](),([ == ++ gafagafa , ⟩⟨=+ gfgafa |)](),([
iar pentru fermioni, anticomutatorii sunt
0)}(),({)}(),({ == ++ gafagafa , ⟩⟨=+ gfgafa |)}(),({
Operatorul anihileaza o particula in starea , iar operatorul creaza o particula
in starea . Starea de vid este o stare fara particule, pentru care
)( fa ⟩f| )( fa+
⟩f| ⟩0|
00|)( =⟩fa
Operatori de creare si anihilare 11
Daca este normalizat, astfel incat ⟩f| 1| =⟩⟨ ff , da numarul de particule in
starea .
)()( fafa+
⟩f|
In teoriile cuantice de camp si in problemele many-body operatorii de creare si
anihilare a starilor cuantice, si schimba valorile proprii ale operatorului numar de
particule
+ia ia
∑=∑= +
iii
ii aanN
cu o unitate, ca si in cazul oscilatorului armonic. In acest caz indicele i reprezinta numere
cuantice care indica starile unei particule din sistem. ),,,( smln
Relatiile de comutare ale operatorilor de creare si anihilare intr-un sistem multi-bozonic sunt
ijijjiji aaaaaa δ=−= +++ ],[ , 0],[],[ ==++jiji aaaa
cu ijδ functiile delta Kronecker,
⟩−=⟩ ,...1,...,,|,...,...,,| 1010 iiii nnnnnnna , = 0, 1, 2, 3, ... in⟩++=⟩+ ,...1,...,,|1,...,...,,| 1010 iiii nnnnnnna
in timp ce pentru fermioni comutatorul este inlocuit de anticomutator
ijijjiji aaaaaa δ=+= +++},{ , 0},{},{ ==++jiji aaaa
⟩−−=⟩ ,...1,...,,|)1(,...,...,,| 1010 iis
ii nnnnnnna i , , = 0, 1 ∑=−
=
1
0
i
jji ns in
⟩+−−=⟩+ ,...1,...,,|1)1(,...,...,,| 1010 iis
ii nnnnnnna i
Deci, schimband intre ei operatori disjuncti (cu ji ≠ ) intr-un produs de operatori de creare si
anihilare schimba semnul produsului intr-un sistem de fermioni, dar nu si intr-un sistem de
bozoni.
Operator 1-particula: , βαβα αβ aaAA +∑= ,ˆ dxxxAxAA )()()(|| *∫=⟩⟨= βααβ ψψβα
Operator 2-particula: γδδγβα
βα γδαβ aaAaaA ⟩⟨∑= ++ ||21ˆ
,,,,
2112212*
1* )()(),()()(|| dxdxxxxxAxxA ∫=⟩⟨ γδβα ψψψψγδαβ
Functii de distributii de ne-echilibru: exemple
In teoria cinetica sistemul este descris prin distributii de
particule
kxkx
dddNtf =),,(
Evolutia sistemului de particule in absenta imprastierii este determinista, cea in prezenta
imprastierii este nedeterminista, iar ecuatiile cinetice descriu evolutia in timp a acestor
distributii sub influenta fortelor externe sau a interactiilor intre particule.
Se presupune in general ca distributiile de particule variaza incet in timp si spatiu. (x,k) sau
(x,p) reprezinta spatiul fazelor in mecanica clasica.
Ecuatia cinetica de transport Boltzmann spune ca, in absenta imprastierilor, probabilitatea se
conserva de-a lungul traiectoriei particulei
determinata de
)(0
0 ττvxx ∫+=t
td , ))(,(
0
0 τττ xFpp ∫+=t
td
adica
),,(),)/(,( tfdttdtdtf krFkvr =+++ h
Exemplu: traiectoria unei particule in spatiul fazelor, in prezenta unui camp electric Ex, este
ilustrata in figura din dreapta.
Functii de distributie de ne-echilibru. Exemple 2
Daca nu sunt aplicate campuri/forte externe si distributia initiala este de ne-echilibru,
imprastierile sunt cele care re-aduc functia de distributie la echilibru dupa un timp suficient de
lung (teoretic, dupa un timp infinit de lung)
In prezenta proceselor de imprastiere, in aproximatia timpului de relaxare,
τ1fff
dtd
tf
−=∇⋅+∇⋅+∂∂
krFrh
Pentru forte aplicate mici solutia in aproximatia timpului de relaxare se poate liniariza.
Exemple de evolutie a functiei de distributie
A) Camp electric extern stationar pe directia x
La aplicarea unui camp electric pe directia x, forta resimtita de un electron este EF e−=
01
1 )()( fef kEkk ∇⋅= −hτ
⎟⎠⎞
⎜⎝⎛ +≅
∂∂
+= xEkekfkfk
ekfkf )()()()()( 000 ττhh
Ek
Deci campul electric Ex deplaseaza suprafata Fermi cu
h/xx Eek τδ −= in directia opusa campului.
Daca imprastierea electronilor/relaxarea este puternica, τ este mic si deplasarea este mica, si
vice-versa.
Functii de distributie de ne-echilibru. Exemple 3
Din relatia de echilibru a probabilitatii de imprastiere in aproximatia timpului de relaxare,
)]'(1)[()',()](1)['(),'( 0000 kkkkkkkk ffPffP −=− P
scrisa sub forma P
]/))'(exp[(]/))(exp[(
)]'(1)[()](1)['(
),'()',(
00
00
TkEETkEE
ffff
PP
BF
BF
−−
=−−
=kk
kkkk
kkkk
rezulta ca, in cazul imprastierii ne-elastice rata de imprastiere pe stari cu energie mai joasa
este mai mare decat pe stari cu energie mai mare.
E posibil sa am curent: sunt mai multi electroni cu +k decat cu –k.
B) Camp electric oscilant
)exp( tie ω−−= EF
Daca functia de distributie nu variaza spatial, ecuatia cinetica Boltzmann se poate scrie
τω 101 )exp( f
Eftie
tf
−=∂∂
⋅−−∂∂ vE
Cautand solutii de tipul )exp()( 101 tiftf ω−= obtinem
vE ⋅∂∂
−−
= eEf
ititf 0
1 1)exp()(
ωτωτ
In acest caz conductivitatea electrica depinde de ω:
Functii de distributie de ne-echilibru. Exemple 4
ωττωσ
imne
−⋅=1
1)(2
in materiale izotrope.
C) Camp magnetic stationar
Din ecuatia de transport Boltzmann avem
feff kBvkk ∇⋅×+= )()()( 0h
τ
Avand in vedere ca vkk ≈⎟⎠⎞
⎜⎝⎛∇=∇
dEdfEf 0
0 in cazul suprafetelor izoenergetice sferice, in
ultimul termen din ecuatia de mai sus f nu mai poate fi inlocuit cu , astfel incat influenta
campului magnetic asupra evolutiei functiei de distributie este mult mai slaba decat in cazul
campului electric. In acest caz,
0f
11 )()( fef kBvk ∇⋅×−=h
τ
In particular, daca , , unde zB || xA || AB ×∇= ,
⎟⎠⎞
⎜⎝⎛∝→⋅=
Ckf
dkdfCf x
xexp1
11
D) Camp electric si magnetic stationar
Din ecuatia de transport Boltzmann avem
00 )()()( feff kBvEkk ∇⋅×++=h
τ
Ca si in cazul campului electric, cautam solutii de tipul ⎟⎠⎞
⎜⎝⎛ += Zkk
h
τeff 0)( , cu Z ce
urmeaza a fi determinat.
In acest caz
Functii de distributie de ne-echilibru. Exemple 5
01 fef kZ ∇⋅=h
τ
si ecuatia cinetica de transport devine
001)( fefefekkk ZEBv ∇⋅=∇⋅+∇⋅×
hhh
τττ
Stim ca pentru )](exp[1
10
FEEf
−+=
β, cu notatiile TkB/1=β si (de fapt
, dar neglijez ultimul termen in campuri magnetice mici)
m/kv h=
me /)( Akv += h
vvk hh )1( 000
0 ffdEdff −−=⎟
⎠⎞
⎜⎝⎛=∇ β
Ca sa obtin , avem, succesiv, 1f
)]()1()1()())(1[(
])1([
000000
00
01
vZvZvZ
vZ
Z
kkk
k
kkk
⋅∇−+−∇⋅+∇⋅−=
⋅−−∇⎟⎠⎞
⎜⎝⎛−=
∇⋅⎟⎠⎞
⎜⎝⎛−∇=∇
ffffffe
ffe
fef
τβ
βτ
τ
hh
h
Primii doi termeni dau contributie nula daca se inmultesc cu Bv × , si presupunand ca
ZvZk mh
≅⋅∇ )(
unde A este potentialul vector al campului magnetic avem
ZBvBv k mffef h
⋅×−=∇⋅× ))(1()( 001 τβ
Ecuatia care trebuie rezolvata este acum
0)1()())(1(
)())(1(
0000
000
=−−+⋅×−=
∇⋅−+⋅×−
vZEZBv
ZEZBv k
hh
h
ffffem
fffem
βτβ
τβ
Functii de distributie de ne-echilibru. Exemple 6
sau
ZBZE ×+=meτ
Z se numeste vector Hall.
Calculand intai EB × si apoi EB ⋅ se obtine Z ca functie de E si B, cu me /τμ −= :
22
2
1)(
Bμμμ
+
⋅+×+=
BEBEBEZ
Curentul se gaseste daca in legea Ohm se inlocuieste E cu Z: Zj σ= , cu μσ ne−=
Daca B || z, si campul electric este in planul (x,y)
zB ˆB= , yxE ˆˆ yx EE +=
avem
221 B
BEEZ yx
xμ
μ
+
−= ,
221 B
BEEZ xy
yμ
μ
+
+=
V
jx
x y
Vxx
Vxy
sau
⎟⎟⎠
⎞⎜⎜⎝
⎛⎟⎟⎠
⎞⎜⎜⎝
⎛ −
+=⎟⎟
⎠
⎞⎜⎜⎝
⎛
y
x
y
x
EE
BB
Bjj
11
1 22 μμ
μσ ,
relatie care poate fi inversata pentru a se obtine matricea rezistivitatii ( σρ /1= )
⎟⎟⎠
⎞⎜⎜⎝
⎛⎟⎟⎠
⎞⎜⎜⎝
⎛−
=⎟⎟⎠
⎞⎜⎜⎝
⎛
y
x
y
x
jj
neBneB
EE
ρρ
)/()/(
In geometria de mai sus, jy trebuie sa fie zero suficient de departe de contacte pentru ca nu
exista surse de electroni pe partile laterale ale probei, a.i. ρ=xx jE / este caderea de tensiune
longitudinala si este caderea de tensiune transversala (Hall), care este
proportionala cu B si invers proportionala cu densitatea de electroni.
)(/ neBjE xy −=
Functii de distributie de ne-echilibru in conductori neomogeni
In acest seminar ne intereseaza conductorii neomogeni, in care fie concentratia de purtatori,
fie temperatura variaza spatial. In acest caz, functia FD se refera la o stare de echilibru local
)](/))(exp[(11),(0 rr
krk TkEE
fBF−+
= ,
(de fapt, in locul energiei Fermi ar trebui sa avem potentialul chimic μ, dar am vazut ca cele
doua notiuni se folosesc cu aceeasi semnificatie in teoria solidului!) si
⎟⎠⎞
⎜⎝⎛
∂∂
−⎟⎠⎞
⎜⎝⎛ ∇
−+∇=∇
EfT
TEEEf F
F0
0 rrr
⎟⎠⎞
⎜⎝⎛∂∂
=⎟⎠⎞
⎜⎝⎛∂∂
∇=∇Ef
EfEf 00
0 vkk h
In cazul cel mai general, in prezenta campului electric si magnetic care produc o forta Lorentz
, la ne-echilibru in aproximatia timpului de relaxare, ecuatia cinetica
Boltzmann pentru electroni cu sarcina
)( BvEF ×+−= e
e− devine
τ10
111 )( f
EfT
TEEEefef
tf
fftf
FF −=⎟
⎠⎞
⎜⎝⎛
∂∂
−⎟⎠⎞
⎜⎝⎛ ∇
−+∇+⋅+∇⋅×+−∇⋅+
∂∂
=∇⋅+∇⋅+∂∂
rrkr
kr
EvBvEv
Fv
h
h
cu , ),(),(),( 10 krkrkr fff += ),(),( 01 krkr ff <<
Pornind de la legatura dintre campul electromagnetic si potential: ,
putem defini un potential electrochimic
Ue rE −∇=−
eEUU Fef /−= astfel incat campul electric efectiv
aplicat este
Fef Ee rEE ∇+= −1
Ecuatia de mai sus este o ecuatie integro-diferentiala neliniara, neliniaritatea provenind
de la termenul de imprastiere. Daca starea sistemului este aproape de starea de echilibru
ecuatia Boltzmann se poate liniariza: termenul neomogen in ecuatia de mai sus implica
gradientul temperaturii si a campului electric efectiv, ceea ce sugereaza ca f1 este proportional
Functii de distributii de ne-echilibru in conductori neomogeni 2
cu acesti parametri. Daca gradientii sunt mici, si f1 este mic, astfel incat gradientul sau este de
ordinul doi ca micime si se poate neglija, obtinandu-se ecuatia Boltzmann liniarizata:
τ10
11 )( f
EfT
TEEEefe
tf F
F −=⎟⎠⎞
⎜⎝⎛
∂∂
−⎟⎠⎞
⎜⎝⎛ ∇
−+∇+⋅+∇⋅×−
∂∂
rrk EvBvh
Exemplu 1
Functia de distributie in gradient de temperatura
In ecuatia de mai sus, considerand problema stationara si E = B = 0, EF = const. Rezulta
⎟⎠⎞
⎜⎝⎛ ∇
−⋅∇=⎟
⎠⎞
⎜⎝⎛∂∂
⎟⎠⎞
⎜⎝⎛ ∇
−⋅= T
TEEf
EfT
TEEf FF
rkrv 00
1h
ττ
In general, kkk k δ⋅∇+≅ 00 )()( fff , de unde rezulta ca si de aceasta data functia de
distributie este deplasata fata de valoarea la echilibru. Dar, spre deosebire de cazul in care am
aplicat doar un camp electric, deplasarea nu mai este aceeasi pentru toate valorile k. Pentru un
sistem in apropiere de echilibru, cu , ))(2/(2))(2/( 2222FFFF kkmkkkmEE −≅−=− hh
TkkmT
kF
Frk ∇−= )(hτδ
kx
ky
T
ceea ce inseamna ca functia de distributie a electronilor
de pe nivelul Fermi nu este afectata de gradientul de
temperatura, functiile de distributie pentru electronii cu
energii mai mari sau mai mici decat EF fiind deplasate
in sensuri opuse: in sensul sau in opozitie cu gradientul
temperaturii.
kvvkkkvkvkvjkk
TdkkdEdf
Tekdfef
Vefe
V FF ∇−∫ ⎟
⎠⎞
⎜⎝⎛−=∫−=∑−=∑−= )()(
4)(
4)(2)(1 0
3
2131
,τ
ππσ
h
Intr-un cristal izotrop cu relatie de dispersie parabolica si gradient de
temperatura dupa x,
mkE 2/22h=
ϕθ cossinkkx = , ϕθ sinsinkk y = , θcoskkz = , ϕθθ dddkkd sin2=k
Functii de distributii de ne-echilibru in conductori neomogeni 3
kk dkkvdEdf
dxdT
Tekj Fx
F )()(4
203
2−∫ ⎟
⎠⎞
⎜⎝⎛−= τ
πh
∫∫∫ −⎟⎠⎞
⎜⎝⎛−=
ππϕϕθθτ
π
2
0
2
0
3403
2cossin)()(
4dddkkkk
dEdf
dxdT
Tek
FF kh
∫ −⎟⎠⎞
⎜⎝⎛−=
max
min
)()(3
403
2 k
kF
F dkkkkdEdf
dxdT
Tek kτπh
Daca definesc coeficientul Seebeck ca
∫ ⎟⎠⎞
⎜⎝⎛−
∫ −⎟⎠⎞
⎜⎝⎛−
−=max
min
max
min
40
402
)(
)()(
k
k
k
kF
F
dkkdEdf
dkkkkdEdf
eTmkS
k
k
τ
τ
unde
∫ ⎟⎠⎞
⎜⎝⎛−=
max
min
4023
22)(
3
k
kdkk
dEdf
me kτπ
σ h
pot scrie curentul datorat gradientului de temperatura ca
TS rj ∇−= σ
Exemplu 2
Derivarea ecuatiei de drift-difuzie
In acest caz B = 0; T = const., sistem stationar. Pentru electroni
( ) ⎟⎠⎞
⎜⎝⎛
∂∂
−∇+⋅=−EfEef F0
1 rEvτ
Deoarece , intr-o prima aproximatie functia de distributie de
neechilibru este deplasata cu
01
0 )/( fEf kv ∇=∂∂ −h
))(/( FEe rEk ∇+= hτδ si densitatea curentului electric se poate
exprima ca
kEvvkkkvkvkvj rkk
dEedEdfedfef
Vefe
V F )()(4
)(4
)(2)(1 03131
,∇+∫ ⎟
⎠⎞
⎜⎝⎛−=∫−=∑−=∑−= τ
ππσ
Functii de distributii de ne-echilibru in conductori neomogeni 4
Intr-un cristal izotrop, la curs am vazut ca integrala dupa primul termen da curentul de drift
EEEj μστ neEm
nedrift ==⟩⟨= )(
2
cu mEe /)( ⟩⟨= τμ mobilitatea in aproximatia timpului de relaxare
Avand in vedere ca intr-un sistem cu T = const. nivelul Fermi nu poate varia decat
datorita unei concentratii neuniforme a electronilor: ndndEE FF rr ∇=∇ )/( , al doilea termen
al curentului, si anume curentul de difuzie, se poate exprima ca
neDndn
dEEmne F
dif rrj ∇=∇⟩⟨= )(τ
unde coeficientul de difuzie se defineste ca
μσdn
dEen
dndE
eD FF == 2
1
In metale sau semiconductori puternic degenerati,
∫+−
=∞
cE BF TkEEdEEmn
1]/)exp[(2)2( 2/1
32
2/3
hπ= 2/3
32
2/3)(
3)2(
cF EEm−
hπ
si deci
μe
EED cF −=
32
cu Ec = 0 pentru metale.
Din contra, pentru semiconductori intrinseci nedegenerati
∫ −−=∞
cEBF dEETkEEmn 2/1
32
2/3]/)(exp[
2)2(hπ
⎟⎠⎞
⎜⎝⎛ −−=
TkEETmk
B
FcB exp4
)2(33
2/3
hππ
si
Functii de distributii de ne-echilibru in conductori neomogeni 5
μeTkD B=
formula cunoscuta ca relatia lui Einstein.
In concluzie, pentru electroni (introduc indicele n)
neDen nnn ∇+= Ej μ ,
o relatie similara existand pentru goluri:
peDep ppp ∇−= Ej μ
Observatie: expresia coeficientului de difuzie depinde de statistica distributiei de purtatori.
Pentru nivele donoare, in care electronii au o functie de distributie
]/)exp[(11)( 1 TkEEg
EfBFdd
d−+
= −
coeficientul de difuzie depinde de temperatura, deoarece concentratia electronilor si relatia
acestora cu nivelul Fermi depinde de temperatura si se determina din conditia de neutralitate
dd NpNn +=+ 0 ,
cu
]/)exp[(1)()( 1
0
TkEEgNdEEENEfN
BFdd
ddddd −+
=∫ −=−
δ
1 2 3
E
T Ts Ti
Ev
Ec
Ed
Eg/2
Egd/2
n
T Ts Ti
intrinsic extrinsic
freeze-out
Functii de distributie de ne-echilibru. Transport ambipolar
Ecuatia cinetica Boltzmann pentru purtatori de sarcina cu sarcina q in conditii stationare se
poate scrie ca
τ1
100 )( ffqfqfff −=∇⋅×+∇⋅+∇⋅=∇⋅+∇⋅ kkrkr BvEvFvhhh
La curs am vazut ca perturbatia functiei de distributie
)](/))(exp[(11),(0 rr
krk TkEE
fBF−+
=
in conditii de neechilibru este de forma
)()()( 00
1 kkk
kχk EfEdEdff <<⋅−=
Vrem sa vedem ce expresie are χ si in ce mod intervine in gasirea coeficientilor de transport
la conductia ambipolara in semiconductori, la care iau parte electronii si golurile. Parametrii
referitori la cele doua tipuri de purtatori de sarcina vor fi indexati cu n si, respectiv, p.
Tinand cont de faptul ca
⎟⎠⎞
⎜⎝⎛
∂∂
−⎟⎠⎞
⎜⎝⎛ ∇
−+∇=∇
EfT
TEEEf F
F0
0 rrr , ⎟⎠⎞
⎜⎝⎛∂∂
=⎟⎠⎞
⎜⎝⎛∂∂
∇=∇Ef
EfEf 00
0 vkk h
)()()( 0001 kχkχkχ kkkk ⋅∇⎟
⎠⎞
⎜⎝⎛∂∂
−⎟⎠⎞
⎜⎝⎛∂∂
∇⋅−=⎟⎠⎞
⎜⎝⎛ ⋅
∂∂
−∇=∇Ef
EfE
Eff ,
primul termen din da rezultat nul cand este inmultit cu 1fk∇ )( Bv × , iar al doilea termen
vχkχχkχkχ kk h⎟⎠⎞
⎜⎝⎛
∂∂⋅+=∇⋅+=⋅∇
EE)()(
astfel incat ecuatia cinetica Boltzmann pentru electronii de conductie, cu si
sarcina devine (am folosit egalitatea
vk )/( hnm=
eq −= vχBχBv )()( ×=× , si notam ) r∇=∇
Functii de distributii de ne-echilibru. Transport ambipolar 2
vχvχBE ⋅∂∂
=⎥⎦⎤
⎢⎣⎡ ×+−∇−∇
−∂∂
nn
nnF
F mEf
EeeET
TEE
Ef
hh00
)(1)(
τ
Din aceasta ecuatie vectorul χ se gaseste ca
⎥⎦⎤
⎢⎣⎡ ×−−∇+∇
−−= )()()()( nF
F
nnn
eeETT
EEm
EE χBχh
h φτ
unde φ este potentialul electrostatic definit prin φ−∇=E .
Similar, pentru goluri cu energie 'E , masa efectiva , sarcina si functie de
distributie la echilibru
pm eq =
Ec
Ev
-EF 0
-Eg E’
1]/)exp[(1)(1)( 0
+−=−=
TkEEEfEf
BFp ,
⎟⎠⎞
⎜⎝⎛∂∂
⎟⎠⎞
⎜⎝⎛ ∇
−−∇=∇
EfT
TEEEf F
Fp0
rrr
cu , din ecuatia cinetica Boltzmann obtinem gFF EEEEE ++=− '
⎥⎦
⎤⎢⎣
⎡×+−∇−∇
++−= )()(
')'()'( pF
gF
ppp
eeETT
EEEm
EE χBχh
h φτ
Pentru a gasi χ din ecuatii de forma χbaχ ×+= , se inmulteste scalar cu b, rezultand
abχb ⋅=⋅ (*)
deoarece , si se substituie apoi χ in partea dreapta cu 0)( =×⋅ χbb χba ×+ :
)( χbbabaχ ××+×+=
Folosind identitatea si (*) din ecuatia de mai sus rezulta 2)()( bχχbbχbb −⋅⋅=××
21)(
b+⋅⋅+×+
=bbaabaχ ,
adica,
Functii de distributii de ne-echilibru. Transport ambipolar 3
22
32
)(1)()()()(
Bnn
nnnnnnnnnn
τατατατα
+
⋅+×+−=
BPBPBPχ h
cu
⎟⎠⎞
⎜⎝⎛ −∇+∇
−= φ
eET
eTEE FF
nP , n
n me
=α
si
22
32
)(1
)()()()(
Bpp
pppppppppp
τα
τατατα
+
⋅+×−=
BPBPBPχ h
cu
⎟⎠⎞
⎜⎝⎛ −∇+∇
++−= φ
eET
eTEEE FFg
p'
P , p
p me
=α
Densitatea curentului electric
kkvkvkvjkk
dfqefVqfq
V ∫=∑=∑= )(4
)(2)(1131
, πσ
Pentru electroni
kkχj kdEf
me
nn
n ∫ ⋅⎟⎠⎞
⎜⎝⎛
∂∂
−−= )(4
03h
π
In coordonate polare cu axa polara paralela cu , pentru care nχ 34sincos
0
22
0
πθθθφππ
=∫∫ dd ,
dkkEf
mej n
nn
4
0
023
∫ ⎟⎠⎞
⎜⎝⎛
∂∂
−−=∞
χπ
h
Intr-un semiconductor pentru care
⎟⎠⎞
⎜⎝⎛ −−=
TkEENn
B
Fcc exp , 3
2/3)2(2h
TkmN Bnc
π=
Functii de distributii de ne-echilibru. Transport ambipolar 4
⎟⎠⎞
⎜⎝⎛ −−=
TkEEEf
B
Fexp)(0 ,
⎟⎠⎞
⎜⎝⎛ −−⎟
⎠⎞
⎜⎝⎛ −−−=⎟
⎠⎞
⎜⎝⎛ −−−=
∂∂
TkEE
TkEE
TkTkEE
TkEEf
B
c
B
Fc
BB
F
Bexpexp1exp1)(0
luand energia de referinta si facand schimbarea de variabila 0=cE TkEz B/= obtinem, in
general,
⟩⟨−= nnen χjh
,
cu ∫ −
∫ −=⟩⟨
∞
∞
0
2/3
0
2/3
)exp(
)exp('
dzzz
dzzzn
n
χχ
Analog, ⟩⟨= ppep χjh
Exemple
1) Curent electric de drift: B = 0, 0=∇=∇ TEF
Eχ nn
nn
n me
me τφτ hh
−=−∇−= )( , EEEj nnnn
n enmne μστ ==⟩⟨=
2
Eχ pp
pp
p me
me τφτ hh
=−∇= )( , EEEj pppp
p epmpe μστ ==⟩⟨=
2
⎪⎩
⎪⎨⎧
+=++
=→+=+=+=
pn
pn
pnpnpn pnpn
pneσσσ
μμμσσμμ EEjjj )()(
2) Curent de difuzie: E = B = 0, 0=∇T
ndn
dEm
Em
Fn
nFn
nn ∇−=∇−= ττ hhχ , neDn
dndE
mne
nF
nn
n ∇=∇⟩⟨= τj
pdp
dEm
Em
Fp
pFp
pp ∇−=∇= ττ hhχ , peDp
dpdE
mpe
pF
pp
p ∇−=∇⟩⟨−= τj
Functii de distributii de ne-echilibru. Transport ambipolar 5
3) Efect Hall: 0=∇=∇ TEF , BE ⊥
EPP =−∇== φpn
22
2
)(1)()(Bnn
nnnnn
τατατα
+
×+−=
EBEχ h , 22
2
)(1
)()(
Bpp
ppppp
τα
τατα
+
×−=
EBEχ h
),0,0( B=B , , )0,,( yx EE=E BEE xy )ˆˆ( yxEB +−=× Pentru electroni:
⎪⎪
⎩
⎪⎪
⎨
⎧
+=
+=
→⎩⎨⎧
+=
−=
2
222
2
21
12
21
)(1)(
)(1
BB
mnea
Bmnea
EaEajEaEaj
nn
nn
n
nn
n
n
yxy
yxx
τατα
τατ
Presupun , si notez . Din 0=yj nx jj =
nnHy
y
x
y
x BjRaa
aEj
jaaaa
aaEE
=+
−=→⎟⎟⎠
⎞⎜⎜⎝
⎛=⎟⎟
⎠
⎞⎜⎜⎝
⎛−+
=⎟⎟⎠
⎞⎜⎜⎝
⎛22
21
2
12
2122
21
01
de unde obtinem expresia constantei Hall in aproximatia campului magnetic slab, pentru care
1<<Bnnτα , caz in care neglijam si in : 22a 2
221 aa +
2
2
22
21
2 1)( ⟩⟨
⟩⟨−=
+−=
n
nnH neBaa
aRττ
In campuri magnetice intense 1>>Bnnτα , n
m
Bnma
τ1
21 = , 12 aBnea >>= ,
neRn
H1
−= .
Pentru cazul in care am conductie ambipolara:
2
2
2
21
)(1)(1 Bmpe
Bmnea
pp
p
pnn
n
n τα
τ
τατ
++
+=
2
22
2
222
)(1)(
)(1)(
BB
mpe
BB
mnea
pp
pp
pnn
nn
n τα
τα
τατα
++
+=
Functii de distributii de ne-echilibru. Transport ambipolar 6
Pentru campuri magnetice slabe, tinand cont ca jjj me /⟩⟨= τμ , jj ejμσ = , j = n,p
2
22
2
2
22
2
22
2
22
22
2
2
)()(11
pn
ppHn
nH
pn
p
pp
n
nn
pp
nn
pp
nn
HRR
pn
pn
e
mpe
mne
mpe
mne
eR
σσ
σσ
μμ
τ
τμ
ττμ
ττ
ττ
+
+=
+
⟩⟨
⟩⟨−
⟩⟨
⟩⟨
−=
⎟⎟⎠
⎞⎜⎜⎝
⎛⟩⟨+⟩⟨
⟩⟨−⟩⟨
−=
Daca timpul de relaxare nu depinde de energie, astfel incat , 2
,2, ⟩⟨≅⟩⟨ pnpn ττ
2
22
)(1
pn
pnH
pn
pne
Rμμ
μμ
+
−−=
expresie care nu se anuleaza nici pentru semiconductorii intrinseci pentru ca mobilitatile
electronilor si golurilor sunt in general diferite.
Pentru campuri magnetice puternice,
Bepna
Bpm
Bnma
pn
nn )(1111
2221 −=<<+=ττ
si
)(1
pneRH
−−= .
Pentru semiconductori intrinseci, cu n = p, constanta Hall se anuleaza in campuri
magnetice puternice. In acest caz este nevoie sa calculam pana la termeni de ordinul : 2a 3−B
⎥⎥⎦
⎤
⎢⎢⎣
⎡−−= 2
22
232
11
pp
nn mm
eBna
ττ,
si
2
22
22
11
11
1
⎟⎟⎠
⎞⎜⎜⎝
⎛+
−
−=
nn
pp
nn
pp
H
mm
mm
enR
ττ
ττ
Functii de distributii de ne-echilibru. Transport ambipolar 7
4) Efect magnetorezistiv: 0=∇=∇ TEF , BE ⊥ .
Efectul magnetorezistiv se refera la modificarea rezistivitatii conductorilor ρ in prezenta
campului magnetic. Clasic, aceasta se datoreaza curbarii traiectoriei electronilor datorata
fortei Lorentz, care duce la scaderea distantei de propagare pe directia campului electric.
Consider,
),0,0( B=B , , )0,0,( xE=E BExyEB ˆ=×
Rezistivitatea in camp magnetic se defineste ca xxH jE /=ρ .
11 −=−=−
=Δ
x
x
x
xH
jE
jE σρρ
ρρρρ
Similar ca in cazul efectului Hall, pentru electroni 0=yj , si din
10
122
21
122
21
1
12
2122
21
−+
=Δ
→+
=→⎟⎟⎠
⎞⎜⎜⎝
⎛=⎟⎟
⎠
⎞⎜⎜⎝
⎛−+
=⎟⎟⎠
⎞⎜⎜⎝
⎛
aaa
aaaE
jj
aaaa
aaEE
xy
x
y
x σρρ
2
21
)(1 Bmnea
nn
n
n τατ
+= , 2
222
)(1)(
BB
mnea
nn
nn
n τατα
+= , ⟩⟨= n
nn m
ne τσ2
Observatie: daca BE ||
EEEχ ∝+
+−= 22
23
)(1)()(
BB
nn
nnnnn
τατατα
h , EEE
χ ∝+
+= 22
23
)(1
)()(
B
B
pp
ppppp
τα
ταταh
si nu exista efect Hall si nici efect magnetorezistiv In campuri magnetice slabe, efectul magnetorezistiv se anuleaza in aproximatia de ordin intai
in B. In aproximatia de ordin doi,
⟩−⟨≅ ])(1[ 22
1 Bmnea nnn
nτατ , ⟩⟨≅⟩−⟨≅ 2
222
22 )(])(1[)( nn
nnnnn
nB
mneBB
mnea τατατα
][ 22322
2
2
2
2232
⟩⟨−⟩⟩⟨⟨⎟⎟⎠
⎞⎜⎜⎝
⎛=
⟩⟨
⟩⟨−⟩⟩⟨⟨⎟⎠
⎞⎜⎝
⎛≅Δ
nnnnnnn
nnn
n mneB
meB τττα
σττττ
ρρ .
Functii de distributii de ne-echilibru. Transport ambipolar 8
efectul magnetorezistiv anulandu-se doar daca timpul de relaxare nu depinde de energie.
In campuri magnetice intense, in care n
m
Bnma
τ1
21 = , 12 aBnea >>= , si
magnetorezistenta nu depinde de campul magnetic aplicat:
11122
1 −⟩⟨=−=Δ
nn
aa
ττσ
ρρ .
In semiconductori cu conductie mixta, in campuri magnetice slabe
⟩−⟨+⟩−⟨≅ ])(1[])(1[ 22
22
1 BmpeB
mnea ppp
pnnn
nταττατ
⟩−⟨+⟩−⟨≅ ])(1[)(])(1[)( 222
222
2 BBmpeBB
mnea pppp
pnnnn
nτατατατα
⎥⎥
⎦
⎤⟩⟩⟨⟨+⟩⟩⟨⟨⎟
⎟⎠
⎞⎜⎜⎝
⎛⎟⎟⎠
⎞⎜⎜⎝
⎛+⟩⟨−⟩⟩⟨⟨⎟
⎟⎠
⎞⎜⎜⎝
⎛+
⎢⎢
⎣
⎡⟩⟨−⟩⟩⟨⟨⎟
⎟⎠
⎞⎜⎜⎝
⎛≅−
+=
Δ
][][
][1
323222
22322
22322
2
2
22
21
1
pnpnpnppn
ppppp
nnnnnn
mpe
mne
mpe
mneB
aaa
ττατταατττα
τττασ
σρρ
In campuri magnetice intense Bepna
Bpm
Bnma
pn
nn )(1111
2221 −=<<+=ττ
1)(
11
1 2222
1 −−
⎟⎟⎠
⎞⎜⎜⎝
⎛+
=−≅Δ
pne
pmnm
aa p
pn
n ττσ
σρρ
Iar magnetorezistenta nu depinde de campul magnetic. Expresia de mai sus nu este valabila
pentru n = p. In acest caz, deoarece 02 =a
111
)(111
1 22
1−
+
⟩⟨+⟩⟨=−
+
=−≅Δ
pp
nn
pp
nn
pp
nn pmnm
mp
mn
eBpmnm
Ba
ττ
ττ
ττ
σσρρ
Functii de distributii de ne-echilibru. Transport ambipolar 9
5) Efectele Nernst si Ettinghausen
Aceste doua efecte, ca si efectul Hall si cel magnetorezistiv, sunt efecte galvanomagnetice.
Acest tip de efecte apar in prezenta simultana a campurilor electric si magnetic. Efectele
Nernst si Ettinghausen sunt insotite de efecte termice care apar datorita faptului ca, din punct
de vedere clasic, forta Lorentz compenseaza actiunea campului Hall doar pentru purtatorii cu
viteza medie. In consecinta, datorita distributiei dupa viteze a purtatorilor, electronii sau
golurile mai rapide vor fi deviate spre fata opusa celei spre care vor fi deviati purtatorii de
sarcina mai lenti. Deoarece purtatorii de sarcina ajung la un echilibru termic cu reteaua
cristalina prin procesele de ciocnire, o parte a conductorului se incalzeste si cealalta se
raceste, aparand un gradient de temperatura pe o directie perpendiculara atat pe campul
cat si pe curentul : ),0,0( B=B )0,0,( xj=j
BjCyT
xE=∂∂
Acesta este efectul Ettinghausen, care este un efect galvanomagnetic transversal.
Efectul Nernst este un efect galvanomagnetic longitudinal cauzat de micsorarea
fluxului de electroni mai rapizi si mai lenti pe directia campului electric datorita fortei
Lorentz, care duce la aparitia unui gradient de temperatura pe directia curentului:
2BjCxT
xN=∂∂
Spre deosebire de efectul Ettinghausen, efectul Nernst nu depinde de orientarea campului
magnetic, ambele efecte depinzand de directia de curgere a curentului electric.