 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3820 | 回復(fù): 14 | |||
| 【懸賞金幣】回答本帖問題，作者半寸日光將贈(zèng)送您 20 個(gè)金幣 | |||
半寸日光新蟲 (小有名氣)
|
[求助]
求助：根據(jù)文獻(xiàn)參數(shù)寫適合lammps中tersoff勢(shì)的Albe格式問題 已有1人參與
|
||
|
求助：請(qǐng)問有人根據(jù)文獻(xiàn)里面的參數(shù)寫過適合lammps中的tersoff勢(shì)？可以幫忙看下我根據(jù)文獻(xiàn)參數(shù)所寫關(guān)于三種元素的tersoff中Albe格式正確嗎？(參考lammps自帶的SiC的Albe)，文獻(xiàn)里面的參數(shù)，自己所仿照寫的W-C-Co體系的勢(shì)函數(shù)如下，感謝! @smutao @oxox6085 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有11人回復(fù)
-
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。
已經(jīng)有14人回復(fù)
-
281電子信息求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料277分求調(diào)劑
已經(jīng)有7人回復(fù)
-
求調(diào)劑院校
已經(jīng)有12人回復(fù)
-
293求調(diào)劑
已經(jīng)有3人回復(fù)
-
化學(xué) 0703求調(diào)劑 總分293 一志愿211
已經(jīng)有10人回復(fù)
-
320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn)
已經(jīng)有7人回復(fù)
-
材料328求調(diào)劑
已經(jīng)有8人回復(fù)
-
304分材料專碩求調(diào)劑
已經(jīng)有12人回復(fù)
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
我現(xiàn)在對(duì)lammps里的這個(gè)Albes勢(shì)格式不確定的地方就是這個(gè)n的值。從Albe勢(shì)的公式里，看不出這個(gè)n對(duì)應(yīng)的是哪個(gè)物理量。lammps程序里附帶的GaN的Albe勢(shì)里，n的值都是設(shè)置1的。我測(cè)試過，將n=0與n=1，其他的都不變的話，兩個(gè)勢(shì)給的結(jié)果基本一樣。這個(gè)我在我的github做過注釋。 從GaN勢(shì)的出處文獻(xiàn)來看，找不到n在公式里對(duì)應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻(xiàn)里（針對(duì)的BNC體系的Albes或Tersoff勢(shì)）的公式來看，n是一個(gè)指數(shù)因子，并也出現(xiàn)在一個(gè)倒數(shù)的分母里，所以n=0似乎是不合理的。這也是讓我困惑的地方。 你參考一下BNC體系的Tersoff勢(shì)： https://www.ctcms.nist.gov/poten ... B-N-C/1/BNC.tersoff 這個(gè)勢(shì)文件里，n并不等于0或1. 從這些文獻(xiàn)和相應(yīng)的勢(shì)文件來看，這個(gè)n的取值應(yīng)根據(jù)文獻(xiàn)里所給的公式，看它們是怎么定的b_{ij}或B_{ij}的表達(dá)式（也就是repulsive potential前面的一個(gè)系數(shù)）。 原本的Tersoff勢(shì)表達(dá)式中b_{ij}表達(dá)式比較復(fù)雜的。而Albe勢(shì)對(duì)這個(gè)b_{ij}簡(jiǎn)化了，簡(jiǎn)化為(1+\chi_{ij})^{-1/2}，這樣n就取為1了。 你可以參考B-N-C體系的Tersoff勢(shì)文件（我提到的那篇PRB文獻(xiàn)）來試試寫W-C-Co三元體系的Albe勢(shì)文件。 與三體相關(guān)的項(xiàng)只是b_ij、\chi_{ij}、g(\theta_{ijk})、cos(theta_{ijk})的項(xiàng)。那篇W-C-Co的勢(shì)參考文件中勢(shì)參數(shù)列表中h對(duì)應(yīng)的就是cos(\theta_{ijk})的值。2\mu對(duì)應(yīng)著lammps的Tersoff勢(shì)參數(shù)格式中的\lambda_{3}。 |
8樓2021-06-10 07:01:24
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
我把這個(gè)W-C二元的勢(shì)文件根據(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小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 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小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 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小時(shí)
- 蟲號(hào): 23577820
- 注冊(cè): 2020-09-21
- 專業(yè): 機(jī)械工程
|
我之前在你的主頁上參考了這個(gè)，照著你的和lammps里面的SiC的Albe一起參考，W-C這個(gè)除了某些n，我給了0，沒給你寫的1，其他基本一樣，這個(gè)n的0和1對(duì)最后的結(jié)果有什么影響嗎？除了W-C，文獻(xiàn)里面給的Co-C，W-Co，Co-Co，整個(gè)W-C-Co三個(gè)元素體系，我也是按照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小時(shí)
- 蟲號(hào): 23577820
- 注冊(cè): 2020-09-21
- 專業(yè): 機(jī)械工程
6樓2021-06-10 01:02:14
半寸日光
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 800.2
- 帖子: 259
- 在線: 40.3小時(shí)
- 蟲號(hào): 23577820
- 注冊(cè): 2020-09-21
- 專業(yè): 機(jī)械工程
|
這個(gè)我對(duì)比了一下和tersoff的Albe格式所用的公式，兩個(gè)一樣，我看了lammps手冊(cè)不是很清楚三元素那種簡(jiǎn)化過的Albe格式可以和二元素那種一樣考慮嗎？比如WCCo，我可以看成WCC中C僅是被Co替換了嗎？其他思考方式還是按照二元素一樣來看嗎？比如說里面僅涉及W-C兩種元素的單獨(dú)拿出來放到一個(gè)文本里面就和體系只有兩元素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小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 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小時(shí)
- 蟲號(hào): 23577820
- 注冊(cè): 2020-09-21
- 專業(yè): 機(jī)械工程
10樓2021-06-10 08:49:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 半寸日光 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研]
材料325求調(diào)劑
 30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 085700資環(huán)求調(diào)劑，初始279，六級(jí)已過，英語能力強(qiáng) +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +15 | February_Feb 2026-03-01 | 17/850 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學(xué)材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|