大家好,我是小马老师。
本文介绍OPLS势函数的设置方法。
lammps力场大体分为两种情况:
一种类似于eam这种势,参数是写到eam力场文件中的,在使用的时候需要下载力场文件,在设置时不需设置力场参数,只需映射原子类型即可。
另一种力场则不需要力场文件,如lj/cut。这种力场大多使用公式描述原子间相互作用力,这些公式内嵌在lammps代码中,在设置力场时,不需要下载力场文件,只需要设定公式中的参数即可。
本文所讲的OPLS属于第二种力场,不需要力场文件,但是需要设置力场参数。
OPLS对应的公式为:
从公式可以看出,OPLS包含了键(bond)、角(angle)、非正常二面角(dihedral或torsion)以及非键接势(non-bonded)。
在lammps中设置OPLS时,这些势都需要单独设置,这就需要找到与之对应的力场类型和参数。
(1)bond
bond势就是谐振势harmonic&