合金体系的势函数除了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.2200
4.4501908328 2.450 1.000 3.000 2.000 3.2000000000 6.360 0.680
1.00 6.300 -3.300 -10.00 1.000 3
'Cu' 'fcc' 12 1 63.5460
5.1548300830 3.830 2.200 6.000 2.200 3.6133156519 3.540 0.940
1.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加油站
