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

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

LSEPK：指定是不是把要計(jì)算的partialchargedensity按每個(gè)k分別寫到各自對(duì)應(yīng)的文件PARCHG.nk.(設(shè)置為.TRUE.)中，還是把它們合并寫到一個(gè)文件中(相當(dāng)于把各個(gè)帶對(duì)應(yīng)的partialchargedensity加起來，設(shè)置為.FALSE.)。默認(rèn)值為.FALSE.。第二步：處理數(shù)據(jù)1〕將計(jì)算結(jié)果中的POSCAR和CHGCAR下載到本地，分別用小程序處理。POSCAR---POStoOKes11.exe〔輸入版本、原子種類數(shù)和原子符號(hào)〕---得到MS可以翻開的結(jié)構(gòu)文件.res文件CHGCAR〔刪除原子符號(hào)那一行〕---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)得到一個(gè)平行于BC面的切面：我們可以對(duì)這個(gè)切面進(jìn)行操作：選定這個(gè)切面，變成黃色，點(diǎn)擊下面的圖標(biāo)AlignLeft/Right：左對(duì)齊/右對(duì)齊AlignVertical:垂直對(duì)齊Alignin/out:對(duì)齊軸AlignWithScreen：對(duì)齊屏幕AlignWithViewYZPlane：對(duì)齊YZ平面AlignHorizontal：