计算化学公社

 找回密码 Forget password
 注册 Register

求助:设置用显卡加速PME计算

查看数: 4852 | 评论数: 5 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-3-8 13:49

正文摘要:

本帖最后由 botu 于 2022-3-8 13:49 编辑 使用CHARMM-GUI生成的膜-蛋白输入文件,想在gmx mdrun命令中设置PME计算也用显卡加速,但是添加上-pme gpu之后,运行会报如下错误: Feature not implemented: Cannot ...

回复 Reply

botu 发表于 Post on 2022-3-12 20:42:27
喵星大佬 发表于 2022-3-8 16:17
你在做能量极小化?
只有md支持pme gpu

是这一步不应该设置pme gpu……去掉之后后面就正常了
botu 发表于 Post on 2022-3-10 20:50:42
喵星大佬 发表于 2022-3-8 16:17
你在做能量极小化?
只有md支持pme gpu

啊是的,极小化那一步就卡住了
(谢谢回复!)
botu 发表于 Post on 2022-3-10 20:49:30

老师好,谢谢回复!!  mdp如下:
1)step6.0_minimization.mdp
define                  = -DPOSRES -DPOSRES_FC_BB=4000.0 -DPOSRES_FC_SC=2000.0 -DPOSRES_FC_LIPID=1000.0 -DDIHRES -DDIHRES_FC=1000.0
integrator              = steep
emtol                   = 1000.0
nsteps                  = 5000
nstlist                 = 10
cutoff-scheme           = Verlet
rlist                   = 1.2
vdwtype                 = Cut-off
vdw-modifier            = Force-switch
rvdw_switch             = 1.0
rvdw                    = 1.2
coulombtype             = PME
rcoulomb                = 1.2
;
constraints             = h-bonds
constraint_algorithm    = LINCS

2)step6.6_equilibration.mdp
define                  = -DPOSRES -DPOSRES_FC_BB=50.0 -DPOSRES_FC_SC=0.0 -DPOSRES_FC_LIPID=0.0 -DDIHRES -DDIHRES_FC=0.0
integrator              = md
dt                      = 0.002
nsteps                  = 250000
nstxtcout               = 5000
nstvout                 = 5000
nstfout                 = 5000
nstcalcenergy           = 100
nstenergy               = 1000
nstlog                  = 1000
;
cutoff-scheme           = Verlet
nstlist                 = 20
rlist                   = 1.2
vdwtype                 = Cut-off
vdw-modifier            = Force-switch
rvdw_switch             = 1.0
rvdw                    = 1.2
coulombtype             = PME
rcoulomb                = 1.2
;
tcoupl                  = berendsen
tc_grps                 = SOLU MEMB SOLV
tau_t                   = 1.0 1.0 1.0
ref_t                   = 294.15 294.15 294.15
;
pcoupl                  = berendsen
pcoupltype              = semiisotropic
tau_p                   = 5.0
compressibility         = 4.5e-5  4.5e-5
ref_p                   = 1.0     1.0
refcoord_scaling        = com
;
constraints             = h-bonds
constraint_algorithm    = LINCS
continuation            = yes
;
nstcomm                 = 100
comm_mode               = linear
comm_grps               = SOLU_MEMB SOLV

3)step7_production.mdp
integrator              = md
dt                      = 0.002
nsteps                  = 500000
nstxout                 = 50000
nstvout                 = 50000
nstfout                 = 50000
nstcalcenergy           = 100
nstenergy               = 1000
nstlog                  = 1000
;
cutoff-scheme           = Verlet
nstlist                 = 20
rlist                   = 1.2
vdwtype                 = Cut-off
vdw-modifier            = Force-switch
rvdw_switch             = 1.0
rvdw                    = 1.2
coulombtype             = PME
rcoulomb                = 1.2
;
tcoupl                  = Nose-Hoover
tc_grps                 = SOLU MEMB SOLV
tau_t                   = 1.0 1.0 1.0
ref_t                   = 294.15 294.15 294.15
;
pcoupl                  = Parrinello-Rahman
pcoupltype              = semiisotropic
tau_p                   = 5.0
compressibility         = 4.5e-5  4.5e-5
ref_p                   = 1.0     1.0
;
constraints             = h-bonds
constraint_algorithm    = LINCS
continuation            = yes
;
nstcomm                 = 100
comm_mode               = linear
comm_grps               = SOLU_MEMB SOLV
sobereva 发表于 Post on 2022-3-8 22:10:50
把mdp贴出来
喵星大佬 发表于 Post on 2022-3-8 16:17:31
你在做能量极小化?
只有md支持pme gpu

评分 Rate

参与人数
Participants 2
eV +6 收起 理由
Reason
botu + 2 谢谢
muuu2333 + 4

查看全部评分 View all ratings

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 14:59 , Processed in 0.197112 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list