Statistical variances of diffusional properties from ab initio molecular dynamics simulations

Abstract

Ab initio molecular dynamics (AIMD) simulation is widely employed in studying diffusion mechanisms and in quantifying diffusional properties of materials. However, AIMD simulations are often limited to a few hundred atoms and a short, sub-nanosecond physical timescale, which leads to models that include only a limited number of diffusion events. As a result, the diffusional properties obtained from AIMD simulations are often plagued by poor statistics. In this paper, we reexamine the process to estimate diffusivity and ionic conductivity from the AIMD simulations and establish the procedure to minimize the fitting errors. In addition, we propose methods for quantifying the statistical variance of the diffusivity and ionic conductivity from the number of diffusion events observed during the AIMD simulation. Since an adequate number of diffusion events must be sampled, AIMD simulations should be sufficiently long and can only be performed on materials with reasonably fast diffusion. We chart the ranges of materials and physical conditions that can be accessible by AIMD simulations in studying diffusional properties. Our work provides the foundation for quantifying the statistical confidence levels of diffusion results from AIMD simulations and for correctly employing this powerful technique.

Computation‐Guided Design of LiTaSiO5, a New Lithium Ionic Conductor with Sphene Structure

“…The total time of AIMD simulations were in the range of 100–600 ps until the diffusivity was converged. The ionic conductivity and their error bars were calculated following established methods in previous studies …”

Computer Simulation of Cathode Materials for Lithium Ion and Lithium Batteries: A Review

“…These calculations often help to clarify the diffusion mechanisms that are difficult to obtain experimentally. [82] For example, no significant Li migration was observed in LiFePO 4 at 600 and 1000 K after 20 ps of AIMD simulations. [18] It suggests that the chemical diffusivity will depends on the local lithium concentration.…”

“…Longer trajectory and higher temperatures can help to capture sufficient number of ion hops to reduce statistical uncertainty. [82] For example, no significant Li migration was observed in LiFePO 4 at 600 and 1000 K after 20 ps of AIMD simulations. [83] At 2000 K, A zigzag diffusion pathway can be identified, as shown in Figure 11.…”

Doping in garnet-type electrolytes: Kinetic and thermodynamic effects from molecular dynamics simulations

“…In particular, it tends to zero for large T , when the simulation is long enough to determine exactly the value of the diffusion coefficient at a fixed configuration Γ D . Sampling noise is the only one which must be considered in simulations of undoped materials [39,58,59], or in an implicit doping setting. In the case of an explicitly doped material under analysis, an additional mean over configurations, as explained in the main text, can be performed:…”