大家好，我是小马老师。
本文继续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加油站
