| 查看: 2479 | 回復(fù): 15 | ||
解子林新蟲 (小有名氣)
|
[求助]
使用berendsen方法控溫控壓問題 已有1人參與
|
|
|
在使用berendsen對含水的MOF進行控溫控壓模擬時,同樣的fix條件,不同的含水量有的可以跑,有的出現(xiàn)溫度升高,體系分解的狀況,哪位大神幫忙給點建議!以下是我的in文件: # ReaxFF potential for MOF units real boundary p p p atom_style charge read_data data.IRMOF1+16H2O-BER.txt mass 1 12.0107 #苯環(huán)上的C mass 2 12.0107 #羧基C mass 3 1.00794 #苯環(huán)上的H mass 4 15.9994 #中心O mass 5 15.9994 #橋O mass 6 65.38 #Zn mass 7 15.9994 #水中的O mass 8 1.00794 #水中的H pair_style reax 6.0 1 1 1.0e-6 pair_coeff * * ffield.reax 1 1 2 3 3 4 3 2 compute reax all pair reax variable eb equal c_reax[1] variable ea equal c_reax[2] variable elp equal c_reax[3] variable emol equal c_reax[4] variable ev equal c_reax[5] variable epen equal c_reax[6] variable ecoa equal c_reax[7] variable ehb equal c_reax[8] variable et equal c_reax[9] variable eco equal c_reax[10] variable ew equal c_reax[11] variable ep equal c_reax[12] variable efi equal c_reax[13] variable eqeq equal c_reax[14] neighbor 2.5 bin neigh_modify every 10 delay 0 check no group IRMOF1 type 1 2 3 4 5 6 group water type 7 8 #minimize 1.0e-4 1.0e-6 100 1000 velocity all create 300.0 825577 dist gaussian units box fix 1 all reax/bonds 1000 bonds.reax fix 2 all nve fix 4 all temp/berendsen 300.0 300.0 10.0 fix 6 all press/berendsen iso 1.0 1.0 1000.0 modulus 1.0e6 run_style verlet dump 1 all custom 1000 trj.IRMOF1+20water-1 id type q x y z #dump_modify 1 scale no dump_modify 1 sort 1 restart 100 restart1 restart2 thermo 1000 thermo_style custom step cpu temp evdwl ecoul epair etotal press vol enthalpy & v_eb v_ea v_elp v_emol v_ev v_epen v_ecoa v_ehb & v_et v_eco v_ew v_ep v_efi v_eqeq thermo_modify flush yes timestep 1 run 2000000 |
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
新蟲 (小有名氣)
新蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
|
新蟲 (小有名氣)
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 小熊joy 2026-03-14 | 7/350 |
|
|
[考研] 301求調(diào)劑 +3 | A_JiXing 2026-03-16 | 3/150 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 2026考研調(diào)劑+本科延邊大學(xué)+山東大學(xué)+生物化學(xué)與分子生物學(xué)+有項目經(jīng)驗 +3 | ccdsscjy 2026-03-10 | 3/150 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|
|
[考研] 085600 材料與化工 295 求調(diào)劑 +10 | dream…… 2026-03-10 | 12/600 |
|