Especially at, 25 • C and 1 atm both models have very close agreement with with experimental density. In Figure , we see that the deviations for both GOPAL and TIP4P-FB are higher for lower temperatures (1 • C-35 • C) as the likely-uncapturable nuclear quantum effects are dominant in that temperature range. For higher temperatures, the deviations become as small as the statistical uncertainties. However, TIP3P-ST consistently overestimates and TIP3P underestimates the heat capacity significantly.