求助MD模拟计算吸附过程中焓变及熵变

Tifo

请问各位老师，我准备使用MD模拟的方法计算离子吸附过程焓变及熵变，已经通过伞状抽样获得离子的PMF曲线，据我查找到的文献上有两种方法计算熵和焓，一种是有限差分法，另一种是将系统势能的平均值看作系统的焓，因为我想得到焓的分解项（电势能、范德华势），因此考虑使用第二种方法但我看到文献中利用系统平均势能计算焓的时候，体系中一般都是只有一个离子或一个肽链，将其从平衡位置向外移动到体相中，但离子在每个位置中取样时间以及受到的约束刚度一般多大都没有提到，而且我的系统里吸附离子有多个，如果我只移动其中一个离子这样获得的结果是否正确？还是说这种方法只能用于只有一个离子的情况
请各位老师不吝赐教

ridermarx

请问通过势能计算焓的方法是哪篇文献呢？最近也在做这方面的工作。谢谢