计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4248|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求助restarint-Lambdas调节不了pull code限制

[复制链接 Copy URL]

5

帖子

0

威望

123

eV
积分
128

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
在跑蛋白质和配体对接的 fep 时,Lambda = 0 时的 vdw 和 库仑力都为0,结果配体会离开口袋。查到 restarint-Lambdas 可以像调节vdw 和库仑力一样调节 pull code, 就在mdp里写了 pull code 限制口袋(pocket)和配体(lig_heavy)的质心距离不变:
        pull                         = yes
        pull_ncoords            = 1         
        pull_ngroups            = 2         
        pull_group1_name        = pocket   
        pull_group2_name        = lig_heavy
        pull_coord1_type        = umbrella
        pull_coord1_geometry    = distance  
        pull_coord1_dim         = Y Y Y
        pull_coord1_groups      = 1 2
        pull_coord1_start       = yes     
        pull_coord1_rate        = 0     
        pull_coord1_k           = 1000     


然后分别设置:

restraint-Lambdas = 1.0 1.0 1.0 .... 1.0

restraint-Lambdas = 0.0 0.0 0.0 .... 0.0
都只跑了Lambda=0的fep,结果ligand都没离开口袋。

请问:
1.正常情况下restraint-Lambdas = 0 或 1 是否表示pull code 生效 或 不生效?
2.是否需要其他关键词才能让restraint-Lambdas有效调节pull code?

本版积分规则 Credits rule

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

GMT+8, 2026-2-23 00:13 , Processed in 0.150853 second(s), 20 queries , Gzip On.

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