The numerial values for the Ds of the cations and anions in ILs are affected by numerous factors. For example, the relative cationic and anionic sizes, geometric shape, the strength of vander Waals and Coulombic interactions between the cation and anion of the corresponding IL are just some of the factors that can affect the Ds of ILs. On top of this, differences in the type of the force fields, equilibration times, length of production runs, and a method for calculating the Ds may have a huge impact on the computed Ds. Therefore, the direct comparison of our computed Ds values with other MD simulations, as well as reported experimental results, is rather difficult. Recent scientific reports have emphasized the effect of the force field on differences in the calculated Ds of ILs. Our own experience from the current work shows that simulation conditions such as the potential truncation at the cutoff radius, finite-size effects due to the use of small simulation cells, numerical imprecision in the integration of the equation of motion, and most importantly, insufficient equilibration of the system before the production phase of a simulation are the most important factors than need to be considered during classical MD simulations . These factors may potentially introduce systematic errors in the MD simulation results, which consequently may give rise to substantial differences in the numerical values of the calculated Ds. Therefore, it requires the careful choice of such simulation parameters in order to eliminate these systematic errors . Most importantly, the choice of the MD simulation method for calculating the Ds is a trade-off between computational costs, the complexity of the IL, the specific property to be computed, the temperature of simulation, and the accuracy required . Because ILs do have large polarizable ions, and thus they require polarizable force fields to represent their behavior. The work of Bhargava and Balasubramanian shows that equilibrium simulations may give Ds values close to experimental values for ILs if a custom designed force field without polarizability is employed.