合金體系的勢函數除了eam勢，還有meam勢。在新版本的lammps中，meam勢類型已經改為meam/c，本文主要介紹meam/c勢的設置方法。
和普通的勢函數文件不同，meam/c勢有兩個勢函數文件：library.meam和**.meam，**表示不同的勢函數名稱。

library.meam通用文件
library.meam類似于參數庫，存儲了多種元素的meam參數，每一種元素占3行，共19個通用參數。
例如，下面一段代碼是ZrCu合金的library.meam文件：

#elt lat z iel atwt #alpha b0 b1 b2 b3 alat esub asub #t0 t1 t2 t3 rhozero ibar 'Zr' 'hcp' 12 1 91.22004.4501908328 2.450 1.000 3.000 2.000 3.2000000000 6.360 0.6801.00 6.300 -3.300 -10.00 1.000 3 'Cu' 'fcc' 12 1 63.54605.1548300830 3.830 2.200 6.000 2.200 3.6133156519 3.540 0.9401.00 2.720 3.040 1.950 1.000 3 …………

library.meam前三行是注釋部分，說明了各行參數的定義，后面分別是Zr和Cu對應的參數。
2.專用meam文件
第二個meam文件存儲合金元素專用的參數，描述合金原子之間的相互作用。
如ZrCu.meam參數如下：

rc = 5.0 delr = 0.1 augt1 = 0 erose_form = 2 ialloy = 2 zbl(1,1) = 0 nn2(1,1) = 1 rho0(1) = 1.000 Ec(1,1) = 6.360 re(1,1) = 3.2000 alpha(1,1) = 4.45019083 repuls(1,1) = 0.00 attrac(1,1) = 0.00 Cmin(1,1,1) = 1.00 Cmax(1,1,1) = 1.44 ……

可以看出，在ZrCu.meam中并沒有出現Zr和Cu元素名稱，而是用序號1、2表示。
在這里，1和2為索引號，并不是in文件中原子類型編號。
3.meam/c設置方法
把library.meam和ZrCu.meam文件保存到in文件同一文件夾。
假設in文件中只有兩種原子，原子類型1為Zr，2為Cu，勢函數設置為：

pair_style meam/c pair_coeff * * library.meam Zr Cu ZrCu.meam Zr Cu

假設原子類型1為Cu，2為Zr，勢函數設置為：

pair_style meam/c pair_coeff * * library.meam Zr Cu ZrCu.meam Cu Zr

假設有四種原子類型，1和2為Cu，3和4為Zr：

pair_style meam/c pair_coeff * * library.meam Zr Cu ZrCu.meam Cu Cu Zr Zr

小結：通過上面例子可以看出，library.meam文件名后緊跟的參數Zr Cu功能是從library.meam中提取Zr和Cu的參數，并將Zr和Cu的索引設為1和2.
ZrCu.meam文件名的參數是將Zr和Cu分別映射到in文件中的原子類型。

---

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

總結

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

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