材料計算學(xué)原理33

1、4v 計算軟件CASTEPv 能帶計算的過程v 晶體的總能量v 幾何優(yōu)化v 能帶結(jié)構(gòu)與能態(tài)密度目錄目錄v 布居分析v 彈性常數(shù)v 熱力學(xué)性質(zhì)v 光學(xué)性質(zhì) 關(guān)于CASTEP CASTEP（Cambridge Sequential Total Energy Package 的縮寫）是特別為固體材料學(xué)而設(shè)計的一個現(xiàn)代的量子力學(xué)基本程序，其使用了密度泛函(DFT)平凡面波贗勢方法，進行第一原理量子力學(xué)計算，以探索如半導(dǎo)體，陶瓷，金屬，礦物和沸石等材料的晶體和表面性質(zhì)。 CASTEP的應(yīng)用包括表面化學(xué)、鍵結(jié)構(gòu)、態(tài)密度和光學(xué)性質(zhì)等研究，CASTEP也可用于研究體系的電荷密度和波函數(shù)的3D形式。此外，CAS

2、TEP可用于有效研究點缺陷(空位，間隙和置換雜質(zhì))和擴展缺陷(如晶界和位錯)的性質(zhì)。適用于固體物理，材料科學(xué)，化學(xué)以及化工領(lǐng)域，可以節(jié)省實驗成本，縮短開發(fā)周期。CASTEP 所采用的一些方法u 超原胞方法 u 自洽電子弛豫方法u 平面波基u 快速傅立葉變換在CASTEP計算中有很多運行步驟，可分為如下幾組： 設(shè)置電子選項：精度設(shè)置、交換-關(guān)聯(lián)函數(shù)的設(shè)置、贗勢的設(shè)置、 截斷能的設(shè)置、K點的設(shè)置（布里淵區(qū)的設(shè)置）等。 結(jié)構(gòu)優(yōu)化任務(wù)的設(shè)置 計算體系性質(zhì)的設(shè)置：在CASTEP 中可以計算體系的性質(zhì)，如能帶結(jié)構(gòu)、態(tài)密度、聚居數(shù)、聲子色散關(guān)系、聲子態(tài)密度、光學(xué)性質(zhì)、應(yīng)力等。 結(jié)果分析CASTEP可以實現(xiàn)如

3、下的功能： 1. 計算體系的總能； 2. 進行結(jié)構(gòu)優(yōu)化； 3. 執(zhí)行動力學(xué)任務(wù)：在設(shè)置的溫度和關(guān)聯(lián)參數(shù)下，研究體系中原子的運動行為； 4. 計算周期體系的彈性常數(shù)； 5. 化學(xué)反應(yīng)的過度態(tài)搜索。 除此之外，計算一些晶體的性質(zhì)，如能帶結(jié)構(gòu)、態(tài)密度、聚居數(shù)分析、聲子色散關(guān)系、聲子態(tài)密度、光學(xué)性質(zhì)、應(yīng)力等。 晶體能帶及晶體物理性質(zhì)的計算過程與分子自洽場計算類似，如圖1.1所示，不同之處在于晶體是一個具有周期性結(jié)構(gòu)的體系，輸入時只能給出一個體積有限的晶體結(jié)構(gòu)模型，需利用周期性邊界條件，才能得到整個晶體的能帶結(jié)構(gòu)。 在計算晶體的物理性質(zhì)之前必須對所設(shè)晶體結(jié)構(gòu)模型進行幾何優(yōu)化，根據(jù)關(guān)于能量、力、應(yīng)力、位移

