Free energy of particle insertion

In free energy simulations involving the creation of a new particle in a condensed medium, there is a well-known ‘origin singularity’ when the particle is first introduced, unless the coupling of the particle with its surroundings is increased very slowly. This singularity is analysed theoretically for a liquid that interacts with the new particle through a soft-sphere or a Lennard-Jones potential, using a virial expansion of the free energy derivative. For a soft-sphere potential u(r) = λα Ar-n (where λ is a coupling constant), we show rigorously that the derivative of the free energy with respect to λ is infinite at the origin if and only if α < n/3. This confirms a previous, approximate, prediction based on scaled particle theory. For a Lennard-Jones potential u(r) = λα Ar -12 + λβ Br -6, such that α ˇ- 2β, the derivative of the free energy is again infinite if and only if α < n/3. If the volume of the solute is proportional to the coupling constant, there is no singularity.

[1]  L. Schwartz.  Méthodes mathématiques pour les sciences physiques , 1961 .

[2]  J. G. Powles,et al.  Non-destructive molecular-dynamics simulation of the chemical potential of a fluid , 1982 .

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

[4]  R. Pierotti,et al.  Aqueous Solutions of Nonpolar Gases1 , 1965 .

[5]  Peter A. Kollman,et al.  A new method for carrying out free energy perturbation calculations: Dynamically modified windows , 1989 .

[6]  William L. Jorgensen,et al.  Efficient computation of absolute free energies of binding by computer simulations. Application to the methane dimer in water , 1988 .

[7]  K. Shing,et al.  The chemical potential from computer simulation , 1981 .

[8]  H. Berendsen,et al.  THERMODYNAMICS OF CAVITY FORMATION IN WATER - A MOLECULAR-DYNAMICS STUDY , 1982 .

[9]  B. Widom,et al.  Some Topics in the Theory of Fluids , 1963 .

[10]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[11]  L. V. Hove,et al.  Radial Distribution Function of a Gas of Hard Spheres and the Superposition Approximation , 1952 .

[12]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[13]  Mihaly Mezei,et al.  A cavity-biased (T, V, μ) Monte Carlo method for the computer simulation of fluids , 1980 .

[14]  J. Prausnitz.  Thermodynamic properties of dense fluid mixtures , 1967 .

[15]  M. Mezei.  Excess free energy of different water models computed by Monte Carlo methods , 1982 .

[16]  Ab initio calculation of the free energy of liquid water , 1978 .

[17]  J. Hermans,et al.  Excess free energy of liquids from molecular dynamics simulations. Application to water models. , 1988, Journal of the American Chemical Society.