Abstract:Most electrochemical reactions should be studied under a grand canonical ensemble condition with a constant potential and/or a constant pH value. Free energy profiles provide key insights into understanding the reaction mechanisms. However, many molecular dynamics (MD)-based theoretical studies for electrochemical reactions did not employ an exact grand canonical ensemble sampling scheme for the free energy calculations, partially due to the issues of discontinuous trajectories induced by the particle-number variations during MD simulations. An alternative statistical sampling approach, the Monte Carlo (MC) method, is naturally appropriate for the open-system simulations if we focus on the thermodynamic properties. An advanced MC scheme, the hybrid Monte Carlo (HMC) method, which can efficiently sample the configurations of a system with large degrees of freedom, however, has limitations in the constrained-sampling applications. In this work, we propose an adjusted constrained HMC method to compute free energy profiles using the thermodynamic integration (TI) method. The key idea of the method for handling the constraint in TI is to integrate the reaction coordinate and sample the rest degrees of freedom by two types of MC schemes, the HMC scheme and the Metropolis algorithm with unbiased trials (M(RT)2-UB). We test the proposed method on three different systems involving two kinds of reaction coordinates, which are the distance between two particles and the difference of particles' distances, and compare the results to those generated by the constrained M(RT)2-UB method serving as benchmarks. We show that our proposed method has the advantages of high sampling efficiency and convenience of implementation, and the accuracy is justified as well. In addition, we show in the third test system that the proposed constrained HMC method can be combined with the path integral method to consider the nuclear quantum effects, indicating a broader application scenario of the sampling method reported in this work.

Pseudo-marginal approximation to the free energy in a micro–macro Markov chain Monte Carlo method

Pseudo-Marginal Approximation to the Free Energy in a Micro-Macro Markov Chain Monte Carlo Method

Computational efficiency study of a micro-macro Markov chain Monte Carlo method for molecular dynamics

Efficiency and Parameter Selection of a micro-macro Markov chain Monte Carlo method for molecular dynamics

A micro-macro Markov chain Monte Carlo method for molecular dynamics using reaction coordinate proposals I: direct reconstruction

A micro-macro Markov chain Monte Carlo method for molecular dynamics using reaction coordinate proposals II: indirect reconstruction

A micro/macro algorithm to accelerate Monte Carlo simulation of stochastic differential equations

Accelerating delayed-acceptance Markov chain Monte Carlo algorithms

Graph-accelerated Markov Chain Monte Carlo using Approximate Samples

Efficient Irreversible Monte Carlo Samplers

Unbiased Markov chain quasi-Monte Carlo for Gibbs samplers

A Vine-Copula Based Adaptive Mcmc Sampler for Efficient Inference of Dynamical Systems

Efficiency of a micro-macro acceleration method for scale-separated stochastic differential equations

A non-asymptotic error analysis for parallel Monte Carlo estimation from many short Markov chains

Microcanonical Hamiltonian Monte Carlo

Unbiased Markov chain Monte Carlo methods with couplings

Accelerating Markov Chain Monte Carlo sampling with diffusion models

A quantitative study on the approximation error and speed-up of the multi-scale MCMC (Monte Carlo Markov chain) method for molecular dynamics

Constrained Hybrid Monte Carlo Sampling Made Simple for Chemical Reaction Simulations

Adaptive, delayed-acceptance MCMC for targets with expensive likelihoods

A practical guide to pseudo-marginal methods for computational inference in systems biology