F6121 Základy fyziky pevných látek — příklady do cvičení 1 Drudeho model volných elektronů 1 1.1 Poissonovo rozdělení............................................ 1 1.2 Jouleho teplo................................................ 1 1.3 Elektrická vodivost kovů.......................................... 1 2 Sommerfeldův model volných elektronů 2 2.1 Nízkorozměrný elektronový plyn ..................................... 2 2.2 Betheho-Sommerfeldův rozvoj ...................................... 2 2.3 Teplotní závislosti v Sommerfeldově modelu............................... 2 2.4 Číselné odhady; 2D plyn volných elektronů............................... 2 2.5 Tepelná vodivost.............................................. 2 3 Krystalová struktura 3 3.1 Kuprátové roviny.............................................. 3 3.2 Operace symetrie 2D krystalových mřížek................................ 3 3.3 Osy rotace v prostorových mřížkách................................... 3 3.4 Hustota diamantu ............................................. 3 3.5 Součinitele zaplnění ............................................ 3 3.6 Hexagonální těsně uspořádaná mřížka.................................. 3 3.7 Dvourozměrné mřížky........................................... 4 3.8 Reciproké mřížky a 1. Brillouinova zóna kubických mřížek....................... 4 3.9 Mezirovinné vzdálenosti a úhly...................................... 4 4 RTG difrakce na krystalech 5 4.1 Strukturní faktor, vyhasínání difrakcí................................... 5 4.2 Difrakční úhly................................................ 5 4.3 Difrakční efekty spojené s konečnou velikostí krystalu ......................... 5 4.4 RTG difrakce na A^Ceo.......................................... 5 5 Elektron v periodickém potenciálu 6 5.1 Jednorozměrný potenciál.......................................... 6 5.2 Model Kroning-Penney.......................................... 6 5.3 Dvourozměrná Fermiho plocha...................................... 7 5.4 Pásové schéma volných elektronů..................................... 7 5.5 Odhad šířky zakázaného pásu - metoda téměř volných elektronů................... 7 5.6 Metoda těsné vazby pro s-pás v fcc mřížce ............................... 7 5.7 Metoda těsné vazby pro p-pásy ve čtvercové mřížce........................... 7 6 Kvaziklasická aproximace 8 6.1 Elektrony v okolí minima pásu...................................... 8 6.2 Oscilace v homogenním elektrostatickém poli.............................. 8 7 Polovodiče 9 7.1 Příměsový stav............................................... 9 7.2 Statistika nositelů náboje v polovodiči typu N.............................. 9 7.3 Hallův jev pro dva typy nositelů ..................................... 9 7.4 Intrinsický polovodič............................................ 9 8 Kmity mřížky v harmonickém přiblížení 10 8.1 Kmity dvouatomového řetězce....................................... 10 8.2 Kmity lineární mřížky s dalekodosahovou interakcí........................... 10 8.3 Konstantní rychlost zvuku......................................... 10 8.4 Měkký fononový mód ........................................... 10 8.5 Rychlost zvuku v křemíku......................................... 10 8.6 Tepelná kapacita jednoduché ID a 2D mřížky.............................. 11 8.7 Hustota stavů akustické fononové větve................................. 11 1 Drudeho model volných elektronů 1.1 Poissonovo rozdělení V Drudeho modelu je pravděpodobnost, že se elektron srazí za elementární časový úsek dí, rovna dí/r. 1. Dokažte, že elektron libovolně vybraný v daný časový okamžik se nesrazil v předchozích t sekundách s pravděpodobností e_í/T. 2. Dokažte, že pravděpodobnost toho, že doba mezi dvěma následujícími srážkami je v intervalu (í,í + dí), jee_*/Tdí/r. 3. Dokažte, že doba od poslední srážky vystředovaná přes všechny elektrony je r. 4. Dokažte, že střední doba mezi dvěma srážkami pro libovolně vybraný elektron je r. 1.2 Jouleho teplo Kus kovu se nachází v homogenním elektrostatickém poli E, teplota kovu je konstantní. Vyberme libovolný elektron z elektronového plynu a předpokládejme, že tento elektron vykonal srážku v čase t = 0 a další srážku v čase í. 1. Dokažte, že střední energie předaná elektronem při druhé uvažované srážce je (eEt)2/2m. 2. Dokažte, že střední energie předaná elektronem při libovolné srážce je (eEr)2/m. 3. Nechť má kus kovu tvar válce s plochou podstavy S a výškou L a nechť je intenzita elektrického pole E rovnoběžná s výškou válce. Z výsledku části 2 odvoďte vztah pro elektrický odpor válce. 4. Najděte tepelný výkon generovaný při průchodu proudu a ověřte, zda v Drudeho modelu platí známý vztah P = RI2. 1.3 Elektrická vodivost kovů 1. Vypočtěte hustotu volných elektronů v mědi, je její hustota pcu = 8960kgm-3 a relativní atomová hmotnost 63.5. 2. Měděným vodičem s příčným průřezem 0.2cm2 prochází proud 1A. Jaká je střední driftová rychlost elektronů? 3. Vypočtěte pohyblivost elektronů v sodíku, je—li jeho specifická vodivost a = 0.23 • 108 Í2_1 m-1 a koncentrace nositelů náboje 2.652 • 1028 m-3. 4. Specifická elektrická vodivost mědi je a = 6 • 107 Í2_1 m_1. Určete relaxační dobu elektronu. 5. Určete střední volnou dráhu vodivostních elektronů v sodíku. Jeho specifická vodivost je a = 0.23 • 108 fi_1 m_1. 1 2 Sommerfeldův model volných elektronů 2.1 Nízkorozměrný elektronový plyn Pro jednorozměrný, dvourozměrný a trojrozměrný plyn volných elektronů najděte: 1. souvislost Uf a £f a hustoty elektronů n (počet elektronů na jednotku délky, plochy resp. objemu) 2. souvislost Iíf a veličiny rs definované jako poloměr koule1 s objemem rovným objemu připadajícímu v elektronovém plynu na jeden elektron 3. energiovou hustotu stavů g(E) Pozn.: Vzájemnou konzistentnost výsledků je možné ověřit vztahem J0 F g(E) áE = n. 2.2 Betheho—Sommerfeldův rozvoj Ukažte, že integrál /0°° H(E)jfd{E) áE je možné aproximovat rozvojem r°°„,„, . ,„, ,„ [» „,„,,„ tt2 „ „,,„„, 7vr4 r°° r^ TT2 7-7T4 J H(E)fFD(E)dE = J H{E)aE+-{kBTfH'{ß) + —{kBTyH"'{ß) + 0 kBT\{ P J 2.3 Teplotní závislosti v Sommerfeldově modelu Pomocí Betheho-Sommerfeldova rozvoje určete teplotní závislost chemického potenciálu /x(T), střední hodnoty hustoty energie u(T) a tepelnou kapacitu 3D elektronového plynu. Předpokládejte přitom, že v uvažovaném intervalu teplot je T/Tp Cla stačí tedy vzít pouze první opravu z příkladu 2.2. 2.4 Číselné odhady; 2D plyn volných elektronů 1. S využitím předchozích výsledků spočtěte Fermiho mez Iíf, Fermiho energii £p, Fermiho rychlost vf, Fermiho teplotu Tp, střední energii elektronu (E) a hustotu energie u v elektronovém plynu s hustotou odpovídající stříbru (n = 5.85 • 1028 m-3). Dále stanovte chemický potenciál a střední hustotu energie při teplotě 300 K. Jaká je tepelná kapacita elektronového plynu? Porovnejte se skutečnou hodnotou a pokuste se vysvětlit případný rozdíl. 2. Určete chemický potenciál dvourozměrného elektronového plynu. Díky příznivému průběhu hustoty stavů není v tomto případě třeba aproximací. 2.5 Tepelná vodivost Spočtěte tepelnou vodivost elektronového plynu v Sommerfeldově modelu. Porovnejte výsledek s tepelnou vodivostí v Drudeho modelu kovu a s tabulkovými hodnotami pro reálné kovy. 1rozumí se zobecněná Z?-dimenzionální koule, která je zadána vztahem \/x\ + ... + x2D < 2 3 Krystalová struktura 3.1 Kuprátové roviny Ve většině vysokoteplotních supravodičů se lze setkat s tzv. kuprátovými rovinami, které jsou tvořeny atomy mědi a kyslíku uspořádanými jako na následujících obrázcích. Atomy mědi jsou znázorněny plnými kroužky, atomy kyslíku prázdnými. 1. Vyznačte bázové vektory, primitivní buňku a atomy báze krystalové mřížky z levého obrázku. 2. Ve skutečnosti leží atomy kyslíku střídavě nad a pod kuprátovou rovinou, což je v pravém obrázku označeno znaménky + a —. Vyznačte bázové vektory, primitivní buňku a atomy báze i v tomto případě. • o • o • o • • ©•©•©• oooo e©©© • o • o • o • •©•©•©• OOOO ©OffiO • o • o • o • • ©•©•©• OOOO OffiO® • o • o • o • •©•©•©• 3.2 Operace symetrie 2D krystalových mřížek Najděte všechny bodové operace symetrie následujících krystalových mřížek a srovnejte je s operacemi symetrie prosté mřížky. F^h džh F^h ?šh Jk 3.3 Osy rotace v prostorových mřížkách Dokažte, že trojrozměrné prostorové mřížky mohou mít pouze 2-, 3-, 4- a 6-četné osy symetrie. 3.4 Hustota diamantu Spočtete hustotu diamantu, víte-li, že jeho mřížkový parametr je a = 3.57 Á a relativní atomová hmotnost uhlíku je 12. 3.5 Součinitele zaplnění Vypočtěte součinitele zaplnění při umístění koulí maximálního poloměru do uzlů prostorové mřížky pro tyto mřížky: prostá kubická, kubická plošně centrovaná, kubická prostorově centrovaná a diamantová. 3.6 Hexagonální těsně uspořádaná mřížka Vypočtěte poměr c/a pro hexagonální těsně uspořádanou mřížku. 3 3.7 Dvourozměrné mřížky Najděte reciproké mřížky a několik prvních Brillouinových zón (alespoň pět) pro dvourozměrnou čtvercovou a hexagonální mřížku. 3.8 Reciproké mřížky a 1. Brillouinova zóna kubických mřížek Najděte reciproké mřížky a první Brillouinovu zónu pro kubické mřížky - prostou, prostorově centrovanou a plošně centrovanou. Porovnejte rozměry první Brillouinovy zóny prosté mřížky s mřížkovým parametrem a = 3 Á a vlnová čísla typická pro viditelné a RTG záření. Najděte souvislost objemu primitivní buňky přímé a reciproké mřížky. Pozn.: Tato souvislost je obecná a lze ji nejsnadněji získat přímo pomocí definice reciproké mřížky. 3.9 Mezirovinne vzdálenosti a úhly Najděte vztahy pro mezirovinne vzdálenosti a mezirovinne úhly pro tyto syngonie: kubická, tetragonální a ortorombická. 4 4 RTG difrakce na krystalech Intenzita RTG záření rozptýleného krystalem je úměrná kvadrátu absolutní hodnoty Fourierova obrazu elektronové hustoty s argumentem rovným rozptylovému vektoru I-\pFJ{Q)\\ Q = Kf-Kl. Zapišme elektronovou hustotu v krystalu ve tvaru N ßei(r) = tt{r)^^pj(r-rj-R) , R j=l kde fž(r) je tvarová funkce krystalu, R značí vektory poloh elementárních buněk (tvoří prostorovou mřížku), j indexuje atomy v elementární buňce, r j jejich polohy v rámci elementární buňky a Pj(r) je nábojová hustota charakteristická pro daný atom. Fourierova transformace dává N pFe?(Q) = y- J2nFT(Q - G)F^G) > F^G) = Ee"ÍG'rj^(G) > MQ = / d3r p,(r)e-ifc- . PB G j=1 J Zde VpB je objem primitivní buňky, Q,FT geometrický faktor, F(G) je strukturní faktor a f j je atomový rozptylový faktor atomu j. 4.1 Strukturní faktor, vyhasínání difrakcí Vypočtěte strukturní faktor difrakce na krystalech s touto strukturou: kubická plošně centrovaná mřížka, kubická prostorově centrovaná mřížka, mřížka se sfaleritovou strukturou a mřížka s diamantovou strukturou. Zjistěte, které difrakce vyhasnou. 4.2 Difrakční úhly Vypočtěte všechny možné difrakční úhly při difrakci záření o vlnové délce 0.1541 nm (charakteristická čára CuKai) na krystalu Si (mřížkový parametr 0.54309 nm). 4.3 Difrakční efekty spojené s konečnou velikostí krystalu Spočítejte Fourierovu transformaci nábojové hustoty malého krystalu s prostou kubickou mřížkou Nx Ny Nz Pei(x,y,z)= Y^ Yl Yl Po(x-Jia,y-J2a,z-j3a) ji = -Nx J2=-Ny js=-Nz a srovnejte výsledek s výrazem pro pFT obsahujícím geometrický faktor. Odpovídá elektronová hustota zadaná v tomto příkladu elektronové hustotě uvažované výše, nebo je zde nějaký rozdíl? Diskutujte o souvislosti konečné velikosti krystalu s divergencí rozptýleného rentgenového záření. 4.4 RTG difrakce na A^Ceo Experimentálně bylo zjištěno, že difrakční pík (200) fcc mřížky fulerenu Cßo (mřížkový parametr a = 14.11 Á) je velmi slabý. Předpokládejte, že nábojová hustota fulerenu je reprezentována nábojem rovnoměrně rozloženým na povrchu koule s poloměrem 3.5 Á. Spočítejte strukturní faktor molekuly Cßo v této aproximaci a s jeho pomocí ukažte, že difrakční pík (200) je mnohem slabší než pík (111). 5 5 Elektron v periodickém potenciálu 5.1 Jednorozměrný potenciál Metodou rozvoje do rovinných vln najděte vlastní energie elektronu v jednodimenzionálním potenciálu s periodou a zadaném funkcí U(x) = —Vq 2_j exP (x — naý jehož Fourierovy složky jsou UG = -VoVir - exp-----— a V 4 G = 2irn Z vlastních energií pro dostatečný počet Blochových vektorů v 1. Brillouinově zóně sestavte pásové schéma. Při numerickém řešení použijte následující hodnoty parametrů: a = 0.5nm, a = 0.1a. Srovnejte výsledky pro Vq = 2 eV a Vq = 10 eV s disperzními relacemi volných elektronů. x [nm] Pozn.: Při srovnávání je výhodné použít energii vztaženou na střední hodnotu potenciálu, tj. E — Ug=o- 5.2 Model Kroning—Penney Vyřešte Schrödingerovu rovnici pro elektron v jednorozměrném periodickém potenciálovém poli, které má tvar U(x) U0 0 a a + b Ukažte, že vlastní hodnoty energie jsou dány rovnicí cos k(a + b) = cos an\ cos bn2-----(-----1-----) sin a«i sin 6^2 2 \Ki K-2 kde «i = V2mĚ a K2 = ^2m(E - UQ) n n Najděte disperzní relaci En(k) pro několik nejnižších pásů. Vyřešte problém numericky pro vhodně zvolené číselné konstanty, např. a = 4Á, 6 = lÁaVb=5 eV. 6 5.3 Dvourozměrná Fermiho plocha Na následujícím obrázku jsou zachyceny energiove pásy pro elektron pohybující se v dvourozměrné analogii potenciálu z úlohy 5.1. S použitím těchto grafů načrtněte Fermiho plochu pro látku s jedním, dvěma a třemi elektrony v primitivní buňce. Pro ilustraci je připojen graf příspěvků jednotlivých pásů do hustoty stavů. Rozsahy energií prvních čtyř pásů jsou 0.00 eV — 2.07 eV, 1.93 eV — 5.23 eV, 3.40 eV — 6.89 eV a 3.46 eV - 8.04 eV. E(k) [eV] 5.4 Pásové schéma volných elektronů Uvažme prostou kubickou mřížku. Sestrojte redukované pásové schéma pro volné elektrony ve směru [100] a [111]. 5.5 Odhad šířky zakázaného pásu — metoda téměř volných elektronů Uvažme dvourozměrnou čtvercovou mřížku s potenciálem U(x,y) = —4ř7ocos(^) cos(^) . Najděte přibližnou velikost šířky zakázaného pásu v bodě M = (f > f )> tj. v rohu první Brillouinovy zóny. 5.6 Metoda těsné vazby pro s—pás v fcc mřížce Odvoďte disperzní relace pásu vycházejícího z s-stavů atomů umístěných v uzlech kubické plošně centrované mřížky. Uvažujte pouze maticové elementy mezi nejbližšími sousedy t = (s|Ař7|s'). Překryv s-orbitalů na sousedních atomech zanedbejte. Výsledek znázorněte graficky obvyklým způsobem, tj. podél lomené čáry L — F — X — K — F. 5.7 Metoda těsné vazby pro p—pásy ve čtvercové mřížce Uvažujme o dvourozměrné čtvercové mřížce s jednoatomovou bází. Najděte disperzní relace pásů odvozených z dvakrát degenerovaných p-orbitalů px a py. Vlnové funkce těchto orbitalů mají tvar ipPx (x, y) = x f(\/x2 + y2) a ipPy(x,y) = y f(\/x2 + y2). Při výpočtu se omezte pouze na maticové elementy mezi nejbližšími sousedy a matici překryvových integrálů aproximujte jednotkovou maticí. Pásové schéma zobrazte podél lomené čáry M — F — X. 7 6 Kvaziklasická aproximace 6.1 Elektrony v okolí minima pásu Pro elektrony v okolí minima pásu platí 2 (m^ 0 0 \ E(k) = E(ko) + —(k-ko)TM-1(k-k0) M~l = 0 m^1 0 , V 0 0 ml1) kde rriT a rriL jsou transverzální a longitudinální efektivní hmotnosti. Ekvienergiové plochy mají tedy tvar rotačních elipsoidů. 1. Ukažte, že cyklotronová frekvence je eB U!c = ------^=^= , 2y/mTmL leží-li homogenní magnetické pole v rovině xy. 2. Vypočtěte elektronovou tepelnou kapacitu. 6.2 Oscilace v homogenním elektrostatickém poli Elektrony vodivostního pásu odvozeného od s-orbitalů atomů v prosté kubické mřížce mají v přiblížení těsné vazby disperzní relaci E(k) = Es —2t [cos^aio) + cos(kyd) + cos^a)] . Najděte časový průběh rychlosti a polohy elektronu v homogenním elektrickém poli E = (Ex, 0, 0), je-li toto pole zapnuto v čase t = 0, kdy se elektron nachází ve stavu s k = (0,0,0). Jaký je příspěvek elektronu do elektrické vodivosti materiálu? 8 7 Polovodiče 7.1 Příměsový stav Polovodič InSb má zakázaný pás o šířce Eg = 0.23eV, statickou permitivitu e = 18 a efektivní hmotnost elektronů mef = 0.15 me. Vypočtěte ionizační energii donoru, poloměr dráhy odpovídající základnímu stavu a minimální koncentraci donoru, při níž se začíná projevovat překrývání elektronových drah sousedních příměsových atomů (vzniká příměsový pás). 7.2 Statistika nositelů náboje v polovodiči typu N V polovodiči je 1013 donoru v cm3, které mají ionizační energii E d = 1 meV a efektivní hmotnost mef = 0.01 me. Žádné akceptorové atomy nejsou přítomny a polovodič je nedegenerovaný, tj. Eg ^> kßT. Odhadněte koncentraci vodivostních elektronů při T = 4 K a hodnotu Hallovy konstanty. 7.3 Hallův jev pro dva typy nositelů Předpokládejte, že koncentrace vodivostních elektronů a děr v polovodiči jsou nap, relaxační doby re a Th a efektivní hmotnosti me a m^. Ukažte, že Hallův koeficient je _ 1 p - nb2 e (p + nb)2 kde b = ue/fih Je poměr pohyblivostí. Při výpočtu zanedbejte členy s B2. 7.4 Intrinsický polovodič Germanium má nepřímý zakázaný pás o šířce 0.67 eV. Ve vodivostním pásu je osm L minim ve tvaru rotačních elipsoidů s efektivními hmotnostmi niT = 1.6 me a ttíl — 0.08 me. Maximum valenčního pásu se nachází v bodě T a vybíhají z něj dvakrát degenerovaný pás těžkých děr s izotropní efektivní hmotností 0.28 me a dvakrát degenerovaný pás lehkých děr s izotropní efektivní hmotností 0.044 me. Vypočtěte intrinsickou koncentraci nositelů náboje při teplotě 300 K. a k b k Pásová struktura germania podle článku Wachs, A. L., Miller, T., Hsieh, T. C, Shapiro, A. P., Chiang. T. C: Phys. Rev. B 32 (1985) 2326 9 8 Kmity mřížky v harmonickém přiblížení 8.1 Kmity dvouatomového řetězce Uvažte normální módy lineárního řetězce, ve kterém jsou hmotnosti atomů střídavě M\ a M2. Silové konstanty interakce mezi nejbližšími atomy jsou rovny / a tyto atomy jsou vzdáleny a/2. Najděte disperzní relace. Ukažte, že pro M\ = M2 se výsledek redukuje na disperzní relaci jednoatomového řetězce. Pozn.: Zajímavý je graf uj/uq, kde Uq = f/fi, pro měnící se a = M\/M2 při konstantní efektivní hmotnosti (i = M1M2/(M1 + M2). 8.2 Kmity lineární mřížky s dalekodosahovou interakcí Předpokládejme jednorozměrný krystal, v němž existuje interakce i mezi dalekými sousedy. Harmonický člen v potenciální energii nechť je tvaru n m>0 Najděte disperzní relaci takového krystalu a její dlouhovlnnou limitu. 8.3 Konstantní rychlost zvuku Jak musí být voleny konstanty Km v předchozí úloze, aby disperzní relace byla čistě lineární, uj = c\k\? 8.4 Měkký fononový mód Uvažte lineární řetězec složený z iontů stejné hmotnosti, ale střídajícího se náboje ±e. Meziatomový potenciál se skládá z krátkodosahové interakce se silovou konstantou C a z elektrostatické (daleko-dosahové) interakce. Ukažte, že elektrostatickou interakci lze popsat silovou konstantou mezi n-tými nejbližšími sousedy Cn_ (~l)we2 27r£o(na)3 kde a je vzdálenost nejbližších sousedů. Najděte disperzní relaci a nakreslete její graf pro vhodně volené parametry. 8.5 Rychlost zvuku v křemíku S použitím následujícího obrázku určete rychlost zvuku v křemíku ve směrech [100], [110] a [111]. r L Wavevector q DOS Fononové disperzní křivky a hustota stavů podle článku Giannozzi, P., de Gironcoli, S., Pavone, P., Baroni, S.: Phys. Rev. B 43 (1991) 7231 10 8.6 Tepelná kapacita jednoduché ID a 2D mřížky Uvažujme o jednoduché jednorozměrné resp. dvourozměrné mřížce s jedním atomem v primitivní buňce, pro jejíž transverzální kmity platí pohybové rovnice müi = K(ui-i-2ui + ui+i) (ID) müij = K(uí-ij + ui+ij + Uíj-i + Uiíj+i - iuíj) (2D) Najděte disperzní relace kmitů mřížky a teplotní závislost jejího specifického tepla při velmi nízkých teplotách. 8.7 Hustota stavů akustické fononové větve Nechť je disperzní relace některé akustické fononové větve dána vztahem u){q) = u)0a \k 2 Numerickým výpočtem zjistěte hustotu stavů od této akustické větve pro dimenzi mřížky D = 1,2,3. Tuto hustotu stavů srovnejte s Debyeovým modelem. 11