lammps教程：官方自带eam_database程序拟合EAM势方法介绍

使用lammps模拟金属体系，最常用的势函数为eam、meam势。

对于大多数金属或者合金来说，可以从网上下载相对应的eam或meam势力场文件。

而对于一些新的合金材料，可能就找不到相对应的力场文件。

解决的方法有两种，一是采用混合势（pair_style hybrid），第二种方法是自己拟合势函数文件。

混合势的写法在前面文章中已经介绍过，本文介绍一种拟合EAM势的方法。

该方法是运行Xiaowang Zhou编写的一个拟合程序，该程序可自动生成所需要的合金势。

这个程序能够拟合的合金元素有：

Cu, Ag, Au, Ni, Pd, Pt, Al, Pb, Fe, Mo, Ta, W, Mg, Co, Ti, Zr。

不在这个列表之内的元素不能使用这个程序。
该程序是一个开源软件，代码收录在官方Ubuntu版lammps源代码tools/eam_database目录内。
具体拟合方法：
（1）编译creat.f文件
在Ubuntu中进入tools/eam_database目录，运行以下代码进行编译：

gfortran create.f

如果提示没有权限，则运行：

sudo gfortran create.f

运行结束后，会生成a.out程序。
（2）编辑EAM.input文件
EAM.input文件包含了合金信息，按照例子中的格式添加或删除合金元素即可。

 &funccard
 atomtype='Al'
 &end
 &funccard
 atomtype='Cu'
 &end
 &funccard
 &end

（3）运行a.out
在命令行输入以下命令即可生成相应的合金势力场文件：

a.out < EAM.input

文件名为合金元素.set，本例中生成的合金是为AlCu.set。
（4）势函数使用方法

pair_style eam/alloy
pair_coeff * * AlCu.set Al Cu

这个程序使用起来比较方便，在找不到相对应势函数的时候，值得一试。

更多lammps教程，扫描关注微信公众号：