Abstract:A robust and simple implementation of the generalized Einstein formulation using single equilibrium molecular dynamics simulation is introduced to compute diffusion and shear viscosity. The unique features underlying this framework are as follows: (1) The use of a simple binary-based method to sample time-dependent transport coefficients results in a uniform distribution of data on a logarithmic time scale. Although we sample “on-the-fly,” the algorithm is readily applicable for post-processing analysis. Overlapping same-length segments are not sampled as they indicate strong correlations. (2) Transport coefficients are estimated using a power law fitting function, a generalization of the standard linear relation, that accurately describes the long-time plateau. (3) The use of a generalized least squares (GLS) fitting estimator to explicitly consider correlations between fitted data points results in a reliable estimate of the statistical uncertainties in a single run. (4) The covariance matrix for the GLS method is estimated analytically using the Wiener process statistics and computed variances. (5) We provide a Python script to perform the fits and automate the procedure to determine the optimal fitting domain. The framework is applied to two fluids, binary hard sphere and a Lennard-Jones near the triple point, and the validity of the single-run estimates is verified against multiple independent runs. The approach should be applicable to other transport coefficients since the diffusive limit is universal to all of them. Given its rigor and simplicity, this methodology can be readily incorporated into standard molecular dynamics packages using on-the-fly or post-processing analysis.

Optimal Estimation Of The Diffusion Coefficient From Noisy Time-Series Measurements

Optimal estimates of diffusion coefficients from molecular dynamics simulations

Local Linear Estimation of Second-Order Diffusion Models

Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation

Maximum likelihood estimates of diffusion coefficients from single-particle tracking experiments

Estimation of the diffusion constant from intermittent trajectories with variable position uncertainties

Modeling of 2D Diffusion Processes Based on Microscopy Data: Parameter Estimation and Practical Identifiability Analysis

Super-resolved anomalous diffusion: deciphering the joint distribution of anomalous exponent and diffusion coefficient

Accurate Estimation of Transport Coefficients Using Model-free Time Correlation Functions in Equilibrium Simulations

Single-image diffusion coefficient measurements of proteins in free solution

A Method to Measure the Diffusion Coefficient in Liquids

Numerical Recovery of the Diffusion Coefficient in Diffusion Equations from Terminal Measurement

An Efficient Non-Iterative Algorithm to Estimate the Apparent Diffusion Coefficient of Water Molecules in the Magnetic Resonance Diffusion Imaging Experiments

Simultaneous Estimation of Relaxation Time and Diffusion Coefficient During the Rapid Transient Mass Transfer

Diffusion constant of a nonspecifically bound protein undergoing curvilinear motion along DNA.

Consistent diffusion matrix estimation from population time series

Statistical Spatially Inhomogeneous Diffusion Inference

Efficient single-run implementation of generalized Einstein relation to compute transport coefficients: A binary-based time sampling

Superstatistical modelling of protein diffusion dynamics in bacteria

A Non-iterative Algorithm to Estimate Water Molecule Diffusion Coefficient in Magnetic Resonance Imaging

Optimizing the diffusion coefficient of overdamped Langevin dynamics