Thermodynamic integration of the free energy along a reaction coordinate in Cartesian coordinates

A generalized formulation of the thermodynamic integration (TI) method for calculating the free energy along a reaction coordinate is derived. Molecular dynamics simulations with a constrained reaction coordinate are used to sample conformations. These are then projected onto conformations with a higher value of the reaction coordinate by means of a vector field. The accompanying change in potential energy plus the divergence of the vector field constitute the derivative of the free energy. Any vector field meeting some simple requirements can be used as the basis of this TI expression. Two classes of vector fields are of particular interest here. The first recovers the conventional TI expression, with its cumbersome dependence on a full set of generalized coordinates. As the free energy is a function of the reaction coordinate only, it should in principle be possible to derive an expression depending exclusively on the definition of the reaction coordinate. This objective is met by the second class of ve...

[1]  M. Fixman.  Classical statistical mechanics of constraints: a theorem and application to polymers. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Martin Zacharias,et al.  Holonomic constraint contributions to free energy differences from thermodynamic integration molecular dynamics simulations , 1992 .

[3]  M. J. Ruiz-Montero,et al.  Efficient schemes to compute diffusive barrier crossing rates , 1997 .

[4]  Jan Hermans,et al.  Change of bond length in free‐energy simulations: Algorithmic improvements, but when is it necessary? , 1994 .

[5]  Harold A. Scheraga,et al.  On the Use of Classical Statistical Mechanics in the Treatment of Polymer Chain Conformation , 1976 .

[6]  Peter A. Kollman,et al.  The overlooked bond‐stretching contribution in free energy perturbation calculations , 1991 .

[7]  P. Krüger,et al.  Free energy as the potential of mean constraint force , 1996 .

[8]  Herman J. C. Berendsen,et al.  Simulation of Water Transport through a Lipid Membrane , 1994 .

[9]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[10]  D. Herschbach,et al.  Molecular Partition Functions in Terms of Local Properties , 1959 .

[11]  Michiel Sprik,et al.  Free energy from constrained molecular dynamics , 1998 .

[12]  M. Karplus,et al.  The Jacobian factor in free energy simulations , 1996 .

[13]  W. Briels,et al.  Molecular dynamics simulations of free energy and conformational transition rates of calix[4]arene in chloroform , 1997 .

[14]  David A. Pearlman,et al.  Determining the contributions of constraints in free energy calculations: Development, characterization, and recommendations , 1993 .

[15]  Giovanni Ciccotti,et al.  Molecular dynamics simulation of rigid molecules , 1986 .

[16]  Wilfred F. van Gunsteren,et al.  Computer Simulation of Biomolecular Systems: Theoretical and Experimental Applications , 1989 .

[17]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[18]  Krzysztof Kuczera.  One‐ and multidimensional conformational free energy simulations , 1996 .

[19]  RESEARCH NOTE Constrained dynamics: compensation of metric tensor effects , 1998 .

[20]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[21]  Giovanni Ciccotti,et al.  Andersen's canonical-ensemble molecular dynamics for molecules with constraints , 1986 .

[22]  M. Fixman,et al.  Simulation of polymer dynamics. I. General theory , 1978 .

[23]  Giovanni Ciccotti,et al.  Sampling of molecular conformations by molecular dynamics techniques , 1993 .

[24]  P. Hänggi,et al.  Reaction-rate theory: fifty years after Kramers , 1990 .

[25]  W. Briels,et al.  THE CALCULATION OF FREE-ENERGY DIFFERENCES BY CONSTRAINED MOLECULAR-DYNAMICS SIMULATIONS , 1998 .

[26]  E. Bright Wilson,et al.  Book Reviews: Molecular Vibrations. The Theory of Infrared and Raman Vibrational Spectra , 1955 .

[27]  David Chandler,et al.  Statistical mechanics of isomerization dynamics in liquids and the transition state approximation , 1978 .

[28]  Jonathan W. Essex,et al.  Generalized alteration of structure and parameters: A new method for free‐energy perturbations in systems containing flexible degrees of freedom , 1995, J. Comput. Chem..

[29]  Ab initio molecular dynamics study of the reaction of water with formaldehyde in sulfuric acid solution. , 1998 .