| 5 | 1/1 | 返回列表 |
| 查看: 1568 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題,作者木棉和我將贈(zèng)送您 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)老師們幫忙解答一下,十分感謝。 |
至尊木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
|
謝謝老師您的解答! 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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +9 | 王曉陽- 2026-03-09 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 材料與化工 一志愿山大 321分 求調(diào)劑 +4 | 每天散步 2026-03-09 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +6 | 小奶蒼 2026-03-07 | 6/300 |
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩,總分302求調(diào)劑 +5 | 五維天空 2026-03-04 | 11/550 |
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
[考博] 26年博士申請 +4 | 科研狗111 2026-03-07 | 4/200 |
|
|
[考研] 301求調(diào)劑 +11 | 朝天椒的雙馬尾 2026-03-05 | 11/550 |
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
[考研] 301求調(diào)劑 +5 | 一二LV 2026-03-07 | 5/250 |
|
|
[考研] 第一志愿上海大學(xué),專業(yè)化學(xué)工程與技術(shù),總分288,求調(diào)劑 +3 | 1829197082 2026-03-07 | 3/150 |
|
|
[考研] 278求調(diào)劑 +5 | Gale1314 2026-03-06 | 5/250 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 7/350 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
[考研] 一志愿清華深研院材料專碩294分,專業(yè)課111分,本科中南大學(xué)材料,有六級(jí),有工作經(jīng)驗(yàn) +3 | H14528 2026-03-04 | 3/150 |
|
|
[考研] 304分材料專碩求調(diào)劑 +11 | qiuzhigril 2026-03-03 | 14/700 |
|
|
[考研] 320材料一志愿西工大專碩英二數(shù)二 有兩年光伏工作經(jīng)驗(yàn) +5 | 錘某人 2026-03-04 | 8/400 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|