Quantum diffusion in liquid para-hydrogen from ring-polymer molecular dynamics
Abstract
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.
Additional Information
© 2005 American Institute of Physics. Received 13 December 2004; accepted 25 February 2005; published 6 May 2005. The authors would like to thank Paul Madden for stimulating discussions, David Chandler for helpful comments on an early draft of this manuscript, and Joel Bowman for providing them with computer time on a linux cluster for their largest (864 molecule) RPMD calculations.Attached Files
Published - MILjcp05d.pdf
Files
Name | Size | Download all |
---|---|---|
md5:83aabe78707f68a61e2fd6db6ea630e9
|
108.0 kB | Preview Download |
Additional details
- Eprint ID
- 11063
- Resolver ID
- CaltechAUTHORS:MILjcp05d
- Created
-
2008-06-25Created from EPrint's datestamp field
- Updated
-
2021-11-08Created from EPrint's last_modified field