材料科學(xué)中的第一原理計(jì)算方法介紹課件

材料科學(xué)中的第一原理計(jì)算方法介紹張平北京應(yīng)用物理與計(jì)算數(shù)學(xué)研究所

Emailaddress:zhang_ping@2006年8月14日目標(biāo)初步掌握計(jì)算半導(dǎo)體（和金屬，絕緣體）及其表面電子結(jié)構(gòu)性質(zhì)的第一原理方法和工具內(nèi)容

晶體結(jié)構(gòu)，固體能帶理論第一原理電子結(jié)構(gòu)計(jì)算程序：VASP

利用VASP程序計(jì)算半導(dǎo)體，金屬及絕緣體的物理性質(zhì)材料表面性質(zhì)的第一原理計(jì)算介紹晶體結(jié)構(gòu)元素周期表列(族):具有類(lèi)似的價(jià)電子結(jié)構(gòu)

行(周期性):從左到右，順次加入一個(gè)電子7?值越大，電負(fù)性越強(qiáng)，即越容易得到電子.AdaptedfromFig.2.7,Callister6e.(Fig.2.7isadaptedfromLinusPauling,TheNatureoftheChemicalBond,3rdedition,Copyright1939and1940,3rdedition.Copyright1960byCornellUniversity.電負(fù)性Mulliken:MM++e-Im,IonizationPotentiale-+MM-Am,ElectronAffinityElectronegativity=(Im+Am)/2He-Ne-Ar

-Kr

-Xe

-Rn

-F

4.0Cl

3.0Br2.8I

2.5At

2.2Li

1.0Na

0.9K

0.8Rb

0.8Cs

0.7Fr

0.7H2.1Be

1.5Mg

1.2Ca

1.0Sr

1.0Ba

0.9Ra

0.9Ti

1.5Cr

1.6Fe

1.8Ni

1.8Zn

1.8As

2.0increasingincreasing(三斜晶系)(單斜晶系)(正交晶系)(四方晶系)(立方晶系)(三角晶系)(六角晶系)平移對(duì)稱(chēng)性7大晶系14種Bravais格子原胞，原胞基矢FccBravais格子的基矢對(duì)上圖所示的晶體結(jié)構(gòu)來(lái)講，由這三個(gè)基矢構(gòu)成的原胞含有一個(gè)原子，其位置在基矢空間中表示為：AAFe,Pb,La,…閃鋅礦結(jié)構(gòu)(GaAs,…)AB對(duì)閃鋅礦結(jié)構(gòu)來(lái)說(shuō)，由這三個(gè)基矢構(gòu)成的原胞含有兩個(gè)原子，其位置在基矢空間中表示為：鹽巖(rocksalt)結(jié)構(gòu)對(duì)rocksalt結(jié)構(gòu)來(lái)說(shuō)，由這三個(gè)基矢構(gòu)成的原胞含有兩個(gè)原子，其位置在基矢空間中表示為：ABhcpBravais格子的基矢Be,Mg,Zn,…對(duì)上圖所示的晶體結(jié)構(gòu)來(lái)講，由這三個(gè)基矢構(gòu)成的原胞含有2個(gè)原子，其位置在基矢空間中表示為：AA鉛鋅礦(wurtzite)結(jié)構(gòu)(GaN,ZnO,…)ABAB對(duì)鉛鋅礦結(jié)構(gòu)來(lái)說(shuō)，由這三個(gè)基矢構(gòu)成的原胞含有四個(gè)原子，其位置在基矢空間中表示為：理想值：u=3/8更復(fù)雜的晶體結(jié)構(gòu)

螢石(CaF2)結(jié)構(gòu)黃鐵礦(FeS2)結(jié)構(gòu)晶面，密勒指數(shù)[hkl]代表晶向(hkl)代表晶面金剛石或閃鋅礦結(jié)構(gòu)中的重要晶面

Bloch定理

能帶：價(jià)帶和導(dǎo)帶能隙費(fèi)米能級(jí)有效質(zhì)量空穴金屬，半導(dǎo)體，絕緣體布洛赫定理近自由電子近似可以定性理解簡(jiǎn)單金屬系統(tǒng)中的電子行為定性理解費(fèi)米面和布里淵區(qū)高對(duì)稱(chēng)點(diǎn)的能譜性質(zhì)定性理解禁帶的物理起源Bloch函數(shù)約化布里淵區(qū)第一布里淵區(qū)每條能帶包含的量子態(tài)數(shù)(k點(diǎn)個(gè)數(shù))等于晶體中的原胞個(gè)數(shù)費(fèi)米面2-D立方格子自由電子設(shè)每個(gè)格點(diǎn)(原子)上平均占據(jù)2個(gè)電子，費(fèi)米球半徑為多少？平移對(duì)稱(chēng)性：Bloch定理：本征態(tài)可寫(xiě)為并且滿(mǎn)足周期性邊界條件下：材料科學(xué)中的第一原理計(jì)算

材料模擬的不同計(jì)算層次

Born-Oppenheimer近似

1)離子自由度和電子自由度的剝離2)Hellmann-Feynman定理第一原理電子結(jié)構(gòu)計(jì)算方法

1)Hartree-Fock方法2)密度泛函理論(DFT)3)

局域密度近似(LDA),廣義梯度近似(GGA)密度泛函(DFT)方法忽略動(dòng)能項(xiàng):Thomas-Fermi近似電子密度寫(xiě)成單電子軌道的函數(shù):Kohn-Sham理論求解單電子軌道的基函數(shù)選取

平面波vs局域軌道贗勢(shì)vs全電子處理求解Kohn-Sham方程總能極小化:Born-Oppenheimer方案vsCar-Parrinello

方案疊代對(duì)角化材料模擬的不同計(jì)算層次從頭計(jì)算技術(shù)

Hartree-Fock

量子化學(xué)中常用密度泛函技術(shù)計(jì)算材料科學(xué)緊束縛計(jì)算技術(shù)力場(chǎng)模擬

分子動(dòng)力學(xué)蒙卡

Born-Oppenheimer近似Hellmann-Feynman定理密度泛函理論Hohenberg-Kohn-Sham定理局域密度近似Kohn-Sham理論基函數(shù)選取交換關(guān)聯(lián)勢(shì)密度泛函理論的局限性密度泛函理論中電子-離子相互作用勢(shì)的處理贗勢(shì)近似

贗勢(shì)的導(dǎo)出不是唯一的，原始的贗勢(shì)方法建立在正交化平面波基礎(chǔ)上，對(duì)一個(gè)由許多原子組成的固體，坐標(biāo)空間根據(jù)波函數(shù)的不同特點(diǎn)可分成兩部分：（1）近原子核區(qū)域，所謂的芯區(qū)。波函數(shù)由緊束縛的芯電子波函數(shù)組成，與近鄰原子的波函數(shù)相互作用很小；(2)其余區(qū)域，價(jià)電子波函數(shù)交疊并相互作用。盡管芯區(qū)的勢(shì)很強(qiáng)地吸引價(jià)電子，但是正交化平面波方法中對(duì)價(jià)態(tài)與芯態(tài)正交的要求而產(chǎn)生的大動(dòng)能，對(duì)價(jià)態(tài)的貢獻(xiàn)如同一個(gè)有效的排斥勢(shì)。兩者之和是價(jià)態(tài)的有效勢(shì)，與核的庫(kù)侖勢(shì)相比，這種有效勢(shì)較弱。

關(guān)