| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1573 | 回復(fù): 4 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者木棉和我將贈(zèng)送您 100 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
木棉和我木蟲(chóng) (職業(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保溫弛豫,模型開(kāi)始計(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)見(jiàn)附件。 另外,MD構(gòu)建初始晶體時(shí),是否需要考慮晶體結(jié)構(gòu)的化學(xué)鍵類(lèi)型?還是只是需要把各類(lèi)原子的初始位置確定,而化學(xué)鍵類(lèi)型主要影響作用力,作用力通過(guò)勢(shì)函數(shù)體現(xiàn)?或者說(shuō),建立MD結(jié)構(gòu)模型時(shí),是不是僅有位置信息,沒(méi)有其他信息,其他信息由勢(shì)函數(shù)決定。 請(qǐng)相關(guān)老師們幫忙解答一下,十分感謝。 |
木蟲(chóng) (職業(yè)作家)
|
謝謝老師您的解答! 1. 也就說(shuō)MD晶體構(gòu)型中,只有位置信息(結(jié)果中,可以根據(jù)原子位置與速度,提取相關(guān)參量);然后原子間作用力的影響因素,如化學(xué)鍵類(lèi)型、電荷等,是直接由勢(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)。 |
至尊木蟲(chóng) (職業(yè)作家)
木蟲(chóng) (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料301分求調(diào)劑 +4 | Liyouyumairs 2026-03-12 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料工程調(diào)劑 +3 | 咪咪空空 2026-03-11 | 3/150 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +5 | mysdl 2026-03-11 | 7/350 |
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽(yáng)- 2026-03-09 | 18/900 |
|
|
[考研] 293求調(diào)劑,一志愿陜師大生物學(xué) +3 | ??????.?.??? 2026-03-09 | 3/150 |
|
|
[考研] 化工0817調(diào)劑 +8 | 燦若星晨 2026-03-10 | 8/400 |
|
|
[考研] 一志愿天津大學(xué),英一數(shù)二305分求調(diào)劑,四六級(jí)已過(guò) +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
[考研] 求調(diào)劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +6 | 呼呼?~+123456 2026-03-08 | 8/400 |
|
|
[考博] 找博導(dǎo) +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
[考研] 085701環(huán)境工程295求調(diào)劑,一志愿A區(qū)211 +8 | 等春來(lái), 2026-03-06 | 9/450 |
|
|
[考研] 320求調(diào)劑 +4 | 魏zy 2026-03-08 | 4/200 |
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研] 085701環(huán)境工程專(zhuān)業(yè),初試305,均過(guò)國(guó)家A區(qū)線 +7 | 卡卡來(lái)了@ 2026-03-07 | 8/400 |
|
|
[考研] 一志愿鄭大071000分?jǐn)?shù)282求調(diào)劑 +3 | 研研顏 2026-03-05 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|