 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3836 | 回復(fù): 14 | |||
| 【懸賞金幣】回答本帖問題，作者半寸日光將贈送您 20 個金幣 | |||
半寸日光新蟲 (小有名氣)
|
[求助]
求助：根據(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客戶端 |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑，307分
已經(jīng)有10人回復(fù)
-
0817化學(xué)工程319求調(diào)劑
已經(jīng)有9人回復(fù)
-
304求調(diào)劑（085602一志愿985）
已經(jīng)有6人回復(fù)
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有20人回復(fù)
-
調(diào)劑的同學(xué)，走過路過，不要錯過
已經(jīng)有16人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有8人回復(fù)
-
物理學(xué)求調(diào)劑
已經(jīng)有4人回復(fù)
-
找博導(dǎo)
已經(jīng)有4人回復(fù)
-
一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學(xué)材料，有六級，有工作經(jīng)驗
已經(jīng)有4人回復(fù)
-
288求調(diào)劑
已經(jīng)有15人回復(fù)
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
我現(xiàn)在對lammps里的這個Albes勢格式不確定的地方就是這個n的值。從Albe勢的公式里，看不出這個n對應(yīng)的是哪個物理量。lammps程序里附帶的GaN的Albe勢里，n的值都是設(shè)置1的。我測試過，將n=0與n=1，其他的都不變的話，兩個勢給的結(jié)果基本一樣。這個我在我的github做過注釋。 從GaN勢的出處文獻(xiàn)來看，找不到n在公式里對應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻(xiàn)里（針對的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. 從這些文獻(xiàn)和相應(yīng)的勢文件來看，這個n的取值應(yīng)根據(jù)文獻(xiàn)里所給的公式，看它們是怎么定的b_{ij}或B_{ij}的表達(dá)式（也就是repulsive potential前面的一個系數(shù)）。 原本的Tersoff勢表達(dá)式中b_{ij}表達(dá)式比較復(fù)雜的。而Albe勢對這個b_{ij}簡化了，簡化為(1+\chi_{ij})^{-1/2}，這樣n就取為1了。 你可以參考B-N-C體系的Tersoff勢文件（我提到的那篇PRB文獻(xiàn)）來試試寫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}。 |
8樓2021-06-10 07:01:24
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
我把這個W-C二元的勢文件根據(jù)這篇文獻(xiàn)整理了的。見https://github.com/houzf/ABOP_lammps 發(fā)自小木蟲Android客戶端 |
2樓2021-06-09 23:35:23
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2021-06-09 23:39:33
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2021-06-09 23:42:38
半寸日光
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 800.2
- 帖子: 259
- 在線: 40.3小時
- 蟲號: 23577820
- 注冊: 2020-09-21
- 專業(yè): 機械工程
|
我之前在你的主頁上參考了這個，照著你的和lammps里面的SiC的Albe一起參考，W-C這個除了某些n，我給了0，沒給你寫的1，其他基本一樣，這個n的0和1對最后的結(jié)果有什么影響嗎？除了W-C，文獻(xiàn)里面給的Co-C，W-Co，Co-Co，整個W-C-Co三個元素體系，我也是按照W-C這只包含兩種元素思考方式來的，比如WCCo，我就根據(jù)WCC這樣的考慮，僅僅當(dāng)做Co替換了C，把它當(dāng)成類似二種元素來思考，這樣思考可以嗎？ 發(fā)自小木蟲Android客戶端 |
5樓2021-06-10 01:01:10
半寸日光
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 800.2
- 帖子: 259
- 在線: 40.3小時
- 蟲號: 23577820
- 注冊: 2020-09-21
- 專業(yè): 機械工程
6樓2021-06-10 01:02:14
半寸日光
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 800.2
- 帖子: 259
- 在線: 40.3小時
- 蟲號: 23577820
- 注冊: 2020-09-21
- 專業(yè): 機械工程
|
這個我對比了一下和tersoff的Albe格式所用的公式，兩個一樣，我看了lammps手冊不是很清楚三元素那種簡化過的Albe格式可以和二元素那種一樣考慮嗎？比如WCCo，我可以看成WCC中C僅是被Co替換了嗎？其他思考方式還是按照二元素一樣來看嗎？比如說里面僅涉及W-C兩種元素的單獨拿出來放到一個文本里面就和體系只有兩元素W-C的tersoff一樣。 發(fā)自小木蟲Android客戶端 |
7樓2021-06-10 01:12:03
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2021-06-10 07:12:40
半寸日光
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 800.2
- 帖子: 259
- 在線: 40.3小時
- 蟲號: 23577820
- 注冊: 2020-09-21
- 專業(yè): 機械工程
10樓2021-06-10 08:49:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 半寸日光 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +9 | 稻妻小編 2026-03-09 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 化工284求調(diào)劑 +12 | 小黑想變白 2026-03-07 | 12/600 |
|
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級已過 +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
|
[考研] 求調(diào)劑，一志愿江南大學(xué)，食品科學(xué)與工程，總分，320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
|
[考研] 286求調(diào)劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
|
[考研] 22408 275分求調(diào)劑 +3 | 宇智波比 2026-03-03 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
|
[考研] 085701環(huán)境工程專業(yè)，初試305，均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
|
[考研] 材料專碩調(diào)劑 +5 | 慕辰123 2026-03-05 | 8/400 |
|
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|
|
[考研] 復(fù)試調(diào)劑 +7 | 呼呼？~+123456 2026-03-05 | 10/500 |
|
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
|
[考研] 274求調(diào)劑 +9 | 一個學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 泡了個椒 2026-03-04 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|