During the optimization, we used residual properties as they collectively nearly isolate the intermolecular interactions while minimizing the effect of the internal quantum degrees of freedom. However, for the final comparison between the simulated properties and the experimental properties, we have to combine the residual properties with the ideal gas properties to compute the complete thermodynamic properties. We have ensured that the simulated residual properties match the experimental residual properties to the extent possible. However, matching the simulated/model ideal gas properties with experimental ideal gas properties, specifically the ones which come from temperature derivatives of Gibbs energy, requires corrections as we use a rigid water model with classical mechanics.