24小时热门版块排行榜    

查看: 1639  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

奋斗1s

金虫 (正式写手)

[求助] Gromacs做steep之后做cg出错

我在做最陡下降法优化mdrun之后想进行一次共轭梯度优化,出现了下面错误:
ERROR: can not do Conjugate Gradients with constraints (14556)
而我的em_cg.mdp中写了constrains = none,请问高手们这个错误时怎么回事呀?

ps:我的输入命令如下:
grompp -f em.mdp -c 1omb_b4em.gro -p 1omb.top -o 1omb_em.tpr

nohup mdrun -s 1omb_em.tpr -o 1omb_em.trr -c 1omb_b4pr.gro -g em.log -e em.edr &

grompp -f em_cg.mdp -c 1omb_b4em.gro -p 1omb.top -o 1omb_em.tpr

nohup mdrun -s 1omb_em.tpr -o 1omb_em.trr -c 1omb_b4pr.gro -g em.log -e em.edr &

我的em.mdp em_cg.mdp分别如下:
em.mdp:
title               =  fws
cpp                 =  /usr/bin/cpp
define              =  -DFLEXIBLE
constraints         =  none
integrator          =  steep
dt                  =  0.002    ; ps !
nsteps              =  500
nstlist             =  10
ns_type             =  grid
rlist               =  1.0
coulombtype         =  PME
rcoulomb            =  1.0
vdwtype             =  cut-off
rvdw                =  1.4
fourierspacing                =  0.12
fourier_nx                =  0
fourier_ny                =  0
fourier_nz                =  0
pme_order                =  4
ewald_rtol                =  1e-5
optimize_fft                =  yes
emtol               =  1000.0
emstep              =  0.01

em_cg.mdp:
title               =  fws
cpp                 =  /usr/bin/cpp
define              =  -DFLEXIBLE
constraints         =  none
integrator          =  cg
dt                  =  0.002    ; ps !
nsteps              =  3000
nstlist             =  10
ns_type             =  grid
rlist               =  1.0
coulombtype         =  PME
rcoulomb            =  1.0
vdwtype             =  cut-off
rvdw                =  1.4
fourierspacing                =  0.12
fourier_nx                =  0
fourier_ny                =  0
fourier_nz                =  0
pme_order                =  4
ewald_rtol                =  1e-5
optimize_fft                =  yes
emtol               =  1000.0
emstep              =  0.01
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋斗1s

金虫 (正式写手)

没有人回复了?怎么老是提醒我有人回复,每次兴匆匆的跑来看,却没有……
6楼2011-05-22 23:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

奋斗1s

金虫 (正式写手)

怎么没有回复呢,自己先顶一下
2楼2011-05-18 09:47:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

奋斗1s

金虫 (正式写手)

怎么没有人回复啊
4楼2011-05-18 19:21:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fantastic

银虫 (著名写手)


【答案】应助回帖

★ ★
奋斗1s(金币+5): thanks a lot 2011-05-21 13:19:56
ghcacj(金币+2): 谢谢 2011-05-21 17:07:20
奋斗1s(金币+1): 谢谢,只做steep会不会有什么不好的影响呢? 2011-05-22 12:09:17
貌似amber也有这个问题,优化做完了steep,再做CG就是会有问题,可能gromacs也有这个问题吧,你可以全部用steep,收敛还是比较容易的
5楼2011-05-21 11:48:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:14 by l0VvVHGBGRLv
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 11:52 by l0VvVHGBGRLv
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 06:11 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 06:10 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见