Jfreda

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時
蟲號: 870275

[交流] mdrun時出錯

mdrun -s aa.tpr -o aa.trr -c aa.gro &
出現(xiàn)了如下的提示錯誤：
step 11: Water molecule starting at atom 180665 can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

step 13: Water molecule starting at atom 180665 can not be settled.
Check for bad contacts and/or reduce the timestep if appropriate.
Wrote pdb files with previous and current coordinates

Stepsize too small, or no change in energy.
Converged to machine precision,
but not to the requested precision Fmax < 1000

Double precision normally gives you higher accuracy.
You might need to increase your constraint accuracy, or turn
off constraints alltogether (set constraints = none in mdp file)

writing lowest energy coordinates.

Steepest Descents converged to machine precision in 31 steps,
but did not reach the requested Fmax < 1000.
Potential Energy  = -3.1653622e+06
Maximum force    =  9.5183238e+05 on atom 180665
Norm of force    =  2.3416899e+03

下面是我的em.mdp文件
em.mdp

er spoel (236)
;    Wed Nov  3 17:12:44 1993
;    Input file
;
cpp                =  /lib/cpp
constraints       =  none
define             =  -DFLEX_SPC

; Parameters describing what to do, when to stop and what to save
integrator       = steep
dt                   = 0.002
emtol             = 1000.0
emstep             = 0.01
nsteps             = 3000

; Parameters describing how to find the neighbors of each atom and how to calculate the interactions
nstlist          = 10
ns_type          = grid
rlist             = 1.0
coulombtype       = pme
rcoulomb          = 1.0
rvdw             = 1.0
fourierspacing    =0.12
pme_order          =4
optimize_fft       =yes

freezegrps    =Protein UNK POP
freezedim       =Y Y Y Y Y Y Y Y Y
我也嘗試了改變em.mdp文件中的一些參數(shù)(只對一個參數(shù)做了修改)：
1.em.mdp文件中的freezegrps改為Protein UNK能量能很好地收斂，下面是它的顯示結果：
Steepest Descents:
Tolerance (Fmax) =  1.00000e+03
Number of steps =       3000

writing lowest energy coordinates.

Steepest Descents converged to Fmax < 1000 in 1678 steps
Potential Energy  = -4.2245720e+06
Maximum force    =  9.5952905e+02 on atom 9032
Norm of force    =  1.1612597e+01

2.em.mdp文件中的emstep改為0.00001（參照別人的回答), 結果顯示為：
Steepest Descents:
Tolerance (Fmax) =  1.00000e+03
Number of steps =       3000

writing lowest energy coordinates.

Steepest Descents converged to Fmax < 1000 in 32 steps
Potential Energy  = -3.0919205e+06
Maximum force    =  9.1417908e+02 on atom 122811
Norm of force    =  1.1197423e+02
這樣的結果應該是有問題的吧，只跑了32步。

3.integrator改為md, 結果為：
Writing final coordinates.

Average load imbalance: 0.2 %
Part of the total run time spent waiting due to load imbalance: 0.1 %

      Parallel run - timing based on wallclock.

            NODE (s) Real (s)    (%)
   Time: 7515.940 7515.940 100.0
                     2h05:15
            (Mnbf/s) (GFlops) (ns/day)  (hour/ns)
Performance:    63.678    3.946    0.115 208.755

gcq#290: "Oh, There Goes Gravity" (Eminem)

請問一下，我應該怎么做才是合理的？謝謝！

回復此樓

» 猜你喜歡

086502化學工程342求調劑已經(jīng)有4人回復
材料求調劑一志愿哈工大324 已經(jīng)有7人回復
328求調劑已經(jīng)有7人回復
340求調劑已經(jīng)有5人回復
一志愿南昌大學324求調劑已經(jīng)有4人回復
求調劑推薦材料 304 已經(jīng)有15人回復
291求調劑已經(jīng)有19人回復
331環(huán)境科學與工程求調劑已經(jīng)有3人回復
289求調劑已經(jīng)有8人回復
315分求調劑已經(jīng)有7人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

Gromacs中，mdrun出錯已經(jīng)有5人回復
在gromacs中運行mdrun命令時出錯已經(jīng)有10人回復
vsap 5.2 環(huán)境和編譯問題已經(jīng)有14人回復
Gromacs做steep之后做cg出錯已經(jīng)有5人回復
編譯vasp出錯已經(jīng)有11人回復
gromacs做全原子模擬時候，突然到MD這一步出錯已經(jīng)有6人回復
【求助】安裝vasp出錯 make: *** [fftmpi_map.o] 錯誤 1【已解決】已經(jīng)有6人回復
【求助】vasp安裝編譯錯誤，請大家?guī)兔聪履抢锍鲥e了。已經(jīng)有6人回復
【求助成功】Vasp 5.2.11編譯出錯已經(jīng)有10人回復
【求助】在編譯時出錯已經(jīng)有6人回復
【求助】PBE0 計算Si、Ge能帶結構已經(jīng)有9人回復

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-04-18 18:28:08

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jawang

木蟲 (正式寫手)

應助: 3 (幼兒園)
金幣: 3533.7
帖子: 376
在線: 214.7小時
蟲號: 838608

★ ★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
Jfreda(金幣+5): 2011-04-18 19:08:54
zh1987hs(金幣+3): 謝謝 2011-04-18 21:33:02

