| 5 | 1/1 | 返回列表 |
| 查看: 1572 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者木棉和我將贈(zèng)送您 100 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
木棉和我木蟲 (職業(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é)鍵類型?還是只是需要把各類原子的初始位置確定,而化學(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)。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們,祝我生日快樂(lè)吧 +8 | myrtle 2026-03-10 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 考研一志愿長(zhǎng)安大學(xué)材料與化工309分請(qǐng)求調(diào)劑 +7 | dtdxzxx 2026-03-06 | 9/450 |
|
|
[考研] 0856材料與化工309分求調(diào)劑 +5 | ZyZy…… 2026-03-10 | 5/250 |
|
|
[考研] 314求調(diào)劑 +4 | 無(wú)懈可擊的巨人 2026-03-12 | 4/200 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評(píng)估A類院;瘜W(xué)工程與工藝,誠(chéng)求調(diào)劑 +4 | 小劉快快上岸 2026-03-11 | 4/200 |
|
|
[考研] 301求調(diào)劑 +5 | Liyouyumairs 2026-03-11 | 5/250 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-11 | 4/200 |
|
|
[考研] 收調(diào)劑 +7 | 調(diào)劑的考研學(xué)生 2026-03-10 | 7/350 |
|
|
[考研] 一志愿天大化工(085600)調(diào)劑總分338 +5 | 蔡大美女 2026-03-09 | 5/250 |
|
|
[考研] 一志愿:武漢理工,材料工程,英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
[考博] 找博導(dǎo) +4 | 小呆呆熊 2026-03-07 | 4/200 |
|
|
[考研] 085600材料與化工 326 求調(diào)劑 +4 | 熱愛(ài)生活ing 2026-03-09 | 4/200 |
|
|
[考研] 310 070300化學(xué)求調(diào)劑 +4 | 撲風(fēng)鈴的貓 2026-03-08 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑,總分309分,一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過(guò)國(guó)家A區(qū)線 +7 | 卡卡來(lái)了@ 2026-03-07 | 8/400 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|