注意使用meam/c在混合势中的小细节

合金体系的势函数除了eam势，还有meam势。在新版本的lammps中，meam势类型已经改为meam/c。和普通的势文件不同，meam/c势有两个函数势文件：library.meam和*.meam，*表示不同的势函数名称。单独在in文件中的具体使用方法请阅读手册，本文详细介绍在混合势中存在多个meam/c势的正确书写格式。

1.介绍meam/c势文件

前已介绍meam/c势有两个函数势文件：library.meam和*.meam，*表示不同的势函数名称。

1.1 library.meam文件内容

下面为NiMn.meam对应的library.meam文件

library.meam前三行是注释部分，说明了各行参数的定义，后面分别是个元素所对应的参数。

1.2 meam文件内容

第二个meam文件存储合金元素专用的参数，描述合金原子之间的相互作用。

如NiMn.meam参数如下（部分）：

1.3 meam/c设置方法

假设in文件中只有两种原子，原子类型1为Mn，2为Ni，势函数设置为：

其中libraray.meam后的元素可和NiMn.meam文件的元素顺序相同，若有的文件不是以元素命名可按其内容中的元素顺序书写。最后的元素顺序按in文件中的顺序书写。

2.meam/c在混合势中——在含有多个meam/c的混合势中

若这样书写会报错ERROR：All pair coeffs are not set。出现这种error，主要是因为同一in文件中使用了两个meam/势，第二个meam/c势会覆盖掉前面的meam/c势，导致Cu的meam/c势参数丢失。

解决方法：

如果在同一in文件中使用多个同种类型的势，为防止覆盖需要对相同类型的势进行编号区分，所以正确写法：

注意：下面有使用几个不同的meam/c势就在pair_style中写几个meam/c。

若只写一个meam/c势就会报错ERROR：Cannot open MEAM library file 见下图所示：