where µ is the dipole moment of the simulated model, µ gas and α gas are the dipole moment and average polarizability of a water molecule in the gas phase, respectively. The last term in equation ( ) contains corrections that account for the vibrational effects of water molecules and non-ideality of the gas phase. These corrections are reported in Ref. for different temperatures. Experimental reference data are taken from Ref. Specific heat capacity . We computed the isobaric heat capacity c p using the enthalpy fluctuation formula, namely: