得到分波電荷密度的詳細(xì)過程

第一步：計(jì)算做一次自洽計(jì)算〔靜態(tài)計(jì)算〕，得到收斂的CHGCAR、CHG、WAVECAR〔注：不計(jì)算分波的話，第一步得到的CHGCAR，就可以經(jīng)過處理得到電荷密度圖，需要自旋極化的情況的話，就在計(jì)算時INCAR中參加ISPIN=2.〕在自洽計(jì)算的根底上，做一次非自洽計(jì)算，讀入文件：第一步得到的：CHGCAR、CHG、WAVECARPOSCAR、POTCAR、KPOINTS、INCAR但是INCAR需要修改：在INCAR中參加參數(shù)：LPARD=.TURE.!表示計(jì)算分波電荷密度，具體分以下三種情況來控制需要計(jì)算的分波情況：IBAND=20212223！表示計(jì)算這四條能帶對應(yīng)的分波電荷密度NBMOD=4！這里代表能帶數(shù)，可以不手動設(shè)置EINT=-0.10.1！表示計(jì)算〔-0.1，0.1〕能量區(qū)間的分波電荷密度圖NBMOD=-2！可以不手動設(shè)置EINT=1！表示計(jì)算1~費(fèi)米能級能量區(qū)間的分波電荷密度圖NBMOD=-3KPUSE=1234！計(jì)算這四個K點(diǎn)對應(yīng)的分波電荷密度以上，1〕和3〕可以同時使用，而2〕一般單獨(dú)使用。此外，還有兩個參數(shù)控制輸出的情況：LSEPB：指定是不是要把計(jì)算的partialchargedensity按每個帶分別寫到各自對應(yīng)的文件PARCHG.nb.(設(shè)置為.TRUE.)中，還是把它們合并寫到一個文件中(相當(dāng)于把各個帶對應(yīng)的partialchargedensity加起來，設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。

LSEPK：指定是不是把要計(jì)算的partialchargedensity按每個k分別寫到各自對應(yīng)的文件PARCHG.nk.(設(shè)置為.TRUE.)中，還是把它們合并寫到一個文件中(相當(dāng)于把各個帶對應(yīng)的partialchargedensity加起來，設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。第二步：處理數(shù)據(jù)1〕將計(jì)算結(jié)果中的POSCAR和CHGCAR下載到本地，分別用小程序處理。POSCAR---POStoOKes11.exe〔輸入版本、原子種類數(shù)和原子符號〕---得到MS可以翻開的結(jié)構(gòu)文件.res文件CHGCAR〔刪除原子符號那一行〕---chgcar2grd111.exe(輸入版本，是否自旋極化)---可以得到.grd文件。2〕將結(jié)構(gòu)文件導(dǎo)入MS導(dǎo)入電荷密度的數(shù)據(jù)，Edit—Insertform—找到文件位置，翻開—OK！看上去沒什么變化，需要點(diǎn)擊鼠標(biāo)右鍵DisplayStyle—field—勾選Visible電荷密度，就在圖中顯示出來，但還需要根據(jù)自己的需要調(diào)整顯示的圖片：1〕右鍵--ColorMaps—Bands可以調(diào)整圖片的細(xì)膩程度—From…to可以調(diào)整區(qū)間的大小，Legend可以在右下角顯示圖例。2〕切片操作，View—Toolbars—Modules顯示切片操作需要的圖標(biāo)得到一個平行于BC面的切面：我們可以對這個切面進(jìn)行操作：選定這個切面，變成黃色，點(diǎn)擊下面的圖標(biāo)AlignLeft/Right：左對齊/右對齊AlignVertical:垂直對齊Alignin/out:對齊軸AlignWithScreen：對齊屏幕AlignWithViewYZPlane：對齊YZ平面AlignHorizontal：