4、的判據(jù)來判斷晶體結(jié)構(gòu)是否為穩(wěn)定結(jié)構(gòu)（總能量最?。Ｈ绻w結(jié)構(gòu)不是穩(wěn)定結(jié)構(gòu)，重新設(shè)置晶格參數(shù)進行計算，直至得到穩(wěn)定的晶體結(jié)構(gòu)。對結(jié)構(gòu)優(yōu)化后的晶體體進行物理性質(zhì)計算，最后輸出計算結(jié)果。 晶體總能量（不包括核的動能部分）可分成兩部分：一部分是原子核與內(nèi)層電子組成的離子實的能量，這部分能量基本上與晶體結(jié)構(gòu)無關(guān)，是一個常數(shù)，贗勢方法中常把總能量中這部分不變的能量設(shè)為零；另一部分是總能量與離子實能量之差，包括離子實與價電子的相互作用、離子實之間的相互作用以及價電子間的相互作用。 在密度泛函理論中，晶體總能量ET 是晶格電子的能量與離子實的排斥能之和，即ET TEextEcoul ExcENNNNEETE

5、rcx Vrcx rr d| r r | r r r rd d21ii）（）（）（）（）（ 對于組成元素確定的體系，可能存在不同的晶體結(jié)構(gòu)。 分別對不同的晶體結(jié)構(gòu)進行總能量計算和Murnaghan方程擬合，得到相應(yīng)的能量極小值E（ N），比較不同晶體結(jié)構(gòu)的能量極小值，便可確定穩(wěn)定的晶體結(jié)構(gòu)。 幾何優(yōu)化是通過調(diào)節(jié)結(jié)構(gòu)模型的幾何參數(shù)來獲得穩(wěn)定結(jié)構(gòu)的過程，其結(jié)果是使模型結(jié)構(gòu)盡可能地接近真實結(jié)構(gòu)。進行幾何優(yōu)化的判據(jù)可以根據(jù)研究的需要而定，一般是幾個判據(jù)組合使用。 進行幾何優(yōu)化常用的判據(jù)有以下幾個： 自洽場收斂判據(jù)。 對給定的結(jié)構(gòu)模型進行自洽場計算時，相繼兩次自恰計計算得到的晶體總能量之差足夠小，即相繼

6、兩次自洽計算的晶體總能量之差小于設(shè)定的最大值。 力判據(jù)。 每個原子所受的晶體內(nèi)作用力（Hellmann-Feynman）足夠小，即單個原子受力小于設(shè)定的最大值。 應(yīng)力判據(jù)。每個結(jié)構(gòu)模型單元中的應(yīng)力足夠小，即應(yīng)力小于設(shè)定的最大值。 位移判據(jù)。 相繼兩次結(jié)構(gòu)參數(shù)變化引起的原子位移的分量足夠小，即原子位移的分量小于設(shè)定的最大值。CASTEP幾何優(yōu)化幾何優(yōu)化 CASTEP幾何優(yōu)化任務(wù)允許改善結(jié)構(gòu)的幾何，獲得穩(wěn)定結(jié)構(gòu)或多晶型物。通過一個迭代過程來完成這項任務(wù)，迭代過程中調(diào)整原子坐標(biāo)和晶胞參數(shù)使結(jié)構(gòu)的總能量最小化。 CASTEP幾何優(yōu)化是基于減小計算力和應(yīng)力的數(shù)量級，直到小于規(guī)定的收斂誤差。也可能給定外部

7、應(yīng)力張量來對拉應(yīng)力，壓應(yīng)力和切應(yīng)力等作用下的體系行為模型化。在這些情況下反復(fù)迭代內(nèi)部 幾何優(yōu)化處理產(chǎn)生的模型結(jié)構(gòu)與真實結(jié)構(gòu)緊密相似。利用CASTEP計算的晶格參數(shù)精度列于右圖。應(yīng)力張量直到與所施加的外部應(yīng)力相等。 第一個是能量的收斂精度，單位為eV/atom ，是體系中每個原子的能量值；第二個是作用在每個原子上的最大力收斂精度，單位為 eV/ ；第三個是最大應(yīng)變收斂精度，單位為 GPa；第四個是最大位移收斂精度，單位為 。這些收斂精度指的兩次迭代求解之間的差，只有當(dāng)某次計算的值與上一次計算的值相比小于設(shè)置的值時，計算才停止。 （1）能帶結(jié)構(gòu)計算 量子力學(xué)計算表明，固體中若有N個原子，由于各原子

