The TIP4P-ST model also predicts ∆H vap slightly higher than the experimental value with a self-polarization correction of 7.03 kJ mol -1 . Notably, the corrected ∆H vap is almost identical to TIP4P-FB, again indicating they are close in the model space. The correction for nuclear quantum effects is relatively small at -0.27 kJ mol -1 at 298 K. The self-diffusion coefficient is another property where TIP3P-ST is different from the other models included in our comparison. The lower self-diffusion coefficient indicates a slightly more structured liquid, with stronger hydrogen bonds needed to reproduce the surface tension. This behavior is also reflected in the radial distribution plot, where the TIP3P-ST curve shows a higher first peak and lower first trough.