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

library.meam通用文件
library.meam類似于參數(shù)庫，存儲了多種元素的meam參數(shù)，每一種元素占3行，共19個(gè)通用參數(shù)。
例如，下面一段代碼是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前三行是注釋部分，說明了各行參數(shù)的定義，后面分別是Zr和Cu對應(yīng)的參數(shù)。
2.專用meam文件
第二個(gè)meam文件存儲合金元素專用的參數(shù)，描述合金原子之間的相互作用。
如ZrCu.meam參數(shù)如下：

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中并沒有出現(xiàn)Zr和Cu元素名稱，而是用序號1、2表示。
在這里，1和2為索引號，并不是in文件中原子類型編號。
3.meam/c設(shè)置方法
把library.meam和ZrCu.meam文件保存到in文件同一文件夾。
假設(shè)in文件中只有兩種原子，原子類型1為Zr，2為Cu，勢函數(shù)設(shè)置為：

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

假設(shè)原子類型1為Cu，2為Zr，勢函數(shù)設(shè)置為：

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

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

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

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

---

更多l(xiāng)ammps教程請關(guān)注微信公眾號：lammps加油站

總結(jié)

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

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