chands 发表于 2022-9-16 16:23
你的三幅图是同一个MD跑出来的吗？第二幅是1000 fs，排在第三幅700 fs后？感觉第二幅图只用了1000 fs附件 ...

是一个md跑的 700fs是一相对完整fes 到后期就变成两截了。还是说因为跑到后期自由能的值超出了色带的范围 引起了这样情况

wwp 发表于 2022-9-16 14:56
MD全过程是没有问题的，可以观察到明显的成断键过程，也符合预期。第三幅图就是在700fs左右的结果，就是 ...

你的三幅图是同一个MD跑出来的吗？第二幅是1000 fs，排在第三幅700 fs后？感觉第二幅图只用了1000 fs附件一截MD计算自由能。

wwp 发表于 2022-9-15 18:28
我这里hillspawn用的参数是 30 0.1 0.1 高斯峰取值很小时，能量的变化就感觉很小。这样他就是采样的时间 ...

第一幅图显然没有充分采样。第二幅图感觉有点怪啊，CV1、CV2应该是连续变化的，为什么图上是两截，你查了MD的全过程了吗？第三幅图应该是最好的，我不知道“收敛没有”，不过可以试一下把Hillspawn的频率调低（即第一个参数调高），或者hill的高度调低，这样采样次数会多一点，看看结果是否一致。

wwp 发表于 2022-9-15 14:26
可是在模拟的前期，我在一些别的地方也能观察到势能的变化，跑到后期就全部变白了，这是怎么回事呢？麻烦 ...

你这个图采样区域的自由能值非常小，在0附件波动，可能有很小的负值，怀疑FES图以最小值为零点，导致没采样区域变成很小的正值。感觉在模拟前期反应尚未发生，如果反应发生了，一个CV会变小同时另一个CV变大。

wwp 发表于 2022-9-15 14:26
可是在模拟的前期，我在一些别的地方也能观察到势能的变化，跑到后期就全部变白了，这是怎么回事呢？麻烦 ...

这个我不知道，是什么体系，你的CV1和CV2是什么东西，还有填色图的标度是什么，体系是怎么跑的？

wwp 发表于 2022-9-15 14:17
嗯嗯，这里的“收敛”这样理解就好理解了。还有一个问题就是：加到把自由能的"坑"填满 这句话。 那这个时 ...

这个大片空白大概率是没采集到（或者正常情况下你的collective variable就不会跑到这些地方），FES是通过加的偏置势算的，如果自由能的“坑”满了，那自由能就是所加的总偏置势的负值，顶多加一个无关紧要的常数。

这时“收敛”不应该理解成体系处于能量极小点（然后就不乱跑了）。一种算自由能的方法是加偏置势一直加到把自由能的"坑"填满，这时就没有“坑”了， 在这个维度上系统当然会到处跑，因为势能剖面是平的了，没有什么限制了。

chands 发表于 2022-9-15 13:23
我的理解是如果你做了1.5 ps和1 ps的结果类似，那就是“收敛”了。

我也这样觉得，就是感觉可能能不能有类似能量收敛那种的硬性标准。还有一个问题就是，在做两个分子化学反应的动力学，整个体系的能量的收敛是不是比较难达到。我的理解是在因为孤立的分子在整个空间里一直在（或者加了限制势）乱跑。

我的理解是如果你做了1.5 ps和1 ps的结果类似，那就是“收敛”了。