Giovannelli Edoardo, Cardini Gianni, Chelli Riccardo
Dipartimento di Chimica, Università di Firenze, Via della Lastruccia 3, I-50019 Sesto Fiorentino, Italy.
Phys Rev E Stat Nonlin Soft Matter Phys. 2015 Oct;92(4):043310. doi: 10.1103/PhysRevE.92.043310. Epub 2015 Oct 30.
Generalized-ensemble simulations, such as replica exchange and serial generalized-ensemble methods, are powerful simulation tools to enhance sampling of free energy landscapes in systems with high energy barriers. In these methods, sampling is enhanced through instantaneous transitions of replicas, i.e., copies of the system, between different ensembles characterized by some control parameter associated with thermodynamical variables (e.g., temperature or pressure) or collective mechanical variables (e.g., interatomic distances or torsional angles). An interesting evolution of these methodologies has been proposed by replacing the conventional instantaneous (trial) switches of replicas with noninstantaneous switches, realized by varying the control parameter in a finite time and accepting the final replica configuration with a Metropolis-like criterion based on the Crooks nonequilibrium work (CNW) theorem. Here we revise these techniques focusing on their correlation with the CNW theorem in the framework of Markovian processes. An outcome of this report is the derivation of the acceptance probability for noninstantaneous switches in serial generalized-ensemble simulations, where we show that explicit knowledge of the time dependence of the weight factors entering such simulations is not necessary. A generalized relationship of the CNW theorem is also provided in terms of the underlying equilibrium probability distribution at a fixed control parameter. Illustrative calculations on a toy model are performed with serial generalized-ensemble simulations, especially focusing on the different behavior of instantaneous and noninstantaneous replica transition schemes.
广义系综模拟,如副本交换和串行广义系综方法,是强大的模拟工具,可增强具有高能垒系统中自由能景观的采样。在这些方法中,通过副本(即系统的拷贝)在不同系综之间的瞬时转变来增强采样,这些系综由与热力学变量(如温度或压力)或集体力学变量(如原子间距离或扭转角)相关的某个控制参数来表征。通过用非瞬时转变取代传统的副本瞬时(尝试)切换,提出了这些方法的一种有趣演变,非瞬时转变是通过在有限时间内改变控制参数并基于克鲁克斯非平衡功(CNW)定理用类似 metropolis 准则接受最终副本构型来实现的。在这里,我们在马尔可夫过程的框架内,着重于这些技术与 CNW 定理的相关性来修订这些技术。本报告的一个成果是推导了串行广义系综模拟中非瞬时切换的接受概率,我们表明在此类模拟中,无需明确了解进入模拟的权重因子的时间依赖性。还根据固定控制参数下的潜在平衡概率分布提供了 CNW 定理的广义关系。使用串行广义系综模拟对一个玩具模型进行了说明性计算,特别关注瞬时和非瞬时副本转变方案的不同行为。