Based on our observations, as described in the supporting information about the choice of optimization function along with the observations from IAPWS95 parameterization we formulated our objective function to the following form, as an optimization in the residual free energy with constraints in physical observables which are derivative properties of the residual free energy.