Quantum diffusion in liquid water from ring polymer molecular dynamics

We have used the ring polymer molecular-dynamics method to study the translational and orientational motions in an extended simple point charge model of liquid water under ambient conditions. We find, in agreement with previous studies, that quantum-mechanical effects increase the self-diffusion coefficient D and decrease the relaxation times around the principal axes of the water molecule by a factor of around 1.5. These results are consistent with a simple Stokes-Einstein picture of the molecular motion and suggest that the main effect of the quantum fluctuations is to decrease the viscosity of the liquid by about a third. We then go on to consider the system-size scaling of the calculated self-diffusion coefficient and show that an appropriate extrapolation to the limit of infinite system size increases D by a further factor of around 1.3 over the value obtained from a simulation of a system containing 216 water molecules. These findings are discussed in light of the widespread use of classical molecular-dynamics simulations of this sort of size to model the dynamics of aqueous systems

Quantum diffusion in liquid para-hydrogen from ring-polymer molecular dynamics

We have used the ring-polymer molecular dynamics method to calculate approximate Kubo-transformed velocity autocorrelation functions and self-diffusion coefficients for low-pressure liquid para-hydrogen at temperatures of 25 and 14 K. The resulting diffusion coefficients are shown to be consistent with experimental shear viscosities and the established finite-size relation D(L)~=D([infinity])–2.837kBT/6pietaL, where kB is the Boltzmann constant, T the absolute temperature, eta the shear viscosity, and L the length of the (cubic) simulation cell. The diffusion coefficients D(L) obtained in simulations with finite system sizes are therefore too small. However, the extrapolation to infinite system size corrects this deficiency and leads to excellent agreement with experimental results. This both demonstrates the influence of system-size effects on quantum mechanical diffusion coefficients and provides further evidence that ring-polymer molecular dynamics is an accurate as well as practical way of including quantum effects in condensed phase molecular dynamics

Comment on "A centroid molecular dynamics study of liquid para hydrogen and ortho deuterium" [J. Chem. Phys. 121, 6412 (2004)]

We show that the two phase points considered in the recent simulations of liquid para hydrogen by Hone and Voth lie in the liquid-vapor coexistence region of a purely classical molecular dynamics simulation. By contrast, their phase point for ortho deuterium was in the one-phase liquid region for both classical and quantum simulations. These observations are used to account for their report that quantum mechanical effects enhance the diffusion in liquid para hydrogen and decrease it in ortho deuterium