lammps教程：lj/cut/coul力场参数设置详解

本文介绍lj/cut/coul/cut以及lj/cut/coul/long力场的区别和用法。
前面一篇文章介绍了LJ势的用法，普通的LJ势只考虑原子之间的吸引力与排斥力，没有考虑原子之间电荷作用。
如果原子带有电荷，用普通的LJ势算出的结果误差较大。
lj/cut/coul系列力场在普通LJ势的基础上增加库伦作用力，根据所用库伦力公式不同分为多种力场，本文只介绍lj/cut/coul/cut以及lj/cut/coul/long两种力场。
1. lj/cut/coul/cut：短程作用
lj/cut/coul/cut力场使用的公式为：

语法为：

lj/cut/coul/cut args = cutoff (cutoff2)
  cutoff = global cutoff for LJ (and Coulombic if only 1 arg) (distance units)
  cutoff2 = global cutoff for Coulombic (optional) (distance units)

lj/cut/coul/cut的截断半径有两个：普通LJ截断半径和库仑力截断半径。
如果只写一个截断半径，则默认cutoff2=cutoff。

#只设置一个截断半径
pair_style lj/cut/coul/cut 10.0
pair_coeff * * 100.0 3.0

#设置两个截断半径
pair_style lj/cut/coul/cut 10.0 8.0
#原子1之间使用截断半径9.0 8.0
pair_coeff 1 1 100.0 3.5 9.0
#原子2之间使用截断半径9.0 9.0
pair_coeff 2 2 100.0 3.5 9.0 9.0

2. lj/cut/coul/long：长程作用
lj/cut/coul/long力场公式与lj/cut/coul/cut类似，但是在计算库仑力时，不仅计算近邻原子之间的库伦力，也计算近邻原子之外的原子对中心原子的库伦作用力。
因此，需要配合kspace_style pppm或者kspace_style ewald命令使用。
lj/cut/coul/long使用方法与lj/cut/coul/long类似：

kspace_style pppm 1.0e-5
pair_style lj/cut/coul/long 10.0
pair_coeff * * 100.0 3.0

在使用lj/cut/coul系列力场时，需要根据库伦力计算公式选择对应的力场，如果是长程库伦作用，务必使用kspace_style命令以免出现模拟错误。

关注微信公众号：lammps加油站，获得更多lammps教程。