计算化学公社

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

[GROMACS] 为什么用gmx wham命令产生的两次伞形采样的边界不同

[复制链接 Copy URL]

75

帖子

0

威望

705

eV
积分
780

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
我用gromacs对体系的某个二面角进行伞形采样,pull code如下:
pull                    = yes
pull-ngroups            = 6
pull-ncoords            = 1
pull-group1-name        = 6
pull-group2-name        = 5
pull-group3-name        = 5
pull-group4-name        = 2
pull-group5-name        = 2
pull-group6-name        = 1
pull-coord1-type        = umbrella      
pull-coord1-geometry    = dihedral      
pull-coord1-init        = -150
pull-coord1-groups      = 1 2 3 4 5 6
pull-coord1-dim         = Y Y Y
pull-coord1-rate        = 0.00
pull-coord1-k           = 200         
pull-coord1-start             = no           
——————————————————————————————————————————————————————————————————————————
用gmx wham -it tpr-files.dat -if pullf-files.dat -o histo -unit kCal进行分析,输出的边界: Determined boundaries to -1629.018311 and 1506.088135,然而边界明明是-160度到-64度。
然而用pullx-files.dat 作为输入文件就对了,这是为什么?pullx.xvg和pullf.xvg两个文件不应该输出的结果是一样的吗?

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 07:13 , Processed in 0.152604 second(s), 20 queries , Gzip On.

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