我把这个W-C二元的势文件根据这篇文献整理了的。见

三元的比较复杂。还没有摸清楚转换关系。

你参考的这篇文献对某个参数的写错了，与原始文献有出了。

引用回帖:

5楼: Originally posted by 半寸日光 at 2021-06-10 02:01:10
我之前在你的主页上参考了这个，照着你的和lammps里面的SiC的Albe一起参考，W-C这个除了某些n，我给了0，没给你写的1，其他基本一样，这个n的0和1对最后的结果有什么影响吗？除了W-C，文献里面给的Co-C，W-Co，Co- ...

我现在对lammps里的这个Albes势格式不确定的地方就是这个n的值。从Albe势的公式里，看不出这个n对应的是哪个物理量。lammps程序里附带的GaN的Albe势里，n的值都是设置1的。我测试过，将n=0与n=1，其他的都不变的话，两个势给的结果基本一样。这个我在我的github做过注释。
从GaN势的出处文献来看，找不到n在公式里对应的变量。 从10.1103/PhysRevB.86.115410文献里（针对的BNC体系的Albes或Tersoff势）的公式来看，n是一个指数因子，并也出现在一个倒数的分母里，所以n=0似乎是不合理的。这也是让我困惑的地方。
你参考一下BNC体系的Tersoff势：
https://www.ctcms.nist.gov/poten ... B-N-C/1/BNC.tersoff
这个势文件里，n并不等于0或1.
从这些文献和相应的势文件来看，这个n的取值应根据文献里所给的公式，看它们是怎么定的b_{ij}或B_{ij}的表达式（也就是repulsive potential前面的一个系数）。
原本的Tersoff势表达式中b_{ij}表达式比较复杂的。而Albe势对这个b_{ij}简化了，简化为(1+\chi_{ij})^{-1/2}，这样n就取为1了。

你可以参考B-N-C体系的Tersoff势文件（我提到的那篇PRB文献）来试试写W-C-Co三元体系的Albe势文件。 与三体相关的项只是b_ij、\chi_{ij}、g(\theta_{ijk})、cos(theta_{ijk})的项。那篇W-C-Co的势参考文件中势参数列表中h对应的就是cos(\theta_{ijk})的值。2\mu对应着lammps的Tersoff势参数格式中的\lambda_{3}
，