24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1353  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂:085600材料与化工,考材科基,总分319 +17 678lucky 2026-03-31 21/1050 2026-04-01 01:40 by 1018329917
[考研] 085602 307分 求调剂 +10 不知道叫什么! 2026-03-26 10/500 2026-03-31 19:53 by Dyhoer
[考研] 省双一流重点一本大学招收调剂 +3 wwwwffffff 2026-03-31 5/250 2026-03-31 19:49 by 曾曾曾!
[考研] 材料科学与工程339求调剂 +3 hyz0119 2026-03-31 3/150 2026-03-31 18:31 by 1939136013狗壮
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 085600 材料与化工 329分求调剂 +20 Mr. Z 2026-03-25 21/1050 2026-03-31 16:53 by Zzxxxs
[考研] 理学07化学 303求调剂 +8 睿08 2026-03-27 8/400 2026-03-31 16:31 by 690616278
[考研] 材料与化工272求调剂 +25 阿斯蒂芬2004 2026-03-28 25/1250 2026-03-31 16:27 by hypershenger
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +12 nnnnnnn5 2026-03-25 12/600 2026-03-31 16:21 by 记事本2026
[考研] 求收留 +8 1943443204 2026-03-28 8/400 2026-03-31 15:00 by -迷了路啊路
[考研] 085600材料与化工调剂 +16 kikiki7 2026-03-30 16/800 2026-03-31 10:03 by 氯化亚硝酰
[考研] 一志愿大连理工大学材料求调剂 +6 Gymno 2026-03-30 6/300 2026-03-31 07:26 by 无际的草原
[考研] 359求调剂 +5 王了个楠 2026-03-25 5/250 2026-03-30 19:36 by 源_2020
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 295求调剂 +5 wei-5 2026-03-26 5/250 2026-03-30 08:34 by 探123
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 352分 化工与材料 +5 海纳百川Ly 2026-03-27 5/250 2026-03-28 03:39 by fmesaito
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
信息提示
请填处理意见