大家好,我是小马老师。

本文介绍lammps中设置Buckingham势的方法。

Buckingham势基本公式为:

lammps教程:Buckingham势设置方法_分子动力学

对应lammps中的pair_style buck命令。

此外,根据是否包含库仑力,又细分了多个类型:buck/coul/cut、buck/coul/long、buck/coul/msm,不同类型对应不同的库仑力计算方法。

Buckingham势设置包含5个参数,前三个参数是必须设置的参数,可从论文中查找不同原子对应的参数值,后两个参数是可选参数,主要设置截断半径。

lammps教程:Buckingham势设置方法_分子动力学_02

设置参数时,按照上图顺序设置对应的参数即可,如:

pair_style buck/coul/long 10.0
pair_coeff 1 1 100.0 1.5 200.0

或者:

pair_style buck/coul/cut 10.0 8.0
pair_coeff * * 100.0 1.5 200.0

如sio2的参数可设置为:

pair_style buck/coul/long 10
pair_coeff 1 1 0 0.1 0 #si-si
pair_coeff 2 2 32025.73059 0.3623188 4035.5875 #o-o
pair_coeff 1 2 415158.1815 0.2052124 3079.455355#si-o
kspace_style ewald 1e-4

论文查到的si的参数大多为0,但是在设置时,第二个参数是分母不能为0,只需要设置A、C为0即可。

公众号:lammps加油站