Coarser properties for the condensed phase liquid such as the OO, OH and HH radial distribution functions and self-diffusion coefficient were accurately obtained using q-AQUA (in NVT calculations) as well as MB-pol (in NPT calculations). These properties are also accurately described by q-AQUA-pol in NPT calculations. Nuclear quantum effects, described using PIMD calculations, using q-AQUA and here using q-AQUA-pol, were found to be noticeable and generally to bring theory closer to experiment. The current results using q-AQUA-pol for the temperature dependence of the density of H 2 O are especially encouraging. They (both classical and PIMD) are significantly closer to experiment than those from MB-pol.