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

Journal of Chemical Physics

Volumn 132, Issue 15, 2010, Pages

  Van Zon, Ramses ^a   Schofield, Jeremy ^a  

^a UNIVERSITY OF TORONTO   (Canada)

Author keywords

[No Author keywords available]

Indexed keywords

A-DENSITY; ANALYTICAL ESTIMATES; ANALYTICAL FIT; BIN SIZE; DIFFUSING PARTICLES; FIRST ARRIVAL; FITTING PROCEDURE; FOURIER MODES; HARRIS METHOD; JACOBIANS; LONG TAIL; NONUNIFORM; PHYSICAL PARAMETERS; PIECE-WISE; POTENTIALS OF MEAN FORCES; PROBABILITY DENSITIES; RADIAL DISTRIBUTION FUNCTIONS; RESAMPLING SCHEME; SAMPLED DATA; STANDARD METHOD; WATER MODELS;

GRAPHIC METHODS; PROBABILITY DENSITY FUNCTION; RANDOM VARIABLES; STATISTICAL METHODS;

DISTRIBUTION FUNCTIONS;

EID: 77951695581     PISSN: 00219606     EISSN: None     Source Type: Journal    
DOI: 10.1063/1.3366523     Document Type: Article

References (31)

1. D. A. McQuarrie, Statistical Mechanics (Harper Collins, New York, 1976).
2. D. Frenkel and B. Smit, Understanding Molecular Simulation: From Algorithms to Applications (Academic, San Diego, 2002).
3. H. J. C. Berendsen, Simulating the Physical World: Hierarchical Modeling from Quantum Mechanics to Fluid Dynamics (Cambridge University Press, Cambridge, 2007).
4. B. A. Berg and R. C. Harris, Comput. Phys. Commun. CPHCBZ 0010-4655 179, 443 (2008). 10.1016/j.cpc.2008.03.010
5. R. van Zon, L. Hernández de la Pea, G. H. Peslherbe, and J. Schofield, Phys. Rev. E PLEEE8 1063-651X 78, 041103 (2008). 10.1103/PhysRevE.78. 041103
6. T. Gericke, P. Würtz, D. Reitz, T. Langen, and H. Ott, Nat. Phys. ZZZZZZ 1745-2473 4, 949 (2008). 10.1038/nphys1102
7. L. Mones, P. Kulhánek, I. Simon, A. Laio, and M. Fuxreiter, J. Phys. Chem. B JPCBFK 1089-5647 113, 7867 (2009). 10.1021/jp9000576
8. Y. Liu and T. Ichiye, J. Phys. Chem. JPCHAX 0022-3654 100, 2723 (1996). 10.1021/jp952324t
9. A. Chandra and T. Ichiye, J. Chem. Phys. JCPSA6 0021-9606 111, 2701 (1999). 10.1063/1.479546
10. M. -L. Tan, J. T. Fisher, A. Chandra, B. R. Brooks, and T. Ichiye, Chem. Phys. Lett. CHPLBC 0009-2614 376, 646 (2003). 10.1016/S0009-2614(03)01044-3
11. C. J. Fennell and J. D. Gezelter, J. Chem. Phys. JCPSA6 0021-9606 120, 9175 (2004). 10.1063/1.1697381
12. R. van Zon and J. Schofield, J. Chem. Phys. JCPSA6 0021-9606 128, 154119 (2008). 10.1063/1.2901173 (Pubitemid 351574621)
13. L. Hernández de la Pea, R. van Zon, J. Schofield, and S. B. Opps, J. Chem. Phys. JCPSA6 0021-9606 126, 074105 (2007). 10.1063/1.2434957
14. L. Hernández de la Pea, R. van Zon, J. Schofield, and S. B. Opps, J. Chem. Phys. JCPSA6 0021-9606 126, 074106 (2007). 10.1063/1.2434959
15. R. van Zon and J. Schofield, J. Comput. Phys. JCTPAH 0021-9991 225, 145 (2007). 10.1016/j.jcp.2006.11.019
16. E. C. Cyr and S. D. Bond, J. Comput. Phys. JCTPAH 0021-9991 225, 714 (2007). 10.1016/j.jcp.2006.12.015
17. A. B. Adib and C. Jarzynski, J. Chem. Phys. JCPSA6 0021-9606 122, 014114 (2005). 10.1063/1.1829631 (Pubitemid 41047757)
18. J. E. Basner and C. Jarzynski, J. Phys. Chem. B JPCBFK 1089-5647 112, 12722 (2008). 10.1021/jp803635e
19. W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in Fortran, The Art of Scientific Computing, 2nd ed. (Cambridge University Press, Cambridge, 1992).
20. Alternatively, one can use other tests such as the Kuiper test, but these are less convenient for the piecewise approach presented later, since they involve several "hard-to-fit" points.
21. M. A. Stephens, J. R. Stat. Soc. Ser. B (Methodol.) JSTBAJ 0035-9246 32, 115 (1970).
22. B. Efron, Ann. Stat. ASTSC7 0090-5364 7, 1 (1979). 10.1214/aos/1176344552
23. H. R. Künsch, Ann. Stat. ASTSC7 0090-5364 17, 1217 (1989). 10.1214/aos/1176347265
24. Another good explanation of the sensitivity of the KS test to values near the median can be found in Ref., Ch. 4, Sec..
25. One could imagine performing the reweighting at the level of KS test instead, but this would modify the q statistic. The problem is that it would then not be known what a reasonable testing criterion for this new q statistic should be, or even whether it allows a simple formulation such as q> qcut.
26. A procedure in which the piecewise Fourier series is restricted to be continuous across the interval boundaries was also considered but found to lead to undesirable oscillations near discontinuities, reminiscent of the Gibbs phenomenon (Ref.).
27. H. Jeffreys and B. S. Jeffreys, Methods of Mathematical Physics, 3rd ed. (Cambridge University Press, Cambridge, 1956).
28. S. Redner, A Guide to First-Passage Processes (Cambridge University Press, Cambridge, 1991).
29. This is the Laplace inverse of Eq. (6.4.8) in Ref..
30. http://www.chem.utoronto.ca/~rzon/Code.html.
31. This analysis was also attempted for the deviations plotted in Fig. but the error bars on the resulting distributions were too large to tell whether they are biased or non-Gaussian.