上期介紹了LJ勢函數參數寫法，本期主要介紹金屬體系勢函數的設置。

金屬原子之間沒有鍵連接，因此，在lammps模擬中，金屬體系的勢函數類型為pair_style，而不是bond_style。

模擬金屬體系時，可以用LJ勢描述金屬原子之間的受力，不過更精確的是嵌入原子勢（EAM），eam勢函數公式為：

eam勢由兩部分組成，在原子對勢（pair）的基礎上添加了電子云密度相關項，比單純的對勢精確度更高。

eam勢函數寫在一個以“.eam”為后綴的文件中，lammps自帶的勢函數包含一部分eam勢文件，也可以到網站下載eam文件，下載eam文件后，保存到in文件所在的文件夾。

常用的勢函數下載網站有：

Interatomic Potentials Repositoryhttp://cst-www.nrl.navy.mil/ccm6/aphttp://enpub.fulton.asu.edu/cms/potentials/main/main.htmhttps://openkim.org

1.單原子體系eam寫法

模擬體系只包含一種金屬原子，勢函數的設置比較簡單，如：

pair_style eam pair_coeff * * Fe.eam

2.合金體系eam寫法

對于合金體系，對應的eam勢函數為eam/alloy或者eam/fs，寫法稍有不同。

pair_style 指明eam合金勢類型，pair_coeff映射原子類型。

單原子eam勢文件名后不需要進行原子類型映射，不用寫原子類型列表。

合金體系必須在pair_coeff語句中勢文件名后面把所有的原子類型全部列出，順序和in文件中原子類型要保持一致。

例如，體系中包含Ni、Al兩種原子，Ni的原子類型分別為1、3、4，Al原子類型為2，eam/alloy類型寫法：

pair_style eam/alloy pair_coeff * * NiAlH.eam.alloy Ni Al Ni Ni

eam/fs類型，Ni原子類型為1、2、3，Al原子類型為4：

pair_style eam/fs pair_coeff * * NiAlH_jea.eam.fs Ni Ni Ni Al

更多lammps教程請關注公眾號：lammps加油站

總結

以上是生活随笔為你收集整理的lammps教程：EAM势函数设置详解的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。