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.

Combine Umbrella Sampling with Integrated Tempering Method for Efficient and Accurate Calculation of Free Energy Changes of Complex Energy Surface

Efficient Free Energy Calculations By Combining Two Complementary Tempering Sampling Methods

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

An Efficient Umbrella Potential for the Accurate Calculation of Free Energies by Molecular Simulation.

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

Free energy calculations of a proton transfer reaction by simulated tempering umbrella sampling first principles molecular dynamics simulations

Buckets Instead of Umbrellas for Enhanced Sampling and Free Energy Calculations

Chapter 10. Free Energy Analysis Algorithms along Transition Paths and Transmembrane Ion Permeation

Combining the Biased and Unbiased Sampling Strategy into One Convenient Free Energy Calculation Method.

Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer

Combining Metadynamics and Integrated Tempering Sampling

Do all roads lead to Rome? Convergence issues in umbrella sampling simulations

A Selective Integrated Tempering Method.

Boosting the Conformational Sampling by Combining Replica Exchange with Solute Tempering and Well-Sliced Metadynamics

Calculation of Collective Variable-Based Pmf by Combining Wham with Umbrella Sampling

Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations

Tutorial on Umbrella Sampling Simulation with a Combined QM/MM Potential: The Potential of Mean Force for an SN2 Reaction in Water

Constrained Hybrid Monte Carlo Sampling Made Simple for Chemical Reaction Simulations

Study on the Application of the Combination of TMD Simulation and Umbrella Sampling in PMF Calculation for Molecular Conformational Transitions

Understanding Free-Energy Perturbation Calculations Through A Model of Harmonic Oscillators: Theory and Implications to Improve the Sampling Efficiency by Molecular Simulation

On The Enhanced Sampling Over Energy Barriers In Molecular Dynamics Simulations