On the analytical representation of free energy profiles with a Morse/long-range model: application to the water dimer.

We investigate the analytical representation of potentials of mean force (pmf) using the Morse/long-range (MLR) potential approach. The MLR method had previously been used to represent potential energy surfaces, and we assess its validity for representing free-energies. The advantage of the approach is that the potential of mean force data only needs to be calculated in the short to medium range region of the reaction coordinate while the long range can be handled analytically. This can result in significant savings in terms of computational effort since one does not need to cover the whole range of the reaction coordinate during simulations. The water dimer with rigid monomers whose interactions are described by the commonly used TIP4P model [W. Jorgensen and J. Madura, Mol. Phys. 56, 1381 (1985)] is used as a test case. We first calculate an "exact" pmf using direct Monte Carlo (MC) integration and term such a calculation as our gold standard (GS). Second, we compare this GS with several MLR fits to the GS to test the validity of the fitting procedure. We then obtain the water dimer pmf using metadynamics simulations in a limited range of the reaction coordinate and show how the MLR treatment allows the accurate generation of the full pmf. We finally calculate the transition state theory rate constant for the water dimer dissociation process using the GS, the GS MLR fits, and the metadynamics MLR fits. Our approach can yield a compact, smooth, and accurate analytical representation of pmf data with reduced computational cost.

[1]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[2]  B. C. Garrett,et al.  Dynamical nucleation theory: Calculation of condensation rate constants for small water clusters , 1999 .

[3]  S. Kerisit,et al.  Transition path sampling of water exchange rates and mechanisms around aqueous ions. , 2009, The Journal of chemical physics.

[4]  B. C. Garrett,et al.  Dynamical benchmarks of the nucleation kinetics of water , 2002 .

[5]  Alexander M. Mebel,et al.  On the calculation of the dissociation rate constant of the water dimer by the ab initio anharmonic RRKM theory , 2009 .

[6]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[7]  D. Truhlar,et al.  Canonical variational theory for enzyme kinetics with the protein mean force and multidimensional quantum mechanical tunneling dynamics. Theory and application to liver alcohol dehydrogenase , 2001 .

[8]  William L. Jorgensen,et al.  Temperature and size dependence for Monte Carlo simulations of TIP4P water , 1985 .

[9]  Calvin Jary,et al.  An accurate analytic potential function for ground-state N2 from a direct-potential-fit analysis of spectroscopic data. , 2006, The Journal of chemical physics.

[10]  R. L. Roy,et al.  Determining Equilibrium Structures and Potential Energy Functions for Diatomic Molecules , 2012 .

[11]  J. Hynes,et al.  Solvent Barriers in Unimolecular Ionizations. 1. Reaction Path Analysis for Alkyl Iodides , 1994 .

[12]  R. L. Roy,et al.  A new potential function form incorporating extended long-range behaviour: application to ground-state Ca2 , 2007 .

[13]  Peter A. Kollman,et al.  Benzene Dimer: A Good Model for π−π Interactions in Proteins? A Comparison between the Benzene and the Toluene Dimers in the Gas Phase and in an Aqueous Solution , 1996 .

[14]  B. Roux,et al.  Computations of standard binding free energies with molecular dynamics simulations. , 2009, The journal of physical chemistry. B.

[15]  D. Xie,et al.  A new four-dimensional ab initio potential energy surface for N2O-He and vibrational band origin shifts for the N2O-He(N) clusters with N = 1-40. , 2012, The Journal of chemical physics.

[16]  Carl C. Haugen,et al.  Long-range damping functions improve the short-range behaviour of ‘MLR’ potential energy functions , 2011 .

[17]  D. Truhlar,et al.  Free Energy Surface, Reaction Paths, and Kinetic Isotope Effect of Short-Chain Acyl-CoA Dehydrogenase , 2003 .

[18]  Bernd A. Berg,et al.  From data to probability densities without histograms , 2007, Comput. Phys. Commun..

[19]  B. Roux.  The calculation of the potential of mean force using computer simulations , 1995 .

[20]  P. Bernath,et al.  Ground state potential energy curve and dissociation energy of MgH. , 2007, The journal of physical chemistry. A.

[21]  A. Pohorille,et al.  Interactions of anesthetics with the water-hexane interface. A molecular dynamics study. , 1997, The journal of physical chemistry. B.

[22]  Donald G. Truhlar,et al.  Hydride transfer catalyzed by xylose isomerase: Mechanism and quantum effects , 2003, J. Comput. Chem..

[23]  The dependence of the potential of mean force on the solvent friction: Consequences for condensed phase activated rate theories , 1993 .

[24]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[25]  Jianzhong Liu,et al.  Comparison of TIP5P water model and TIP4P water model in cluster nucleation kinetics study through umbrella sampling and free energy perturbation , 2006 .

[26]  Yi Ming,et al.  Free energy perturbation study of water dimer dissociation kinetics. , 2004, The Journal of chemical physics.

[27]  A. Pohorille,et al.  Insights into the recognition and association of transmembrane α-helices. The free energy of α-helix dimerization in glycophorin A , 2005 .

[28]  A. Laio,et al.  Equilibrium free energies from nonequilibrium metadynamics. , 2006, Physical Review Letters.

[29]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[30]  F. Calvo,et al.  Sampling along reaction coordinates with the Wang-Landau method , 2002, cond-mat/0205428.

[31]  R. Le Roy,et al.  Analytic three-dimensional 'MLR' potential energy surface for CO(2)-He, and its predicted microwave and infrared spectra. , 2008, Physical chemistry chemical physics : PCCP.

[32]  Eric C. Cyr,et al.  Using the method of weighted residuals to compute potentials of mean force , 2007, J. Comput. Phys..

[33]  N. Blinov,et al.  Path-integral Monte Carlo simulation of nu3 vibrational shifts for CO2 in (He)n clusters critically tests the He-CO2 potential energy surface. , 2009, The Journal of chemical physics.

[34]  P. Roy,et al.  Analytic Morse/long-range potential energy surfaces and predicted infrared spectra for CO2-H2. , 2010, The Journal of chemical physics.

[35]  N. Dattani,et al.  Accurate analytic potentials for Li(2)(X (1)Sigma(g) (+)) and Li(2)(A (1)Sigma(u) (+)) from 2 to 90 A, and the radiative lifetime of Li(2p). , 2009, The Journal of chemical physics.

[36]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[37]  B. C. Garrett,et al.  Variational transition state theory evaluation of the rate constant for proton transfer in a polar solvent , 2001 .

[38]  A. McKellar,et al.  Molecular superfluid: nonclassical rotations in doped para-hydrogen clusters. , 2010, Physical review letters.

[39]  R. L. Roy,et al.  Uncertainty, Sensitivity, Convergence, and Rounding in Performing and Reporting Least-Squares Fits. , 1998 .

[40]  Jeremy Schofield,et al.  Constructing smooth potentials of mean force, radial distribution functions, and probability densities from sampled data. , 2009, The Journal of chemical physics.