| 5 | 1/1 | 返回列表 |
| 查看: 1562 | 回復(fù): 4 | ||
| 【懸賞金幣】回答本帖問題,作者木棉和我將贈(zèng)送您 100 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
木棉和我木蟲 (職業(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文件及初始晶體請(qǐng)見附件。 另外,MD構(gòu)建初始晶體時(shí),是否需要考慮晶體結(jié)構(gòu)的化學(xué)鍵類型?還是只是需要把各類原子的初始位置確定,而化學(xué)鍵類型主要影響作用力,作用力通過勢函數(shù)體現(xiàn)?或者說,建立MD結(jié)構(gòu)模型時(shí),是不是僅有位置信息,沒有其他信息,其他信息由勢函數(shù)決定。 請(qǐng)相關(guān)老師們幫忙解答一下,十分感謝。 |
至尊木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
木蟲 (職業(yè)作家)
|
謝謝老師您的解答! 1. 也就說MD晶體構(gòu)型中,只有位置信息(結(jié)果中,可以根據(jù)原子位置與速度,提取相關(guān)參量);然后原子間作用力的影響因素,如化學(xué)鍵類型、電荷等,是直接由勢函數(shù)體現(xiàn),對(duì)嗎? 2. 這個(gè)模型是這樣的,就是根據(jù)文獻(xiàn)中的勢函數(shù)計(jì)算出的NbC晶格常數(shù),我建立了一個(gè)完美NbC晶體(NaCl結(jié)構(gòu))。我想確定一下,這個(gè)勢函數(shù)下NbC的晶格常數(shù)是不是如文獻(xiàn)中所述,于是簡單的對(duì)完美NbC晶體做了NPT保溫,結(jié)果體系“爆炸”了。我也嘗試過對(duì)這個(gè)NbC晶體進(jìn)行cg法的能量最小化,但最小化好像對(duì)整個(gè)晶體原子位置毫無影響,因?yàn)楸旧斫⒌腘bC晶體就是完美的Nb:C=1:1的NaCl結(jié)構(gòu)。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩,總分302求調(diào)劑 +4 | 五維天空 2026-03-04 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 0307化學(xué)求調(diào)劑 +4 | 0ok0k 2026-03-07 | 4/200 |
|
|
[考研] 誠求調(diào)劑,323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
[考研] 一志愿清華深研院材料專碩294分,專業(yè)課111分,本科中南大學(xué)材料,有六級(jí),有工作經(jīng)驗(yàn) +3 | H14528 2026-03-04 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
[考研] 中國科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑,有科研和獲獎(jiǎng)經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
[考研] 一志愿鄭州大學(xué),學(xué)碩,物理化學(xué), 333求調(diào)劑 +5 | 李魔女斗篷 2026-03-04 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 化工專碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|