Author(s)
Thomas W. Rosch, Paul Patrone
Abstract
Despite more than 40 years of active research in computational condensed-matter physics, state-of-the-art approaches for simulating radial distribution functions (RDFs) g(r) still rely on binning pair-separations into a histogram. Such representations suffer from undesirable properties, including high uncertainty and slow convergence. Moreover, such problems go undetected by the metrics often used to assess RDFs. To address these issues, we propose (I) a spectral Monte Carlo (SMC) method that yields g(r) as an analytical series expansion; and (II) a Sobolev-type norm that better assesses quality of RDFs by quantifying their fluctuations. Using the latter, we show that, as compared to histograms, SMC reduces both the noise in g(r) and the number of pair separations needed by orders of magnitude. Moreover, SMC yields simple, differentiable formulas for the RDF, which are useful for tasks such as force-field calibration via iterative Boltzmann inversion.
Citation
Physical Review Letters
Citation
Rosch, T.
and Patrone, P.
(2017),
Efficient estimation of simulated radial distribution functions via spectral Monte Carlo, Physical Review Letters, [online], https://tsapps.nist.gov/publication/get_pdf.cfm?pub_id=920798 (Accessed April 30, 2026)
Additional citation formats
Issues
If you have any questions about this publication or are having problems accessing it, please contact [email protected].