方法"/>
lammps教程:Buckingham势设置方法
大家好,我是小马老师。
本文介绍lammps中设置Buckingham势的方法。
Buckingham势基本公式为:
对应lammps中的pair_style buck命令。
此外,根据是否包含库仑力,又细分了多个类型:buck/coul/cut、buck/coul/long、buck/coul/msm,不同类型对应不同的库仑力计算方法。
Buckingham势设置包含5个参数,前三个参数是必须设置的参数,可从论文中查找不同原子对应的参数值,后两个参数是可选参数,主要设置截断半径。
设置参数时,按照上图顺序设置对应的参数即可,如:
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加油站
更多推荐
lammps教程:Buckingham势设置方法
发布评论