Steepest Descents converged to Fmax < 1000 in 32 steps
Potential Energy  = -3.0919205e+06
Maximum force    =  9.1417908e+02 on atom 122811
Norm of force    =  1.1197423e+02
這樣的結果應該是有問題的吧，只跑了32步。

雖然只跑了32步，不過你的力已經(jīng)基本收斂了，應該沒有什么問題了。應該不用freegroup,如果要限制原子的話，可以嘗試使用define.

贊一下(2人)

回復此樓

2樓2011-04-18 18:42:24

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

Jfreda

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時
蟲號: 870275

引用回帖:

Originally posted by jawang at 2011-04-18 18:42:24:
Steepest Descents converged to Fmax < 1000 in 32 steps
Potential Energy  = -3.0919205e+06
Maximum force    =  9.1417908e+02 on atom 122811
Norm of force    =  1.1197423e+02
這樣的結果應該是 ...

好的，謝謝，但是我想問一下freegroup和define有什么區(qū)別嗎？我的目的是在優(yōu)化的過程中逐步放開限制。

贊一下

回復此樓

3樓2011-04-18 19:10:46

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

jawang

木蟲 (正式寫手)

應助: 3 (幼兒園)
金幣: 3533.7
帖子: 376
在線: 214.7小時
蟲號: 838608

★ ★ ★
小木蟲(金幣+0.5):給個紅包，謝謝回帖
御劍江湖(金幣+2): 謝謝 2011-04-19 10:02:25
Jfreda(金幣+5): 2011-04-19 13:29:55

定義了freegrp就是讓這個原子徹底不動了（當然它可以定義不同方向），而使用define是使這個原子多了一個額外的力，使其保持在原來的位置�？茨阆胍_到的效果吧，如果使用freegrp本人覺得會需要你定義的這些原子位置已經(jīng)平衡了，要不然的話，你還是讓人家動動，調整下距離比較好點，手冊上說對壓力耦合還有影響

贊一下(2人)

回復此樓

4樓2011-04-19 09:07:48

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

Jfreda

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時
蟲號: 870275

怎么又沒人回答了？

贊一下

回復此樓

5樓2011-04-19 10:53:16

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

Jfreda

鐵蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時
蟲號: 870275

引用回帖:

Originally posted by jawang at 2011-04-19 09:07:48:
定義了freegrp就是讓這個原子徹底不動了（當然它可以定義不同方向），而使用define是使這個原子多了一個額外的力，使其保持在原來的位置�？茨阆胍_到的效果吧，如果使用freegrp本人覺得會需要你定義的這些原子位 ...

好的，非常感謝！

回復此樓

6樓2011-04-19 13:30:13

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 085602 化工專碩 338分求調劑 +10	路癡小琪 2026-03-27	10/500	2026-03-28 03:36 by fmesaito
[考研] 070300求調劑306分 +3	26要上岸 2026-03-27	3/150	2026-03-27 17:57 by arrow8852
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調劑 +10	幸運的醬醬 2026-03-22	12/600	2026-03-27 16:08 by muchong357
[考研] 085600材料與化工調劑 +10	A-哆啦Z夢 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[考研] 考研化學308分求調劑 +10	你好明天你好 2026-03-23	12/600	2026-03-27 14:43 by shangxh
[考研] 一志愿太原理工安全工程300分，求調劑 +4	0857求調劑. 2026-03-24	4/200	2026-03-27 14:03 by 靜靜靜靜靜靜靜?/a>
[考研] 316求調劑 +5	Pigcasso 2026-03-24	5/250	2026-03-27 12:10 by zhshch
[考研] 321求調劑 +4	璞玉~~ 2026-03-25	4/200	2026-03-27 11:41 by sanrepian
[考研] 考研調劑 +9	小蠟新筆 2026-03-26	9/450	2026-03-27 11:10 by 不吃魚的貓
[考研] 求調劑一志愿本科北科大化學 343 +6	13831862839 2026-03-24	7/350	2026-03-26 22:57 by 不吃魚的貓
[考研] 349求調劑 +4	李木子啊哈哈 2026-03-25	4/200	2026-03-26 22:49 by fmesaito
[考研] 281求調劑 +3	亞克西good 2026-03-26	5/250	2026-03-26 19:48 by 不吃魚的貓
[考研] 材料考研求調劑 +3	Dendel 2026-03-23	6/300	2026-03-26 17:51 by fmesaito
[考研] 085601求調劑總分293英一數(shù)二 +4	鋼鐵大炮 2026-03-24	4/200	2026-03-26 16:28 by dick_runner
[考研] 材料與化工328分調劑 +6	。，。，。，。i 2026-03-23	6/300	2026-03-25 22:30 by 418490947
[考研] 機械學碩總分317求調劑�。。�！ +4	Acaciad 2026-03-25	4/200	2026-03-25 19:59 by hanserlol
[考研] 293求調劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] B區(qū)考研調劑 +4	yqdszhdap－ 2026-03-22	5/250	2026-03-25 08:51 by baoball
[考研] 300分，材料，求調劑，英一數(shù)二 +5	超贊的 2026-03-24	5/250	2026-03-24 21:07 by 星空星月
[考研] 307求調劑 +3	余意卿 2026-03-21	6/300	2026-03-24 15:03 by 余意卿