Issue 11, 2021

Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water

Abstract

We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T ≥ 273 K) and supercooled (210 ≤ T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid–liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T ≈ 230–235 K. No maximum in CP(T) could be detected down to T ≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.

Graphical abstract: Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water

Supplementary files

Article information

Article type
Paper
Submitted
14 Aug 2020
Accepted
02 Mar 2021
First published
17 Mar 2021

Phys. Chem. Chem. Phys., 2021,23, 6914-6928

Author version available

Nuclear quantum effects on the thermodynamic, structural, and dynamical properties of water

A. Eltareb, G. E. Lopez and N. Giovambattista, Phys. Chem. Chem. Phys., 2021, 23, 6914 DOI: 10.1039/D0CP04325G

To request permission to reproduce material from this article, please go to the Copyright Clearance Center request page.

If you are an author contributing to an RSC publication, you do not need to request permission provided correct acknowledgement is given.

If you are the author of this article, you do not need to request permission to reproduce figures and diagrams provided correct acknowledgement is given. If you want to reproduce the whole article in a third-party publication (excluding your thesis/dissertation for which permission is not required) please go to the Copyright Clearance Center request page.

Read more about how to correctly acknowledge RSC content.

Social activity

Spotlight

Advertisements