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势所需要的参数，这些参数均可在相关参考文献中找到：