上期介绍了LJ势函数参数写法,本期主要介绍金属体系势函数的设置。
金属原子之间没有键连接,因此,在lammps模拟中,金属体系的势函数类型为pair_style,而不是bond_style。
模拟金属体系时,可以用LJ势描述金属原子之间的受力,不过更精确的是嵌入原子势(EAM),eam势函数公式为:
eam势由两部分组成,在原子对势(pair)的基础上添加了电子云密度相关项,比单纯的对势精确度更高。
eam势函数写在一个以“.eam”为后缀的文件中,lammps自带的势函数包含一部分eam势文件,也可以到网站下载eam文件,下载eam文件后,保存到in文件所在的文件夹。
常用的势函数下载网站有:
Interatomic Potentials Repository