In molecular dynamics simulations, the desired accuracy of an estimator is limited by the integrated autocorrelation time τ int because the statistical error decreases roughly as the variance times (2τ int /t) 1/2 with the simulation time t. From the results above, this implies that simulations with the TIP3P-F model can result in smaller estimation errors for the same simulation lengths. Equivalently, "independent" samples can be generated within fewer molecular dynamics steps. Indeed, the error bars calculated here with block averages (which do not require the autocorrelation time explicitly), tend to be smaller for the systems simulated with the TIP3P-F model over the same total time.