| 5 | 1/1 | 返回列表 |
| 查看: 1560 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者木棉和我將贈(zèng)送您 100 個(gè)金幣 | ||
木棉和我木蟲 (職業(yè)作家)
|
[求助]
求助:建立的初始簡(jiǎn)單晶體體系(NaCl結(jié)構(gòu)),npt弛豫時(shí)即發(fā)生體系“爆炸” 已有1人參與
|
|
|
如題,建立了一個(gè)NbC完美晶體,尺寸與取向分別為:10a*10a*10a; X-[100], Y-[010], and Z-[001]。 準(zhǔn)備簡(jiǎn)單計(jì)算一下相關(guān)勢(shì)函數(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ì)函數(shù)、in文件及初始晶體請(qǐng)見附件。 另外,MD構(gòu)建初始晶體時(shí),是否需要考慮晶體結(jié)構(gòu)的化學(xué)鍵類型?還是只是需要把各類原子的初始位置確定,而化學(xué)鍵類型主要影響作用力,作用力通過(guò)勢(shì)函數(shù)體現(xiàn)?或者說(shuō),建立MD結(jié)構(gòu)模型時(shí),是不是僅有位置信息,沒(méi)有其他信息,其他信息由勢(shì)函數(shù)決定。 請(qǐng)相關(guān)老師們幫忙解答一下,十分感謝。 |
至尊木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
|
謝謝老師您的解答! 1. 也就說(shuō)MD晶體構(gòu)型中,只有位置信息(結(jié)果中,可以根據(jù)原子位置與速度,提取相關(guān)參量);然后原子間作用力的影響因素,如化學(xué)鍵類型、電荷等,是直接由勢(shì)函數(shù)體現(xiàn),對(duì)嗎? 2. 這個(gè)模型是這樣的,就是根據(jù)文獻(xiàn)中的勢(shì)函數(shù)計(jì)算出的NbC晶格常數(shù),我建立了一個(gè)完美NbC晶體(NaCl結(jié)構(gòu))。我想確定一下,這個(gè)勢(shì)函數(shù)下NbC的晶格常數(shù)是不是如文獻(xiàn)中所述,于是簡(jiǎn)單的對(duì)完美NbC晶體做了NPT保溫,結(jié)果體系“爆炸”了。我也嘗試過(guò)對(duì)這個(gè)NbC晶體進(jìn)行cg法的能量最小化,但最小化好像對(duì)整個(gè)晶體原子位置毫無(wú)影響,因?yàn)楸旧斫⒌腘bC晶體就是完美的Nb:C=1:1的NaCl結(jié)構(gòu)。 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
程晴之 2026-03-06 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩,總分302求調(diào)劑 +4 | 五維天空 2026-03-04 | 8/400 |
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
[考研] 290求調(diào)劑 +5 | Bananaiy 2026-03-04 | 5/250 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
[考研] 282求調(diào)劑 +5 | 2103240126 2026-03-02 | 8/400 |
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
[考研] 武理材料工程302調(diào)劑 +3 | Doleres 2026-03-03 | 6/300 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考博] 26申博 +5 | north, 2026-02-28 | 5/250 |
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 1 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[基金申請(qǐng)]
剛錄用,沒(méi)有期刊號(hào),但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-02-28 | 3/150 |
|