lammps模拟中tersoff势设置方法

大家好，我是小马老师。
本文继续lammps基础知识介绍，主要讲解tersoff势的设置方法。
tersoff势是一种非键接（non-bond）势，在SiC、GaAs、GaN等体系中用的较多。
tersoff势参数保存在一个文本文件中，通常以".tersoff"为后缀名，因此，在lammps中不需要设置tersoff势的具体参数，仅需指定对应的原子类型即可。
下面介绍不同情况下，tersoff势设置方法。
1.体系只包含tersoff势
假设体系只包含Si、C两种原子，对应的原子类型分别为：Si（type 1）、Si（type 2）、C（type 3）。
tersoff势写法：

pair_style tersoff
pair_coeff * * SiC.tersoff Si Si C

第一行pair_style指定势函数类型为tersoff。
第二行pair_coeff映射原子类型，pair_coeff命令后面必须为“”，不能写具体的原子类型（如1 1）。
“”后面为tersoff势文件名称，最后一部分为原子列表。
“Si Si C”表示前两种原子类型为Si，第3种原子为C，lammps在积分运算时会自动根据这个映射关系到“SiC.tersoff”文件找出原子间的参数。
这部分顺序必须与data文件或者体系模型中的原子类型相对应，否则会计算出错。
如体系中C原子类型为2，Si原子类型分别为1和3，则代码要改为：

pair_coeff * * SiC.tersoff Si C Si 

2.混合势写法
如果体系中还包含使用其它势的原子，则需要使用混合势写法。
以SiC和Cu体系为例，体系中的原子类型分别为：Si（type 1）、C（type 2）、Cu（type 3），对应的写法为：

pair_style  hybrid tersoff eam lj/cut 10.0
pair_coeff  * * tersoff SiC.tersoff Si C NULL
pair_coeff  3 3 Cu.eam
pair_coeff  1 3 lj/cut 0.0034 3.12 
pair_coeff  2 3 lj/cut 0.053  2.98

pair_style命令指出需要设置的三种势函数类型。
第2行表示前2种原子（Si、C）使用tersoff势，第三种原子位置为NULL表示该原子（Cu）不使用SiC势。
第3行设置Cu为eam势。
第4、5行设置Si-Cu、C-Cu为lj/cut势，参数仅为演示。
tersoff势设置方法就介绍到这里，如有问题或者需要其它势的介绍，可在文后留言。

公众号案例代码以及学习交流，请加QQ群：754749935

扫描关注微信公众号：lammps加油站。