| 5 | 1/1 | 返回列表 |
| 查看: 4421 | 回復(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) (小有名氣)

木蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽(yáng)CF 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +4 | shxchem 2026-03-20 | 5/250 |
|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
[考研] 287求調(diào)劑 +7 | 晨昏線(xiàn)與星海 2026-03-19 | 8/400 |
|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專(zhuān)材355 2026-03-19 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|