24小时热门版块排行榜    

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

金虫 (正式写手)

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

奋斗1s

金虫 (正式写手)

再顶一下
3楼2011-05-18 09:47:14
已阅   回复此楼   关注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的回帖

奋斗1s

金虫 (正式写手)

没有人回复了?怎么老是提醒我有人回复,每次兴匆匆的跑来看,却没有……
6楼2011-05-22 23:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 奋斗1s 的主题更新
信息提示
请填处理意见