24小时热门版块排行榜    

查看: 1319  |  回复: 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

金虫 (正式写手)

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

奋斗1s

金虫 (正式写手)

怎么没有回复呢,自己先顶一下
2楼2011-05-18 09:47:02
已阅   回复此楼   关注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的回帖

奋斗1s

金虫 (正式写手)

没有人回复了?怎么老是提醒我有人回复,每次兴匆匆的跑来看,却没有……
6楼2011-05-22 23:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见