注意使用meam/c在混合勢中的小細(xì)節(jié)
合金體系的勢函數(shù)除了eam勢,還有meam勢。在新版本的lammps中,meam勢類型已經(jīng)改為meam/c。和普通的勢文件不同,meam/c勢有兩個函數(shù)勢文件:library.meam和*.meam,*表示不同的勢函數(shù)名稱。單獨在in文件中的具體使用方法請閱讀手冊,本文詳細(xì)介紹在混合勢中存在多個meam/c勢的正確書寫格式。
1.介紹meam/c勢文件
前已介紹meam/c勢有兩個函數(shù)勢文件:library.meam和*.meam,*表示不同的勢函數(shù)名稱。
1.1 library.meam文件內(nèi)容
下面為NiMn.meam對應(yīng)的library.meam文件
library.meam前三行是注釋部分,說明了各行參數(shù)的定義,后面分別是個元素所對應(yīng)的參數(shù)。
1.2 meam文件內(nèi)容
第二個meam文件存儲合金元素專用的參數(shù),描述合金原子之間的相互作用。
如NiMn.meam參數(shù)如下(部分):
1.3 meam/c設(shè)置方法
假設(shè)in文件中只有兩種原子,原子類型1為Mn,2為Ni,勢函數(shù)設(shè)置為:
其中libraray.meam后的元素可和NiMn.meam文件的元素順序相同,若有的文件不是以元素命名可按其內(nèi)容中的元素順序書寫。最后的元素順序按in文件中的順序書寫。
2.meam/c在混合勢中——在含有多個meam/c的混合勢中
若這樣書寫會報錯ERROR:All pair coeffs are not set。出現(xiàn)這種error,主要是因為同一in文件中使用了兩個meam/勢,第二個meam/c勢會覆蓋掉前面的meam/c勢,導(dǎo)致Cu的meam/c勢參數(shù)丟失。
解決方法:
如果在同一in文件中使用多個同種類型的勢,為防止覆蓋需要對相同類型的勢進(jìn)行編號區(qū)分,所以正確寫法:
注意:下面有使用幾個不同的meam/c勢就在pair_style中寫幾個meam/c。
若只寫一個meam/c勢就會報錯ERROR:Cannot open MEAM library file 見下圖所示:
