使用lammps模拟金属体系,最常用的势函数为eam、meam势。
对于大多数金属或者合金来说,可以从网上下载相对应的eam或meam势力场文件。
而对于一些新的合金材料,可能就找不到相对应的力场文件。
解决的方法有两种,一是采用混合势(pair_style hybrid),第二种方法是自己拟合势函数文件。
混合势的写法在前面文章中已经介绍过,本文介绍一种拟合EAM势的方法。
该方法是运行Xiaowang Zhou编写的一个拟合程序,该程序可自动生成所需要的合金势。
这个程序能够拟合的合金元素有:
Cu, Ag, Au, Ni, Pd