where N is the number of water molecules, r i (t) the position of each oxygen atom i at time t, and the brackets indicate an average over the NVE ensemble. Improved statistics were gathered across multiple lag times and time origins according to the default parameters of the Diffusion Coefficient Tool plugin 186 for VMD to finally obtain D L from the average slope of the mean-squared displacement (MSD).