 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1563 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題，作者木棉和我將贈送您 100 個(gè)金幣 | ||
木棉和我木蟲 (職業(yè)作家)
|
[求助]
求助：建立的初始簡單晶體體系(NaCl結(jié)構(gòu))，npt弛豫時(shí)即發(fā)生體系“爆炸” 已有1人參與
|
|
|
如題，建立了一個(gè)NbC完美晶體，尺寸與取向分別為：10a*10a*10a; X-[100], Y-[010], and Z-[001]。 準(zhǔn)備簡單計(jì)算一下相關(guān)勢函數(shù)下的晶格常數(shù)，即在5K下進(jìn)行npt保溫弛豫，模型開始計(jì)算則發(fā)生“爆炸”。具體Lammps運(yùn)行腳本如下： units metal boundary p p p atom_style atomic read_data nbc_c10-10-10.lmp #perfect NbC crystal (NaCl structure)；由Atomsk建立. pair_style meam/c pair_coeff * * library.meam Nb Fe C NbFeC.meam Nb C timestep 0.001 thermo 100 thermo_style custom step temp pe ke etotal epair press vol lx ly lz restart 100000 fcc-Fe.restart neighbor 0.5 bin neigh_modify every 2 delay 10 check yes dump 1 all custom 100000 dump.lattice-nbc-5K.* id type x y z velocity all create 5 123456 dist gaussian fix 1 all npt temp 5 5 0.1 iso 0 0 0.1 run 500000 unfix 1 undump 1 相關(guān)勢函數(shù)、in文件及初始晶體請見附件。 另外，MD構(gòu)建初始晶體時(shí)，是否需要考慮晶體結(jié)構(gòu)的化學(xué)鍵類型？還是只是需要把各類原子的初始位置確定，而化學(xué)鍵類型主要影響作用力，作用力通過勢函數(shù)體現(xiàn)？或者說，建立MD結(jié)構(gòu)模型時(shí)，是不是僅有位置信息，沒有其他信息，其他信息由勢函數(shù)決定。 請相關(guān)老師們幫忙解答一下，十分感謝。 |
回復(fù)此樓» 猜你喜歡
26申博自薦
已經(jīng)有8人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有15人回復(fù)
-
我嘗試在小木蟲上招計(jì)算機(jī)專業(yè)的調(diào)劑學(xué)生
已經(jīng)有6人回復(fù)
-
售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急
已經(jīng)有5人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有4人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有3人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有3人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有3人回復(fù)
-
一志愿山東大學(xué)105500藥學(xué)專碩，總分302求調(diào)劑
已經(jīng)有9人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有5人回復(fù)
valenhou001
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2021-12-11 06:34:54
3樓2021-12-12 00:26:14
木棉和我
木蟲 (職業(yè)作家)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 8996.1
- 散金: 2312
- 紅花: 3
- 帖子: 3299
- 在線: 676.2小時(shí)
- 蟲號: 1852823
- 注冊: 2012-06-08
- 性別: GG
- 專業(yè): 冶金物理化學(xué)與冶金原理
4樓2021-12-13 08:34:26
木棉和我
木蟲 (職業(yè)作家)
- 應(yīng)助: 5 (幼兒園)
- 金幣: 8996.1
- 散金: 2312
- 紅花: 3
- 帖子: 3299
- 在線: 676.2小時(shí)
- 蟲號: 1852823
- 注冊: 2012-06-08
- 性別: GG
- 專業(yè): 冶金物理化學(xué)與冶金原理
|
謝謝老師您的解答！ 1. 也就說MD晶體構(gòu)型中，只有位置信息(結(jié)果中，可以根據(jù)原子位置與速度，提取相關(guān)參量)；然后原子間作用力的影響因素，如化學(xué)鍵類型、電荷等，是直接由勢函數(shù)體現(xiàn)，對嗎？ 2. 這個(gè)模型是這樣的，就是根據(jù)文獻(xiàn)中的勢函數(shù)計(jì)算出的NbC晶格常數(shù)，我建立了一個(gè)完美NbC晶體(NaCl結(jié)構(gòu))。我想確定一下，這個(gè)勢函數(shù)下NbC的晶格常數(shù)是不是如文獻(xiàn)中所述，于是簡單的對完美NbC晶體做了NPT保溫，結(jié)果體系“爆炸”了。我也嘗試過對這個(gè)NbC晶體進(jìn)行cg法的能量最小化，但最小化好像對整個(gè)晶體原子位置毫無影響，因?yàn)楸旧斫⒌腘bC晶體就是完美的Nb：C=1:1的NaCl結(jié)構(gòu)。 |
5樓2021-12-13 08:47:32
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 木棉和我 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[找工作] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-06 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 273求調(diào)劑 +5 | 星星111222 2026-03-02 | 7/350 |
|
|
|
[考研] 一志愿211，0860總分286食品工程專業(yè)求調(diào)劑 +6 | 林林Winnie 2026-03-05 | 6/300 |
|
|
|
[考研] 085600材料與化工，一志愿廣州985，求調(diào)劑 +10 | qqyyaill 2026-03-05 | 10/500 |
|
|
|
[考研]
 材料工程調(diào)劑
+6
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-03-05 | 4/200 |
|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
|
[考研] 0856材料與化工，270求調(diào)劑 +17 | YXCT 2026-03-01 | 20/1000 |
|
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 085601 材料工程 320 +6 | 和樂瑤 2026-03-03 | 6/300 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
|
[考研] 291分工科求調(diào)劑 +9 | science餓餓 2026-03-01 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|