腳本分享

在通過vaspkit使用能量應(yīng)變方法計(jì)算材料的彈性常數(shù)時(shí)，往往是依據(jù)vaspkit識別輸入文件的晶體結(jié)構(gòu)，判斷對稱性，依據(jù)晶系、空間群等施加應(yīng)變，生成計(jì)算任務(wù)文件。

對于摻雜體系來講，由于會對原胞進(jìn)行擴(kuò)胞并摻雜替換原子，破壞對稱性并改變點(diǎn)群等，vaspkit會將其識別為其他空間群，計(jì)算彈性常數(shù)的任務(wù)生成也會受到影響。

這里主要介紹如何更改輸入文件，使得vaspkit可以接受指定對稱性即指定的空間群進(jìn)行彈性常數(shù)的計(jì)算任務(wù)的生成。

筆者使用一個(gè)空間群為Fm-3m的面心立方結(jié)構(gòu)計(jì)算彈性常數(shù)時(shí)，使用vaspkit生成計(jì)算文件的同時(shí)會將vaspkit識別的晶體對稱信息保存到文件SYMMETRY中，內(nèi)容如下，其中225為空間群Fm-3m的編號。

# Read the symmetry of structure from the SYMMETRY.in file if it exists.?

?225 ? ? ? ? ?#Space group number of the input structure

所需要計(jì)算的彈性模量主要為C11，C12和C44。

而基于原胞擴(kuò)包并摻雜之后結(jié)構(gòu)對稱性變?yōu)榱薈m。所生成的計(jì)算文件夾任務(wù)數(shù)量顯著增多，具體為：

vaspkit同時(shí)也生成了其對稱性信息的文件保存在SYMMETRY中。為了在計(jì)算其彈性常數(shù)的過程中施加和未擴(kuò)胞結(jié)構(gòu)同樣方向的應(yīng)變，可將原空間群為Fm-3m的SYMMETRY文件更改為SYMMETRY.in，并重新運(yùn)行vaspkit以獲得基于指定空間群的彈性常數(shù)的計(jì)算。

愿有所成

IEchoQ

引喻失義? ?妄自菲薄