 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3822 | 回復(fù): 14 | ||
| 【懸賞金幣】回答本帖問題，作者半寸日光將贈(zèng)送您 20 個(gè)金幣 | ||
半寸日光新蟲 (小有名氣)
|
[求助]
求助：根據(jù)文獻(xiàn)參數(shù)寫適合lammps中tersoff勢的Albe格式問題 已有1人參與
|
|
|
求助：請(qǐng)問有人根據(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ù)此樓» 猜你喜歡
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有9人回復(fù)
-
328求調(diào)劑
已經(jīng)有3人回復(fù)
-
歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有21人回復(fù)
-
材料與化工304求調(diào)劑
已經(jīng)有6人回復(fù)
-
282求調(diào)劑
已經(jīng)有3人回復(fù)
-
347求調(diào)劑
已經(jīng)有8人回復(fù)
-
成果系統(tǒng)訪問量大，請(qǐng)15分鐘后再嘗試。由此給您造成的不便，敬請(qǐng)諒解。
已經(jīng)有16人回復(fù)
-
080500材料科學(xué)與工程
已經(jīng)有10人回復(fù)
-
材料類考研調(diào)劑
已經(jīng)有8人回復(fù)
-
化學(xué)專業(yè)調(diào)劑
已經(jīng)有3人回復(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勢格式不確定的地方就是這個(gè)n的值。從Albe勢的公式里，看不出這個(gè)n對(duì)應(yīng)的是哪個(gè)物理量。lammps程序里附帶的GaN的Albe勢里，n的值都是設(shè)置1的。我測試過，將n=0與n=1，其他的都不變的話，兩個(gè)勢給的結(jié)果基本一樣。這個(gè)我在我的github做過注釋。 從GaN勢的出處文獻(xiàn)來看，找不到n在公式里對(duì)應(yīng)的變量。 從10.1103/PhysRevB.86.115410文獻(xiàn)里（針對(duì)的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勢對(duì)這個(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對(duì)應(yīng)的就是cos(\theta_{ijk})的值。2\mu對(duì)應(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小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
我把這個(gè)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小時(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è)不是很清楚三元素那種簡化過的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ā)表 | |
|---|---|---|---|---|
|
[考研] 328求調(diào)劑 +3 | vuzhdkfjkx 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 302材料工程求調(diào)劑 +7 | Doleres 2026-03-01 | 8/400 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +3 | T123 tt 2026-03-04 | 3/150 |
|
|
|
[考研] 070300化學(xué) 280 一志愿太原理工 求調(diào)劑 +8 | 拾玖壹 2026-03-04 | 8/400 |
|
|
|
[考研] 376求調(diào)劑 +3 | 王浩然sam 2026-03-04 | 3/150 |
|
|
|
[考研] 材料復(fù)試調(diào)劑 +7 | 學(xué)材料的點(diǎn) 2026-03-01 | 8/400 |
|
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
|
[考博] 26申博-目前4篇SCI一作 +3 | chen_2024 2026-03-02 | 3/150 |
|
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
|
[考研] 主題 +3 | realstar2006 2026-02-27 | 3/150 |
|
|
|
[考研] 0856材料工程，初試313調(diào)劑 +7 | 賣個(gè)關(guān)子吧 2026-03-03 | 7/350 |
|
|
|
[考研] 085700資環(huán)求調(diào)劑，初始279，六級(jí)已過，英語能力強(qiáng) +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
|
[考研] 【2026 碩士調(diào)劑】課題組 招收調(diào)劑生 +3 | 考研版棒棒 2026-03-02 | 5/250 |
|
|
|
[考研] 295求調(diào)劑。一志愿報(bào)考鄭州大學(xué)化學(xué)工藝學(xué)碩，總分295分 +8 | yl1 2026-03-02 | 9/450 |
|
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
|
[考研] 299求調(diào)劑 +3 | Y墨明棋妙Y 2026-02-28 | 5/250 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|