Abstract:Molecular‐dynamics (MD) simulations constitute an important technique for studying mass transport in the solid state. In this study, we investigated how different simulation parameters (e.g., time and cell size) determine the statistical error in the tracer diffusion coefficient obtained from an MD simulation. In addition, we distill the characteristics of the uncertainty into a mathematical expression that captures the interrelated influences of the simulation parameters. Analysis of the mean squared displacement of species k , rk2 , as a function of simulation time t constitutes a powerful method for extracting, from a molecular‐dynamics (MD) simulation, the tracer diffusion coefficient, Dk* . The statistical error in Dk* is seldom considered, and when it is done, the error is generally underestimated. In this study, we examined the statistics of rk2t curves generated by solid‐state diffusion by means of kinetic Monte Carlo sampling. Our results indicate that the statistical error in Dk* depends, in a strongly interrelated way, on the simulation time, the cell size, and the number of relevant point defects in the simulation cell. Reducing our results to one key quantity—the number of k particles that have jumped at least once—we derive a closed‐form expression for the relative uncertainty in Dk* . We confirm the accuracy of our expression through comparisons with self‐generated MD diffusion data. With the expression, we formulate a set of simple rules that encourage the efficient use of computational resources for MD simulations.

The need for accuracy and smoothness in numerical simulations

Probabilistic Richardson Extrapolation

Richardson Extrapolation-Based Verification Method of Scientific Calculation Program Without the Oracles: A Case Study

Numerical accuracy of the derivative-expansion-based functional renormalization group

Truncation Error Calculation Based on Richardson Extrapolation for Variable-Step Collaborative Simulation

Explicit and implicit error inhibiting schemes with post-processing

A posteriori error estimates for the Richards equation

Efficient simulation of complex Ginzburg--Landau equations using high-order exponential-type methods

Study of Numerical Errors in Direct Numerical Simulation and Large Eddy Simulation

Evaluating a double integral using Euler's method and Richardson extrapolation

On the reliability of simulations of collisional stellar systems

Probabilistic evolution of the error of numerical method for linear stochastic differential equation

A posteriori error analysis of round-off errors in the numerical solution of ordinary differential equations

Improved Accuracy for Trotter Simulations Using Chebyshev Interpolation

Numerical Issues for a Non-autonomous Logistic Model

Extrapolation of Nystrom Solution for Two-Dimensional Nonlinear Fredholm Integral Equations

Exponential Asymptotics Using Numerical Rational Approximation In Linear Differential Equations

Convergence analysis of Richardson extrapolation for a quasilinear singularly perturbed problem with an integral boundary condition on an adaptive grid

A general expression for the statistical error in a diffusion coefficient obtained from a solid‐state molecular‐dynamics simulation