White Ronald P, Meirovitch Hagai
Department of Computational Biology, University of Pittsburgh School of Medicine, Biomedical Science Tower, Pittsburgh, Pennsylvania 15261, USA.
J Chem Phys. 2005 Dec 1;123(21):214908. doi: 10.1063/1.2132285.
Hypothetical scanning Monte Carlo (HSMC) is a method for calculating the absolute entropy S and free energy F from a given MC trajectory developed recently and applied to liquid argon, TIP3P water, and peptides. In this paper HSMC is extended to random coil polymers by applying it to self-avoiding walks on a square lattice--a simple but difficult model due to strong excluded volume interactions. With HSMC the probability of a given chain is obtained as a product of transition probabilities calculated for each bond by MC simulations and a counting formula. This probability is exact in the sense that it is based on all the interactions of the system and the only approximation is due to finite sampling. The method provides rigorous upper and lower bounds for F, which can be obtained from a very small sample and even from a single chain conformation. HSMC is independent of existing techniques and thus constitutes an independent research tool. The HSMC results are compared to those obtained by other methods, and its application to complex lattice chain models is discussed; we emphasize its ability to treat any type of boundary conditions for which a reference state (with known free energy) might be difficult to define for a thermodynamic integration process. Finally, we stress that the capability of HSMC to extract the absolute entropy from a given sample is important for studying relaxation processes, such as protein folding.
假设扫描蒙特卡罗(HSMC)是一种最近开发的从给定的蒙特卡罗轨迹计算绝对熵S和自由能F的方法,并已应用于液态氩、TIP3P水和肽。在本文中,通过将HSMC应用于正方形晶格上的自回避行走,将其扩展到无规卷曲聚合物——这是一个由于强排斥体积相互作用而简单但困难的模型。利用HSMC,给定链的概率通过蒙特卡罗模拟计算每个键的转移概率和一个计数公式得到。这种概率在基于系统的所有相互作用的意义上是精确的,唯一的近似是由于有限采样。该方法为F提供了严格的上下界,可从非常小的样本甚至单个链构象中获得。HSMC独立于现有技术,因此构成了一个独立的研究工具。将HSMC的结果与其他方法获得的结果进行了比较,并讨论了其在复杂晶格链模型中的应用;我们强调了它处理任何类型边界条件的能力,对于这些边界条件,在热力学积分过程中可能难以定义参考状态(具有已知自由能)。最后,我们强调HSMC从给定样本中提取绝对熵的能力对于研究诸如蛋白质折叠等弛豫过程很重要。