The low mass of the hydrogen atom makes nuclear quantum effects (NQEs) significant when simulating water properties . For example, tunneling effects can affect the formation and breaking of hydrogen bonds and influence the dynamics . The results presented in this work (Fig. ) should therefore be interpreted in light of the fact that NQEs are absent in CPMD or BOMD dynamics with a classical propagation of nuclei. As an illustration, taking into account NQEs with revPBE-D3 revealed that its good agreement with water properties using classical nuclei is due to a fortuitous cancellation of errors, where the neglect of exact exchange compensates for the neglect of quantum nuclei . Advanced path integral molecular dynamics (PIMD) methods are necessary for quantum-mechanical treatment of nuclei, particularly when comparing high-level electronic structure calculations with experimental results . However, this comes at the cost of approximately two orders of magnitude more computational expense than simulations where the nuclei are treated classically. As a result, it has been common practice to mimic NQEs by performing classical (nuclei) MD at elevated temperatures increased by around 30 K . While this ad hoc technique was found to provide reasonable accuracy for RDFs, it often fails to correctly reproduce the dynamical properties that become too fast compared to proper NQEs . Alternatively, recent advances have enabled the acceleration of PIMD dynamics, especially with the help of ML potentials that infer DFT energies and forces at a much reduced cost .