Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation

Fuente: arXiv
Salvato in:
Dettagli Bibliografici
Autori principali: McCluskey, Andrew R., Coles, Samuel W., Morgan, Benjamin J.
Natura: Preprint
Pubblicazione: 2023
Soggetti:
Accesso online:
Tags: Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
_version_ 1866917178854342656
author McCluskey, Andrew R.
Coles, Samuel W.
Morgan, Benjamin J.
author_facet McCluskey, Andrew R.
Coles, Samuel W.
Morgan, Benjamin J.
contents Self-diffusion coefficients, $D^*$, 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^*$. An optimal scheme for estimating $D^*$ 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^*$ and an accurate estimate of the associated statistical uncertainty.
format Preprint
id arxiv_https___arxiv_org_abs_2305_18244
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation
McCluskey, Andrew R.
Coles, Samuel W.
Morgan, Benjamin J.
Statistical Mechanics
Materials Science
Self-diffusion coefficients, $D^*$, 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^*$. An optimal scheme for estimating $D^*$ 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^*$ and an accurate estimate of the associated statistical uncertainty.
title Accurate Estimation of Diffusion Coefficients and their Uncertainties from Computer Simulation
topic Statistical Mechanics
Materials Science
url https://arxiv.org/abs/2305.18244