 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 4417 | 回復(fù): 7 | ||
minbad木蟲 (正式寫手)
|
[求助]
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)能逐漸變的很大，即原先鍵連的原子分開了，所以找不到。 不知道為什么它們會(huì)分開�。� 懇請(qǐng)大俠幫忙 |
回復(fù)此樓» 猜你喜歡
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦
已經(jīng)有8人回復(fù)
-
本人考085602 化學(xué)工程 專碩
已經(jīng)有21人回復(fù)
-
265求調(diào)劑
已經(jīng)有12人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中南化學(xué)337求調(diào)劑
已經(jīng)有6人回復(fù)
-
350求調(diào)劑
已經(jīng)有3人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有22人回復(fù)
-
287求調(diào)劑
已經(jīng)有7人回復(fù)
-
329求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 請(qǐng)問(wèn)lammps初始溫度過(guò)高的問(wèn)題
已經(jīng)有4人回復(fù)
- 關(guān)于lammps fix npt/deform做拉伸時(shí)溫度的變化
已經(jīng)有6人回復(fù)
- LAMMPS溫度控制----如何達(dá)到平衡態(tài)
已經(jīng)有4人回復(fù)
- LAMMPS溫度－步長(zhǎng)曲線震蕩
已經(jīng)有11人回復(fù)
- 在LAMMPS中怎么控制溫度變化速率
已經(jīng)有4人回復(fù)
ktwyw
木蟲 (小有名氣)
- 模擬EPI: 1
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2727.6
- 散金: 500
- 紅花: 41
- 帖子: 162
- 在線: 114.6小時(shí)
- 蟲號(hào): 143500
- 注冊(cè): 2005-12-22
- 專業(yè): 高分子物理與高分子物理化
2樓2014-12-17 18:40:13
minbad
木蟲 (正式寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1209.6
- 散金: 20
- 紅花: 1
- 帖子: 609
- 在線: 80.1小時(shí)
- 蟲號(hào): 235223
- 注冊(cè): 2006-04-01
- 性別: GG
- 專業(yè): 化學(xué)動(dòng)力學(xué)
3樓2014-12-17 19:03:48
ktwyw
木蟲 (小有名氣)
- 模擬EPI: 1
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2727.6
- 散金: 500
- 紅花: 41
- 帖子: 162
- 在線: 114.6小時(shí)
- 蟲號(hào): 143500
- 注冊(cè): 2005-12-22
- 專業(yè): 高分子物理與高分子物理化
4樓2014-12-17 19:11:22
minbad
木蟲 (正式寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 1209.6
- 散金: 20
- 紅花: 1
- 帖子: 609
- 在線: 80.1小時(shí)
- 蟲號(hào): 235223
- 注冊(cè): 2006-04-01
- 性別: GG
- 專業(yè): 化學(xué)動(dòng)力學(xué)
5樓2014-12-17 19:26:04
ktwyw
木蟲 (小有名氣)
- 模擬EPI: 1
- 應(yīng)助: 11 (小學(xué)生)
- 金幣: 2727.6
- 散金: 500
- 紅花: 41
- 帖子: 162
- 在線: 114.6小時(shí)
- 蟲號(hào): 143500
- 注冊(cè): 2005-12-22
- 專業(yè): 高分子物理與高分子物理化
6樓2014-12-17 19:55:19
7樓2014-12-18 19:09:47
hustli2014
新蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 39.2
- 散金: 82
- 紅花: 2
- 帖子: 84
- 在線: 134.5小時(shí)
- 蟲號(hào): 3093527
- 注冊(cè): 2014-03-27
- 性別: GG
- 專業(yè): 高分子材料的加工與成型
8樓2016-07-18 15:54:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 minbad 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +7 | Ncdx123456 2026-03-19 | 7/350 |
|
|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3 | lemonzzn 2026-03-17 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|