Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation

Abstract

Self-diffusion coefficients, $D^{\ast }$, are routinely estimated from molecular dynamics simulations by fitting a linear model to the observed mean-squared displacements (MSDs) of mobile species. MSDs derived from simulation exhibit statistical noise that causes uncertainty in the resulting estimate of $D^{\ast }$. An optimal scheme for estimating $D^{\ast }$ minimises this uncertainty, i.e., it will have high statistical efficiency, and also gives an accurate estimate of the uncertainty itself. We present a scheme for estimating $\D$ from a single simulation trajectory with high statistical efficiency and accurately estimating the uncertainty in the predicted value. The statistical distribution of MSDs observable from a given simulation is modelled as a multivariate normal distribution using an analytical covariance matrix for an equivalent system of freely diffusing particles, which we parameterise from the available simulation data. We use Bayesian regression to sample the distribution of linear models that are compatible with this multivariate normal distribution, to obtain a statistically efficient estimate of $D^{\ast }$ and an accurate estimate of the associated statistical uncertainty.

References in corpus (5)

• emcee v3: A Python ensemble sampling toolkit for affine-invariant MCMC
• Power spectral density of a single Brownian trajectory: What one can and cannot learn from it
• Sparse Cyclic Excitations Explain the Low Ionic Conductivity of Stoichiometric Li${}_7$La${}_3$Zr${}_2$O${}_{12}$
• Optimal estimates of diffusion coefficients from molecular dynamics simulations
• Relationships Between Atomic Diffusion Mechanisms and Ensemble Transport Coefficients in Crystalline Polymorphs

Cited by in corpus (4)

• Flow Matching for Accelerated Simulation of Atomic Transport in Crystalline Materials
• Improved mean squared displacement analysis for anomalous single particle trajectories
• A Residence-Time Approach for Determining Position-Dependent Diffusivities from Biased Molecular Simulations
• VacHopPy: A Python package for vacancy hopping analysis based on molecular dynamics simulations