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.

Importance Sampling for Equilibrium Properties in Non-Adiabatic Systems: Application to Anomalous Interfacial Charge Transfer

Dynamical computation of the density of states and Bayes factors using nonequilibrium importance sampling

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

Monte Carlo simulation of electrostatic interactions in inhomogeneous dielectric media: Correct sampling for the local lattice simulation algorithm

Efficient Sampling for Machine Learning Electron Density and Its Response in Real Space

On the use of QM/MM Surface Hopping simulations to understand thermally-activated rare event nonadiabatic transitions in the condensed phase

Equilibrium sampling by re-weighting non-equilibrium simulation trajectories

Efficient Phase Diagram Sampling by Active Learning

Towards accurate thermodynamics from random energy sampling

Training Transferable Interatomic Neural Network Potentials for Reactive Chemistry: Improved Chemical Space Sampling

Importance Sampling within Configuration Space Integration for Adsorbate Thermophysical Properties: A Case Study for CH3/Ni(111)

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation

Constrained Hybrid Monte Carlo Sampling Made Simple for Chemical Reaction Simulations

An Integrate-over-temperature Approach for Enhanced Sampling.

Mixed Quantum-Classical Approaches to Calculating Charge Transfer Rate Constants: Applications to Realistic Systems

Two Improved Electronegativity Equalization Methods for Charge Distribution in Large Scale Non-Uniform System

Machine-learning-accelerated simulations to enable automatic surface reconstruction

Density Fluctuations, Solvation Thermodynamics and Coexistence Curves in Grand Canonical Molecular Dynamics Simulations

Molecular free energies, rates, and mechanisms from data-efficient path sampling simulations

Improving the Efficiency of Non-equilibrium Sampling in the Aqueous Environment via Implicit-Solvent Simulations.

Using Overlap and Funnel Sampling to Obtain Accurate Free Energies from Nonequilibrium Work Measurements.