求助：根據(jù)文獻參數(shù)寫適合lammps中tersoff勢的Albe格式問題
求助：請問有人根據(jù)文獻里面的參數(shù)寫過適合lammps中的tersoff勢？可以幫忙看下我根據(jù)文獻參數(shù)所寫關(guān)于三種元素的tersoff中Albe格式正確嗎？(參考lammps自帶的SiC的Albe)，文獻里面的參數(shù)，自己所仿照寫的W-C-Co體系的勢函數(shù)如下，感謝!
返回小木蟲查看更多
今日熱帖
猜你喜歡
板塊導(dǎo)航
- 網(wǎng)絡(luò)生活
- 資源共享
- 化學(xué)化工
- 專業(yè)學(xué)科
- 科研生活
- 生物醫(yī)藥
- 材料
- 計算模擬
- 學(xué)術(shù)交流
- 出國留學(xué)
- 注冊執(zhí)考
24小時熱帖
換一批
- 下載小木蟲APP
- 與700萬科研達人隨時交流
-
二維碼
-
IOS
-
安卓
京公網(wǎng)安備 11010802022153號
我把這個W-C二元的勢文件根據(jù)這篇文獻整理了的。見
三元的比較復(fù)雜。還沒有摸清楚轉(zhuǎn)換關(guān)系。
你參考的這篇文獻對某個參數(shù)的寫錯了，與原始文獻有出了。
我現(xiàn)在對lammps里的這個Albes勢格式不確定的地方就是這個n的值。從Albe勢的公式里，看不出這個n對應(yīng)的是哪個物理量。lammps程序里附帶的GaN的Albe勢里，n的值都是設(shè)置1的。我測試過，將n=0與n=1，其他的都不變的話，兩個勢給的結(jié)果基本一樣。這個我在我的github做過注釋。
從GaN勢的出處文獻來看，找不到n在公式里對應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻里（針對的BNC體系的Albes或Tersoff勢）的公式來看，n是一個指數(shù)因子，并也出現(xiàn)在一個倒數(shù)的分母里，所以n=0似乎是不合理的。這也是讓我困惑的地方。
你參考一下BNC體系的Tersoff勢：
https://www.ctcms.nist.gov/poten ... B-N-C/1/BNC.tersoff
這個勢文件里，n并不等于0或1.
從這些文獻和相應(yīng)的勢文件來看，這個n的取值應(yīng)根據(jù)文獻里所給的公式，看它們是怎么定的b_{ij}或B_{ij}的表達式（也就是repulsive potential前面的一個系數(shù)）。
原本的Tersoff勢表達式中b_{ij}表達式比較復(fù)雜的。而Albe勢對這個b_{ij}簡化了，簡化為(1+\chi_{ij})^{-1/2}，這樣n就取為1了。
你可以參考B-N-C體系的Tersoff勢文件（我提到的那篇PRB文獻）來試試寫W-C-Co三元體系的Albe勢文件。 與三體相關(guān)的項只是b_ij、\chi_{ij}、g(\theta_{ijk})、cos(theta_{ijk})的項。那篇W-C-Co的勢參考文件中勢參數(shù)列表中h對應(yīng)的就是cos(\theta_{ijk})的值。2\mu對應(yīng)著lammps的Tersoff勢參數(shù)格式中的\lambda_{3}
，