| 5 | 1/1 | 返回列表 |
| 查看: 3831 | 回復(fù): 14 | |||
| 【懸賞金幣】回答本帖問題,作者半寸日光將贈(zèng)送您 20 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
半寸日光新蟲 (小有名氣)
|
[求助]
求助:根據(jù)文獻(xiàn)參數(shù)寫適合lammps中tersoff勢的Albe格式問題 已有1人參與
|
||
|
求助:請問有人根據(jù)文獻(xiàn)里面的參數(shù)寫過適合lammps中的tersoff勢?可以幫忙看下我根據(jù)文獻(xiàn)參數(shù)所寫關(guān)于三種元素的tersoff中Albe格式正確嗎?(參考lammps自帶的SiC的Albe),文獻(xiàn)里面的參數(shù),自己所仿照寫的W-C-Co體系的勢函數(shù)如下,感謝! @smutao @oxox6085 發(fā)自小木蟲Android客戶端 |
至尊木蟲 (職業(yè)作家)
|
我現(xiàn)在對lammps里的這個(gè)Albes勢格式不確定的地方就是這個(gè)n的值。從Albe勢的公式里,看不出這個(gè)n對應(yīng)的是哪個(gè)物理量。lammps程序里附帶的GaN的Albe勢里,n的值都是設(shè)置1的。我測試過,將n=0與n=1,其他的都不變的話,兩個(gè)勢給的結(jié)果基本一樣。這個(gè)我在我的github做過注釋。 從GaN勢的出處文獻(xiàn)來看,找不到n在公式里對應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻(xiàn)里(針對的BNC體系的Albes或Tersoff勢)的公式來看,n是一個(gè)指數(shù)因子,并也出現(xiàn)在一個(gè)倒數(shù)的分母里,所以n=0似乎是不合理的。這也是讓我困惑的地方。 你參考一下BNC體系的Tersoff勢: https://www.ctcms.nist.gov/poten ... B-N-C/1/BNC.tersoff 這個(gè)勢文件里,n并不等于0或1. 從這些文獻(xiàn)和相應(yīng)的勢文件來看,這個(gè)n的取值應(yīng)根據(jù)文獻(xiàn)里所給的公式,看它們是怎么定的b_{ij}或B_{ij}的表達(dá)式(也就是repulsive potential前面的一個(gè)系數(shù))。 原本的Tersoff勢表達(dá)式中b_{ij}表達(dá)式比較復(fù)雜的。而Albe勢對這個(gè)b_{ij}簡化了,簡化為(1+\chi_{ij})^{-1/2},這樣n就取為1了。 你可以參考B-N-C體系的Tersoff勢文件(我提到的那篇PRB文獻(xiàn))來試試寫W-C-Co三元體系的Albe勢文件。 與三體相關(guān)的項(xiàng)只是b_ij、\chi_{ij}、g(\theta_{ijk})、cos(theta_{ijk})的項(xiàng)。那篇W-C-Co的勢參考文件中勢參數(shù)列表中h對應(yīng)的就是cos(\theta_{ijk})的值。2\mu對應(yīng)著lammps的Tersoff勢參數(shù)格式中的\lambda_{3}。 |
至尊木蟲 (職業(yè)作家)
|
我把這個(gè)W-C二元的勢文件根據(jù)這篇文獻(xiàn)整理了的。見https://github.com/houzf/ABOP_lammps 發(fā)自小木蟲Android客戶端 |
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑,不管什么專業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +7 | cgyzqwn 2026-03-01 | 14/700 |
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 7/350 |
|
|
[考研] 085600材料與化工,一志愿廣州985,求調(diào)劑 +10 | qqyyaill 2026-03-05 | 10/500 |
|
|
[考研] 085602高分子方向求調(diào)劑 +7 | tlgudy 2026-03-04 | 7/350 |
|
|
[考研] 085600,一志愿鄭州大學(xué),280分求調(diào)劑 +7 | Wuqi725 2026-03-05 | 7/350 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 334求調(diào)劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 泡了個(gè)椒 2026-03-04 | 4/200 |
|
|
[考研] 材料類考研調(diào)劑 +7 | gemmgemm 2026-03-01 | 8/400 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
[考研] 環(huán)境工程專碩307 一志愿211 四六級已過 求調(diào)劑 +5 | ccc! 2026-03-03 | 6/300 |
|
|
[考研] 320材料與化工,求調(diào)劑 +6 | 鶴遨予卿 2026-03-04 | 8/400 |
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研]
材料270求調(diào)劑
6+6
|
Eiiiio 2026-03-01 | 11/550 |
|