On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water

Abstract

The effect of the applied trajectory length on the convergence of the static dielectric constant and the self-diffusion coefficient were examined for the SPC/E water model in the NVT ensemble with different system size at 293 K. Very long simulation times of 6-8 ns were employed in order to track the convergence of these properties. Temperature dependence and isotope effects, via using D${}_2$O instead of H${}_2$O, were also investigated. A simulation for the polarizable SWM4-DP model was also carried out to compare the effect of different potential models. Radial distribution functions and the neutron weighted structure factor were also calculated; they were found to be insensitive to changing the system size in the range of 216-16000 molecules. On the other hand, the static dielectric constant and the diffusion coefficient are rather sensitive to the applied trajectory length, system size and the method of calculation. These latter properties are therefore not appropriate for assessing, and distinguishing between, potential models of water. It is clearly shown that trajectories shorter than about 6 ns are not sufficient for a sufficiently accurate determination of the dielectric constant of this water model.

References in corpus (1)

• What ice can teach us about water interactions: a critical comparison of the performance of different water models

Cited by in corpus (10)

• Polar nanoregions in water - a study of the dielectric properties of TIP4P/2005, TIP4P/2005f and TTM3F
• A comparison of classical interatomic potentials applied to highly concentrated aqueous lithium chloride solutions
• Dielectric constant of supercritical water in a large pressure-temperature range
• Temperature dependence of the microscopic structure and density anomaly of the SPC/E and TIP4P-Ew water models. Molecular dynamics simulation results
• Performance Evaluation of the Zero-Multipole Summation Method in Modern Molecular Dynamics Software
• Uncertainties in the Static Dielectric Constants computed from Molecular Dynamics Simulations
• On the composition dependence of thermodynamic, dynamic and dielectric properties of water-methanol model mixtures. Molecular dynamics simulation results
• Comparison of the TIP4P-2005, SWM4-DP and BK3 interaction potentials of liquid water with respect to their consistency with neutron and X-ray diffraction data of pure water
• High-Accuracy Physical Property Prediction for Organics via Molecular Representation Learning: Bridging Data to Discovery
• Static Dielectric Permittivity Profiles and Coarse-graining Approaches for Water in Graphene Slit Pores