In molecular-dynamics simulations of fluids, the Einstein-Helfand (EH) and Green-Kubo (GK) relationships are frequently used to compute a variety of transport coefficients, including diffusion coefficients. These relationships are formally valid in the limit of infinite sampling time: The error in the estimate of a transport coefficient (relative to an infinitely long simulation) asymptotically approaches zero as more dynamics are simulated and recorded. In practice, of course, one can only simulate a finite number of particles for a finite amount of time. In this work, we show that in this pre-asymptotic regime, an approach for estimating diffusion coefficients based upon excess entropy scaling (EES) achieves a significantly lower error than either EH or GK relationships at fixed online sampling time. This approach requires access only to structural information at the level of the radial distribution function (RDF). We further demonstrate that the use of a recently developed RDF mollification scheme significantly reduces the amount of sampling time needed to converge to the long-time value of the diffusion coefficient. We also demonstrate favorable sample-to-sample variances in the diffusion coefficient estimate obtained using EES as compared to those obtained using EH and GK.
Read full abstract