Estimating Statistical Distributions Using an Integral Identity

Cheng Zhang,Jianpeng Ma
DOI: https://doi.org/10.1063/1.4721638
2012-01-01
Abstract:We present an identity for an unbiased estimate of a general statisticaldistribution. The identity computes the distribution density from dividing ahistogram sum over a local window by a correction factor from a mean-forceintegral, and the mean force can be evaluated as a configuration average. Weshow that the optimal window size is roughly the inverse of the localmean-force fluctuation. The new identity offers a more robust and preciseestimate than a previous one by Adib and Jarzynski [J. Chem. Phys. 122, 014114,(2005)]. It also allows a straightforward generalization to an arbitraryensemble and a joint distribution of multiple variables. Particularly we derivea mean-force enhanced version of the weighted histogram analysis method (WHAM).The method can be used to improve distributions computed from molecularsimulations. We illustrate the use in computing a potential energydistribution, a volume distribution in a constant-pressure ensemble, a radialdistribution function and a joint distribution of amino acid backbone dihedralangles.
What problem does this paper attempt to address?