您好同学，请问你最后用什么办法解决了反应坐标描述这个运动的呢？

如果计算资源充足的话，可以考虑进行大量unbiased MD，并构建马尔可夫模型来研究这两个loop区之间的距离变化

azero 发表于 2020-5-17 08:57
谢谢付老师的解答~
感觉这体系eABF跑很久都难以收敛

如果反应坐标选择的不好，不管哪种方法都不会收敛，只不过是表现形式不一样

用SMD，多半还是正着拉和反着拉的结果千差万别
用US的话结果和每个小窗口内的初始结构会有很大的关系，你给不同的初始结构结果完全不同

你这个问题的话更适合选择全空间增强采样方法（权重系综或者广义系综），比如aMD或者REMD，然后观察轨迹是否出现了你想看到的自发变化

fhh2626 发表于 2020-5-16 23:27
因为你的自由能计算没有收敛

或者换句话说，你两次模拟的采样是不相同的（都没有达到平衡），如果你仔细 ...

谢谢付老师的解答~
感觉这体系eABF跑很久都难以收敛

我只是想简单判断这两段loop区靠近和远离这两种情况能否自发发生。
1. 能否用SMD分别把两种情况分别跑一个模拟。
2. 或者用伞形采样？伞形采样跟SMD应该不是同一个东西，两者区别弄得不是很清楚。
3. 或者有没更适合的模拟方法？

因为你的自由能计算没有收敛

或者换句话说，你两次模拟的采样是不相同的（都没有达到平衡），如果你仔细观察两次模拟的轨迹的话，会发现它们在反应坐标值相同时采样到了不同的结构，也就是说你选择的反应坐标可能无法充分描述这个运动