lammps教程：ljcutcoul⼒场参数设置详解
本⽂介绍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.03.0
#设置两个截断半径
pair_style lj/cut/coul/cut 10.08.0
#原⼦1之间使⽤截断半径9.08.0
pair_coeff 11100.03.59.0香港著名男演员
#原⼦2之间使⽤截断半径9.09.0
pair_coeff 22100.03.59.09.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
hatsukoi
pair_coeff **100.03.0
在使⽤lj/cut/coul系列⼒场时，需要根据库伦⼒计算公式选择对应的⼒场，如果是长程库伦作⽤，务必使⽤kspace_style命令以免出现模拟错误。
关注微信公众号：lammps加油站，获得更多lammps教程。
