引用回帖:

8楼: Originally posted by valenhou001 at 2021-06-10 08:01:24
我现在对lammps里的这个Albes势格式不确定的地方就是这个n的值。从Albe势的公式里，看不出这个n对应的是哪个物理量。lammps程序里附带的GaN的Albe势里，n的值都是设置1的。我测试过，将n=0与n=1，其他的都不变的话 ...

补充h对应的是cos(\theta_0)，还需要注意h的正负号。Albe势，在很多文献里写的是[h+cos(\theta_{ijk})]^2。Tersoff势在很文献里写成的是[cos(\theta_{ijk}-h]^2或[h-cos(\theta_{ijk}]^2。 lammps主页上对Tersoff势表达式写的是[cos(\theta_{ijk}-h]^2，这里h=cos(\theta_0)。

引用回帖:

10楼: Originally posted by 半寸日光 at 2021-06-10 09:49:28
感谢！
...

再补充一下：写Albe势的格式或Tersoff1势格式时：
以X， Y,  Z三个不同的元素为例，当然了，也可以类推到它们俩俩为为相同的元素。

X  Y  Z 的其中所有与二体势相关的参数(n, beta, lambda2, B, labmda1, A)都与X-Y对的列表值一样。而与三体势相关的参数(m, gamma, lambda3，c, d， costheta0, R, D)的值需要要考虑到与元素C相关的设置，或者说这些三体势相关的参数的值是与X-Z对的参数列表值给出的，而二体势相关的参数有X-Y对的参数列表值给出。
A,B, labmda1, lambda2的值是基于参数r0, D0, S, beta的值按公式计算出来的：
lambda1=beta*sqrt(2.0*S)
lambda2=beta*sqrt(2.0/S)
A= D0/(S-1.0) * exp(lambda1*r0)
B= S*D0/(S-1.0)*exp(lambda2*r0)
，

引用回帖:

11楼: Originally posted by valenhou001 at 2021-06-10 09:55:25
再补充一下：写Albe势的格式或Tersoff1势格式时：
以X， Y,  Z三个不同的元素为例，当然了，也可以类推到它们俩俩为为相同的元素。

X  Y  Z 的其中所有与二体势相关的参数(n, beta, lambda2, B, labmda1, A)都 ...

我是按这个方式来写的X-Y两种不同元素的Albe势的文件的。按公式应该是可以类推到X-Y-Z三种不同元素的。主要原因是，很多文献上，包括Albe势函数的文章中，有的将将three-body项相关的参数中的索引(ik}简写为了{i}。

您好，γ, c, d, and cosθ0 等参数有拟合公式吗？

您好，请问您写出修改过的wc-co的tersoff势函数了吗，可以分享一下吗。（有偿也行）