引用回帖:

8樓: Originally posted by valenhou001 at 2021-06-10 08:01:24
我現(xiàn)在對lammps里的這個(gè)Albes勢格式不確定的地方就是這個(gè)n的值。從Albe勢的公式里，看不出這個(gè)n對應(yīng)的是哪個(gè)物理量。lammps程序里附帶的GaN的Albe勢里，n的值都是設(shè)置1的。我測試過，將n=0與n=1，其他的都不變的話 ...

補(bǔ)充h對應(yīng)的是cos(\theta_0)，還需要注意h的正負(fù)號(hào)。Albe勢，在很多文獻(xiàn)里寫的是[h+cos(\theta_{ijk})]^2。Tersoff勢在很文獻(xiàn)里寫成的是[cos(\theta_{ijk}-h]^2或[h-cos(\theta_{ijk}]^2。 lammps主頁上對Tersoff勢表達(dá)式寫的是[cos(\theta_{ijk}-h]^2，這里h=cos(\theta_0)。

引用回帖:

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

再補(bǔ)充一下：寫Albe勢的格式或Tersoff1勢格式時(shí)：
以X， Y,  Z三個(gè)不同的元素為例，當(dāng)然了，也可以類推到它們倆倆為為相同的元素。

X  Y  Z 的其中所有與二體勢相關(guān)的參數(shù)(n, beta, lambda2, B, labmda1, A)都與X-Y對的列表值一樣。而與三體勢相關(guān)的參數(shù)(m, gamma, lambda3，c, d， costheta0, R, D)的值需要要考慮到與元素C相關(guān)的設(shè)置，或者說這些三體勢相關(guān)的參數(shù)的值是與X-Z對的參數(shù)列表值給出的，而二體勢相關(guān)的參數(shù)有X-Y對的參數(shù)列表值給出。
A,B, labmda1, lambda2的值是基于參數(shù)r0, D0, S, beta的值按公式計(jì)算出來的：
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
再補(bǔ)充一下：寫Albe勢的格式或Tersoff1勢格式時(shí)：
以X， Y,  Z三個(gè)不同的元素為例，當(dāng)然了，也可以類推到它們倆倆為為相同的元素。

X  Y  Z 的其中所有與二體勢相關(guān)的參數(shù)(n, beta, lambda2, B, labmda1, A)都 ...

我是按這個(gè)方式來寫的X-Y兩種不同元素的Albe勢的文件的。按公式應(yīng)該是可以類推到X-Y-Z三種不同元素的。主要原因是，很多文獻(xiàn)上，包括Albe勢函數(shù)的文章中，有的將將three-body項(xiàng)相關(guān)的參數(shù)中的索引(ik}簡寫為了{(lán)i}。

您好，γ, c, d, and cosθ0 等參數(shù)有擬合公式嗎？

您好，請問您寫出修改過的wc-co的tersoff勢函數(shù)了嗎，可以分享一下嗎。（有償也行）