8、間的相互作用，對應(yīng)于原來孤立原子的每一個能級，變成了N條靠得很近的能級，稱為能帶。固體中的電子能級固體中的電子能級有什么特點？有什么特點？ 能帶的寬度記作E，數(shù)量級為EeV。 若N1023，則能帶中兩能級的間距約10-23eV。 能帶結(jié)構(gòu)就是晶體電子的能量E與波矢k之間的關(guān)系曲線。能帶結(jié)構(gòu)分析在各個領(lǐng)域的第一原理計算工作中用得非常普遍。首先可以看出這個體系是金屬、半導(dǎo)體還是絕緣體。對于本征半導(dǎo)體，還可以看出是直接能隙還是間接能隙：如果導(dǎo)帶的最低點和價帶的最高點在同一個k點處，則為直接能隙，否則為間接能隙。 目前的計算大多采用超單胞（supercell）的形式，在一個單胞里有幾十個原子以及上百個

9、電子，所以得到的能帶圖往往在遠(yuǎn)低于費米能級處非常平坦，也非常密集。但是，我們主要關(guān)心的還是費米能級附近的能帶形狀。 能帶的寬窄在能帶的分析中占據(jù)很重要的位置。能帶越寬，也即在能帶圖中的起伏越大，說明處于這個帶中的電子有效質(zhì)量越小、非局域（non-local）的程度越大、組成這條能帶的原子軌道擴展性越強。 一條比較窄的能帶表明對應(yīng)于這條能帶的本征態(tài)主要是由局域于某個格點的原子軌道組成，這條帶上的電子局域性非常強，有效質(zhì)量相對較大。 能態(tài)密度反映了EdE這個能量范圍內(nèi)能級數(shù)的多少，EdE這個能量范圍內(nèi)軌道（能級數(shù)）越多越密集，態(tài)的密度越大。（2）能態(tài)密度的計算能態(tài)密度（DOS）的分類： 總態(tài)密度

10、分波態(tài)密度（PODS）1.局域態(tài)密度（LDOS） 兩個原子的原子軌道組合以后，得到兩個分子軌道，在周期邊界條件下，這兩個分子軌道形成兩個能帶，根據(jù)能帶的寬度和斜率，可以得到態(tài)密度的近似圖。BN 的總的能態(tài)密度 對電子電荷在各組分原子之間的分布情況進行計算，稱為布居分析。有多種布居分析方法，其中被廣泛采用的布居分析方法是Mulliken 布居分析。布居分析可以給出原子上、原子軌道上、兩原子間的電子電荷分布，依次稱為原子布居、軌道布居、鍵布居。 布居分析為原子間的成鍵提供了一個客觀判據(jù)，鍵布居的值高表明鍵是共價的，鍵布居的值低表示鍵是一種離子相互作用。 有效離子價來也可以用來評價鍵的離子性，若這個

11、值為零，則表明該鍵是完全的離子鍵，若這個值大于零，則表明該鍵的共價成分增加。 材料的彈性常數(shù)描述了它對所加應(yīng)力的響應(yīng)，或者說，彈性常數(shù)描述了為維持一個給定的形變所需的應(yīng)力。應(yīng)力和應(yīng)變均為二階對稱張量，可分別用i 和j （i、 j ，2， ，）來表示，則線彈性常數(shù)可表示為一個的對稱矩陣Cij。 由晶體總能量ET ，可算出彈性常數(shù)Cij。利用計算得到的彈性常數(shù)Cij，可以計算體彈性模量、泊松系數(shù)等性質(zhì)。在計算彈性常數(shù)時，對能量的計算精度要求很高，因此，k的取值不少于15 5 15。 使用CASTEP來計算彈性常數(shù)和其他的力學(xué)性能。首先我們要優(yōu)化BN立方晶體的結(jié)構(gòu)，然后計算它的彈性常數(shù)。 1 優(yōu)化B

