lammps教程：EAM势参数设置详解

合金体系的势函数除了eam势，还有meam势。在新版本的lammps中，meam势类型已经改为meam/c，本文主要介绍meam/c势的设置方法。
和普通的势函数文件不同，meam/c势有两个势函数文件：library.meam和**.meam，**表示不同的势函数名称。

1. 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加油站