| 5 | 1/1 | 返回列表 |
| 查看: 3889 | 回復: 14 | |||
| 【懸賞金幣】回答本帖問題,作者半寸日光將贈送您 20 個金幣 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
半寸日光新蟲 (小有名氣)
|
[求助]
求助:根據(jù)文獻參數(shù)寫適合lammps中tersoff勢的Albe格式問題 已有1人參與
|
||
|
求助:請問有人根據(jù)文獻里面的參數(shù)寫過適合lammps中的tersoff勢?可以幫忙看下我根據(jù)文獻參數(shù)所寫關(guān)于三種元素的tersoff中Albe格式正確嗎?(參考lammps自帶的SiC的Albe),文獻里面的參數(shù),自己所仿照寫的W-C-Co體系的勢函數(shù)如下,感謝! @smutao @oxox6085 發(fā)自小木蟲Android客戶端 |
至尊木蟲 (職業(yè)作家)
|
我現(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}表達式比較復雜的。而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}。 |
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學芳草園] 伙伴們,祝我生日快樂吧 +16 | myrtle 2026-03-10 | 25/1250 |
|
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +3 | 為學666 2026-03-16 | 3/150 |
|
|
[基金申請] 國自科面上基金字體 +5 | iwuli 2026-03-12 | 6/300 |
|
|
[碩博家園] 深圳大學碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 8/400 |
|
|
[考研] 一志愿985,本科211,0817化學工程與技術(shù)319求調(diào)劑 +4 | Liwangman 2026-03-15 | 4/200 |
|
|
[考研] 化學調(diào)劑0703 +7 | 啊我我的 2026-03-11 | 7/350 |
|
|
[考博] 歡迎申博同學聯(lián)系 +3 | 天道酬勤2026686 2026-03-10 | 7/350 |
|
|
[基金申請] 面上和青基一樣限30頁不合理 +5 | wowsunflower 2026-03-10 | 7/350 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-09 | 4/200 |
|
|
[考研] 求調(diào)劑 +6 | yfihxh 2026-03-09 | 6/300 |
|
|
[考研] 265求調(diào)劑 +9 | 小木蟲085600 2026-03-09 | 12/600 |
|
|
[基金申請]
有必要更換申報口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 0703化學調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|