计算得到的吸附自由能比吸附能弱，请问是否正常？

k64_cc

a9471163 发表于 2020-11-4 18:44
谢谢您的提醒~
请问PMF曲线上两个state之间的差值如何求出来，我目前想的办法是把state1和state2的分子 ...

不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential Energy就行。另外，水和溶质的相互作用能也会变化，中和一部分AB间作用能的变化，所以只分析这AB间相互作用也不是很充分。

饭得一口一口吃，你得先确定H主导了自由能变化，再分析什么因素主导了H变化。对于<H>的计算，你就在CV曲线上选两个点，用非常强的restraint把CV控制在那个点附近，然后跑三组长时间的模拟，算每组的平均势能，并用三组数据算variance。
Umbrella Sampling的原理是打格子数sample。你可以用一个模拟给多个格子提供数据。但是算PMF上某点平均势能的时候，你需要要尽量让sample处于你选定的那个格子里，所以要用一个很强的restraint来避免CV的波动。也是由于同个原因，US时的数据不能拿来直接用。

a9471163

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

我明白了，您讲得很清楚，非常感谢老师

a9471163

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

我想用您说的这个方法计算熵变，请问有没有用过这种方法的参考文献可以引用

k64_cc

a9471163 发表于 2020-11-6 15:36
我想用您说的这个方法计算熵变，请问有没有用过这种方法的参考文献可以引用

随便整一篇我老板从前的工作：

https://pubs.acs.org/doi/abs/10.1021/jacs.5b07232

Tifo

k64_cc 发表于 2020-11-5 11:58
不行，能量收敛比PMF慢，别想rerun了。而且G=H-TS的H是总势能，不是分解的，不用定义组，看总的Potential ...

老师您好，我也用和您说的一样的方法计算系统的总势能，就是在多个吸附离子中选择一个使用restraint将其控制在CV路径上的某一点3ns，再将获得的系统势能平均，但这样最后获得的势能随距离的变化曲线波动很大，与该离子的自由能变化对不上，请问是还有其他措施我没做吗？

k64_cc

Tifo 发表于 2022-1-10 21:10
老师您好，我也用和您说的一样的方法计算系统的总势能，就是在多个吸附离子中选择一个使用restraint将其 ...

Variance大是一码事，和自由能曲线对不上是另一码事。

3 ns模拟，variance大是很正常的。一般来说，平均势能的收敛需要更长的时间。我之前做水的PMF，已经很好弛豫了，还是需要几十ns的轨迹才能确保它收敛。如果收敛之后势能曲线还是和自由能对不上，那就说明是熵驱动的，不是焓驱动。这俩玩意也不一定非要对上。

Tifo

k64_cc 发表于 2022-1-12 13:17
Variance大是一码事，和自由能曲线对不上是另一码事。

3 ns模拟，variance大是很正常的。一般来说，平 ...

感谢老师回复，我还有个问题想请教一下。就是我现在以一个离子为探针，把它从吸附稳定位置拉到体相位置，那等这个探针离子到体相位置的时候，其他离子难道不会重新吸附到之前探针离子所在的位点上吗，那这时候这个系统的总势能为什么还会和我拉那个探针之前的时候不一样？