计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 分子模拟 (Molecular Modeling) › 请教一下伞形采样窗口覆盖不全问题

请教一下伞形采样窗口覆盖不全问题

查看数: 1198 | 评论数: 6 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

liuxdhs

发布时间: 2023-4-15 22:20

正文摘要:

我使用gmx模拟两条互锁的分子链沿着y轴拉伸过程，并想进行伞形采样计算其PMF。两条分子链被拉伸分别为上部与底部的碳原子。互锁分子在拉伸过程明显可以看出来构象发生明显的变化，如图是smd的mdp与伞形采样过程 ...

回复 Reply

liuxdhs 发表于 Post on 2023-4-16 22:07:49

Entropy.S.I 发表于 2023-4-16 21:42
你该不会给要计算PMF的分子加限制势（define = -DPOSRES）了吧？？

您好，我可能知道什么原因了。我freeze错对象了，我重新跑一次试一下。非常感谢。

liuxdhs 发表于 Post on 2023-4-16 21:51:01

Entropy.S.I 发表于 2023-4-16 21:42
你该不会给要计算PMF的分子加限制势（define = -DPOSRES）了吧？？

我没有加入这个条件。我是使用smd产生的轨迹做PMF
后续做了一个500ps的eq，将top跟btm两组拉伸的原子固定，rate设置为0。再做了一个3nm的采样md。这两个md都在nvt下计算。

同时，我查看了不同frame计算出来的pullf.xvg，都是定值，不知道这个会不会影响呢？(不同frame是定值，但不相同)

Entropy.S.I 发表于 Post on 2023-4-16 21:42:26

liuxdhs 发表于 2023-4-16 21:06
您好，我将k从1500减小到100甚至是10，还是出现这种情况？请问这是什么问题呢？
同时，将k减少要将跑的 ...

你该不会给要计算PMF的分子加限制势（define = -DPOSRES）了吧？？

liuxdhs 发表于 Post on 2023-4-16 21:06:50

Entropy.S.I 发表于 2023-4-15 22:35
减小k

您好，我将k从1500减小到100甚至是10，还是出现这种情况？请问这是什么问题呢？
同时，将k减少要将跑的时间延长才可以拉伸到当k为1500时的距离。如下图。

Entropy.S.I 发表于 Post on 2023-4-15 22:35:31

减小k

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register