| 5 | 1/1 | 返回列表 |
| 查看: 4422 | 回復(fù): 7 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
minbad木蟲(chóng) (正式寫(xiě)手)
|
[求助]
Lammps fixnve 為什么溫度也一直升高? 已有2人參與
|
|
|
一個(gè)有機(jī)分子和一個(gè)水,進(jìn)行測(cè)試模擬,溫度一致升高, 輸入in文件 —————————————————————————————————— 1 ## Initialization 2 units real #default lj 3 atom_style full #default atomic 4 pair_style lj/cut/coul/cut 10.0 5 pair_modify mix geometric 6 bond_style harmonic 7 angle_style harmonic 8 dihedral_style opls 9 10 read_data data.glucose1wat 11 12 neighbor 2 bin 13 neigh_modify delay 0 every 1 check yes 14 15 timestep 0.5 # default 1fs in real units 16 17 group all molecule 1:2 18 19 thermo 20 20 thermo_style multi 21 22 dump min all custom 500 min.lammpstrj id mol type x y z ix iy iz 23 dump_modify min sort id 24 minimize 1.0e-8 1.0e-8 500 20000 25 undump min 26 27 dump 2 all custom 20 nveMD1.lammpstrj id mol type x y z ix iy iz 28 dump_modify 2 sort id 29 fix fixnve all nve #temp 300.0 300.0 100 30 run 20000 31 undump 2 32 unfix fixnve —————————————————————————————————— 最后就出現(xiàn)bonds atoms missing: ................. ---------------- Step 800 ----- CPU = 0.0614 (sec) ---------------- TotEng = 37433.3818 KinEng = 28019.5768 Temp = 361537.9752 PotEng = 9413.8050 E_bond = 8593.2434 E_angle = 694.9632 E_dihed = 30.0571 E_impro = 0.0000 E_vdwl = 198.2975 E_coul = -102.7561 E_long = 0.0000 Press = 42398.0431 ---------------- Step 820 ----- CPU = 0.0632 (sec) ---------------- TotEng = 173543.1558 KinEng = 50823.3926 Temp = 655776.7303 PotEng = 122719.7632 E_bond = 121031.4794 E_angle = 1585.6246 E_dihed = 35.0079 E_impro = 0.0000 E_vdwl = 186.2945 E_coul = -118.6432 E_long = 0.0000 Press = -262854.6486 ---------------- Step 840 ----- CPU = 0.0648 (sec) ---------------- TotEng = 173729.2625 KinEng = 78851.3657 Temp = 1017423.0440 PotEng = 94877.8968 E_bond = 93040.0890 E_angle = 1838.2278 E_dihed = 42.4633 E_impro = 0.0000 E_vdwl = 15.7058 E_coul = -58.5890 E_long = 0.0000 Press = -110094.8168 ERROR on proc 4: Bond atoms 4 5 missing on proc 4 at step 852 (../neigh_bond.cpp:65) application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4 經(jīng)查,發(fā)現(xiàn)動(dòng)能逐漸變的很大,即原先鍵連的原子分開(kāi)了,所以找不到。 不知道為什么它們會(huì)分開(kāi)。 懇請(qǐng)大俠幫忙 |
木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (小有名氣)

木蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (小有名氣)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 機(jī)械專(zhuān)碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +6 | 煙火先于春 2026-03-17 | 6/300 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
[考研] 一志愿中南化學(xué)(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 290求調(diào)劑 +7 | ^O^乜 2026-03-19 | 7/350 |
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學(xué)289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|