计算化学公社

标题: 求助gromacs伞形采样数据缺失 [打印本页]

作者
Author: zhangzh 时间: 2020-7-11 00:17
标题: 求助gromacs伞形采样数据缺失
本帖最后由 zhangzh 于 2020-7-11 00:21 编辑

我在伞形采样过程中，用pull命令生成了一条轨迹，然后抽出了前110帧的结构：

gmx trjconv -f pull.xtc -s pull.tpr -o conf.gro -b 0 -e 110 -sep

复制代码

然后每5帧为一个窗口，对22个窗口进行100ps的npt平衡和300ps的md平衡。将平衡之后得到的pullx.xvg和tpr文件放到wham中分析：

gmx wham -it tpr.dat -if pullx.dat -o -hist -unit kJ

复制代码

最后得到的直方图和PMF曲线均有数据缺失，提示无法读取某些位置的数据。可是我检查了xvg文件中均是有数据的，请问是哪一步出错了呢？应该如何解决
附件中是错误的直方图，PMF曲线，输入文件和报错信息
(, 下载次数 Times of downloads: 29) (, 下载次数 Times of downloads: 31)
(, 下载次数 Times of downloads: 14) (, 下载次数 Times of downloads: 21)

(, 下载次数 Times of downloads: 29) (, 下载次数 Times of downloads: 23)

(, 下载次数 Times of downloads: 30) (, 下载次数 Times of downloads: 17)

作者
Author: zhangzh 时间: 2020-7-15 01:18
补充一下问题，我的模拟的时间，步长等参考环糊精体系又试了一遍，npt和md过程正常输出了数据，可是还是出现了同样的问题，为什么会有一段根本没有数据呢

作者
Author: wyb 时间: 2020-7-15 09:32
把所有输出数据在一起做出histgram图，然后看你的数据是否有问题

作者
Author: zhangzh 时间: 2020-7-17 00:39

wyb 发表于 2020-7-15 09:32
把所有输出数据在一起做出histgram图，然后看你的数据是否有问题

这个就是所有的输出一起做的

作者
Author: liuyuje714 时间: 2020-7-17 08:03
你的窗口选取就有问题，不是10 15 20 25这样选的，一般是等距离选取的，比如0.2 nm

作者
Author: wyb 时间: 2020-7-17 08:39
看看那几个windows的轨迹，是否配体在指定的位置，看你的输出文件，应该跑掉了，如果是跑掉了，考虑增大constraints以确保在该位置。

作者
Author: yuanen 时间: 2022-3-21 20:54

wyb 发表于 2020-7-17 08:39
看看那几个windows的轨迹，是否配体在指定的位置，看你的输出文件，应该跑掉了，如果是跑掉了，考虑增大con ...

我想请问，增大constraints是在楼主的npt_um.mdp和ummd.mdp中的pull_cood1_k这个参数吗？这个参数与pull.mdp要保持一致吗？我最近在学这方面的，还请劳烦指点

作者
Author: Frozen-Penguin 时间: 2022-3-22 02:02

zhangzh 发表于 2020-7-17 00:39
这个就是所有的输出一起做的

在xmgrace界面里，data->import->选择直方图xvg->读取格式改成Nxy，这样应该可以看到很多条线，直接用xmgrace打开只能看到一条线。

作者
Author: thor 时间: 2023-2-27 10:38
你这个问题最终解决了吗

作者
Author: 熊小熊 时间: 2023-4-16 18:37

yuanen 发表于 2022-3-21 20:54
我想请问，增大constraints是在楼主的npt_um.mdp和ummd.mdp中的pull_cood1_k这个参数吗？这个参数与pull. ...

您的这个问题解决了吗？

作者
Author: yuanen 时间: 2023-4-23 16:00

熊小熊发表于 2023-4-16 18:37
您的这个问题解决了吗？

可以通过增加pull_cood1_k的大小使运动组保持在你的采样窗口范围内，但是需要增加采样窗口以保证数据点的连续

欢迎光临计算化学公社 (http://bbs.keinsci.com/)