24小时热门版块排行榜    

查看: 1331  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +4 Liwangman 2026-03-15 4/200 2026-03-16 08:01 by wang_dand
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
信息提示
请填处理意见