Constructing smooth potentials of mean force, radial distribution functions, and probability densities from sampled data.

In this paper a method of obtaining smooth analytical estimates of probability densities, radial distribution functions, and potentials of mean force from sampled data in a statistically controlled fashion is presented. The approach is general and can be applied to any density of a single random variable. The method outlined here avoids the use of histograms, which require the specification of a physical parameter (bin size) and tend to give noisy results. The technique is an extension of the Berg-Harris method [B. A. Berg and R. C. Harris, Comput. Phys. Commun. 179, 443 (2008)], which is typically inaccurate for radial distribution functions and potentials of mean force due to a nonuniform Jacobian factor. In addition, the standard method often requires a large number of Fourier modes to represent radial distribution functions, which tends to lead to oscillatory fits. It is shown that the issues of poor sampling due to a Jacobian factor can be resolved using a biased resampling scheme, while the requirement of a large number of Fourier modes is mitigated through an automated piecewise construction approach. The method is demonstrated by analyzing the radial distribution functions in an energy-discretized water model. In addition, the fitting procedure is illustrated on three more applications for which the original Berg-Harris method is not suitable, namely, a random variable with a discontinuous probability density, a density with long tails, and the distribution of the first arrival times of a diffusing particle to a sphere, which has both long tails and short-time structure. In all cases, the resampled, piecewise analytical fit outperforms the histogram and the original Berg-Harris method.

[1]  M. Stephens.  Use of the Kolmogorov-Smirnov, Cramer-Von Mises and Related Statistics without Extensive Tables , 1970 .

[2]  Discontinuous molecular dynamics for rigid bodies: applications. , 2007, The Journal of chemical physics.

[3]  Christopher J Fennell,et al.  On the structural and transport properties of the soft sticky dipole and related single-point water models. , 2004, The Journal of chemical physics.

[4]  Amalendu Chandra,et al.  Dynamical properties of the soft sticky dipole model of water: Molecular dynamics simulations , 1999 .

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

[6]  B. Brooks,et al.  A temperature of maximum density in soft sticky dipole water , 2003 .

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

[8]  Toshiko Ichiye,et al.  Soft Sticky Dipole Potential for Liquid Water: A New Model , 1996 .

[9]  Ramses van Zon,et al.  Numerical implementation of the exact dynamics of free rigid bodies , 2006, J. Comput. Phys..

[10]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[11]  Philip M. Morse,et al.  Methods of Mathematical Physics , 1947, The Mathematical Gazette.

[12]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[13]  Berend Smit,et al.  Molecular Dynamics Simulations , 2002 .

[14]  P. Würtz,et al.  High-resolution scanning electron microscopy of an ultracold quantum gas , 2008 .

[15]  Gilles H. Peslherbe,et al.  Quantum free-energy differences from nonequilibrium path integrals. I. Methods and numerical application. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  S. Redner.  A guide to first-passage processes , 2001 .

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

[18]  Christopher Jarzynski,et al.  Binless estimation of the potential of mean force. , 2008, The journal of physical chemistry. B.

[19]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[20]  Alessandro Laio,et al.  The energy gap as a universal reaction coordinate for the simulation of chemical reactions. , 2009, The journal of physical chemistry. B.

[21]  Jeremy Schofield,et al.  Discontinuous molecular dynamics for semiflexible and rigid bodies. , 2006, The Journal of chemical physics.

[22]  H. Berendsen.  Simulating the Physical World: Hierarchical Modeling from Quantum Mechanics to Fluid Dynamics , 2007 .

[23]  Unbiased estimators for spatial distribution functions of classical fluids. , 2004, The Journal of chemical physics.

[24]  H. Künsch.  The Jackknife and the Bootstrap for General Stationary Observations , 1989 .

[25]  Jeremy Schofield,et al.  Event-driven dynamics of rigid bodies interacting via discretized potentials. , 2007, The Journal of chemical physics.

[26]  B. Efron.  Bootstrap Methods: Another Look at the Jackknife , 1979 .