In the case of hybrid functionals, there are only marginal improvements. Results using zero damping are generally unsatisfactory, as the density dierence between ice and water is too large, precluding the observation of an isobar density maximum resembling experimental results. Although BJ damping yields more reasonable results, they are not markedly superior to those using the semi-local functional RPBE-D3/zd. Empirically, we only achieve excellent agreement with the experiment when we combine an MLP with zero damping and BJ damping. We conclude that the state of density functionals for water is not satisfactory, and much more development eort will be required to resolve this issue.