We introduce a generalised micro-macro Markov chain Monte Carlo (mM-MCMC) method with pseudo-marginal approximation to the free energy, that is able to accelerate sampling of the microscopic Gibbs distributions when there is a time-scale separation between the macroscopic dynamics of a reaction coordinate and the remaining microscopic degrees of freedom. The mM-MCMC method attains this efficiency by iterating four steps: i) Propose a new value of the reaction coordinate; ii) Accept or reject the macroscopic sample; iii) Run a biased simulation that creates a microscopic molecular instance that lies close to the newly sampled macroscopic reaction coordinate value; iv) Microscopic accept/reject step for the new microscopic sample. In the present paper, we eliminate the main computational bottleneck of earlier versions of this method: the necessity to have an accurate approximation of the free energy. We show that introduction of a pseudo-marginal approximation significantly reduces the computational cost of the microscopic accept/reject step, while still providing unbiased samples. We illustrate the method's behaviour on several molecular systems with low-dimensional reaction coordinates.
翻译:Translated Abstract:
我们引入了一种广义的微观-宏观马尔科夫链蒙特卡罗(mM-MCMC)方法,并采用伪边缘近似方法逼近自由能。当反应坐标的宏观动力学和剩余微观自由度之间存在时间尺度分离时,该方法能够加速对微观吉布斯分布的采样。mM-MCMC方法通过迭代四个步骤来实现这种效率:i)提出反应坐标的新值; ii)接受或拒绝宏观样本; iii)运行一个有偏置的模拟,创建一个微观分子实例,该实例接近新采样的宏观反应坐标值; iv)新微观样本的微观接受/拒绝步骤。在本文中,我们消除了该方法早期版本的主要计算瓶颈:准确逼近自由能的必要性。我们表明,伪边缘近似的引入显著降低了微观接受/拒绝步骤的计算成本,同时仍然提供无偏的样本。我们在具有低维反应坐标的几个分子系统上说明了该方法的行为。