Figure shows performance over the entire temperature and pressure range examined, we can see that GOPAL has the lowest RMSE in molar volume V , density ρ, isothermal compressibility (κ), and volume expansivity α over the T and P range. GOPAL is the only model which has RMSE in the ρ(T, P ) surface within the upper limit of experimental uncertainty of 0.02%. In comparison, OPC has the lowest RMSE for internal energy U over the T and P surface, though GOPAL, SPC, SPC-E and OPC3 have errors less than 0.1 kJ/mol. Somewhat surprisingly, TIP3P has the lowest RMSE in entropy surface, despite being worse in most other properties. As TIP3P parameterization did not use entropy in its objective function, it is not entirely clear how TIP3P model has the lowest RMSE in entropy over all T and P . GOPAL and all other water models except TIP3P-ST have deviations of ∼ 1.5 J/mol/K in the entropy surface. SPC has the lowest RMSE for Gibbs energy surface. While developing the objective function, we showed that some parameter combinations can minimize errors in Gibbs energy very well but are not able to keep errors low for other property surfaces. SPC model clearly shows this behavior. The RMSE for GOPAL, TIP4P-FB, TIP4P-ST and OPC3 are below 1.5 kJ/mol/K in G surface. GOPAL, SPC, SPC-E, TIP4P, OPC, OPC3, TIP4P-ST, and TIP4P-FB have ∼ 2 J/mol/K RMSE in the heat capacity. The RMSEs for the enthalpy are identical to the internal energy surface as the errors in the molar volume estimates are small and errors in internal energies dominate.