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

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

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

eam勢(shì)由兩部分組成，在原子對(duì)勢(shì)（pair）的基礎(chǔ)上添加了電子云密度相關(guān)項(xiàng)，比單純的對(duì)勢(shì)精確度更高。

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

常用的勢(shì)函數(shù)下載網(wǎng)站有：

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

1.單原子體系eam寫法

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

pair_style eam pair_coeff * * Fe.eam

2.合金體系eam寫法

對(duì)于合金體系，對(duì)應(yīng)的eam勢(shì)函數(shù)為eam/alloy或者eam/fs，寫法稍有不同。

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

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

合金體系必須在pair_coeff語句中勢(shì)文件名后面把所有的原子類型全部列出，順序和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

更多l(xiāng)ammps教程請(qǐng)關(guān)注公眾號(hào)：lammps加油站

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。