LAMMPS学习系列（5）

定义原子间相互作用势（一）

无论是力场还是原子间相互作用势，其作用都是告诉LAMMPS如何去计算原子间的相互作用力，而且都是人为拟合得到，LAMMPS所支持的原子相互作用势以及力场类型有很多：L-JEAMMEAMREAXFFAIREBO … 虽然LAMMPS能支持的原子相互作用势以及力场类型有很多，但是经验势毕竟需要靠人为拟合，所以还有很多体系（由于找不到合适的力场和相互作用势）是无法用LAMMPS去进行模拟的，这也是为什么DFT用的人更多。 至于如何拟合相互作用势，推荐一篇良心好文章：

此经验势是由湖南大学物理学院陈阳春博士花费数年时间拟合得出~ 所以一般想看自己的体系是否能用LAMMPS模拟时，一定要先多查文献，看看有没有合适的力场或者相互作用势，如果找不到，那么最好暂时放下这个想法，或者转去想其它方法（例如DFT），当然了，如果实力允许的话，可以考虑自行拟合一个，顺便也能造福他人~ 抛开力场不谈（因为本人基本不用力场），LAMMPS的in文件里面定义原子间相互作用势涉及到两个命令：pair_stylepair_coeffpair_style命令表示告诉LAMMPS相互作用势的类型pair_coeff命令给出势函数中的参数或者数值列表例子1. L-J势pair_style lj/cut 2.5pair_coeff * * 1 1 其中* * 表示考虑任意的两个原子间的相互作用，2.5, 1和1为L-J势所需要的参数，这些参数均可在相关参考文献中找到：

2. EAM势pair_style eam/fspair_coeff * * ~/Ti.eam.fs Ti其中，~/表示路径，Ti表示元素类型，W_Re.eam.fs为势文件，里面是一堆数值列表：具体含义有兴趣的同学可以详看手册。一般来说势文件的后缀便是它的类型，例如W_He.eam.alloy,则其势类型便是eam/alloy:pair_style eam/alloypair_coeff * * ~/W_He.eam.alloy W He注：对于含多元的势文件，如果只用到其中一部分元素，则其它的元素就不需要写出来，例如只用到W，则He就不必写，而且W，He的书写顺序要与定义的原子类型序号对应（W为类型1，He为类型2）。至于如何得到这个文件，三个方法：
   直接发邮件找作者要（推荐），一般作者都愿意给，毕竟多一个人用，多一个引用~LAMMPS软件包中有个potentials文件夹，里面有一些常用的势文件和力场文件由于不同的经验势文章作者的表达不尽相同，所以暂时还没能有一个统一的程序去帮助生成这种列表，只有自己根据手册说明以及文章内公式表达，写程序生成，这种方法可参考LAMMPS内置工具：
   _{/lammps/tools/eam_database}/lammps/tools/eam_generate