12、N立方晶體的結(jié)構(gòu) 2 計算BN的彈性常數(shù)（1）優(yōu)化BN立方晶體的結(jié)構(gòu)首先應(yīng)導(dǎo)入BN結(jié)構(gòu)。 在菜單欄中選擇File | Import，從structures/semiconductors中選中BN.xsd。BN的晶體結(jié)構(gòu)如右：現(xiàn)在設(shè)置幾何優(yōu)化。 從工具欄中選擇CASTEP 工具，然后選擇Calculation或從菜單欄中選擇Modules | CASTEP | Calculation。 CASTEP Calculation對話框如下： 在Setup標(biāo)簽中，把Task設(shè)置為Geometry Optimization，把Quality 設(shè)置為Fine，并且把Functional設(shè)置為GGA and

13、PW91。 選擇Electronic標(biāo)簽，按下More.按鈕以得到CASTEP Electronic Options對話框。把Derived grid的設(shè)置從Standard改為Fine。關(guān)閉CASTEP Electronic Options對話框。 選擇Job Control標(biāo)簽，選中你想要運行CASTEP工作的Gateway。 按下CASTEP Calculation對話框中的Run按鈕。優(yōu)化之后，此結(jié)構(gòu)的晶胞參數(shù)應(yīng)為a=b=c=3.60556。現(xiàn)在我們可以繼續(xù)計算優(yōu)化結(jié)構(gòu)的彈性常數(shù)。 （2）計算BN的彈性常數(shù)選擇CASTEP Calculation對話框中的Setup標(biāo)簽。從Task的下

14、拉清單中選擇Elastic Constants。按下More.按鈕。 把Number of steps for each strain從4增加到6。關(guān)閉CASTEP Elastic Constants對話框。按下CASTEP Calculation對話框中的Run按鈕。 CASTEP的彈性常數(shù)計算任務(wù)的結(jié)果以一批.castep輸出文件的形式給出。這些文件中的每一個文件都代表確定的晶胞在假設(shè)的應(yīng)變模式和應(yīng)變振幅下的幾何優(yōu)化運行結(jié)果。對這些文件的習(xí)慣命名方式為：seedname_cij_m_n。對于給定的模式來說，m代表當(dāng)前的應(yīng)變模式，n代表當(dāng)前的應(yīng)變振幅。 CASTEP可以使用這些結(jié)果來分析每一

15、個運行計算出來的壓力張量和產(chǎn)生一個有關(guān)彈性屬性信息的文件。 CASTEP Elastic Constants對話框如右： 在結(jié)果文件夾中創(chuàng)建了一個新的文檔BN Elastic Constants.txt。 此文檔中的信息包括輸入的應(yīng)變摘要和計算出的應(yīng)力；每一種應(yīng)變模式的線性適配結(jié)果包括適配的質(zhì)量；在計算出的應(yīng)力和給定對稱性的彈性常數(shù)之間的對應(yīng)；彈性常數(shù)Cij和彈性柔量Sij的表格以及最后推知的屬性如體積模量和于之相反的可壓縮性楊氏模量以及三唯方向上的Poisson比率和建造一個物質(zhì)為各向同性媒介的模型所需要的Lame 常數(shù)。 從工具欄中選擇CASTEP 工具，然后選擇Analysis或者從菜單欄中選擇Modules | CASTEP | Analysis。 從屬性清單中選擇Elastic constants，從BN的彈性常數(shù)計算工作中得到的結(jié)果文件應(yīng)自動顯示在Results file選框中。按下Calculate按鈕。 8、熱力學(xué)性質(zhì) 對體系熱力學(xué)性質(zhì)的描述基于聲子，聲子是晶格振動的能量子。聲子的角頻率與波矢q的函數(shù)關(guān)系（q）稱為聲子譜或色散關(guān)系。 利用第一性原理計算聲子譜（q）的方法有兩種：超胞法和線性響應(yīng)法。 由聲子譜（