 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3817 | 回復(fù): 14 | |||
| 【懸賞金幣】回答本帖問題，作者半寸日光將贈(zèng)送您 20 個(gè)金幣 | |||
半寸日光新蟲 (小有名氣)
|
[求助]
求助：根據(jù)文獻(xiàn)參數(shù)寫適合lammps中tersoff勢(shì)的Albe格式問題 已有1人參與
|
||
|
求助：請(qǐng)問有人根據(jù)文獻(xiàn)里面的參數(shù)寫過(guò)適合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ù)此樓» 猜你喜歡
322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷
已經(jīng)有3人回復(fù)
-
320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn)
已經(jīng)有3人回復(fù)
-
一志愿鄭州大學(xué)，學(xué)碩，物理化學(xué)， 333求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中石油，化學(xué)工程，301求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué)
已經(jīng)有3人回復(fù)
-
320材料與化工，求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑
已經(jīng)有5人回復(fù)
-
286
已經(jīng)有3人回復(fù)
-
070300化學(xué) 280 一志愿太原理工 求調(diào)劑
已經(jīng)有5人回復(fù)
-
085600材料與化工（高分子）290分求調(diào)劑
已經(jīng)有4人回復(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è)試過(guò)，將n=0與n=1，其他的都不變的話，兩個(gè)勢(shì)給的結(jié)果基本一樣。這個(gè)我在我的github做過(guò)注釋。 從GaN勢(shì)的出處文獻(xiàn)來(lái)看，找不到n在公式里對(duì)應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻(xiàn)里（針對(duì)的BNC體系的Albes或Tersoff勢(shì)）的公式來(lái)看，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ì)文件來(lái)看，這個(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)）來(lái)試試寫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ī)械工程
|
我之前在你的主頁(yè)上參考了這個(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這只包含兩種元素思考方式來(lái)的，比如WCCo，我就根據(jù)WCC這樣的考慮，僅僅當(dāng)做Co替換了C，把它當(dāng)成類似二種元素來(lái)思考，這樣思考可以嗎？ 發(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)化過(guò)的Albe格式可以和二元素那種一樣考慮嗎？比如WCCo，我可以看成WCC中C僅是被Co替換了嗎？其他思考方式還是按照二元素一樣來(lái)看嗎？比如說(shuō)里面僅涉及W-C兩種元素的單獨(dú)拿出來(lái)放到一個(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ā)表 | |
|---|---|---|---|---|
|
[考研] 320材料與化工，求調(diào)劑 +5 | 鶴遨予卿 2026-03-04 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 286 +3 | ksncj 2026-03-04 | 3/150 |
|
|
|
[考研] 學(xué)碩材料275調(diào)劑 +9 | 路三三 2026-03-03 | 9/450 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +6 | xiiiia 2026-03-04 | 6/300 |
|
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|
|
|
[考博] 26申博 +5 | north， 2026-02-28 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[基金申請(qǐng)] 請(qǐng)問大家，研究風(fēng)險(xiǎn)與應(yīng)對(duì)措施那里， 大家都怎么寫呢 ？ +3 | cauasen 2026-03-02 | 3/150 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
|
[考研] 材料類求調(diào)劑 +11 | wana_kiko 2026-02-28 | 14/700 |
|
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|