We evaluated the surface tension using the test-area method, for TIP3P at 298 K, 1.0 atm, using various van der Waals cutoff distances. Figure shows that as the cutoff distance increases, the simulated surface tension continues to increase even at the distance of 18 Å. To estimate the surface tension at infinite cutoff distance, we performed a linear extrapolation on the surface tension against the inverse of the cutoff value, as shown in Figure . The final value of 53.87 mJ m -2 obtained from the intercept of the linear extrapolation is about 4 mJ m -2 higher than the value computed with cutoff at 15 Å. This indicates that our previously shown surface tension evaluated with the cutoff distance of 15 Å may be underestimated by a small amount. Since our simulated value in Figure were below the experimental value, a larger cutoff brings the theoretical surface tension into a closer agreement. In Table , the surface tension extrapolated to infinite cutoff are reported as γ ∞ . The TIP3P-ST model achieved the best agreement of within 1 mJ m -2 to the experiment. However, since the extrapolation is computationally time-consuming, we decided to use the 15 Å cutoff distance through out our parameterization. As a result, the weight of the surface tension property