Å. In addition to the peak positions, other regions of the OO RDF also show better agreement with experiment, especially the fact that the distribution becomes more delocalized in the region of [3,4] Å. This is consistent with our previous finding in the NVT simulations with original q-AQUA potential. Similar behavior is seen in the and HH RDFs which further verifies the accuracy of q-AQUA-pol potential and the importance of nuclear quantum effects in describing the covalent OH bonds in liquid water, as well as the H-bonded and non-H-bonded solvation shells.