Jfreda

應(yīng)助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時(shí)
蟲(chóng)號(hào): 870275

[交流] mdrun時(shí)出錯(cuò)

mdrun -s aa.tpr -o aa.trr -c aa.gro &
出現(xiàn)了如下的提示錯(cuò)誤：
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ù)(只對(duì)一個(gè)參數(shù)做了修改)：
1.em.mdp文件中的freezegrps改為Protein UNK能量能很好地收斂，下面是它的顯示結(jié)果：
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（參照別人的回答), 結(jié)果顯示為：
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
這樣的結(jié)果應(yīng)該是有問(wèn)題的吧，只跑了32步。

3.integrator改為md, 結(jié)果為：
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)

請(qǐng)問(wèn)一下，我應(yīng)該怎么做才是合理的？謝謝！

回復(fù)此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有12人回復(fù)
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！已經(jīng)有3人回復(fù)
329求調(diào)劑已經(jīng)有5人回復(fù)
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù) 已經(jīng)有4人回復(fù)
材料學(xué)碩318求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑已經(jīng)有6人回復(fù)
081700化工學(xué)碩調(diào)劑已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

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

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

查看全部散金貼

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

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

jawang

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 3 (幼兒園)
金幣: 3533.7
帖子: 376
在線: 214.7小時(shí)
蟲(chóng)號(hào): 838608

★ ★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖
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
這樣的結(jié)果應(yīng)該是有問(wèn)題的吧，只跑了32步。

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

贊一下(2人)

回復(fù)此樓

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

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

Jfreda

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時(shí)
蟲(chóng)號(hào): 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
這樣的結(jié)果應(yīng)該是 ...

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

贊一下

回復(fù)此樓

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

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

jawang

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 3 (幼兒園)
金幣: 3533.7
帖子: 376
在線: 214.7小時(shí)
蟲(chóng)號(hào): 838608

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

定義了freegrp就是讓這個(gè)原子徹底不動(dòng)了（當(dāng)然它可以定義不同方向），而使用define是使這個(gè)原子多了一個(gè)額外的力，使其保持在原來(lái)的位置。看你想要達(dá)到的效果吧，如果使用freegrp本人覺(jué)得會(huì)需要你定義的這些原子位置已經(jīng)平衡了，要不然的話，你還是讓人家動(dòng)動(dòng)，調(diào)整下距離比較好點(diǎn)，手冊(cè)上說(shuō)對(duì)壓力耦合還有影響

贊一下(2人)

回復(fù)此樓

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

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

Jfreda

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時(shí)
蟲(chóng)號(hào): 870275

怎么又沒(méi)人回答了？

贊一下

回復(fù)此樓

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

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

Jfreda

鐵蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 21.8
帖子: 91
在線: 31.5小時(shí)
蟲(chóng)號(hào): 870275

引用回帖:

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

好的，非常感謝！

回復(fù)此樓

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

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

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

[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	3/150	2026-03-19 23:58 by 23Postgrad
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +18	rare12345 2026-03-18	18/900	2026-03-19 22:36 by 學(xué)員8dgXkO
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +5	Doleres 2026-03-19	5/250	2026-03-19 20:14 by 制度的
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3	葵梓衛(wèi)隊(duì) 2026-03-18	5/250	2026-03-19 19:35 by 給你你注意休息
[考研] 321求調(diào)劑 +8	何潤(rùn)采123 2026-03-18	10/500	2026-03-19 16:46 by 何潤(rùn)采123
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +3	@taotao 2026-03-19	3/150	2026-03-19 14:07 by peike
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 304求調(diào)劑 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 工科材料085601 279求調(diào)劑 +6	困于星晨 2026-03-17	6/300	2026-03-18 10:21 by kkcoco25
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么