| 5 | 1/1 | 返回列表 |
| 查看: 1571 | 回復: 4 | ||
| 【懸賞金幣】回答本帖問題,作者木棉和我將贈送您 100 個金幣 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
木棉和我木蟲 (職業(yè)作家)
|
[求助]
求助:建立的初始簡單晶體體系(NaCl結(jié)構(gòu)),npt弛豫時即發(fā)生體系“爆炸” 已有1人參與
|
|
|
如題,建立了一個NbC完美晶體,尺寸與取向分別為:10a*10a*10a; X-[100], Y-[010], and Z-[001]。 準備簡單計算一下相關(guān)勢函數(shù)下的晶格常數(shù),即在5K下進行npt保溫弛豫,模型開始計算則發(fā)生“爆炸”。具體Lammps運行腳本如下: 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)建初始晶體時,是否需要考慮晶體結(jié)構(gòu)的化學鍵類型?還是只是需要把各類原子的初始位置確定,而化學鍵類型主要影響作用力,作用力通過勢函數(shù)體現(xiàn)?或者說,建立MD結(jié)構(gòu)模型時,是不是僅有位置信息,沒有其他信息,其他信息由勢函數(shù)決定。 請相關(guān)老師們幫忙解答一下,十分感謝。 |
至尊木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
|
謝謝老師您的解答! 1. 也就說MD晶體構(gòu)型中,只有位置信息(結(jié)果中,可以根據(jù)原子位置與速度,提取相關(guān)參量);然后原子間作用力的影響因素,如化學鍵類型、電荷等,是直接由勢函數(shù)體現(xiàn),對嗎? 2. 這個模型是這樣的,就是根據(jù)文獻中的勢函數(shù)計算出的NbC晶格常數(shù),我建立了一個完美NbC晶體(NaCl結(jié)構(gòu))。我想確定一下,這個勢函數(shù)下NbC的晶格常數(shù)是不是如文獻中所述,于是簡單的對完美NbC晶體做了NPT保溫,結(jié)果體系“爆炸”了。我也嘗試過對這個NbC晶體進行cg法的能量最小化,但最小化好像對整個晶體原子位置毫無影響,因為本身建立的NbC晶體就是完美的Nb:C=1:1的NaCl結(jié)構(gòu)。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703求調(diào)劑 +6 | jtyq001 2026-03-10 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +5 | Liyouyumairs 2026-03-11 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | zchqwer 2026-03-10 | 6/300 |
|
|
[考研] 一志愿河海大學085900土木水利專碩279求調(diào)劑不挑專業(yè) +3 | SunWwWwWw 2026-03-10 | 3/150 |
|
|
[考研] 材料工程307,求調(diào)劑 +7 | 我要燃燒你的夢 2026-03-08 | 7/350 |
|
|
[考研] 420求調(diào)劑 +3 | 莫向外求11 2026-03-10 | 3/150 |
|
|
[考研] 化工0817調(diào)劑 +8 | 燦若星晨 2026-03-10 | 8/400 |
|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +9 | dream…… 2026-03-10 | 9/450 |
|
|
[考研] 研究生招生 +3 | 徐海濤11 2026-03-10 | 6/300 |
|
|
[考研] 求調(diào)劑 一志愿蘇州大學,0856化工323分 | 本科應化 | 有專利/競賽/科研助手經(jīng)歷 | +7 | 橙子cyx 2026-03-06 | 9/450 |
|
|
[考研] 考研求調(diào)劑,工科,總分304 +4 | 矮子不想長高了 2026-03-08 | 7/350 |
|
|
[考研] 理學287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 求調(diào)劑,不管什么專業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | 呼呼?~+123456 2026-03-06 | 4/200 |
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風了呢 2026-03-05 | 9/450 |
|
|
[考研]
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|